Skip to main content
5 events
when toggle format what by license comment
Feb 13, 2016 at 8:21 comment added Federico Poloni Yes, QZ uses no inverses, it is backward stable, and it is implemented in most scientific computing software (Matlab, R, scipy...). The only problem with it is that it is not backward stable in a way that preserves symmetry, so the computed eigenvalues may have a small imaginary part instead of being real (but I think that if you have formed the product $\Sigma_{XX}$ you have already lost hope of preserving realness). Maybe you can fix this, too, with a generalized SVD, if I understand your notation correctly. Check Section 8.7 of Golub-Van Loan, Matrix computations, 4th ed.
Feb 12, 2016 at 23:11 comment added user41037 The issue is also that both matrices $\Sigma_{XX}$ and $\Sigma_{XY}\Sigma_{YX}$ are ill-conditioned.
Feb 12, 2016 at 21:55 comment added user41037 It is a problem with large matrices but they are not sparse. And thank you for QZ, i did not know this method. You mean that with QZ method, one can solved this generalized eigenvalue problem whithout adding the regularization parameter ?
Feb 12, 2016 at 21:13 comment added Federico Poloni Is this a problem with large and sparse matrices? Because otherwise you shouldn't be doing inverses to compute eigenvalues, but rather use a backward stable method such as QZ, which avoids the issue completely.
Feb 12, 2016 at 18:55 history asked user41037 CC BY-SA 3.0