src/directsolveroperator.h

Go to the documentation of this file.

//
// C++ Interface: directsolveroperator
//
// Description: Epetra_Operator wrapper for direct solvers (e.g. Amesos solvers)
//
//
// Author: Roman Geus <roman.geus@psi.ch>, (C) 2004
//
// Copyright: See COPYING file that comes with this distribution
//
//

#ifndef DIRECTSOLVEROPERATOR_H
#define DIRECTSOLVEROPERATOR_H

#include <Epetra_MultiVector.h>
#include <Epetra_LinearProblem.h>
#include <Teuchos_ParameterList.hpp>
#include <Amesos.h>

class Epetra_MultiVector;
class Epetra_Map;
class Epetra_Comm;

class DirectSolverOperator : public Epetra_Operator {
    public:
        DirectSolverOperator(Epetra_RowMatrix& A, Teuchos::ParameterList& params);

        ~DirectSolverOperator();

        int SetUseTranspose(bool use_transpose);

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

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

        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;

    private:
        const Epetra_Map& domain_map_;
        const Epetra_Map& range_map_;
        const Epetra_Comm& comm_;
        bool use_transpose_;
        Epetra_LinearProblem* linear_problem_;
        Amesos_BaseSolver* amesos_solver_;
};

#endif

---