Changeset 252f5c for src/level/level_operator_fas.cpp

Timestamp:

May 2, 2013, 11:58:02 AM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

092b175

Parents:

da7c98

Message:

Added tolerance parameter to Helper::AssertVectorsEqual.

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));
+    Helper::AssertVectorsEqual(sol_c_undist.GetSpatialPos(*iter_c), sol_f.GetSpatialPos(*iter_f), 4.0 * 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));
+    Helper::AssertVectorsEqual(rhs_c_undist.GetSpatialPos(*iter_c), (*res_grid).GetSpatialPos(*iter_f), 4.0 * 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));
+    Helper::AssertVectorsEqual(sol_c.GetSpatialPos(*iter_c), sol_f_undist.GetSpatialPos(*iter_f), 4.0 * 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;