Reading mesh and functions from XDMF Timeseries

8 months ago by
Can I use the XDMFFile Class to get the mesh and functions of an earlier written time series?
Minimum not working example with dolfin 2017.2.0.dev0 on debian:

#creating and saving mesh
from fenics import *

XDMFFile = XDMFFile("temp.xdmf")
mesh = BoxMesh(Point(0, 0, 0), Point(1, 1, 1),5,5,5)

V=FunctionSpace(mesh, 'P', 1)
u = Function(V)

#write function with timestamp.
XDMFFile.write(u, float(0))

and after resetting / new console:

from fenics import *
mesh = Mesh()
XDMFFILE = XDMFFile("temp.xdmf")

Thanks for your help


The Error message:

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
*** -------------------------------------------------------------------------
*** Error:   Unable to recognise cell type.
*** Reason:  Unknown value "".
*** Where:   This error was encountered inside XDMFFile.cpp.
*** Process: 0
*** DOLFIN version: 2017.2.0.dev0
*** Git changeset:  8a9b150311981abdcb72e808bd5b1657f2384aa4
*** -------------------------------------------------------------------------
Community: FEniCS Project

1 Answer

8 months ago by
The version of XDMF in the current interface is not suitable for restoring function values. You should look at the "HDF5File" interface instead.
So I have to save it twice, once for visualization and once for later acces? (I'll use the TimeSeries class and the XDMF File now)
written 8 months ago by Baltasar  
Unfortunately, that is the case - you must save as HDF5 and XDMF. The XDMF format is really only for visualisation, at present. We have a project currently to improve the situation.
written 8 months ago by Chris Richardson  
Please login to add an answer/comment or follow this question.

Similar posts:
Search »