src/eigsolv/HierarchicalBasisPrec.h

Go to the documentation of this file.

#ifndef HIERARCHICALBASISPREC_H
#define HIERARCHICALBASISPREC_H

#include "Epetra_Operator.h"

class Epetra_LinearProblem;
class Amesos_BaseSolver;
class BlockPCGSolver;
class AztecOO;
namespace Teuchos {
    class ParameterList;
}

#include <ml_include.h>
#include <ml_epetra_preconditioner.h>

#include "Teuchos_ParameterList.hpp"

class HierarchicalBasisPrec : public virtual Epetra_Operator {

public:

    HierarchicalBasisPrec(Epetra_CrsMatrix* Block11,
                          Epetra_CrsMatrix* TMatrix,
                          Epetra_CrsMatrix* NodeMatrix,
                          Epetra_CrsMatrix* Block12,
                          Epetra_CrsMatrix* Block22,
                          const Epetra_Map& DMap,
                          const Epetra_Map& RMap,
                          Teuchos::ParameterList& params);

    ~HierarchicalBasisPrec();
  


    int SetUseTranspose(bool UseTranspose);
  


    int Apply(const Epetra_MultiVector& X, Epetra_MultiVector& Y) const;


    int ApplyInverse(const Epetra_MultiVector& X, Epetra_MultiVector& Y) const;

    /* Returns the quantity \f$ \| A \|_\infty\f$ such that
       \f[\| A \|_\infty = \max_{1\lei\lem} \sum_{j=1}^n |a_{ij}| \f].

       \warning This method must not be called unless HasNormInf() returns true.
    */ 
    double NormInf() const;
  

    const char * Label() const;

    bool UseTranspose() const;

    bool HasNormInf() const;

    const Epetra_Comm & Comm() const;

    const Epetra_Map & OperatorDomainMap() const;

    const Epetra_Map & OperatorRangeMap() const;


protected:

    Epetra_CrsMatrix* K11;
    Epetra_CrsMatrix* TMatrix11;
    Epetra_CrsMatrix* NodeMatrix11;
    Epetra_CrsMatrix* K12;
    Epetra_CrsMatrix* K22; 
    const Epetra_Map& DomainMap;
    const Epetra_Map& RangeMap;

    Amesos_BaseSolver* Abase;
    Epetra_LinearProblem* ProblemL;
    Teuchos::ParameterList* ParamList;

    AztecOO* solver22_;
    Epetra_Operator *operator22_;
    Epetra_Operator *operator11_;

    const Epetra_Comm& _Comm;

    bool _useTranspose;
    mutable bool isFactorized;

private:
    Teuchos::ParameterList& params_;

};

#endif /* HIERARCHICALBASISPREC_H */

---