[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
Trilinos/packages/nox/examples/lapack/LOCA_Brusselator/BrusselatorContinuation.C
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
and
// 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);
-Eric
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,
Greeting,
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.
Thanks,
Regards,
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