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

Reply via email to