solve.h

Dolfin provide a high-performance linear algebra library

// Copyright (C) 2007 Anders Logg.
// Licensed under the GNU LGPL Version 2.1.
//
// First added:  2007-04-30
// Last changed: 2007-04-30

#ifndef __SOLVE_H
#define __SOLVE_H

#include <dolfin/constants.h>

namespace dolfin
{  

  class Matrix;
  class Vector;

  /// Solve linear system Ax = b using a direct method (LU factorization).
  /// Note that iterative methods (preconditioned Krylov methods including
  /// GMRES) are also available through the KrylovSolver interface.

  /// Solve linear system Ax = b
  void solve(const Matrix& A, Vector& x, const Vector& b);
  
  /// Solve linear system Ax = b
  //void solve(const PETScKrylovMatrix& A, PETScVector& x, const PETScVector& b);
  
}

#endif