creating _facets.xdmf mesh file

[teresak]

When I create an .xdmf file from an .xml and try to read it back in, for the _facets.xdmf I get the following error.

*** -------------------------------------------------------------------------
*** Error:   Unable to read mesh from XDMF/H5 files.
*** Reason:  XML parsing error. XDMF file should contain only one mesh/dataset.
*** Where:   This error was encountered inside XDMFxml.cpp.
*** Process: 0
*** 
*** DOLFIN version: 1.6.0dev
*** Git changeset:  5ec6384529c117e06d2a65900e78759ada9761ab
*** -------------------------------------------------------------------------

Example test code (error created at mesh_facets2):

mesh=Mesh("Sinusoid.xml")
mesh_facets = MeshFunction("size_t", mesh, "Sinusoid_facet_region.xml")
File("Sinusoid.xdmf") << mesh
File("Sinusoid_facet_region.xdmf") << mesh_facets
mesh2 = Mesh("Sinusoid.xdmf")
mesh_facets2 = MeshFunction("size_t", mesh2, "Sinusoid_facet_region.xdmf")

[Chris Richardson]

Can you ask this on the Q&A forum (fenicsproject.org/qa)? I think it is more of a question about usage, although maybe there could be a discussion about improving the user interface.

[Chris Richardson]

Actually, this does look like a bug in XDMFFile.

[Chris Richardson]

I think this is fixed in the latest master. Can you check?

[Chris Richardson]

Works OK with latest dev version

[teresak]

Yes, works perfectly now. Thank you!

[Anders Logg (Chalmers)]

Removing milestone: 1.7 (automated comment)