On Mon, Nov 09, 2009 at 09:30:44PM +0000, Garth N. Wells wrote: > 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.
Since it requires information from two different cells that intersect along the "facet". And those cells come from two different meshes. -- Anders > 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). > > >
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list DOLFIN-dev@fenics.org http://www.fenics.org/mailman/listinfo/dolfin-dev