src/test/radial.c
Check refinement of radial metric.
#include "radial.h"
#include "run.h"
int main() {
init_grid (4);
run();
}
event init (i = 0) {
refine (sq(x - 0.5) + sq(y - 0.5) < sq(0.3) && level < 3);
foreach() {
// fprintf (stderr, "%g %g %g %g\n", x, y, cm[], r*dtheta/L0);
assert (fabs(cm[] - r*dtheta/L0) < 1e-12);
}
foreach_face (x) {
// fprintf (stderr, "%g %g %g %g\n", x, y, fm.x[], r*dtheta/L0);
assert (fabs(fm.x[] - r*dtheta/L0) < 1e-12);
}
}