dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

Is the constructor correct in uBlasMatrix? It just looks like

 template <class Mat> 
 uBlasMatrix<Mat>::uBlasMatrix(const uint M, const uint N) : assembled(true)
 { 
   // Do nothing 
 }

The size does not seem to be set correctly. Should there be a call to init()?

/Anders

On Thu, Jul 06, 2006 at 05:13:35PM +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:   2037:c7974f4aff20c967f84e6fd6e7d6812e3a746f96
> tag:         tip
> user:        "Garth N. Wells <g.n.wells@xxxxxxxxxx>"
> date:        Thu Jul 06 17:12:19 2006 +0200
> files:       src/kernel/la/Makefile.am src/kernel/la/Makefile.in src/kernel/la/dolfin/DenseMatrix.h src/kernel/la/dolfin/uBlasDenseMatrix.h src/kernel/la/dolfin/uBlasMatrix.h src/kernel/la/dolfin/uBlasSparseMatrix.h src/kernel/la/uBlasDenseMatrix.cpp src/kernel/la/uBlasSparseMatrix.cpp src/test/main.cpp
> description:
> Add constructors to uBlasDenseMatrix and uBlasSparseMatrix.
> 
> I've removed the compilation flag -DNDEBUG when using --enable-debug, so uBlas operations will be extremely slow when using --enable-debug.
> 
> 
> changeset:   2036:8db0b44a707ccf43aaea7f01ecc1aab3a0dbe5f7
> 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.
> 
> 
> -------------------------------------------------------
> 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