2020
DOI: 10.48550/arxiv.2011.12573
|View full text |Cite
Preprint
|
Sign up to set email alerts
|

On a fast and nearly division-free algorithm for the characteristic polynomial

Abstract: We review the Preparata-Sarwate algorithm, a simple O(n 3.5 ) method for computing the characteristic polynomial, determinant and adjugate of an n×n matrix using only ring operations together with exact divisions by small integers. The algorithm is a baby-step giant-step version of the more well-known Faddeev-Leverrier algorithm. We make a few comments about the algorithm and evaluate its performance empirically.

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
1
1
1

Citation Types

0
5
0

Year Published

2021
2021
2024
2024

Publication Types

Select...
2
1

Relationship

0
3

Authors

Journals

citations
Cited by 3 publications
(5 citation statements)
references
References 22 publications
0
5
0
Order By: Relevance
“…Berkowitz algorithm [49] is faster), it is a rather simple and general way to solve the inverse of a polynomial matrix problem. Despite the recursive nature of FLA, it can be easily modified to carry out the N matrix multiplications in parallel [46,47,50,51].…”
Section: Discussionmentioning
confidence: 99%
See 1 more Smart Citation
“…Berkowitz algorithm [49] is faster), it is a rather simple and general way to solve the inverse of a polynomial matrix problem. Despite the recursive nature of FLA, it can be easily modified to carry out the N matrix multiplications in parallel [46,47,50,51].…”
Section: Discussionmentioning
confidence: 99%
“…For larger matrix dimensions or higher degree polynomials that the ones analyzed in this paper, FLA might suffer from numerical instability in the computation of the polynomial coefficients due to accumulated errors in the trace in Eq. (10) and from the recursive nature of the successive polynomial coefficients [45,46]. However in Ref.…”
Section: Comparison With Recursive Approachesmentioning
confidence: 99%
“…Berkowitz algorithm [41] is faster), it is a rather simple and general way to solve the inverse of a polynomial matrix problem. Despite the recursive nature of FLA, it can be easily modified to carry out the N matrix multiplications in parallel [38,39,[42][43][44].…”
Section: Discussionmentioning
confidence: 99%
“…For larger matrix dimensions or higher degree polynomials that the ones analyzed in this paper, FLA might suffer from numerical instability in the computation of the polynomial coefficients due to accumulated errors in the trace in Eq. (10) and from the recursive nature of the successive polynomial coefficients [37,38]. However in Ref.…”
Section: ĝRcmentioning
confidence: 99%
“…This algorithm is poorly known and has been rediscovered several times (see e.g. [Joh20]). We adapt this algorithm for secure MPC using classical techniques.…”
Section: Our Contributionsmentioning
confidence: 99%