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

[narendiran kumar]

Dear Eric,

It works.


Thanks for your reply.  While compiling It's showing only 4 largest real
componet out off 200 eigen values.  Is there any possibility to store all
the eigen values into data file like continuation data file
(ChanContinuation.dat).  If it is possible then can you please give an idea.

-Narendiran

On Mon, Jul 11, 2011 at 8:16 PM, Phipps, Eric T <etphipp at sandia.gov> wrote:

>  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>
> Date: Mon, 11 Jul 2011 14:42:15 +0530
> To: <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
> 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/2eeda762/attachment.html