dolfin team mailing list archive

[Jake Ostien <tostien@xxxxxxxxx>]

Assembly for interior facets was broken if you were using subdomains.  
Attached is a patch, but its really just a small change.


diff -r 9672255be1b5 dolfin/fem/Assembler.cpp
--- a/dolfin/fem/Assembler.cpp	Wed May 28 10:11:25 2008 +0100
+++ b/dolfin/fem/Assembler.cpp	Wed May 28 13:16:14 2008 -0400
@@ -332,7 +332,7 @@ void Assembler::assembleInteriorFacets(G
     }
 
     // Tabulate exterior interior facet tensor on macro element
-    ufc.interior_facet_integrals[0]->tabulate_tensor(ufc.macro_A, ufc.macro_w, ufc.cell0, ufc.cell1, facet0, facet1);
+    integral->tabulate_tensor(ufc.macro_A, ufc.macro_w, ufc.cell0, ufc.cell1, facet0, facet1);
 
     // Add entries to global tensor
     A.add(ufc.macro_A, ufc.macro_local_dimensions, ufc.macro_dofs);