← Back to team overview

dolfin team mailing list archive

Re: [HG DOLFIN] Revert to old assembly in LinearPDE due to bug in symmetric assembly.

 

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

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
>




Follow ups

References