Loading mesh twice fails

If I load the same mesh from a file twice, the second mesh is not correct anymore (cell connectivity appears to be random). I tested this on Ubuntu 14.04 with the Fenics 1.5 release. See also this question. Minimal example:

from dolfin import *
filename = "mesh.xml"
File(filename) << UnitSquareMesh(2,2)

for i in range(2):
    mesh = Mesh(filename)
    plot(mesh)

interactive()

Comments (4)