dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

I can't create a uBlasSparseMatrix of given size M x N:

    uBlasSparseMatrix A(M, N);

Do you need to add the constructors to uBlasSparseMatrix and
uBlasDenseMatrix?

/Anders


On Thu, Jul 06, 2006 at 01:44:29PM +0200, DOLFIN wrote:
> 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
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev