Add explict check cell-to-face and face-to-cell in simple grid, with 1x2x1 children per parent

Author: aritorto

opm/grid/cpgrid/LgrFaultHelpers.hpp

@@ -104,9 +104,10 @@ struct FieldVectorLess {
     bool operator()(const Dune::FieldVector<double,3>& v,
                     const Dune::FieldVector<double,3>& w) const
     {
+        double tolerance = 1e-8;
         for (int i = 0; i < 3; ++i) {
-            if (v[i] < w[i]) return true; 
-            if (v[i] > w[i]) return false; 
+            if (v[i] < w[i] - tolerance) return true; 
+            if (v[i] > w[i] + tolerance) return false; 
         }
         return false;
     }

opm/grid/cpgrid/LgrHelpers.cpp

@@ -1393,7 +1393,7 @@ void populateRefinedCells(const Dune::cpgrid::CpGridData& current_data,
             if (!Opm::Lgr::areClose(w,v)) {
                 continue;
             }
-            else { std::cout<< v[0]<< " " << v[1] << " " << v[2] <<std::endl;
+            else {
                 return lastCellAndCorner;
             }
         }
@@ -1553,7 +1553,7 @@ void populateLeafGridFaces(const Dune::cpgrid::CpGridData& current_data,
             if (!Opm::Lgr::areClose(w,v)) {
                 continue;
             }
-            else { std::cout<< v[0]<< " " << v[1] << " " << v[2] <<std::endl;
+            else { 
                 return lastCellAndCorner;
             }
         }
@@ -1667,7 +1667,7 @@ void populateLeafGridCells(const Dune::cpgrid::CpGridData& current_data,
             if (!Opm::Lgr::areClose(w,v)) {
                 continue;
             }
-            else { std::cout<< v[0]<< " " << v[1] << " " << v[2] <<std::endl;
+            else { 
                 return lastCellAndCorner;
             }
         }