GaussSeidel.h

Go to the documentation of this file.

// Copyright (C) 2004 Johan Hoffman and Anders Logg.
// Licensed under the GNU GPL Version 2.

#ifndef __GAUSS_SEIDEL_H
#define __GAUSS_SEIDEL_H

#include <dolfin/constants.h>
#include <dolfin/Preconditioner.h>
#include <dolfin/LinearSolver.h>

namespace dolfin
{
  
  class Matrix;
  class Vector;

    
  class GaussSeidel : public Preconditioner, public LinearSolver
  {
  public:

    GaussSeidel(const Matrix& A, real tol, unsigned int maxiter);

    ~GaussSeidel();

    void solve(Vector& x, const Vector& b);

    static void solve(const Matrix& A, Vector& x, const Vector& b,
                      real tol, unsigned int maxiter);

  private:
    
    // Perform one Gauss-Seidel iteration
    void iteration(const Matrix& A, Vector& x, const Vector& b);
    
    // The matrix
    const Matrix& A;

    // Tolerance
    real tol;

    // Maximum number of iterations
    unsigned int maxiter;

  };

}

#endif

Documentation automatically generated with Doxygen on 10 Sep 2004.