Time series IO with XDMFFile.write_checkpoint and append=False

Issue #1019 resolved
Michal Habera created an issue

This PR https://bitbucket.org/fenics-project/dolfin/pull-requests/454/apply-to-dolfinx-append-option-for/diff introduced an issue.

MWE

from dolfin import *

mesh = UnitDiscMesh.create(MPI.comm_world, 20, 1, 2)

S = VectorFunctionSpace(mesh, "CG", 1)
s = interpolate(Expression(("x[0]","x[1]"), degree=1), S)

output = XDMFFile("gradient.xdmf")
for i in range(10):
    output.write_checkpoint(s, "output", float(i))

gives

*** -------------------------------------------------------------------------
*** Error:   Unable to complete call to function open_file().
*** Reason:  Assertion file_id != HDF5_FAIL failed.

All is fine when write_checkpoint is in append=True mode.

Comments (1)

  1. Log in to comment