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_ */