precond_ifpack.h

Go to the documentation of this file.

// This file is part of HermesCommon
//
// Copyright (c) 2009 hp-FEM group at the University of Nevada, Reno (UNR).
// Email: hpfem-group@unr.edu, home page: http://www.hpfem.org/.
//
// Hermes2D is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published
// by the Free Software Foundation; either version 2 of the License,
// or (at your option) any later version.
//
// Hermes2D is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with Hermes2D; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
#ifndef __HERMES_COMMON_PRECOND_IFPACK_H_
#define __HERMES_COMMON_PRECOND_IFPACK_H_
#include "../../config.h"
#ifdef HAVE_IFPACK

#include "solvers/precond.h"
#include "solvers/interfaces/epetra.h"
#include <Ifpack_Preconditioner.h>

namespace Hermes
{
  namespace Preconditioners
  {
    using namespace Hermes::Solvers;
    template <typename Scalar>
    class HERMES_API IfpackPrecond : public EpetraPrecond < Scalar >
    {
    public:
      IfpackPrecond(const char *cls, const char *type = "Jacobi");
      IfpackPrecond(const char *cls, const char *type, int overlap);
      IfpackPrecond(Ifpack_Preconditioner *ipc);
      virtual ~IfpackPrecond();

      void set_param(const char *name, const char *value);
      void set_param(const char *name, int value);
      void set_param(const char *name, double value);
    protected:

      virtual Epetra_Operator *get_obj() { return prec; }

      virtual void create(Matrix<Scalar> *mat);
      virtual void destroy() { }
      virtual void compute();
      virtual void recompute() { destroy(); compute(); }

      // Epetra_Operator interface
      virtual int ApplyInverse(const Epetra_MultiVector &r, Epetra_MultiVector &z) const;
      virtual const Epetra_Comm &Comm() const;
      virtual const Epetra_Map &OperatorDomainMap() const;
      virtual const Epetra_Map &OperatorRangeMap() const;

      void create_point_relax(EpetraMatrix<Scalar> *a, const char *name);
      void create_block_relax(EpetraMatrix<Scalar> *a, const char *name);
      void create_add_schwartz(EpetraMatrix<Scalar> *a, const char *name, int overlap);
      int initialize();
      void apply_params();
      Ifpack_Preconditioner *prec;
      Teuchos::ParameterList ilist;
      EpetraMatrix<Scalar> *mat;
      unsigned owner : 1;
      // class of the preconditioner
      const char *cls;
      const char *type;
      int overlap;

      friend class AztecOOSolver < Scalar > ;
    };
  }
}
#endif
#endif