dolfin team mailing list archive

[dolfin@xxxxxxxxxx (DOLFIN)]

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

changeset:   2036:8db0b44a707ccf43aaea7f01ecc1aab3a0dbe5f7
tag:         tip
user:        "Garth N. Wells <g.n.wells@xxxxxxxxxx>"
date:        Thu Jul 06 13:25:34 2006 +0200
files:       configure configure.ac src/kernel/la/DenseMatrix.cpp src/kernel/la/Makefile.am src/kernel/la/Makefile.in src/kernel/la/dolfin/DenseMatrix.h src/kernel/la/dolfin/DenseVector.h src/kernel/la/dolfin/Makefile.am src/kernel/la/dolfin/Makefile.in src/kernel/la/dolfin/dolfin_la.h src/kernel/la/dolfin/uBlasDenseMatrix.h src/kernel/la/dolfin/uBlasMatrix.h src/kernel/la/dolfin/uBlasSparseMatrix.h src/kernel/la/dolfin/ublas.h src/kernel/la/uBlasDenseMatrix.cpp src/kernel/la/uBlasKrylovSolver.cpp src/kernel/la/uBlasSparseMatrix.cpp src/test/main.cpp
description:
Add templated class uBlasMatrix. uBlas matrix classes should be derived from this class.


changeset:   2035:cd58b3f34efd8831c44c358f934cd5a0c5ad22e3
user:        "Anders Logg <logg@xxxxxxxxx>"
date:        Wed Jul 05 11:49:34 2006 +0200
files:       src/bench/ode/reaction/bench.log src/bench/ode/reaction/main.cpp src/demo/ode/courtemanche/main.cpp src/demo/ode/lorenz/main.cpp src/demo/ode/stiff/TestProblem1.h src/demo/ode/stiff/TestProblem2.h src/demo/ode/stiff/TestProblem3.h src/demo/ode/stiff/TestProblem4.h src/demo/ode/stiff/TestProblem5.h src/demo/ode/stiff/TestProblem6.h src/demo/ode/stiff/TestProblem7.h src/demo/ode/stiff/TestProblem8.h src/demo/ode/stiff/TestProblem9.h src/kernel/ode/Dependencies.cpp src/kernel/ode/MonoAdaptiveJacobian.cpp src/kernel/ode/MonoAdaptiveNewtonSolver.cpp src/kernel/ode/MonoAdaptiveTimeSlab.cpp src/kernel/ode/MultiAdaptiveFixedPointSolver.cpp src/kernel/ode/MultiAdaptiveJacobian.cpp src/kernel/ode/MultiAdaptiveTimeSlab.cpp src/kernel/ode/ODE.cpp src/kernel/ode/TimeSlab.cpp src/kernel/ode/dolfin/Dependencies.h src/kernel/ode/dolfin/MonoAdaptiveJacobian.h src/kernel/ode/dolfin/MonoAdaptiveNewtonSolver.h src/kernel/ode/dolfin/MonoAdaptiveTimeSlab.h src/kernel/ode/dolfin/MultiAdaptiveTimeSlab.h src/kernel/ode/dolfin/ODE.h src/kernel/ode/dolfin/TimeSlab.h
description:
Update interface for ODE solvers to use uBlas vectors.
Only the interface has changed. The internal representation
and solvers are still PETSc vectors. The solvers have been
tested on a number of the demo problems, but not all demos
have been ported to the new interface. A lot of editing is
needed. The speed seems to remain the same.

The new interface is either

    void f(const DenseVector& u, real t, DenseVector& y)

for the mono-adaptive solvers or

    real f(const DenseVector& u, real t, uint i)

for the multi-adaptive solvers.

The next step is to use uBlas vectors and solvers internally.


changeset:   2034:a5328ce357611ce6c8df269c208fd4827844d22d
user:        "Anders Logg <logg@xxxxxxxxx>"
date:        Tue Jul 04 19:51:27 2006 +0200
files:       src/bench/ode/reaction/bench.log src/kernel/ode/MonoAdaptiveFixedPointSolver.cpp src/kernel/ode/MonoAdaptiveNewtonSolver.cpp src/kernel/ode/MonoAdaptiveTimeSlab.cpp src/kernel/ode/MultiAdaptiveTimeSlab.cpp src/kernel/ode/TimeSlab.cpp src/kernel/ode/dolfin/MonoAdaptiveTimeSlab.h src/kernel/ode/dolfin/MultiAdaptiveTimeSlab.h src/kernel/ode/dolfin/TimeSlab.h
description:
Initial work to prepare for porting of ODE solvers to uBlas.
The solvers now rely less on working with offsets on large
arrays. Data is now copied before and after the call to f().
The extra cost seems to be neglibible.


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