dolfin team mailing list archive

[DOLFIN <dolfin@xxxxxxxxxx>]

One or more new changesets pushed to the primary DOLFIN repository.
A short summary of the last three changesets is included below.

changeset:   1796:7040f666a592834cf0537031095d76def5022a9a
tag:         tip
user:        "Anders Logg <logg@xxxxxxxxx>"
date:        Wed Mar 22 17:35:00 2006 -0500
files:       ChangeLog src/demo/pde/convection-diffusion/main.cpp
description:
Update convection-diffusion demo:

 - Fix bug in boundary conditions for pressure (same as in Stokes demos)
 - Use direct solver to speed up computation


changeset:   1795:2bae12e08f59ee393180b256387e5479e2203a87
user:        "Garth N. Wells <g.n.wells@xxxxxxxxxx>"
date:        Wed Mar 22 17:03:04 2006 +0100
files:       ChangeLog
description:
Update ChangeLog.


changeset:   1794:236275c67493735a92352a3735c6e525313ae723
user:        "Garth N. Wells <g.n.wells@xxxxxxxxxx>"
date:        Wed Mar 22 16:41:19 2006 +0100
files:       src/kernel/la/KrylovSolver.cpp src/kernel/la/Matrix.cpp src/kernel/la/dolfin/Matrix.h src/kernel/nls/NewtonSolver.cpp src/kernel/nls/dolfin/NewtonSolver.h src/kernel/parameter/dolfin/DefaultParameters.h src/kernel/pde/LinearPDE.cpp
description:
Allow PETSc matrix type to be set through the Matrix constructor.

The type of PETSc matrix determines which LU solver is used when doing a direct solve. Currently implemented  matrix types are:
  Matrix::default_solver (PETSc MATSEQAIJS)
  Matrix::spooles        (PETSC MATSEQAIJSPOOLES)
  Matrix::superlu        (PETSc MATSUPERLU)
  Matrix::umfpack        (PETSC MATUMFPACK) 

By default, PDE.solve() will use Matrix::umfpack if Umfpack has been compiled with PETSc, otherwise it will use the default PETSc LU solver.


Benchmark (for the linear solve only) for the problem in src/demo/pde/elasticity
  
  default_solver: 48.4s
  spooles:        3.45s
  superlu:        5.77s
  umfpack:        1.68s


-------------------------------------------------------
For more details, visit http://www.fenics.org/hg/dolfin