dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

kent-and@xxxxxxxxx wrote:
> > One or more new changesets pushed to the primary dolfin repository.
> > A short summary of the last three changesets is included below.
> >
> > changeset:   4730:6548c25c33352492c9279c035509a139caab323b
> > tag:         tip
> > user:        "Garth N. Wells <gnw20@xxxxxxxxx>"
> > date:        Tue Sep 09 13:08:58 2008 +0100
> > files:       dolfin/fem/Assembler.cpp dolfin/pde/LinearPDE.cpp
> > description:
> > Revert to old assembly in LinearPDE due to bug in symmetric assembly.
> >
> > There is a problem with exterior facets in the symmetric assembly. The
> > code needs to be broken up to make debugging easier.
>
> Agree that it should be broken up. But it will take some effort.
> Do you have an example where the bug is apparent ?

I was using a RT0 element (not one of the demos).  It should be 
reproducible by adding a non-zero Dirichlet bc (given be a Function) to 
the mixed-poisson demo.

I had a quick look, but I couldn't find the problem, so I thought it 
better to wait until the code is broken up.

Garth


> Kent
>
>
> > changeset:   4729:a84cf2b5cda6170367b30e1c6bcffa00bc4f36d1
> > parent:      4727:04152c9160962363af3a996773ced092caed339d
> > parent:      4728:a5ad46cb91f08dac998bd2a506f6c7e6aa4d0da1
> > user:        Anders Logg <logg@xxxxxxxxx>
> > date:        Mon Sep 08 00:03:12 2008 +0200
> > files:
> > description:
> > Merge patch from Shawn Walker
> >
> >
> > changeset:   4728:a5ad46cb91f08dac998bd2a506f6c7e6aa4d0da1
> > parent:      4703:0e3198c7a287f5d7e8604f0f26fa8cc958fbc5b4
> > user:        Shawn Walker <walker@xxxxxxxxxxxx>
> > date:        Wed Sep 03 10:56:23 2008 -0400
> > files:       dolfin/function/DiscreteFunction.h
> > dolfin/function/Function.cpp dolfin/function/Function.h
> > dolfin/io/XMLMesh.cpp dolfin/io/XMLMesh.h dolfin/io/XMLObject.cpp
> > dolfin/io/XMLObject.h dolfin/io/XMLVector.h dolfin/mesh/MeshEditor.cpp
> > dolfin/mesh/MeshEditor.h dolfin/mesh/MeshGeometry.cpp
> > dolfin/mesh/MeshGeometry.h
> > description:
> > This patch adds in support for reading in an additional 'coordinates' tag
> > in the mesh.xml file, and uses this information to create a Function*
> > called mesh_coordinates that is contained in the MeshGeometry class.  This
> > function stores higher order mesh coordinate data that comes from assuming
> > the mesh elements are mapped using a finite element function higher order
> > than linear.  Another boolean array is created inside MeshGeometry that
> > indicates which elements are ACTUALLY curved, and which are just straight
> > (standard).  This is completely backwards compatible with previous mesh
> > formats.
> >
> > ----------------------------------------------------------------------
> > 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