Re: [Libmesh-users] Complex EigenVector from Real Eigensystem

2014-09-09 Thread David Knezevic
The EPS object is private, so I think the easiest thing here would be to modify SlepcEigenSolver to suit your needs. The beauty of open source software :) You could change SlepcEigenSolver so that it can access Vi, like I said in my previous email (e.g. a bool to pick Vr or Vi). Or you could ju

Re: [Libmesh-users] Complex EigenVector from Real Eigensystem

2014-09-09 Thread Harshad Sahasrabudhe
Hi David, Thanks for your reply. I'm working on a larger software which needs LibMesh configured with real numbers. Does this mean I won't be able to get complex eigen vectors? Is it possible to get the EPS object from SlepcEigenSolver so that I can directly call SLEPc? Thanks, Harshad On Tue,

Re: [Libmesh-users] Complex EigenVector from Real Eigensystem

2014-09-09 Thread David Knezevic
This is a good point. SLEPc does return both real and imaginary parts, e.g. see here: http://www.grycap.upv.es/slepc/documentation/current/docs/manualpages/EPS/EPSGetEigenpair.html Note that it says "If PETSc is configured with complex scalars the eigenvalue is stored directly in eigr (eigi is s

[Libmesh-users] Complex EigenVector from Real Eigensystem

2014-09-09 Thread Harshad Sahasrabudhe
Hi, I'm using LibMesh compiled with SLEPc and PETSc in real mode. I'm using EigenSystem class to solve an eigen problem with real matrix coefficients. The eigen vectors of such a system can be complex. How do I get the complex eigenvectors? As far as I have looked into LibMesh source code, I have