parallel dump/restore

[Stephane Popinet]

Dear all,

With the latest release it is possible to dump and restore the state of
a simulation also for parallel (MPI) runs. See examples in:

http://basilisk.fr/src/examples/bubble.c
http://basilisk.fr/src/examples/tsunami.c

enjoy,

Stephane

[jose.lopez...@gmail.com]

Hi all, 

I had problems restoring a (axi) simulation from a dump file. Using the -catch -events flags
I saw the problem was located at,

  advection_term            src/navier-stokes/centered.h:245

more precisely at bcg.h at the tangential components

#if dimension > 1
      double vn = uf.y[i,0]/fm.y[i,0] + uf.y[i,1]/fm.y[i,1];

The problem is caused because by fm.y at ghost "face?" which is not initialized. So after restoring you get the
"Floating Point Exception". Applying BC to fm as below solves the problem ,

if (!restore (file = "dump")) {
....
}
else
    boundary((scalar *){fm});

Cheers

Jose

[Robson Nazareth]

Dear all,

Is it possible to use dump/restore with mask?

I'm trying but I cannot open the restore file.

In the init event, when I use:

if (!restore (file = "snapshot"))

it goes inside the if statement.

Regards,

Robson

[Robson Nazareth]

My mistake, the dump file wasn't in the same folder than the executable file.

Now it reads the dump file but I'm getting the following error:

restore(): error: expecting 'flags'

I found this error message in output.h inside the restore function, but I have no clue what is the problem.

It's strange that I got no error about flags when I write the dump file.

Anyone has some suggestion about how to sort this out?

Regards,

Robson

[Hiroumi Tani]

Hi Robson,

To be honest, I'm not sure, but dump/restore are not compatible with solid, i.e. mask().
http://basilisk.fr/src/COMPATIBILITY

Best,
Hiroumi

2016年12月6日火曜日 1時04分46秒 UTC+9 Robson Nazareth:

[Robson Nazareth]

Hi Hiroumi,

Thanks for your reply. I saw that compatibility chart, I was not sure if mask() and solid were considered the same thing.

I tried to run without mask() and I got the same problem. So, I'm still wondering about what is this error about flags.

Regards,

Robson

[Ning Wang]

Dear everyone,

During axisymmetric calculation (AXI.H), segmentation fault (signal 11) appears when using RESTORE(). I tried to use "Boundary ((SCALAR*) {fm})" and the error has not been solved. Do I have other methods to solve this error?

Regrads,

Ning Wang

2016年10月27日木曜日 18:53:18 UTC+8 jose.lopez...@gmail.com: