on computing the jordan structure of totally nonnegative ... · on computing the jordan structure...
TRANSCRIPT
![Page 1: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/1.jpg)
On computing the Jordan structure of Totally NonnegativeMatrices with high relative accuracy
Teresa Laudadio (IAC–CNR Bari),Nicola Mastronardi (IAC–CNR Bari),
Paul Van Dooren ( UCLouvain)
Due giorni di Algebra Lineare NumericaRoma, 18–19 febbraio, 2019
![Page 2: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/2.jpg)
On computing the Jordan structure of Totally NonnegativeMatrices with high relative accuracy
The revenge of LR method
Teresa Laudadio (IAC–CNR Bari),Nicola Mastronardi (IAC–CNR Bari),
Paul Van Dooren ( UCLouvain)
Due giorni di Algebra Lineare NumericaRoma, 18–19 febbraio, 2019
![Page 3: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/3.jpg)
Table of contents
I Eigenvalue computation: LR vs QR method
I Totally Nonnegative (TN) matrices
I Computation of the Jordan form of TN matrices
I Numerical Examples
![Page 4: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/4.jpg)
Eigenvalue computation: LR vs QR method
“... the QR algorithm is the ubiquitous tool for computingeigenvalues of dense matrices1. Its predecessor, the LR algorithm,is now largely forgotten and rarely taught to students. What wewish to say here is that, from an intellectual viewpoint, it was theLR algorithm that made the seminal contribution, and QR is animproved stable version of LR.”
1M. H. Gutknecht and B. N. Parlett, From qd to LR, or, how were the qdand LR algorithms discovered?, IMA Journal of Numerical Analysis, vol. 31,no. 3, pp. 741-754, 2011.
![Page 5: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/5.jpg)
Eigenvalue computation: LR vs QR method
“... the QR algorithm is the ubiquitous tool for computingeigenvalues of dense matrices1. Its predecessor, the LR algorithm,is now largely forgotten and rarely taught to students. What wewish to say here is that, from an intellectual viewpoint, it was theLR algorithm that made the seminal contribution, and QR is animproved stable version of LR.”
1M. H. Gutknecht and B. N. Parlett, From qd to LR, or, how were the qdand LR algorithms discovered?, IMA Journal of Numerical Analysis, vol. 31,no. 3, pp. 741-754, 2011.
![Page 6: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/6.jpg)
Eigenvalue computation: the LR method
I A1 = A ∈ Cn×n
for i=1,2,. . .LiRi = Ai
Ai+1 = RiLiend
I Li lower triangular, Li (j , j) = 1, j = 1, . . . , n,Ri upper triangular.
I If the factorization LiRi = Ai exists for any i , the matrices ofthe sequence {Ai}i∈N are similar to A.
I Theorem2. Let Fk = L1L2 · · · Lk . If the sequence {Fk}k∈Nconverges as k →∞, the the sequence {Ak}k∈N convergestoward an upper triangular matrix (symmetric case)
2H. R. Schwarz, H. Rutishauser, E. Stiefel, Numerical analysis of symmetricmatrices, Englewood Cliffs, N.J., Prentice-Hall, 1973
![Page 7: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/7.jpg)
Eigenvalue computation: the LR method
I A1 = A ∈ Cn×n
for i=1,2,. . .LiRi = Ai
Ai+1 = RiLiend
I Li lower triangular, Li (j , j) = 1, j = 1, . . . , n,Ri upper triangular.
I If the factorization LiRi = Ai exists for any i , the matrices ofthe sequence {Ai}i∈N are similar to A.
I Theorem2. Let Fk = L1L2 · · · Lk . If the sequence {Fk}k∈Nconverges as k →∞, the the sequence {Ak}k∈N convergestoward an upper triangular matrix (symmetric case)
2H. R. Schwarz, H. Rutishauser, E. Stiefel, Numerical analysis of symmetricmatrices, Englewood Cliffs, N.J., Prentice-Hall, 1973
![Page 8: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/8.jpg)
Eigenvalue computation: the LR method
I A1 = A ∈ Cn×n
for i=1,2,. . .LiRi = Ai
Ai+1 = RiLiend
I Li lower triangular, Li (j , j) = 1, j = 1, . . . , n,Ri upper triangular.
I If the factorization LiRi = Ai exists for any i , the matrices ofthe sequence {Ai}i∈N are similar to A.
I Theorem2. Let Fk = L1L2 · · · Lk . If the sequence {Fk}k∈Nconverges as k →∞, the the sequence {Ak}k∈N convergestoward an upper triangular matrix (symmetric case)
2H. R. Schwarz, H. Rutishauser, E. Stiefel, Numerical analysis of symmetricmatrices, Englewood Cliffs, N.J., Prentice-Hall, 1973
![Page 9: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/9.jpg)
Eigenvalue computation: the LR method
I A1 = A ∈ Cn×n
for i=1,2,. . .LiRi = Ai
Ai+1 = RiLiend
I Li lower triangular, Li (j , j) = 1, j = 1, . . . , n,Ri upper triangular.
I If the factorization LiRi = Ai exists for any i , the matrices ofthe sequence {Ai}i∈N are similar to A.
I Theorem2. Let Fk = L1L2 · · · Lk . If the sequence {Fk}k∈Nconverges as k →∞, the the sequence {Ak}k∈N convergestoward an upper triangular matrix (symmetric case)
2H. R. Schwarz, H. Rutishauser, E. Stiefel, Numerical analysis of symmetricmatrices, Englewood Cliffs, N.J., Prentice-Hall, 1973
![Page 10: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/10.jpg)
Eigenvalue computation: the QR method
I A1 = A ∈ Rn×n
for i=1,2,. . .QiRi = Ai
Ai+1 = RiQi
end
I Qi unitary, Ri upper triangular.
I The matrices of the sequence {Ai}i∈N are similar to A.
I Teorema3. Sia A ∈ Cn×n con autovalori in modulo tuttidistinti: |λ1| > |λ2| > · · · > |λn| > 0. Sia X la matrice degliautovettori di A, tale che A = XΛX−1, e si supponga cheX−1 ammetta la fattorizzazione LU. Allora esistono matrici difase Sk tali che
limk→∞
SHk RkSk−1 = lim
k→∞SHk−1AkSk−1 = T , lim
k→∞SHk QkSk−1 = I .
3D. Bini, M. Capovani, O. Menchi, Metodi Numerici per l’Algebra Lineare,Zanichelli, 1988
![Page 11: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/11.jpg)
Eigenvalue computation: the QR method
I A1 = A ∈ Rn×n
for i=1,2,. . .QiRi = Ai
Ai+1 = RiQi
end
I Qi unitary, Ri upper triangular.
I The matrices of the sequence {Ai}i∈N are similar to A.
I Teorema3. Sia A ∈ Cn×n con autovalori in modulo tuttidistinti: |λ1| > |λ2| > · · · > |λn| > 0. Sia X la matrice degliautovettori di A, tale che A = XΛX−1, e si supponga cheX−1 ammetta la fattorizzazione LU. Allora esistono matrici difase Sk tali che
limk→∞
SHk RkSk−1 = lim
k→∞SHk−1AkSk−1 = T , lim
k→∞SHk QkSk−1 = I .
3D. Bini, M. Capovani, O. Menchi, Metodi Numerici per l’Algebra Lineare,Zanichelli, 1988
![Page 12: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/12.jpg)
Eigenvalue computation: the QR method
I A1 = A ∈ Rn×n
for i=1,2,. . .QiRi = Ai
Ai+1 = RiQi
end
I Qi unitary, Ri upper triangular.
I The matrices of the sequence {Ai}i∈N are similar to A.
I Teorema3. Sia A ∈ Cn×n con autovalori in modulo tuttidistinti: |λ1| > |λ2| > · · · > |λn| > 0. Sia X la matrice degliautovettori di A, tale che A = XΛX−1, e si supponga cheX−1 ammetta la fattorizzazione LU. Allora esistono matrici difase Sk tali che
limk→∞
SHk RkSk−1 = lim
k→∞SHk−1AkSk−1 = T , lim
k→∞SHk QkSk−1 = I .
3D. Bini, M. Capovani, O. Menchi, Metodi Numerici per l’Algebra Lineare,Zanichelli, 1988
![Page 13: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/13.jpg)
Eigenvalue computation: the QR method
I A1 = A ∈ Rn×n
for i=1,2,. . .QiRi = Ai
Ai+1 = RiQi
end
I Qi unitary, Ri upper triangular.
I The matrices of the sequence {Ai}i∈N are similar to A.
I Teorema3. Sia A ∈ Cn×n con autovalori in modulo tuttidistinti: |λ1| > |λ2| > · · · > |λn| > 0. Sia X la matrice degliautovettori di A, tale che A = XΛX−1, e si supponga cheX−1 ammetta la fattorizzazione LU. Allora esistono matrici difase Sk tali che
limk→∞
SHk RkSk−1 = lim
k→∞SHk−1AkSk−1 = T , lim
k→∞SHk QkSk−1 = I .
3D. Bini, M. Capovani, O. Menchi, Metodi Numerici per l’Algebra Lineare,Zanichelli, 1988
![Page 14: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/14.jpg)
LR vs QR method: stability and convergence
I The QR method is backward stable.
I The speed of convergence of both methods is linear.
I Shift strategies speed up the methods (quadraticconvergence) (cubic in the Hermitian case).
I Drawback: high relative accuracy can be lost in computingthe eigenvalues4.
I Recommendation : avoid subtraction between number withsame sign!
4J. Demmel,W. Kahan, Accurate Singular Values of Bidiagonal Matrices,1990, SIAM J. Sci. and Stat. Comput., 11(5), . 873-912.
![Page 15: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/15.jpg)
LR vs QR method: stability and convergence
I The QR method is backward stable.
I The speed of convergence of both methods is linear.
I Shift strategies speed up the methods (quadraticconvergence) (cubic in the Hermitian case).
I Drawback: high relative accuracy can be lost in computingthe eigenvalues4.
I Recommendation : avoid subtraction between number withsame sign!
4J. Demmel,W. Kahan, Accurate Singular Values of Bidiagonal Matrices,1990, SIAM J. Sci. and Stat. Comput., 11(5), . 873-912.
![Page 16: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/16.jpg)
LR vs QR method: stability and convergence
I The QR method is backward stable.
I The speed of convergence of both methods is linear.
I Shift strategies speed up the methods (quadraticconvergence) (cubic in the Hermitian case).
I Drawback: high relative accuracy can be lost in computingthe eigenvalues4.
I Recommendation : avoid subtraction between number withsame sign!
4J. Demmel,W. Kahan, Accurate Singular Values of Bidiagonal Matrices,1990, SIAM J. Sci. and Stat. Comput., 11(5), . 873-912.
![Page 17: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/17.jpg)
LR vs QR method: stability and convergence
I The QR method is backward stable.
I The speed of convergence of both methods is linear.
I Shift strategies speed up the methods (quadraticconvergence) (cubic in the Hermitian case).
I Drawback: high relative accuracy can be lost in computingthe eigenvalues4.
I Recommendation : avoid subtraction between number withsame sign!
4J. Demmel,W. Kahan, Accurate Singular Values of Bidiagonal Matrices,1990, SIAM J. Sci. and Stat. Comput., 11(5), . 873-912.
![Page 18: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/18.jpg)
LR vs QR method: stability and convergence
I The QR method is backward stable.
I The speed of convergence of both methods is linear.
I Shift strategies speed up the methods (quadraticconvergence) (cubic in the Hermitian case).
I Drawback: high relative accuracy can be lost in computingthe eigenvalues4.
I Recommendation : avoid subtraction between number withsame sign!
4J. Demmel,W. Kahan, Accurate Singular Values of Bidiagonal Matrices,1990, SIAM J. Sci. and Stat. Comput., 11(5), . 873-912.
![Page 19: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/19.jpg)
Condition of an eigenvalue λ
The “condition” of the eigenvalue λ of A ∈ Cn×n is defined as
s(λ) =1
|yHx |
with x and y the right and left normalized eigenvectors of Aassociated to λ. “Roughly speaking, this means that O(ε)perturbation in A can induce ε/s(λ) changes in an eigenvalue5”
5G.H. Golub, C.F. Van Loan, Matrix Computations, 4th ed., The JohnsHopkins University Press Baltimore, 2013.
![Page 20: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/20.jpg)
Condition of an eigenvalue λ
The “condition” of the eigenvalue λ of A ∈ Cn×n is defined as
s(λ) =1
|yHx |
with x and y the right and left normalized eigenvectors of Aassociated to λ. “Roughly speaking, this means that O(ε)perturbation in A can induce ε/s(λ) changes in an eigenvalue5”
5G.H. Golub, C.F. Van Loan, Matrix Computations, 4th ed., The JohnsHopkins University Press Baltimore, 2013.
![Page 21: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/21.jpg)
Example: Bessel’s polynomials
Eigenvalues6 of a generalized Bessel matrix of order 20, B(4.5,2)20 :
(◦), exact by Mathematica,(·) approximations by Matlab with the matrix,
(+), approximations by Matlab with the matrix transpose.6L. Pasquini, Accurate computation of the zeros of the generalized Bessel polynomials, Numer. Math., 2000,
86, pp. 507-538.C. Ferreira, B. Parlett, F.M. Dopico, Sensitivity of eigenvalues of an unsymmetric tridiagonal matrix, Numer.Math., 2012, 122(3), pp. 527-555.
![Page 22: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/22.jpg)
Condition of an eigenvalue λ. Hermitian caseTheorem7 Let A and A = A + E be Hermitian. Let theeigenvalues of A, A and E be
λ1 ≥ · · · ≥ λn, λ1 ≥ · · · ≥ λn, ε1 ≥ · · · ≥ εn.
Thenλ1 + εn ≤ λi ≤ λi + ε1, i = 1, . . . , n.
Consequently,
|λi − λi | ≤ ‖E‖2, i = 1, . . . , n. (1)
In the terminology of perturbation theory, (1) says that theeigenvalues of a Hermitain matrix are perfectly conditioned ...
unfortunately this does not imply that all the eigenvalues aredetermined to high relative accuracy. In fact, the relative error inλi is
|λi − λi ||λi |
≤ ‖E‖2|λi |
.
7G. W. Stewart, Matrix Algorithms: Vol. 2, Eigensystems, SIAM,Philadelphia, 2001
![Page 23: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/23.jpg)
Condition of an eigenvalue λ. Hermitian caseTheorem7 Let A and A = A + E be Hermitian. Let theeigenvalues of A, A and E be
λ1 ≥ · · · ≥ λn, λ1 ≥ · · · ≥ λn, ε1 ≥ · · · ≥ εn.
Thenλ1 + εn ≤ λi ≤ λi + ε1, i = 1, . . . , n.
Consequently,
|λi − λi | ≤ ‖E‖2, i = 1, . . . , n. (1)
In the terminology of perturbation theory, (1) says that theeigenvalues of a Hermitain matrix are perfectly conditioned ...unfortunately this does not imply that all the eigenvalues aredetermined to high relative accuracy. In fact, the relative error inλi is
|λi − λi ||λi |
≤ ‖E‖2|λi |
.
7G. W. Stewart, Matrix Algorithms: Vol. 2, Eigensystems, SIAM,Philadelphia, 2001
![Page 24: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/24.jpg)
Condition of an eigenvalue λ. Hermitian caseTheorem7 Let A and A = A + E be Hermitian. Let theeigenvalues of A, A and E be
λ1 ≥ · · · ≥ λn, λ1 ≥ · · · ≥ λn, ε1 ≥ · · · ≥ εn.
Thenλ1 + εn ≤ λi ≤ λi + ε1, i = 1, . . . , n.
Consequently,
|λi − λi | ≤ ‖E‖2, i = 1, . . . , n. (1)
In the terminology of perturbation theory, (1) says that theeigenvalues of a Hermitain matrix are perfectly conditioned ...unfortunately this does not imply that all the eigenvalues aredetermined to high relative accuracy. In fact, the relative error inλi is
|λi − λi ||λi |
≤ ‖E‖2|λi |
.
7G. W. Stewart, Matrix Algorithms: Vol. 2, Eigensystems, SIAM,Philadelphia, 2001
![Page 25: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/25.jpg)
Totally Nonnegative Matrices: eigenvalue computation
Revenge of the LR method
A modification of LR method based on the Neville methodcomputes the eigenvalues of a totally nonnegative matrix A withhigh relative accuracy supposed A is factored as a product ofparticular bidiagonal matrices8
8P. Koev,Accurate eigenvalues and SVDs of totally nonnegative matrices,SIAM J. Matrix Anal. Appl., 27(1), pp. 1–23, 2005.
![Page 26: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/26.jpg)
Totally Nonnegative Matrices: eigenvalue computation
Revenge of the LR method
A modification of LR method based on the Neville methodcomputes the eigenvalues of a totally nonnegative matrix A withhigh relative accuracy supposed A is factored as a product ofparticular bidiagonal matrices8
8P. Koev,Accurate eigenvalues and SVDs of totally nonnegative matrices,SIAM J. Matrix Anal. Appl., 27(1), pp. 1–23, 2005.
![Page 27: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/27.jpg)
Totally Nonnegative Matrices
I A matrix is totally positive if each of its minors is positive.
I A matrix is totally nonnegative if each of its minors isnonnegative
I C. Cryer. Some properties of totally positive matrices. Linear Algebra Appl., 15,1–25, 1976.
I Shaun M. Fallat and Charles R. Johnson. Totally nonnegative matrices. Princeton Series in AppliedMathematics. Princeton University Press, Princeton, NJ, 2011.
I M. Gasca and J. M. Pena. Total positivity and Neville elimination. Linear Algebra Appl., 165,25–44, 1992.
I P. Koev, Accurate computations with totally nonnegative matrices. SIAM J. Matrix Anal. Appl.,29,731–751, 2007.
I P. Koev, Accurate Eigenvalues and Exact Zero Jordan Blocks of Totally Nonnegative Matrices, 2018.
![Page 28: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/28.jpg)
Totally Nonnegative Matrices
I A matrix is totally positive if each of its minors is positive.
I A matrix is totally nonnegative if each of its minors isnonnegative
I C. Cryer. Some properties of totally positive matrices. Linear Algebra Appl., 15,1–25, 1976.
I Shaun M. Fallat and Charles R. Johnson. Totally nonnegative matrices. Princeton Series in AppliedMathematics. Princeton University Press, Princeton, NJ, 2011.
I M. Gasca and J. M. Pena. Total positivity and Neville elimination. Linear Algebra Appl., 165,25–44, 1992.
I P. Koev, Accurate computations with totally nonnegative matrices. SIAM J. Matrix Anal. Appl.,29,731–751, 2007.
I P. Koev, Accurate Eigenvalues and Exact Zero Jordan Blocks of Totally Nonnegative Matrices, 2018.
![Page 29: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/29.jpg)
Totally nonnegative matrices
A totally nonnegative (TN) matrix can be factorized as theproduct of lower and upper bidiagonal (BD) TN matricesExample
1 2 41 3 91 4 16
=
111 1
11 1
1 1
××
11
2
1 21 3
1
11 2
1
= L2L1DU1U2
![Page 30: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/30.jpg)
Totally nonnegative matrices
A totally nonnegative (TN) matrix can be factorized as theproduct of lower and upper bidiagonal (BD) TN matricesExample
1 2 41 3 91 4 16
=
111 1
11 1
1 1
××
11
2
1 21 3
1
11 2
1
= L2L1DU1U2
![Page 31: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/31.jpg)
Storage of the entries of the bidiagonal matrices
⇓
![Page 32: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/32.jpg)
Storage of the entries of the bidiagonal matrices
![Page 33: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/33.jpg)
Storage of the entries of the bidiagonal matrices
![Page 34: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/34.jpg)
Storage of the entries of the bidiagonal matrices
![Page 35: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/35.jpg)
Storage of the entries of the bidiagonal matrices
![Page 36: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/36.jpg)
Storage of the entries of the bidiagonal matrices
![Page 37: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/37.jpg)
Storage of the entries of the bidiagonal matrices
![Page 38: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/38.jpg)
Storage of the entries of the bidiagonal matrices
![Page 39: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/39.jpg)
Storage of the entries of the bidiagonal matrices
![Page 40: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/40.jpg)
Storage of the entries of the bidiagonal matrices
![Page 41: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/41.jpg)
Examples of Totally nonnegative matricesPascal matrix9 :
A =
1 1 1 1 1 11 2 3 4 5 61 3 6 10 15 211 4 10 20 35 561 5 15 35 70 1261 6 21 56 126 252
Bidiagonal representation :
B =
1 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 1
9L. Aceto, D. Trigiante, The Matrices of Pascal and Other Greats, The
American Mathematical Monthly,108, 2001, pp. 232–245.
![Page 42: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/42.jpg)
Examples of Totally nonnegative matricesPascal matrix9 :
A =
1 1 1 1 1 11 2 3 4 5 61 3 6 10 15 211 4 10 20 35 561 5 15 35 70 1261 6 21 56 126 252
Bidiagonal representation :
B =
1 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 1
9L. Aceto, D. Trigiante, The Matrices of Pascal and Other Greats, The
American Mathematical Monthly,108, 2001, pp. 232–245.
![Page 43: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/43.jpg)
Example: Pascal matrices
46994838541.8026473129040614264498618434813808464 4.699483854180264e + 10 4.699483854180265e + 10554956959.8629315601218316654368527363981179735 5.549569598629313e + 08 5.549569598629270e + 0813982719.8768657579641932785550920189007397791 1.398271987686575e + 07 1.398271987686483e + 07
566741.0592514715502709682155302349922279150 5.667410592514719e + 05 5.667410592512758e + 0533061.3738535312656510304726125038782408082 3.306137385353126e + 04 3.306137385339706e + 042624.5205773818364234350310755556790308814 2.624520577381835e + 03 2.624520577264530e + 03275.3624643918710918790678389797005985188 2.753624643918711e + 02 2.753624643088091e + 0237.7838289328834004091152152207125557078 3.778382893288337e + 01 3.778382888204030e + 016.8607164608962072080859087566867594124 6.860716460896203e + 00 6.860716432387068e + 001.7487615421998684648886862656955835553 1.748761542199869e + 00 1.748761527349773e + 000.5718332522008930167280076975865210505 5.718332522008932e − 01 5.718332416239250e − 010.1457573717992378878694930414593258721 1.457573717992378e − 01 1.457573629317858e − 010.0264663489181133903334588757511631763 2.646634891811337e − 02 2.646634328298966e − 020.0036315770277857839944323164609232541 3.631577027785781e − 03 3.631574328494989e − 030.0003810219697334504586080098453478379 3.810219697334519e − 04 3.810212348803231e − 040.0000302467769316002152635002881173181 3.024677693160011e − 05 3.024699622580157e − 050.0000017644742403537149182711011678706 1.764474240353710e − 06 1.764921336162508e − 06
0.00000007151684427680540072441234011313 7.151684427680527e − 08 7.183315146869231e − 080.00000000180194154200172445069224455960 1.801941542001719e − 09 1.948297884304045e − 090.00000000002127893256001899627761173655 2.127893256001896e − 11 7.542297033530567e − 11
Exact eigenvalues (first column), Koev’s method (second column),QR method (eig of matlab) (third column).
![Page 44: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/44.jpg)
Examples of Totally nonnegative matricesVariance–covariance (semi–separable) matrix of order-statistics fora random sample of size k from a population with exponentialdensity10
A =
1k2
1k2
1k2 · · · 1
k2
1k2
∑2j=1
1(k−j+1)2
∑2j=1
1(k−j+1)2
· · ·∑2
j=11
(k−j+1)2
1k2
∑2j=1
1(k−j+1)2
∑3j=1
1(k−j+1)2
· · ·∑3
j=11
(k−j+1)2
......
.... . .
...1k2
∑2j=1
1(k−j+1)2
∑3j=1
1(k−j+1)2
· · ·∑k
j=11
(k−j+1)2
Bidiagonal representation :
B =
1k2 1 · · · 1 11 1
(k−1)2...
. . .
1 122
1 112
10
F.A. Graybill, Matrices with Applications in Statistics, 2nd ed., Brooks/Cole, 1983.
![Page 45: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/45.jpg)
Examples of Totally nonnegative matricesVariance–covariance (semi–separable) matrix of order-statistics fora random sample of size k from a population with exponentialdensity10
A =
1k2
1k2
1k2 · · · 1
k2
1k2
∑2j=1
1(k−j+1)2
∑2j=1
1(k−j+1)2
· · ·∑2
j=11
(k−j+1)2
1k2
∑2j=1
1(k−j+1)2
∑3j=1
1(k−j+1)2
· · ·∑3
j=11
(k−j+1)2
......
.... . .
...1k2
∑2j=1
1(k−j+1)2
∑3j=1
1(k−j+1)2
· · ·∑k
j=11
(k−j+1)2
Bidiagonal representation :
B =
1k2 1 · · · 1 11 1
(k−1)2...
. . .
1 122
1 112
10
F.A. Graybill, Matrices with Applications in Statistics, 2nd ed., Brooks/Cole, 1983.
![Page 46: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/46.jpg)
Examples of Totally nonnegative matricesA correlation (Toeplitz semi–separable) matrix arising in the studyof Time Series and Markov chains is11
A =
1 ρ ρ2 · · · ρn−2 ρn−1
ρ 1 ρ ρ2. . . ρn−2
ρ2 ρ 1. . .
. . ....
.... . .
. . .. . .
. . . ρ2
ρn−2. . .
. . . ρ 1 ρρn−1 ρn−2 · · · ρ2 ρ 1
Bidiagonal representation :
B =
1− ρ2 ρ ρ · · · ρρ 1− ρ2ρ 1− ρ2...
. . .
ρ 1− ρ2
11M.J.R. Healy, Matrices for Statistics, Oxford University press, 2000.
![Page 47: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/47.jpg)
Examples of Totally nonnegative matricesA correlation (Toeplitz semi–separable) matrix arising in the studyof Time Series and Markov chains is11
A =
1 ρ ρ2 · · · ρn−2 ρn−1
ρ 1 ρ ρ2. . . ρn−2
ρ2 ρ 1. . .
. . ....
.... . .
. . .. . .
. . . ρ2
ρn−2. . .
. . . ρ 1 ρρn−1 ρn−2 · · · ρ2 ρ 1
Bidiagonal representation :
B =
1− ρ2 ρ ρ · · · ρρ 1− ρ2ρ 1− ρ2...
. . .
ρ 1− ρ2
11M.J.R. Healy, Matrices for Statistics, Oxford University press, 2000.
![Page 48: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/48.jpg)
Examples of Totally nonnegative matrices
Vandermonde matrices with ordered knots
V =
1 1 · · · 1 1
x1 x2... xn−1 xn
x21 x22... x2n−1 x2n
......
......
...
xn−21 xn−22
... xn−2n−1 xn−2n
xn−11 xn−12
... xn−1n−1 xn−1n
x1 ≥ x2 ≥ · · · ≥ xn−1 ≥ xn.
![Page 49: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/49.jpg)
Totally nonnegative matrices
Jordan blocks corresponding to zero eigenvalues
I n − rank(A) = # zero eigenvalues
I rank(Ai−1)− rank(Ai ) = # of Jordan blocks of size i
I A2 is TN (as a product of TN) and its BD is a TN-preservingop, thus BD accurate
I need to form BD of A2; . . . ,An, a potential O(n4) algorithm
![Page 50: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/50.jpg)
Totally nonnegative matrices
Sketch of t Koev’s algorithm
r0 = n; i = 1;Ai = A;ri = rank(Ai );while i < n and ri−1 − ri > 0
i = i + 1;Ai = Ai−1 × A;ri = rank(Ai );
end
![Page 51: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/51.jpg)
Totally nonnegative matrices
The proposed algorithm is based on the reduction of the TNmatrix A to a similar Hessenberg one via Neville eliminationsupposed A factored as a product of bidiagonal matrices.
![Page 52: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/52.jpg)
Reduction of a TN matrix to Hessenberg form 1
![Page 53: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/53.jpg)
Reduction of a TN matrix to Hessenberg form 2
![Page 54: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/54.jpg)
Reduction of a TN matrix to Hessenberg form 3
![Page 55: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/55.jpg)
Reduction of a TN matrix to Hessenberg form 4
![Page 56: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/56.jpg)
Reduction of a TN matrix to Hessenberg form 5
![Page 57: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/57.jpg)
Reduction of a TN matrix to Hessenberg form 6
![Page 58: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/58.jpg)
Reduction of a TN matrix to Hessenberg form 7
![Page 59: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/59.jpg)
Reduction of a TN matrix to Hessenberg form 8
![Page 60: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/60.jpg)
Reduction of a TN matrix to Hessenberg form 9
![Page 61: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/61.jpg)
Reduction of a TN matrix to Hessenberg form 10
![Page 62: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/62.jpg)
Reduction of a TN matrix to Hessenberg form 11
![Page 63: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/63.jpg)
Reduction of a TN matrix to Hessenberg form 12
![Page 64: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/64.jpg)
Reduction of a TN matrix to Hessenberg form 13
![Page 65: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/65.jpg)
Reduction of a TN matrix to Hessenberg form 14
![Page 66: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/66.jpg)
Reduction of a TN matrix to Hessenberg form 15
![Page 67: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/67.jpg)
Reduction of a TN matrix to Hessenberg form 16
![Page 68: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/68.jpg)
Reduction of a TN matrix to Hessenberg form 17
![Page 69: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/69.jpg)
Reduction of a TN matrix to Hessenberg form 18
![Page 70: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/70.jpg)
Reduction of a TN matrix to Hessenberg form 19
![Page 71: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/71.jpg)
Reduction of a TN matrix to Hessenberg form 20
![Page 72: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/72.jpg)
Reduction of a TN matrix to Hessenberg form 21
![Page 73: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/73.jpg)
Reduction of a TN matrix to Hessenberg form 22
![Page 74: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/74.jpg)
Reduction of a TN matrix to Hessenberg form 23
![Page 75: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/75.jpg)
Reduction of a TN matrix to Hessenberg form 24
![Page 76: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/76.jpg)
Reduction of a TN matrix to Hessenberg form 25
![Page 77: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/77.jpg)
Reduction of a TN matrix to Hessenberg form 26
![Page 78: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/78.jpg)
Reduction of a TN matrix to Hessenberg form 27
![Page 79: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/79.jpg)
Reduction of a TN matrix to Hessenberg form 28
![Page 80: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/80.jpg)
Reduction of a TN matrix to Hessenberg form 29
![Page 81: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/81.jpg)
Reduction of a TN matrix to Hessenberg form 30
![Page 82: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/82.jpg)
Reduction of a TN matrix to Hessenberg form 31
![Page 83: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/83.jpg)
Reduction of a TN matrix to Hessenberg form 32
![Page 84: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/84.jpg)
Reduction of a TN matrix to Hessenberg form 33
![Page 85: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/85.jpg)
Reduction of a TN matrix to Hessenberg form 34
![Page 86: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/86.jpg)
Reduction of a TN matrix to Hessenberg form 35
![Page 87: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/87.jpg)
Reduction of a TN matrix to Hessenberg form 36
![Page 88: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/88.jpg)
Reduction of a TN matrix to Hessenberg form 37
![Page 89: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/89.jpg)
Reduction of a TN matrix to Hessenberg form 37
![Page 90: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/90.jpg)
Numerical exampleThe lower Hessenberg matrix
Hn =
1 1 0 · · · 0
1 1 1. . .
.
.
.
.
.
.
.
.
.. . .
. . . 01 1 · · · 1 11 1 · · · 1 1
n×n
is irreducible and totally nonnegative with rank(Hn) = n − 1.Hn has
⌈n2
⌉distinct positive eigenvalues and exactly one Jordan
block of size⌊n2
⌋corresponding to the eigenvalue zero.
nz = 2090 5 10 15 20
0
2
4
6
8
10
12
14
16
18
20
Sparsity structure of the computed H20
![Page 91: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/91.jpg)
Numerical example
3.918985947228994e + 00 3.918985947228999e + 00 + 0.000000000000000e + 00i3.682507065662362e + 00 3.682507065662364e + 00 + 0.000000000000000e + 00i3.309721467890570e + 00 3.309721467890570e + 00 + 0.000000000000000e + 00i2.830830026003772e + 00 2.830830026003772e + 00 + 0.000000000000000e + 00i2.284629676546570e + 00 2.284629676546571e + 00 + 0.000000000000000e + 00i1.715370323453430e + 00 1.715370323453430e + 00 + 0.000000000000000e + 00i1.169169973996227e + 00 1.169169973996176e + 00 + 0.000000000000000e + 00i6.902785321094297e − 01 6.902785321142250e − 01 + 0.000000000000000e + 00i3.174929343376376e − 01 3.174929292099200e − 01 + 0.000000000000000e + 00i8.101405277100522e − 02 8.199240254858542e − 02 + 0.000000000000000e + 00i0.000000000000000e + 00 5.116195905098703e − 02 + 2.135157259497554e − 02i0.000000000000000e + 00 5.116195905098703e − 02 − 2.135157259497554e − 02i0.000000000000000e + 00 2.537136779504252e − 02 + 4.389486455178316e − 02i0.000000000000000e + 00 2.537136779504252e − 02 − 4.389486455178316e − 02i0.000000000000000e + 00 −4.297088165307477e − 03 + 4.770120243833010e − 02i0.000000000000000e + 00 −4.297088165307477e − 03 − 4.770120243833010e − 02i0.000000000000000e + 00 −2.931310764823774e − 02 + 3.568142654713092e − 02i0.000000000000000e + 00 −2.931310764823774e − 02 − 3.568142654713092e − 02i0.000000000000000e + 00 −4.341230335978553e − 02 + 1.311956655234002e − 02i0.000000000000000e + 00 −4.341230335978553e − 02 − 1.311956655234002e − 02i
Computation of the eigenvalues of H20 : proposed method (firstcolumn), eig of matlab (second column).
![Page 92: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/92.jpg)
Announcement 1
11th Workshop SDS2020STRUCTURAL DYNAMICAL SYSTEMS:
Computational AspectsPorto Giardino Resort
Capitolo (Monopoli) - Bari, ItalyJune 9-12, 2020
Local organizing committeeNicoletta del Buono et al.
![Page 93: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/93.jpg)
Announcement 2
Householder Symposium 2020Hotel SIERRA SILVANA, Selva di Fasano (Br), Italy
June 14–19, 2020Local organizing committee
Nicola Mastronardi (chair) Dario A. BiniFasma Diele Teresa LaudadioCarmela Marangi Beatrice MeiniStefano Serra Capizzano Valeria Simoncini
![Page 94: On computing the Jordan structure of Totally Nonnegative ... · On computing the Jordan structure of Totally Nonnegative Matrices with high relative accuracy The revenge of LR method](https://reader030.vdocument.in/reader030/viewer/2022041102/5edbf247ad6a402d666665c1/html5/thumbnails/94.jpg)
Announcement 3
Sparse DaysSt Girons, France
either the week before or the week after the HouseholderSymposium
Local organizing committeeIain Duff et al.