On Wed, Feb 4, 2015 at 10:21 AM, Garth N. Wells <[email protected]> wrote: > 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 . . . .
Yes, it is easy to switch to MPICH with hashdist. I will try that and if it works for me on both Linux and OS X, we might consider switching to MPICH in the fenics-install.sh script. Johannes _______________________________________________ fenics-support mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics-support
