src/eigsolv/BlockDiagOperator.h

Go to the documentation of this file.

#ifndef BLOCKDIAGOPERATOR_H
#define BLOCKDIAGOPERATOR_H

#include "Epetra_Operator.h"
#include "Epetra_CrsMatrix.h"
#include "Epetra_MultiVector.h"


class BlockDiagOperator : public virtual Epetra_Operator {

public:

    BlockDiagOperator::BlockDiagOperator(Epetra_CrsMatrix* A);


    ~BlockDiagOperator();
  


    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* A;
    int _NumMyRows;
    int* _diag;
    double* _temp;
    bool _useTranspose;

};

#endif /* BLOCKDIAGOPERATOR_H */

---