Here’s the thread I mentioned http://www.open-mpi.org/community/lists/users/2014/05/24288.php <http://www.open-mpi.org/community/lists/users/2014/05/24288.php>
Problem boils down to MPI file open, this should crash as well: import mpi4py.MPI as MPI comm = MPI.COMM_WORLD MPI.File.Open(comm, filename) M > On 04 Feb 2015, at 18:03, Johannes Ring <[email protected]> wrote: > > On Wed, Feb 4, 2015 at 10:01 AM, Mikael Mortensen > <[email protected] <mailto:[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. > > Johannes
_______________________________________________ fenics-support mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics-support
