This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

TimeSeriesHDF5 reader: default MPI_Comm?

+2 votes

At the beginning of my code, I would like to read in an HDF5 file previously created with Dolfin:

tst = TimeSeriesHDF5(comm, temperature_file)

Apparently, an (MPI) communicator needs to be provided to the reader. What's the default?

(MPI_Comm() isn't working.)

asked Jan 28, 2014 by nschloe FEniCS User (7,120 points)

1 Answer

+3 votes

There's recently been many changes to the MPI interface. I think what you want to do is

comm = mpi_comm_world()
tst = TimeSeriesHDF5(comm, temperature_file)
answered Jan 28, 2014 by mikael-mortensen FEniCS Expert (29,340 points)

This

from dolfin import *
comm = mpi_comm_world()
tst = TimeSeriesHDF5(comm, 'lossless_T')

gives

*** The MPI_Info_create() function was called before MPI_INIT was invoked.
*** This is disallowed by the MPI standard.
*** Your MPI job will now abort.
[fuji:31914] Abort before MPI_INIT completed successfully; not able to guarantee that all other processes were killed!

on my machine.

commented Jan 29, 2014 by nschloe FEniCS User (7,120 points)

No problems for me. Are you trying to read an old file and getting the error on the tst line? Perhaps the file is no longer compatible with your dolfin version?

commented Jan 29, 2014 by mikael-mortensen FEniCS Expert (29,340 points)

I'm still getting this error with a recent build on two different machines with data generated from the latest sources.

commented Feb 7, 2014 by nschloe FEniCS User (7,120 points)
...