[Trilinos-Users] [EXTERNAL] Trilinos question
Thornquist, Heidi K
hkthorn at sandia.gov
Mon Jun 25 04:11:05 MDT 2012
There are several examples in Anasazi that show how to solve the generalized eigenvalue problem depending on whether your problem is symmetric (positive definite) or non-symmetric when K and M are Epetra_CrsMatrix objects. You can find them in the anasazi/epetra/example/<solver_type> directories. Any file with "Gen" in its name is an example that solves a generalized eigenvalue. Regarding an example that constructs an eigenvalue problem with Epetra_FECrsMatrix objects, we do not have that. You can look in at Epetra to find examples on how to construct matrices of this type. However, since an Epetra_FECrsMatrix is also an Epetra_Object, any of the aforementioned examples will work for those objects.
Anasazi does not provide an interface to ARPACK. The files in the anasazi/test/ARPACKExamples directory were used to perform a comparison of Anasazi's block Krylov-Schur algorithm with ARPACK, which was published in the ACM TOMS paper:
Hope this helps.
From: trilinos-users-bounces at software.sandia.gov [trilinos-users-bounces at software.sandia.gov] on behalf of Cesar Fernandez Casanova [ceferca at cam.upv.es]
Sent: Friday, June 22, 2012 8:09 AM
To: trilinos-users at software.sandia.gov
Subject: [EXTERNAL] [Trilinos-Users] Trilinos question
I am trying to compute the typical eigenvalue problem det( K - s * M ) = 0, with FECrsMatrix K and M.
where K is the stiffness matrix and M the mass matrix.
I am wondering if you could give an example code to make this computation of this exact problem.
I have been looking the looking the Anasazi examples, but I din't find any similar.
The Housholder algorithm would also ok.
Can you explain me how to work with Arpack? I have copied the Arpack_Operators.hpp & MyMultivector.hpp in the include folder but I get some errors.
Is there any procedure?
Trilinos-Users mailing list
Trilinos-Users at software.sandia.gov
More information about the Trilinos-Users