dolfin team mailing list archive

["Garth N. Wells" <g.n.wells@xxxxxxxxxx>]

On Thu, 2006-07-06 at 19:26 +0200, Anders Logg wrote:
> 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()?
> 

Yes. Should be fixed now.

> /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/Mu
 lt
>  iAdaptiveTimeSlab.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
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev