I found I have to use SubMesh. I thought since the submesh demo said it was broken that it wouldn't work. Apparently, it does.
Gennadiy On 22 April 2013 19:20, Gennadiy Rishkin <[email protected]> wrote: > Hi > > How do I write out a marked subdomain to a mesh/meshfunction for a) > plotting and b) to a mesh file. for example: > > // Create mesh functions over the cell facets > MeshFunction<std::size_t> sub_domains(mesh, mesh.topology().dim()); > > // Mark all facets as sub domain 3 > sub_domains = 3; > > // Mark hole as sub domain 0 > Hole hole; > hole.mark(sub_domains, 0); > > I'd like to write the new subdomain 'hole' to a mesh, e.g., hole.xml.gz > and also be able to plot it using: plot(hole_mesh); > > Thanks, > > Gennadiy > > > >
_______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
