[Trilinos-Users] How to interface LOCA with Anasazi solver

Phipps, Eric T etphipp at sandia.gov
Mon Jul 11 08:46:03 MDT 2011

Currently the Anasazi eigensolver package is only interfaced to LOCA using the Epetra matrix/vector support, not the LAPACK support (which is what ChanContinuation.C uses). In this case you should use the LAPACK eigensolver.  An example is in


The relevant lines are:

    // Create Eigensolver list using LAPACK eigensolver
    stepperList.set("Compute Eigenvalues",true);
    Teuchos::ParameterList& eigenList = stepperList.sublist("Eigensolver");
    eigenList.set("Method", "DGGEV");     // LAPACK eigensolver
    eigenList.set("NEV", 10);             // Show 10 eigenvalues
    eigenList.set("Sorting Order", "LR"); // Sort by largest real component


    // Create LAPACK Factory (necessary for LAPACK eigensolver)
    Teuchos::RCP<LOCA::LAPACK::Factory> lapackFactory =
      Teuchos::rcp(new LOCA::LAPACK::Factory);

    // Create global data object
    Teuchos::RCP<LOCA::GlobalData> globalData =
      LOCA::createGlobalData(paramList, lapackFactory);


From: narendiran kumar <narendirankumar at gmail.com<mailto:narendirankumar at gmail.com>>
Date: Mon, 11 Jul 2011 14:42:15 +0530
To: <trilinos-users at software.sandia.gov<mailto:trilinos-users at software.sandia.gov>>
Subject: [Trilinos-Users] How to interface LOCA with Anasazi solver

Dear Sir/Madam,


I have used LOCA continuation for my model and i got continuation plot also.  Thanks for that effective example program (ChanContinuation).  My question is that i got steady state solution by varying single parameter, meanwhile i have to find out eigen value for each and every parameter values.  Already i have written jacobin matrix.  I have little bit confession with how to call "Anasazi Block Eigensolver Package"<http://trilinos.sandia.gov/packages/anasazi/index.html> in that same program.   And also i have seen in ChanContinuation example line "stepperList.set<http://trilinos.sandia.gov/packages/docs/dev/packages/teuchos/doc/html/classTeuchos_1_1ParameterList.html#a1e3073f9cfc4ba92bcd6c3f73bd865e2>("Compute Eigenvalues",false);", I made it as true and i tried to make interface with eigen solver in that same program.

I am beginner in c/c++ as well as trilinos package.  Can anyone please give me an idea or tutorial for that.



K. Narendiran,
Research Scholar,
Department of Chemical Engineering,
IITBombay, India.
_______________________________________________ Trilinos-Users mailing list Trilinos-Users at software.sandia.gov<mailto:Trilinos-Users at software.sandia.gov> http://software.sandia.gov/mailman/listinfo/trilinos-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20110711/3d2e3cf5/attachment.html 

More information about the Trilinos-Users mailing list