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 <[email protected]>" >> 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 <[email protected]>" >> 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 <[email protected]>" >> 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 >> [email protected] >> http://www.fenics.org/mailman/listinfo/dolfin-dev > _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
