dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Thu, Jul 06, 2006 at 07:54:02PM +0200, Garth N. Wells wrote:
> 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.

Works fine, thanks.

/Anders

> > /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
>