left_evectors

template<typename Type >

DenseTensor<2, Scalar> syten::DenseEigensolver::DenseEigensolverBase< Type >::left_evectors

Matrix of left eigenvectors, each row contains one complex-conjugated vector.

Specifically, conj(left_evectors)·M is the same as D·conj(left_evectors) or in code:

DenseTensor<2> mat({10,10}); makeRandom(mat);
auto solver = DenseEigensolver::DenseEigensolver<SDef>(transpose(mat));
DenseTensor<2> dia({10, 10});
for(Index i(0); i != 10; ++i) { dia[{i,i}] = solver.evalues[i]; }
auto lvecs = conj(solver.left_evectors);
 
auto lm = prodD<1>(lvecs, mat, {-1, 1}, {1, -2});
auto dl = prodD<1>(dia, lvecs, {-1, 1}, {1, -2}); // dl == lm