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