Changeset 4087cb for src/level/level_operator_fas.cpp

Timestamp:

Oct 15, 2013, 2:55:12 PM (12 years ago)

Author:

Frederik Heber <heber@…>

Children:

fb3b4c

Parents:

b4ac80

git-author:

Frederik Heber <heber@…> (08/30/13 20:45:35)

git-committer:

Frederik Heber <heber@…> (10/15/13 14:55:12)

Message:

Increased numerical tolerance of assert in LevelOperatorFAS::Restrict().

• modified VMG::Vector::IsComponentWise..(), e.g. ..LessOr..() to accept equality within interval of width tol.

File:

• src/level/level_operator_fas.cpp (modified) (3 diffs)

src/level/level_operator_fas.cpp

@@ -81 +81 @@
 
   for (iter_f=bounds_f.Begin(), iter_c=bounds_c.Begin(); iter_c!=bounds_c.End(); iter_f+=2, ++iter_c) {
-    Helper::AssertVectorsEqual(sol_c_undist.GetSpatialPos(*iter_c), sol_f.GetSpatialPos(*iter_f), 4.0 * std::numeric_limits<vmg_float>::epsilon());
+    Helper::AssertVectorsEqual(sol_c_undist.GetSpatialPos(*iter_c), sol_f.GetSpatialPos(*iter_f), 1e6 * std::numeric_limits<vmg_float>::epsilon());
     sol_c_undist(*iter_c) = op_sol.Apply(sol_f, *iter_f);
   }
@@ -88 +88 @@
 
   for (iter_f=bounds_f.Begin(), iter_c=bounds_c.Begin(); iter_c!=bounds_c.End(); iter_f+=2, ++iter_c) {
-    Helper::AssertVectorsEqual(rhs_c_undist.GetSpatialPos(*iter_c), (*res_grid).GetSpatialPos(*iter_f), 4.0 * std::numeric_limits<vmg_float>::epsilon());
+    Helper::AssertVectorsEqual(rhs_c_undist.GetSpatialPos(*iter_c), (*res_grid).GetSpatialPos(*iter_f), 1e6 * std::numeric_limits<vmg_float>::epsilon());
     rhs_c_undist(*iter_c) = op_res.Apply(*res_grid, *iter_f) + prefactor * op_pde.Apply(sol_c_undist, *iter_c);
   }
@@ -134 +134 @@
 
   for (iter_f=bounds_f.Begin(), iter_c=bounds_c.Begin(); iter_c!=bounds_c.End(); iter_f+=2, ++iter_c) {
-    Helper::AssertVectorsEqual(sol_c.GetSpatialPos(*iter_c), sol_f_undist.GetSpatialPos(*iter_f), 4.0 * std::numeric_limits<vmg_float>::epsilon());
+    Helper::AssertVectorsEqual(sol_c.GetSpatialPos(*iter_c), sol_f_undist.GetSpatialPos(*iter_f), 1e6 * std::numeric_limits<vmg_float>::epsilon());
     val = sol_c.GetVal(*iter_c) - sol_old->GetVal(*iter_c);
     sol_f_undist(*iter_f) += op.GetDiag() * val;