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.

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

Reply via email to