]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix CylindricalManifold::get_new_point
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Tue, 5 Mar 2019 17:41:11 +0000 (18:41 +0100)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Tue, 5 Mar 2019 17:44:10 +0000 (18:44 +0100)
source/grid/manifold_lib.cc

index c08bff2f32e5c4a08e57b99d1266aa69039b42b1..eb0be274ff17419fccf391a3153ff102863458bc 100644 (file)
@@ -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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.