Fix CylindricalManifold::get_new_point

diff --git a/source/grid/manifold_lib.cc b/source/grid/manifold_lib.cc
index c08bff2f32e5c4a08e57b99d1266aa69039b42b1..eb0be274ff17419fccf391a3153ff102863458bc 100644 (file)
--- a/source/grid/manifold_lib.cc
+++ b/source/grid/manifold_lib.cc
@@ -1049,7 +1049,7 @@ CylindricalManifold<dim, spacedim>::get_new_point(
   const double lambda = middle * direction;
 
   if ((middle - direction * lambda).square() < tolerance * average_length)
-    return Point<spacedim>() + direction * lambda;
+    return point_on_axis + direction * lambda;
   else // If not, using the ChartManifold should yield valid results.
     return ChartManifold<dim, spacedim, 3>::get_new_point(surrounding_points,
                                                           weights);