]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Adjust tolerance in push_forward 8191/head
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Thu, 16 May 2019 19:56:20 +0000 (21:56 +0200)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Thu, 16 May 2019 19:56:20 +0000 (21:56 +0200)
source/grid/manifold_lib.cc

index c61028bc39ccb6f00bd98f6c04f601d69403df32..75138e0dbd624b353b3b1d773cd977fb585edef1 100644 (file)
@@ -1942,9 +1942,9 @@ TransfiniteInterpolationManifold<dim, spacedim>::push_forward(
   AssertDimension(cell->level(), level_coarse);
 
   // check that the point is in the unit cell which is the current chart
-  // Tolerance 1e-6 chosen that the method also works with
-  // SphericalManifold
-  Assert(GeometryInfo<dim>::is_inside_unit_cell(chart_point, 1e-6),
+  // Tolerance 5e-4 chosen that the method also works with manifolds
+  // that have some discretization error like SphericalManifold
+  Assert(GeometryInfo<dim>::is_inside_unit_cell(chart_point, 5e-4),
          ExcMessage("chart_point is not in unit interval"));
 
   return compute_transfinite_interpolation(*cell,

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.