Merge branch 'yplus-length-scale' into 'main'

Y+ length scale See merge request fun3d-developers/refine!964
nasa · Oct 19, 2022 · a5bc19a · a5bc19a
2 parents 1cf54cc + 83d820a
commit a5bc19a
Show file tree

Hide file tree

Showing 4 changed files with 431 additions and 7 deletions.
diff --git a/src/ref_interp_test.c b/src/ref_interp_test.c
@@ -1340,6 +1340,72 @@ int main(int argc, char *argv[]) {
     RSS(ref_grid_free(from), "free");
   }
 
+  { /* twod tri bricks */
+    REF_GRID from, to;
+    char file[] = "ref_interp_test_tri.meshb";
+    REF_INTERP ref_interp;
+    REF_DBL max_error, min_bary;
+
+    if (ref_mpi_once(ref_mpi)) {
+      REF_INT dim = 11;
+      RSS(ref_fixture_twod_brick_grid(&from, ref_mpi, dim), "brick");
+      RSS(ref_export_by_extension(from, file), "export");
+      RSS(ref_grid_free(from), "free");
+    }
+    RSS(ref_part_by_extension(&from, ref_mpi, file), "import");
+    RSS(ref_part_by_extension(&to, ref_mpi, file), "import");
+    if (ref_mpi_once(ref_mpi)) REIS(0, remove(file), "test clean up");
+
+    RSS(ref_interp_create(&ref_interp, from, to), "make interp");
+    RSS(ref_interp_locate(ref_interp), "map");
+    REIS(4, ref_interp->n_geom, "geom missing");
+    REIS(0, ref_interp->n_geom_fail, "geom fail");
+    if (!ref_mpi_para(ref_mpi)) {
+      REIS(117, ref_interp->n_walk, "walk count");
+      REIS(0, ref_interp->n_tree, "tree count");
+    }
+    RSS(ref_interp_min_bary(ref_interp, &min_bary), "min bary");
+    RAS(-0.00001 < min_bary, "large extrapolation");
+    RSS(ref_interp_max_error(ref_interp, &max_error), "err");
+    RAS(7.0e-16 > max_error, "large interp error");
+    RSS(ref_interp_free(ref_interp), "interp free");
+
+    RSS(ref_grid_free(to), "free");
+    RSS(ref_grid_free(from), "free");
+  }
+
+{ /* twod tri brick (same) */
+    REF_GRID ref_grid;
+    char file[] = "ref_interp_test_tri_same.meshb";
+    REF_INTERP ref_interp;
+    REF_DBL max_error, min_bary;
+
+    if (ref_mpi_once(ref_mpi)) {
+      REF_INT dim = 11;
+      RSS(ref_fixture_twod_brick_grid(&ref_grid, ref_mpi, dim), "brick");
+      RSS(ref_export_by_extension(ref_grid, file), "export");
+      RSS(ref_grid_free(ref_grid), "free");
+    }
+    RSS(ref_part_by_extension(&ref_grid, ref_mpi, file), "import");
+    if (ref_mpi_once(ref_mpi)) REIS(0, remove(file), "test clean up");
+
+    RSS(ref_interp_create(&ref_interp, ref_grid, ref_grid), "make interp");
+    RSS(ref_interp_locate(ref_interp), "map");
+    REIS(4, ref_interp->n_geom, "geom missing");
+    REIS(0, ref_interp->n_geom_fail, "geom fail");
+    if (!ref_mpi_para(ref_mpi)) {
+      REIS(117, ref_interp->n_walk, "walk count");
+      REIS(0, ref_interp->n_tree, "tree count");
+    }
+    RSS(ref_interp_min_bary(ref_interp, &min_bary), "min bary");
+    RAS(-0.00001 < min_bary, "large extrapolation");
+    RSS(ref_interp_max_error(ref_interp, &max_error), "err");
+    RAS(7.0e-16 > max_error, "large interp error");
+    RSS(ref_interp_free(ref_interp), "interp free");
+
+    RSS(ref_grid_free(ref_grid), "free");
+  }
+
   { /* integrate scalar */
     char grid[] = "ref_interp_test_scalar.meshb";
     REF_GRID ref_grid;