> On 04 Feb 2015, at 17:53, Johannes Ring <[email protected]> wrote:
> 
> Hi,
> 
> Running the following script works fine in serial but fails to run in 
> parallel:
> 
> from dolfin import *
> 
> filename = "mesh3d.h5"
> 
> # Write to file
> mesh0 = UnitCubeMesh(10, 10, 10)
> mesh_file = HDF5File(mesh0.mpi_comm(), filename, "w")
> mesh_file.write(mesh0, "/my_mesh")
> mesh_file.close()
> 
> # Read from file
> mesh1 = Mesh()
> mesh_file = HDF5File(mesh0.mpi_comm(), filename, "r")
> mesh_file.read(mesh1, "/my_mesh", False)
> mesh_file.close()
> 
> The error message is attached. Any idea what the problem could be?

I think the MPI version is too old. I had the same problem with h5py some time 
ago and I made it go away with a more recent MPI version. I found some advice 
on the net saying that the threading support in the openMPI 1.6 series what not 
very good and the advice was to upgrade. Cannot remember where I found the 
advice, but I followed it and it worked.

Mikael

> 
> Johannes
> <error.log>_______________________________________________
> fenics-support mailing list
> [email protected]
> http://fenicsproject.org/mailman/listinfo/fenics-support

_______________________________________________
fenics-support mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics-support

Reply via email to