Most of the time the mesh is available in your form via a Function or otherwise. In a few cases, most notably 1*dx to compute the volume of the domain, but possibly other cases as well, the mesh could previously be passed as an argument to assemble. Forms without a domain are now illegal. Here's how to do it:
assemble(1*dx(mesh)) # volume of entire mesh assemble(1*dx(0, domain=mesh, subdomain_data=cellfunction)) # volume of subdomain 0 Martin
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
