aztecoo_solver.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://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_AZTECOO_SOLVER_H_
#define __HERMES_COMMON_AZTECOO_SOLVER_H_
#include "../config.h"
#ifdef HAVE_AZTECOO
#include "epetra.h"
#include "linear_matrix_solver.h"
#include "precond_ifpack.h"
#include <AztecOO.h>

namespace Hermes
{
  namespace Solvers
  {
    template <typename Scalar>
    class HERMES_API AztecOOSolver : public IterSolver<Scalar>
    {
    public:

      void set_solver(const char *solver);
      void set_tolerance(double tol);
      void set_max_iters(int iters);

      virtual void set_precond(const char *name);

      AztecOOSolver(EpetraMatrix<Scalar> *m, EpetraVector<Scalar> *rhs);
      virtual ~AztecOOSolver();
      virtual bool solve();
      virtual int get_matrix_size();
    protected:
      virtual int get_num_iters();
      virtual double get_residual();

      virtual void set_precond(Precond<Scalar> *pc);

      void set_option(int option, int value);

      void set_param(int param, double value);

      AztecOO aztec;    
      EpetraMatrix<Scalar> *m;
      EpetraVector<Scalar> *rhs;

      Precond<Scalar> *pc;

      template<typename T> friend LinearMatrixSolver<T>* create_linear_solver(Matrix<T>* matrix, Vector<T>* rhs);
    };
  }
}
#endif
#endif