sandbox/bugs/axiRestore.c
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.
init_grid(128);
run();
}
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.
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
Note
This bug is not appening if we want to dump and restore the simulation in main(), even with all the include library