On Wed, Feb 4, 2015 at 9:03 AM, Johannes Ring <[email protected]> wrote: > On Wed, Feb 4, 2015 at 10:01 AM, Mikael Mortensen > <[email protected]> wrote: >> >>> 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. > > Ok, great. I will try a newer version of OpenMPI. >
If you have the choice, I recommend MPICH over OpenMPI. Most MPI problems will then go away . . . . Garth > Johannes > _______________________________________________ > 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
