dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

I can't get this to work, so I'll you take care of it when you have
time.

/Anders


On Fri, Jul 07, 2006 at 02:52:43PM +0200, Anders Logg wrote:
> Looks good. I'm not touching the linear algebra currently, so feel
> free to fix it. Otherwise, I'll get to it in a while.
> 
> /Anders
> 
> 
> On Fri, Jul 07, 2006 at 02:50:18PM +0200, Garth N. Wells wrote:
> > On Fri, 2006-07-07 at 14:32 +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:   2047:5d8d7a8843ca9ca8673d5e36f0a20acd7c08696a
> > > tag:         tip
> > > user:        "Anders Logg <logg@xxxxxxxxx>"
> > > date:        Fri Jul 07 14:30:58 2006 +0200
> > > files:       src/demo/ode/courtemanche/main.cpp src/kernel/la/Makefile.am src/kernel/la/Makefile.in src/kernel/la/dolfin/uBlasKrylovMatrix.h src/kernel/la/dolfin/uBlasMatrix.h src/kernel/la/uBlasKrylovMatrix.cpp src/kernel/ode/MonoAdaptiveNewtonSolver.cpp
> > > description:
> > > Implement disp() for uBlasKrylovMatrix and redirect output for
> > > uBlasSparseMatrix through the log system.
> > > 
> > 
> > This can be made a bit simpler by using the matrix iterators. For sparse
> > matrices, they iterate over the non-zero terms only. Something like
> > 
> >   typedef Mat::const_iterator1 it1;
> >   typedef Mat::const_iterator2 it2;
> > 
> >   for (it1 i1 = this->begin1(); i1 != this->end1(); ++i1) 
> >   {
> >     // code for start of row
> >     // i1.index() should give the row number 
> > 
> >     for (it2 i2 = i1.begin(); i2 != i1.end(); ++i2)
> >     {
> >   	// code for printing each term
> > 	// *i2 is the value
> >         // i2.index1() and i2.index2() are the indexes 
> >     }
> >   }
> > 
> > should work.
> > 
> > Garth
> > 
> > > 
> > > changeset:   2046:0fb144578ec060f0f2c108e89f177988e2c9d802
> > > user:        "Anders Logg <logg@xxxxxxxxx>"
> > > date:        Fri Jul 07 13:42:25 2006 +0200
> > > files:       src/kernel/ode/HomotopyODE.cpp src/kernel/ode/dolfin/HomotopyODE.h
> > > description:
> > > Make ODE solvers independent of PETSc. The ODE solvers now compile
> > > with --disable-petsc.
> > > 
> > > 
> > > changeset:   2045:61891eb4c072b4ca71b8be1facabc64f0c41a333
> > > user:        "Anders Logg <logg@xxxxxxxxx>"
> > > date:        Fri Jul 07 12:18:21 2006 +0200
> > > files:       TODO src/kernel/ode/ComplexODE.cpp src/kernel/ode/Dependencies.cpp src/kernel/ode/Homotopy.cpp src/kernel/ode/HomotopyJacobian.cpp src/kernel/ode/HomotopyODE.cpp src/kernel/ode/MultiAdaptiveTimeSlab.cpp src/kernel/ode/TimeStepper.cpp src/kernel/ode/dolfin/ComplexODE.h src/kernel/ode/dolfin/Homotopy.h src/kernel/ode/dolfin/HomotopyJacobian.h
> > > description:
> > > Fixes for homotopy solver. Some of the demos don't work since a
> > > direct solver is needed and is currently not supported by the ODE solvers.
> > > 
> > > 
> > > -------------------------------------------------------
> > > 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
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev