MueLu_AztecEpetraOperator.hpp

Go to the documentation of this file.

 
#ifndef PACKAGES_MUELU_ADAPTERS_AZTECOO_MUELU_AZTECEPETRAOPERATOR_HPP_
#define PACKAGES_MUELU_ADAPTERS_AZTECOO_MUELU_AZTECEPETRAOPERATOR_HPP_
 
#include <Epetra_Operator.h>
 
#include "Xpetra_Operator.hpp"
 
#if defined(HAVE_MUELU_SERIAL) and defined(HAVE_MUELU_EPETRA)
 
namespace MueLu {

  class AztecEpetraOperator : public Epetra_Operator {
    typedef double              SC;
    typedef int                 LO;
    typedef int                 GO;
    typedef Xpetra::EpetraNode  NO;
 
    typedef Xpetra::Map<LO,GO,NO>                        Map;
    typedef Xpetra::EpetraMapT<GO,NO>                    EpetraMap;
    typedef Xpetra::Operator<SC,LO,GO,NO>                Operator;
 
  public:



    AztecEpetraOperator(const Teuchos::RCP<Operator>& Op) : xOp_(Op) { }

    virtual ~AztecEpetraOperator() { }

 
    int SetUseTranspose(bool /* UseTransposeBool */) { return -1; }




    int Apply(const Epetra_MultiVector& /* X */, Epetra_MultiVector& /* Y */) const { return -1; }


    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 { return 0; }



    const char * Label() const { return "MueLu::AztecEpetraOperator"; }

    bool UseTranspose() const { return false; }

    bool HasNormInf() const { return 0; }

    const Epetra_Comm & Comm() const;

    const Epetra_Map & OperatorDomainMap() const;

    const Epetra_Map & OperatorRangeMap() const;




    Teuchos::RCP<Operator> GetOperator() const { return xOp_; }

 
 
  private:
 
    Teuchos::RCP<Operator> xOp_;
 
  };
 
} // namespace
 
#endif
 
#endif /* PACKAGES_MUELU_ADAPTERS_AZTECOO_MUELU_AZTECEPETRAOPERATOR_HPP_ */