[petsc-users] Slepc - Non Hermitian eigenvalue problem.

Jose E. Roman jroman at dsic.upv.es
Sat Jan 11 12:02:37 CST 2014


El 11/01/2014, a las 18:40, James A Charles escribió:

> Hello,
> 
> I am trying to solve a non-Hermitian eigenvalue problem. Unless I'm mistaken, the eigenpairs that come out are not orthogonal. Is there an option to make them orthonormal? I just output the file to MATLAB format and do dot(V(:,1),V(:,2) and that does not come to 0.
> Am I missing something?
> 
> 
> Thanks,
> James

Eigenvector are not mutually orthogonal in general. This is the case for Hermitian problems, but not in the non-Hermitian case. Schur vectors are always orthogonal, and these can be retrived with EPSGetInvariantSubspace().

Jose



More information about the petsc-users mailing list