adds twod brick interp test

src/ref_interp_test.c

@@ -1340,6 +1340,40 @@ 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");
+  }
+
   { /* integrate scalar */
     char grid[] = "ref_interp_test_scalar.meshb";
     REF_GRID ref_grid;