Enable DUNE grid tests

For this to work it is additionally needed that
* Entities return the correct number of faces (to count intersections in the tests)
* The refinement does not throw. Instead, it returns false when it cannot make a refinement.

Author: SoilRos

opm/grid/cpgrid/CpGridData.cpp

@@ -1577,8 +1577,8 @@ void CpGridData::distributeGlobalGrid(CpGrid& grid,
 
     // Compute partition type for points
     computePointPartitionType();
-   
-    computeCommunicationInterfaces(noExistingPoints);   
+
+    computeCommunicationInterfaces(noExistingPoints);
 #else // #if HAVE_MPI
     static_cast<void>(grid);
     static_cast<void>(view_data);
@@ -2065,9 +2065,9 @@ int CpGridData::sharedFaceTag(const std::vector<std::array<int,3>>& startIJK_2Pa
     assert(endIJK_2Patches.size() == 2);
 
     int faceTag = -1; // 0 represents I_FACE, 1 J_FACE, and 2 K_FACE. Use -1 for no sharing face case.
-     
+
     if (patchesShareFace(startIJK_2Patches, endIJK_2Patches)) {
-        
+
         const auto& detectSharing = [](const std::vector<int>& faceIdxs, const std::vector<int>& otherFaceIdxs){
             bool faceIsShared = false;
             for (const auto& face : faceIdxs) {
@@ -2080,7 +2080,7 @@ int CpGridData::sharedFaceTag(const std::vector<std::array<int,3>>& startIJK_2Pa
             }
             return faceIsShared; // should be false here
         };
-     
+
         const auto& [iFalse, iTrue, jFalse, jTrue, kFalse, kTrue] = this->getBoundaryPatchFaces(startIJK_2Patches[0], endIJK_2Patches[0]);
         const auto& [iFalseOther, iTrueOther, jFalseOther, jTrueOther, kFalseOther, kTrueOther] =
             this->getBoundaryPatchFaces(startIJK_2Patches[1], endIJK_2Patches[1]);
@@ -2668,11 +2668,11 @@ CpGridData::refinePatch(const std::array<int,3>& cells_per_dim, const std::array
 bool CpGridData::mark(int refCount, const cpgrid::Entity<0>& element)
 {
     if (refCount == -1) {
-        OPM_THROW(std::logic_error, "Coarsening is not supported yet.");
+        return false; // Coarsening is not supported yet.
     }
     // Prevent refinement if the cell has a non-neighbor connection (NNC).
     if (hasNNCs({element.index()}) && (refCount == 1)) {
-        OPM_THROW(std::logic_error, "Refinement of cells with face representing an NNC is not supported yet.");
+        return false; // Refinement of cells with face representing an NNC is not supported yet
     }
     assert((refCount == 0) || (refCount == 1)); // Do nothing (0), Refine (1), Coarsen (-1) not supported yet.
     if (mark_.empty()) {
@@ -2699,7 +2699,7 @@ bool CpGridData::preAdapt()
             if (local_empty)
                 mark_.resize(size(0));
         }
-       
+
         // Detect the maximum mark across processes, and rewrite
         // the local entry in mark_, i.e.,
         // mark_[ element.index() ] = max{ local marks in processes where this element belongs to}.

opm/grid/cpgrid/Entity.hpp

@@ -382,6 +382,12 @@ unsigned int Entity<codim>::subEntities ( const unsigned int cc ) const
         if ( cc == 3 ) { // Get number of corners of the element.
             return 8;
         }
+        if ( cc == 2 ) { // Get number of lines of the element.
+            return 12;
+        }
+        if ( cc == 1 ) { // Get number of faces of the element.
+            return 6;
+        }
     }
     return 0;
 }

tests/test_cpgrid.cpp

@@ -91,7 +91,6 @@ template <class Grid>
 void testGrid(Grid& grid, const std::string& name, const size_t nElem, const size_t nVertices)
 {
     typedef typename Grid::LeafGridView GridView;
-    /*
 
     try {
       gridcheck( grid );
@@ -100,7 +99,7 @@ void testGrid(Grid& grid, const std::string& name, const size_t nElem, const siz
     {
       std::cerr << "Warning: " << e.what() << std::endl;
     }
-*/
+
     std::cout << name << std::endl;
 
     testGridIteration( grid.leafGridView(), nElem );
@@ -166,7 +165,7 @@ int main(int argc, char** argv )
 
     grid.processEclipseFormat(&ecl_grid, nullptr, false, false, false);
     testGrid( grid, "CpGrid_ecl", 8, 27 );
-    
+
     const auto& grid_leafView = grid.leafGridView();
     Dune::CartesianIndexMapper<Grid> grid_cartMapper =  Dune::CartesianIndexMapper<Grid>(grid);
     for (const auto& element: Dune::elements(grid_leafView)){
@@ -183,8 +182,8 @@ int main(int argc, char** argv )
         }
         std::cout << " " << '\n';
     }
-    
-    
+
+
 #endif
 
     std::stringstream dgfFile;