LGRs and Faults, computing intersection of segments approach

Author: aritorto

CMakeLists_files.cmake

@@ -37,6 +37,7 @@ list(APPEND MAIN_SOURCE_FILES
   opm/grid/cpgrid/GridHelpers.cpp
   opm/grid/cpgrid/Iterators.cpp
   opm/grid/cpgrid/Indexsets.cpp
+  opm/grid/cpgrid/LgrFaultHelpers.cpp
   opm/grid/cpgrid/LgrHelpers.cpp
   opm/grid/cpgrid/LgrOutputHelpers.cpp
   opm/grid/cpgrid/NestedRefinementUtilities.cpp
@@ -112,6 +113,7 @@ list(APPEND TEST_SOURCE_FILES
   tests/cpgrid/lgr/addLgrs_in_allActiveCartesianGrid_test.cpp
   tests/cpgrid/lgr/addLgrsOnDistributedGrid_test.cpp
   tests/cpgrid/lgr/autoRefine_test.cpp
+  tests/cpgrid/lgr/compute_new_vertices_of_overlapping_intersections_test.cpp
   tests/cpgrid/lgr/global_refine_test.cpp
   tests/cpgrid/lgr/id_entity_entityrep_test.cpp
   tests/cpgrid/lgr/level_and_grid_cartesianIndexMappers_test.cpp
@@ -211,6 +213,7 @@ list(APPEND PUBLIC_HEADER_FILES
   opm/grid/cpgrid/Indexsets.hpp
   opm/grid/cpgrid/Intersection.hpp
   opm/grid/cpgrid/Iterators.hpp
+  opm/grid/cpgrid/LgrFaultHelpers.hpp
   opm/grid/cpgrid/LgrHelpers.hpp
   opm/grid/cpgrid/LgrOutputHelpers.hpp
   opm/grid/cpgrid/ElementMarkHandle.hpp

opm/grid/cpgrid/CpGrid.cpp

@@ -2004,8 +2004,15 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
         // In entry 'level cell index', we store 'leafview cell index', or -1 when the cell vanished.
         preAdapt_level_to_leaf_cells_vec[preAdaptLevel].resize(data[preAdaptLevel]->size(0), -1);
     }
+
+    std::vector<Opm::Lgr::CellRefinementBoundaryInfo> cellRefinements(currentLeafData().size(0)); // if a cell has been refined, then
+    // its cell-refinement info is stored in gridCellRefinements[ element.index() ]
+
+    std::vector<Opm::Lgr::BoundaryInfo> boundaryInfoCellRefs(currentLeafData().size(0)); // if a cell has been refined, then
+
+    bool withoutFaults = true;
     // Ignore marked aquifer cells/connections in globalRefine()/adapt(); throw in addLgrsUpdateLeafView()/autoRef().
-    Opm::Lgr::filterMarkedAquiferCellsAndConnections(*this, throwOnFailure); 
+    Opm::Lgr::filterMarkedAquiferCellsAndConnections(*this, throwOnFailure);
     Opm::Lgr::refineAndProvideMarkedRefinedRelations( *this,
                                                       /* Marked elements parameters */
                                                       markedElem_to_itsLgr,
@@ -2025,7 +2032,10 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
                                                       cell_count,
                                                       preAdapt_level_to_leaf_cells_vec,
                                                       /* Additional parameters */
-                                                      cells_per_dim_vec);
+                                                      cells_per_dim_vec,
+                                                      cellRefinements,
+                                                      boundaryInfoCellRefs,
+                                                      withoutFaults);
 
 #if HAVE_MPI
     auto global_markedElem_count = comm().sum(markedElem_count);
@@ -2095,7 +2105,10 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
                                              assignRefinedLevel,
                                              cornerInMarkedElemWithEquivRefinedCorner,
                                              faceInMarkedElemAndRefinedFaces,
-                                             cells_per_dim_vec);
+                                             cells_per_dim_vec,
+                                             cellRefinements,
+                                             boundaryInfoCellRefs,
+                                             withoutFaults);
 
     // --- Adapted corners and PreAdapt corners relations ---
     std::map<std::array<int,2>,int>           elemLgrAndElemLgrCorner_to_adaptedCorner;
@@ -2114,7 +2127,10 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
                                       cornerInMarkedElemWithEquivRefinedCorner,
                                       vanishedRefinedCorner_to_itsLastAppearance,
                                       faceInMarkedElemAndRefinedFaces,
-                                      cells_per_dim_vec);
+                                      cells_per_dim_vec,
+                                      cellRefinements,
+                                      boundaryInfoCellRefs,
+                                      withoutFaults);
 
     // FACES
     // Stablish relationships between PreAdapt faces and refined or adapted ones ---
@@ -2131,22 +2147,21 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
                                            markedElem_to_itsLgr,
                                            assignRefinedLevel,
                                            faceInMarkedElemAndRefinedFaces,
-                                           cells_per_dim_vec);
+                                           cellRefinements);
 
     // --- Adapted faces and PreAdapt faces relations ---
     std::map< std::array<int,2>, int >           elemLgrAndElemLgrFace_to_adaptedFace;
     std::unordered_map< int, std::array<int,2> > adaptedFace_to_elemLgrAndElemLgrFace;
     // Integer to count adapted faces (mixed between faces from pre-refined-leaf faces not involved in LGRs and new-born refined faces).
     int face_count = 0;
     Opm::Lgr::identifyLeafGridFaces(currentLeafData(),
-                                    preAdaptMaxLevel,
                                     elemLgrAndElemLgrFace_to_adaptedFace,
                                     adaptedFace_to_elemLgrAndElemLgrFace,
                                     face_count,
                                     markedElem_to_itsLgr,
                                     assignRefinedLevel,
                                     faceInMarkedElemAndRefinedFaces,
-                                    cells_per_dim_vec);
+                                    cellRefinements);
 
     // Set refined level grids geometries
     // --- Refined corners  ---
@@ -2187,7 +2202,9 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
                                    preAdaptMaxLevel,
                                    markedElemAndEquivRefinedCorn_to_corner,
                                    cornerInMarkedElemWithEquivRefinedCorner,
-                                   cells_per_dim_vec);
+                                   cells_per_dim_vec,
+                                   cellRefinements,
+                                   withoutFaults);
 
     // Update leaf grid geometries
     // --- Adapted corners ---
@@ -2211,7 +2228,8 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
                                     markedElemAndEquivRefinedCorn_to_corner,
                                     cornerInMarkedElemWithEquivRefinedCorner,
                                     cells_per_dim_vec,
-                                    preAdaptMaxLevel);
+                                    preAdaptMaxLevel,
+                                    cellRefinements);
     // --- Adapted cells ---
     Opm::Lgr::populateLeafGridCells(currentLeafData(),
                                     adapted_cells,
@@ -2230,7 +2248,9 @@ bool CpGrid::refineAndUpdateGrid(bool throwOnFailure,
                                     markedElemAndEquivRefinedCorn_to_corner,
                                     cornerInMarkedElemWithEquivRefinedCorner,
                                     cells_per_dim_vec,
-                                    preAdaptMaxLevel);
+                                    preAdaptMaxLevel,
+                                    cellRefinements,
+                                    withoutFaults);
 
     for (int level = 0; level < levels; ++level) {
         const int refinedLevelGridIdx = level + preAdaptMaxLevel +1;