dolfin team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>]

Quoting Anders Logg <logg@xxxxxxxxx>:

> On Mon, Nov 17, 2008 at 02:46:58PM +0100, Kristian Oelgaard wrote:
> > Quoting Anders Logg <logg@xxxxxxxxx>:
> > 
> > > On Mon, Nov 17, 2008 at 02:25:37PM +0100, Kristian Oelgaard wrote:
> > > > Quoting Anders Logg <logg@xxxxxxxxx>:
> > > > 
> > > > > On Fri, Nov 14, 2008 at 09:24:37PM +0100, 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:   5158:8f5618bc84be42a36517969786370689478b2fe2
> > > > > > tag:         tip
> > > > > > user:        Anders Logg <logg@xxxxxxxxx>
> > > > > > date:        Fri Nov 14 21:24:28 2008 +0100
> > > > > > files:       dolfin/mesh/Mesh.cpp sandbox/misc/main.cpp
> > > > > test/unit/function/cpp/test.cpp test/unit/mesh/cpp/test.cpp
> > > > > > description:
> > > > > > Fix bug in init: only init when needed. Unit tests work again now.
> > > > > 
> > > > > I think we're back on track again. The mesh-ordering/init strategy
> is
> > > > > worked out (at least I think it looks good), and there is no call to
> > > > > mesh.order() in any of the demos.
> > > > 
> > > > The function
> > > > 
> > > > void PointCell::order(Cell& cell) const
> > > > {
> > > >   error("PointCell::order() not defined.");
> > > > }
> > > > 
> > > > in dolfin/mesh/PointCell.cpp kills the assembly over facet integrals.
> > > > Should it throw an error or should it simply do nothing?
> > > 
> > > I think it should do nothing.
> > 
> > OK, I think so too.
> > 
> > > But why is this not detected when running the demos?
> > 
> > Looks like the poisson1D demo has changed to only include *dx, I'll take a
> look
> > and fix it.
> 
> Good.

Should be fixed now. BTW, I ran the test/memory/test.py and there seem to be
leaks a plenty...

Kristian

> -- 
> Anders
>