gsTrilinosEigenProblem.h

Go to the documentation of this file.

 
#pragma once
 
// FD Trilinos
class Epetra_MultiVector;
class Epetra_Operator;
 
namespace gismo
{
 
namespace trilinos
{
 
namespace solver
{
 
 
enum AnasaziMethod { BlockDavidson = 1, LOBPCG = 2, BlockKrylovSchur = 3 };
 
class EigenProblemPrivate;
 
class GISMO_EXPORT EigenProblem
{
public:
    typedef Epetra_MultiVector MV;
    typedef Epetra_Operator OP;
 
public:
 
    explicit EigenProblem(const SparseMatrix & A,
                          const AnasaziMethod & method = BlockKrylovSchur);
    // BlockDavidson
 
    ~EigenProblem();
 
    void solve() const;
 
 
protected:
        EigenProblemPrivate * my;
};
 
 
}// namespace solver
}// namespace trilinos
}// namespace gismo