Schur decomposition via IMSL and Matlab

Schur decomposition via IMSL and Matlab

Post by Abdul Qa » Tue, 23 Aug 2005 08:20:07


I'm trying to reconcile differences between IMSL MP F90's and Matlab's
Schur [U,T] decomposition of a real matrix of moderate size whose spectral
radius and departure from normality are equally large (~10^14) and whose
eigenvalues are observed values lying along the diagonal. In both cases the
norm of the residual in the decomposition is finite (~10^-1).

IMSL's lin_eig_gen correctly computes the eigenvalues. However, they differ
significantly (>> 1ulp) from the diagonal entries in IMSL's T, especially
as they decrease in magnitude. In contrast, Matlab's diagonal entries in T
are indistinguishable from the eigenvalues. In both cases, U is apparently
unitary. Why this difference?
 
 
 

1. schur decomposition and mathematica

2. Does anybody have a Matlab wrapper for IMSL and GSL?

Hi all,

While IMSL and GSL's numerical integration and ODE/PDE functions are faster
and sometimes "smarter" than Matlab's functions, I still like the fast
visualization and prototyping in Matlab.

Are there wrappers for IMSL and GSL for Matlab?

I mainly use numerical integration and ODE/PDE's, and FFTs... etc.

Thanks a lot!

3. Wavelet packet decomposition (sdspdwt) issue in Matlab/Simulink

4. Need Matlab code for Wavelet packet image decomposition

5. Energy preservation is not satisfied for wavelet packet decomposition in MATLAB

6. MATLAB codes for Empirical Mode Decomposition

7. Schur Complement

8. what is schur cohn stability test?

9. Schur Complements

10. Schur-Cohn stability test for digital filters

11. Lehmer-Schur pseudocode?

12. Schur-Hadamand Product

13. Need help sorting out ANDs and ORs in parameter query

14. 3 ANDS + 1 OR

15. IF's & AND's