sandbox/Antoonvh/MPI_and_embed.c

Another assertion by refine_embed_linear

Apart from this example, an assertion can be triggered when running this script in parallel with 4 threads…

It does not trigger when

• Running with a quadtree
• Running with 1,2,3 or 5 threads
• doing init_grid(N) and refine() in the main() function (no call to run())
• When using adapt_wavelet() iteratively to refine.

It seems the refine() function allows non-proper trees?

#include "grid/octree.h" 
#include "embed.h"
#include "run.h"

scalar s[];
int maxlevel = 7;

double M_top = 5;
double base_height = 10;
#define MOUNTAIN (base_height + M_top*exp(-((sq(x + L0/4.)/sq(10.)) + sq(z)/sq(30.))) - y)

int main() {
  s.refine = refine_embed_linear;
  L0 = 100;
  X0 = Z0 = -L0/2.;
  run();
}

event init (t = 0) {
  refine (fabs(MOUNTAIN) < 2.5 && level < maxlevel);
}