dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #07895
[HG DOLFIN] Some fixes in uBlasKrylovSolver
One or more new changesets pushed to the primary dolfin repository.
A short summary of the last three changesets is included below.
changeset: 4192:40698451fd5958be0f97312cfdb3d812ebb23c49
tag: tip
user: "Anders Logg <logg@xxxxxxxxx>"
date: Thu May 15 20:59:58 2008 +0200
files: dolfin/la/SingularSolver.cpp dolfin/la/uBlasKrylovSolver.cpp dolfin/la/uBlasKrylovSolver.h
description:
Some fixes in uBlasKrylovSolver
changeset: 4191:8585a81f729f5267636f9d7f768457ebe1b9aa11
user: "Anders Logg <logg@xxxxxxxxx>"
date: Thu May 15 18:02:44 2008 +0200
files: dolfin/fem/SparsityPatternBuilder.h dolfin/la/GenericSparsityPattern.h dolfin/la/SingularSolver.cpp dolfin/la/SingularSolver.h dolfin/la/SparsityPattern.cpp dolfin/la/SparsityPattern.h dolfin/la/dolfin_la.h dolfin/swig/dolfin_docstrings.i dolfin/swig/dolfin_headers.i sandbox/la/hypre/main.cpp site-packages/dolfin/__init__.py
description:
Add new class SingularSolver:
This class provides a linear solver for singular linear systems
Ax = b where A has a one-dimensional null-space (kernel). This
may happen for example when solving Poisson's equation with
pure Neumann boundary conditions.
The solver attempts to create an extended non-singular system
by adding the constraint [1, 1, 1, ...]^T x = 0.
If an optional mass matrix M is supplied, the solver attempts
to create an extended non-singular system by adding the
constraint m^T x = 0 where m is the lumped mass matrix. This
corresponds to setting the average (integral) of the finite
element function with coefficients x to zero.
The solver makes not attempt to check that the null-space is
indeed one-dimensional. It is also assumed that the system
Ax = b retains its sparsity pattern between calls to solve().
changeset: 4190:2f32e4e3748d435fbdf0c2b1fae1f13834d42fa0
user: "Anders Logg <logg@xxxxxxxxx>"
date: Thu May 15 15:11:58 2008 +0200
files: dolfin/la/PETScKrylovSolver.cpp
description:
Add PCSetFromOptions after initializing Hypre. Enables setting
command-line options to Hypre. We really need to clean up the
mess that is PETScKrylovSolver.cpp!
----------------------------------------------------------------------
For more details, visit http://www.fenics.org/hg/dolfin