KrylovSolver.hpp

Go to the documentation of this file.

//  This file is part of ff3d - http://www.freefem.org/ff3d
//  Copyright (C) 2001, 2002, 2003 Stéphane Del Pino

//  This program 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, or (at your option)
//  any later version.

//  This program 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 this program; if not, write to the Free Software Foundation,
//  Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  

//  $Id: KrylovSolver.hpp,v 1.3 2008/02/13 21:25:27 delpinux Exp $

#ifndef KRYLOV_SOLVER_HPP
#define KRYLOV_SOLVER_HPP

#include <ReferenceCounting.hpp>

#include <KrylovSolverOptions.hpp>
#include <GetParameter.hpp>

#include <Structured3DMeshShape.hpp>

#include <Vector.hpp>

class BaseMatrix;
class BaseVector;

class DegreeOfFreedomSet;
class Problem;

class KrylovSolver
{
private:
  const BaseMatrix& __A;        
  const BaseVector& __b;        
  const DegreeOfFreedomSet&
  __degreeOfFreedomSet;         
  KrylovSolverOptions::Type
  __type;                       
  KrylovSolverOptions::PreconditionerType
  __pType;                      
  GetParameter<KrylovSolverOptions>
  __options;                    
public:
  KrylovSolver(const BaseMatrix& A,
               const BaseVector& b,
               const DegreeOfFreedomSet& degreeOfFreedomSet)
    : __A(A),
      __b(b),
      __degreeOfFreedomSet(degreeOfFreedomSet)
  {
    __type = __options.value().type();
    __pType = __options.value().precond();
  }

  void solve(const Problem& problem, ReferenceCounting<Vector<real_t> > u);
};

#endif // KRYLOVSOLVER_HPP

---