I see that facet_normals has been added as an argument to custom_integral::tabulate_tensor. If you actually want to use FacetNormal with custom integrals, I think the better way is to have custom_integral and custom_facet_integral, with the latter taking 'int facet' just like exterior_facet_integral, referring to the first cell. That will allow any facet related geometric quantity to be computed just as for other facet integrals.
On a related note, the num_cells metadata hack should be made to fit with development of domain representation in ufl. Martin
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
