dolfin team mailing list archive

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

On Nov 9 2009, Anders Logg wrote:

> Incidentally, we (myself, Andre Massgin and Mats Larson) discussed a
> generalization of assembly over interior facets last week that also
> covers the case of intersections between non-matching meshes.
>
> The central concept in assemble_interior_facets is that of a facet and
> its two neighboring cells. 

If the 'facet' is interior to a cell, why not just assemble it as part of 
the cell integral? We're doing this already for surfaces which pass through 
cells.

Garth

Our current implementation assumes that
> this facet is a facets for both cells in the pair. But a
> generalization would include assembly over triples
>
>  (facet, cell0, cell1)
>
> where facet is contained in the closure of each of the two cells, but
> it may or may not be a facet for the two cells. It can be part of a
> facet or part of the interior. If we can implement this (and generate
> the appropriate code), we would be able to cover quite a few
> interesting cases (like Nitsche on overlapping non-matching meshes).
>
> --
> Anders
>
>
> On Wed, Nov 04, 2009 at 11:14:40AM +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:   7413:4a1d877ec9a4
> > tag:         tip
> > user:        "Garth N. Wells <gnw20@xxxxxxxxx>"
> > date:        Wed Nov 04 10:07:45 2009 +0000
> > files:       dolfin/fem/Assembler.cpp dolfin/mesh/Facet.h
> > description:
> > First simple steps towards parallel assembly of interior facets.
> >
> >
> > changeset:   7412:a34b1f81d1d2
> > user:        "Johan Hake <hake@xxxxxxxxx>"
> > date:        Sat Oct 31 23:29:25 2009 +0100
> > files:       test/unit/function/python/test.py
> > description:
> > Fix unit test for Function
> >
> >
> > changeset:   7411:84931010048c
> > user:        "Johan Hake <hake@xxxxxxxxx>"
> > date:        Fri Oct 30 11:02:28 2009 +0100
> > files:       dolfin/swig/function_pre.i
> > description:
> > Fix a bug in typemap used for the Constant interface.
> >   - We can now set negative values!
> >   - Thanks to Aron for reporting
> >
> > ----------------------------------------------------------------------
> > 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