[Trilinos-Users] Eigenvalues and constraint interface

[Veltz Romain]

Dear all,

I am studying an equation F(x)=0 with continuous symmetries and I want to compute the eigenvalues of the jacobian. More precisely, I am a performing numerical continuation with the factory LOCA::Epetra and my interface is coupled to a constraint interface LOCA::MultiContinuation::ConstraintInterfaceMVDX in order to remove the undetermined phase associated with the continuous symmetry. 

Then I compute the eigenvalues with anasazi. My problem is I still have the zero eigenvalues associated with the undetermined phases: they should be removed by using my constraint problem. It is as if the constraint interface were not used by anasazi. Also, I put some cout in my interface  LOCA::MultiContinuation::ConstraintInterfaceMVDX in order to check if it is used by anasazi but it does not seem to be the case.

Hence, I would like to know if the anasazi methods are linked to the constrained interface LOCA::MultiContinuation::ConstraintInterfaceMVDX in LOCA. Did I miss something ?

Thank you for your help,

Romain.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20130606/e5e3231a/attachment.html