Hello,

I am trying to use Eigen in a new project and was using the SelfAdjointEigenSolver class. It seems that it does not work with complex Hermitian matrices (which are self-adjoint). Is there a way to make it work? By definition SelfAdjoint operators includes Hermitian complex matrices.

My code works with ComplexEigenSolver but I am sure can be made faster if it could use the selfadjointeigensolver.

Any ideas?

Kaveh
PS. LAPACK has CHEEV function which does the above. I might actually consider using it for this portion of my code!

SelfAdjointEigenSolver definitely works on complex matrices. It's covered by tests, see tests/eigensolver.cpp in Eigen2 and tests/eigensolver_selfadjoint.cpp in Eigen3.

Then perhaps, I am not coding it right (it's a template based defined function). Thanks for the reply!