Bug in the dump-restore function

    This bug occur in axisymetric case. If we try to restore the simulation, we get a floating point exception. This doesn’t occur if we don’t have the library two-phase.h. It also doesn’t occur if we don’t include the library “axi.h”

    #include "axi.h" 
    #include "navier-stokes/centered.h" 
    #include "two-phase.h"
    int main () {

    we are just creating a grid for the simulation, no more.


    We initialise the simulation, without any specific geometry. If we had dump a file, we are restoring the dump file. In the other case, we are just outputing a string in the log file.

    event init(t=0){
      if (!restore (file = "dump"))
        fprintf(stderr, "No restoration\n");

    If we uncomment the two following line, we solve the floating point exception link to the restoration of the axisymetric case.

      //  else
      //    boundary((scalar *){fm});

    We dump the simulation

    event snapshot(i=1; i+=10; i<=100) {
      dump (file = "dump");


    This bug is not appening if we want to dump and restore the simulation in main(), even with all the include library