Snapshot when using MPI

[Guangnian Ji]

Dear all,

I am experiencing a strange issue with dump() while running a 2D atomisation simulation using MPI.

I used the code in atomisation.c to save the snapshot (dump file):

event snapshot (t = 0.1; t += 0.1; t <= 3.8) {
  char name[80];
  sprintf (name, "snapshot-%g", t);
  scalar pid[];
  foreach()
    pid[] = fmod(pid()*(npe() + 37), npe());
  dump (name);
}

When I run my simulation using 1 core, the snapshot files are created correctly and I can visualize them using bview2D. However, when I run the same code using MPI appear to contain only zeros when opened in bview2D.

Has anyone encountered this discrepancy between serial and MPI dump outputs? Could this be related to how Rank 0 gathers data from other ranks before writing the header?

Any suggestions or insights would be greatly appreciated.

Best regards,

Guangnian

[alienor...@gmail.com]

Dear Guangnian,

Can you restart your simulation using the dump files generated with MPI?

If so, I am experiencing the same problem but unfortunately could not find a solution either:

I can visualize the simulation on the fly when compiled in parallel. I can also save dumps and restart from them. Yet, when using dump2D to open a previously generated dump, it's seems like the dump is empty or contains only zeros (which is obviously not true since I can restart simulations from it).

The statistics show the correct fields' names but with all the values equal to 0.

I experience the same problem with dump3D.

I tried reinstalling everything with no success. If anyone has a clue I'll be happy to hear about it as well.

Happy New Year to all of you!

Best,

Aliénor