Joint computation of the left and right eigenvectors using the nonsymmetric_lanczos_tridiagonalization works but can rapidly run into floating point accuracy problems. This mainly comes from the necessity of satisfying the bi-orthogonality condition. Possibly remedy is actually the two-sided Arnoldi decomposition, see here.
Joint computation of the left and right eigenvectors using the
nonsymmetric_lanczos_tridiagonalization
works but can rapidly run into floating point accuracy problems. This mainly comes from the necessity of satisfying the bi-orthogonality condition. Possibly remedy is actually the two-sided Arnoldi decomposition, see here.