Attaching mesh functions to the form in legacy syntax looks like this:
a = f*ds[facetfunction](1)
or this:
dss = ds[facetfunction]
b = f*dss(1) + g*dss(2)
I would now recommend using the new syntax:
a = f*ds(1, subdomain_data=facetfunction)
dss = ds(subdomain_data=facetfunction)
b = f*dss(1) + g*dss(2)
and passing the mesh as well is a good habit:
dss = ds(domain=mesh, subdomain_data=facetfunction)
b = f*dss(1) + g*dss(2)
because this will work even if f is an expression that doesn't contain the
mesh.
Martin
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
