From ef2709b309e64c7bde65e472075d73af0706ecb9 Mon Sep 17 00:00:00 2001 From: Rene Gassmoeller Date: Fri, 8 Dec 2017 10:17:52 -0700 Subject: [PATCH] Fix 2D case. --- source/grid/manifold_lib.cc | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/source/grid/manifold_lib.cc b/source/grid/manifold_lib.cc index ff3b239193..3c79f1c550 100644 --- a/source/grid/manifold_lib.cc +++ b/source/grid/manifold_lib.cc @@ -474,14 +474,14 @@ get_new_points (const ArrayView> &surrounding_points, // If not in 3D, just use the implementation from PolarManifold // after we verified that the candidate is not the center. if (spacedim<3) - { - new_points[row] = polar_manifold.get_new_point(surrounding_points, ArrayView(&weights[row*weight_columns],weight_columns)); - accurate_point_was_found[row] = true; - continue; - } - + new_points[row] = polar_manifold.get_new_point(surrounding_points, ArrayView(&weights[row*weight_columns],weight_columns)); } + // In this case, we treated the case that the candidate is the center and obtained + // the new locations from the PolarManifold object otherwise. + if (spacedim<3) + return; + // If all the points are close to each other, we expect the estimate to // be good enough. This tolerance was chosen such that the first iteration // for a at least three time refined HyperShell mesh with radii .5 and 1. -- 2.39.5