[Trilinos-Users] How to Save Eigen Data in LOCA

[Phipps, Eric T]

To save eigenvalues, you need to create a class derived from LOCA::SaveEigenData::AbstractStrategy

http://trilinos.sandia.gov/packages/docs/dev/packages/nox/doc/html/classLOCA_1_1SaveEigenData_1_1AbstractStrategy.html

The implementation of the save() method should do whatever you need to do in order to save the eigenvectors/eigenvalues to a file.  There are then a couple of ways to register your class with LOCA.  The simplest is to pass an RCP to an object of this class in the parameter list:

Teuchos::RCP<LOCA::SaveEigenData::AbstractStrategy> saveEigData =
  Teuchos::rcp(new MySaveEigDataStrategy);

eigenList.set("Save Eigen Data Method", "User-Defined");
eigenList.set("User-Defined Save Eigen Data Name", "My Strategy");
eigenList.set("My Strategy", saveEigData);

where MySaveEigDataStrategy is the name of your class.  Most likely additional constructor arguments will be needed when you instantiate the object.  Unfortunately there don't appear to be any examples in LOCA demonstrating this.

-Eric

From: narendiran kumar <narendirankumar at gmail.com<mailto:narendirankumar at gmail.com>>
Date: Thu, 14 Jul 2011 22:25:11 +0530
To: <trilinos-users at software.sandia.gov<mailto:trilinos-users at software.sandia.gov>>, Eric Phipps <etphipp at sandia.gov<mailto:etphipp at sandia.gov>>
Subject: How to Save Eigen Data in LOCA

Dear Sir/Madam,

Greeting,

I am using ChanContinuation example for continuation (Bifurcation analysis) technique.  Along with this, I have tried to calculate eigen values by adding following line in that program(ChanContinuation.C).

------------------------------------------------------------------------------------------------------------------------
    ofstream outFile1("Eigenvalue.dat");
    outFile1.setf(ios::scientific, ios::floatfield);
    outFile1.precision(18);
                   :

    stepperList.set("Compute Eigenvalues",true);
    Teuchos::ParameterList& eigenList = stepperList.sublist("Eigensolver");
    eigenList.set("Method", "DGGEV");     // LAPACK eigensolver
    eigenList.set("Num Eigenvalues", n);   // Show "n" nof eigenvalues
------------------------------------------------------------------------------------------------------------------------

After this It started to compute eigenvalues for all the value of bifurcation parameter.  Now the question is, How to store all this eigenvalues (for all bifurcation parameter value) like ChanContinuation.dat format.  I have search "LOCA_SaveEigenData_DefaultStrategy.H", "LOCA_SaveEigenData_DefaultStrategy.H" and "LOCA_SaveEigenData_Factory.H" header files, Is it will help to store the data.


Thanks in Advance,

Regards,
K. Narendiran,
Research Scholar,
Department of Chemical Engineering,
Indian Institute of Technology Bombay,
Mumbai, India.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20110715/20a59690/attachment.html