src/eigsolv/Projector.h

Go to the documentation of this file.

#ifndef PROJECTOR_H
#define PROJECTOR_H

#include "Epetra_Operator.h"


class Projector : public virtual Epetra_Operator {

public:

    Projector(Epetra_Operator* Y, const Epetra_Operator* H, Epetra_Operator* C);


    ~Projector();
  


    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_Operator* _Y;
    const Epetra_Operator* _H;
    Epetra_Operator* _C;

    bool _useTranspose;

    Epetra_Vector* _vector1;
    Epetra_Vector* _vector2;
    Epetra_Vector* _vector3;
    const Epetra_Comm& _Comm;

private:

};

#endif /* PROJECTOR_H */

---