remapycon: added call to yac_get_cell_bounding_circle_reg_quad().

src/remap_conserv.cc

@@ -756,7 +756,13 @@ remapConservWeights(RemapSearch &rsearch, RemapVars &rv)
       if ( useCellsearch )
         {
           struct bounding_circle bnd_circle;
-          yac_get_cell_bounding_circle(tgt_grid_cell[ompthID], &bnd_circle);
+          double *xyz = tgt_grid_cell[ompthID].coordinates_xyz;
+          if (tgt_remap_grid_type == REMAP_GRID_TYPE_REG2D)
+            yac_get_cell_bounding_circle_reg_quad(xyz, xyz+3, xyz+6, xyz+9, &bnd_circle);
+          else if (tgt_num_cell_corners==3)
+            yac_get_cell_bounding_circle_unstruct_triangle(xyz, xyz+3, xyz+6, &bnd_circle);
+          else
+            yac_get_cell_bounding_circle(tgt_grid_cell[ompthID], &bnd_circle);
 
           struct dep_list result_list;
           yac_do_bnd_circle_search((struct grid_search *)rsearch.yacSearch, &bnd_circle, 1, &result_list);