← Back to team overview

dolfin team mailing list archive

[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