← Back to team overview

dolfin team mailing list archive

[HG] Update interface for ODE solvers to use uBlas vectors.

 

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

changeset:   2035:cd58b3f34efd8831c44c358f934cd5a0c5ad22e3
tag:         tip
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.


changeset:   2033:f803317cded1e2f71fb98ddc7701070a8d87232e
user:        "Anders Logg <logg@xxxxxxxxx>"
date:        Tue Jul 04 14:59:30 2006 +0200
files:       src/bench/ode/reaction/bench.log src/demo/ode/homotopy/simple/main.cpp src/demo/ode/stiff/run-tests.sh src/kernel/ode/Homotopy.cpp
description:
Fig bug in homotopy solver: old names used for parameters


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