fix: balance compilation time and memory usage for MFD kernels

src/cmake/GeosxMacros.cmake

@@ -262,8 +262,8 @@ function(generateKernels)
   json_list_of_pairs_to_two_lists( "${constantsJSON}" constantSymbols constantValues )
 
   set( combinatoricSymbolList "" )
-  string(LENGTH ${combinationsJSON} comboStrLen)
-  if(comboStrLen GREATER 0)
+  string(JSON comboLength LENGTH ${combinationsJSON})
+  if(comboLength GREATER 0)
     set(symbolListLists "")
     foreach(symbol IN LISTS symbolList)
       string(JSON symbolArray${symbol} GET ${combinationsJSON} ${symbol})

src/coreComponents/physicsSolvers/fluidFlow/CMakeLists.txt

@@ -162,7 +162,6 @@ set( fluidFlowSolvers_sources
      kernels/compositional/AquiferBCKernel.cpp
      kernels/compositional/CFLKernel.cpp
      kernels/compositional/FluxComputeKernelBase.cpp
-     kernels/compositional/CompositionalMultiphaseHybridFVMKernelsInstantiations.cpp     
      wells/CompositionalMultiphaseWell.cpp
      wells/kernels/CompositionalMultiphaseWellKernels.cpp
      wells/SinglePhaseWell.cpp
@@ -178,7 +177,11 @@ set( kernelTemplateFileList "" )
 
 # Keep only the templates that are unrelated to hybrid flux/dirichlet
 list( APPEND kernelTemplateFileList
-      ReactiveCompositionalMultiphaseOBLKernels.cpp.template )
+      ReactiveCompositionalMultiphaseOBLKernels.cpp.template
+      CompositionalMultiphaseHybridFVMKernels_flux_grouped.cpp.template
+      CompositionalMultiphaseHybridFVMKernels_dirichlet_grouped.cpp.template
+      CompositionalMultiphaseHybridFVMKernels_bcprops.cpp.template
+      SinglePhaseHybridFVMKernels_grouped.cpp.template )
 
 foreach( kernelTemplateFile ${kernelTemplateFileList} )
   get_filename_component( jsonKey ${kernelTemplateFile} NAME_WE )
@@ -222,9 +225,6 @@ blt_add_library( NAME       fluidFlowSolvers
                  SHARED     ${GEOS_BUILD_SHARED_LIBS}
                )
 
-# Always enable the manual instantiation compile flag
-target_compile_definitions( fluidFlowSolvers PUBLIC GEOS_ENABLE_MANUAL_HYBRID_FVM_INST )
-
 target_include_directories( fluidFlowSolvers PUBLIC ${CMAKE_SOURCE_DIR}/coreComponents )
 
 install( TARGETS fluidFlowSolvers LIBRARY DESTINATION ${CMAKE_INSTALL_PREFIX}/lib )

src/coreComponents/physicsSolvers/fluidFlow/CompositionalMultiphaseHybridFVMKernels_bcprops.cpp.template

@@ -0,0 +1,205 @@
+/*
+ * ------------------------------------------------------------------------------------------------------------
+ * SPDX-License-Identifier: LGPL-2.1-only
+ *
+ * Copyright (c) 2016-2024 Lawrence Livermore National Security LLC
+ * Copyright (c) 2018-2024 TotalEnergies
+ * Copyright (c) 2018-2024 The Board of Trustees of the Leland Stanford Junior University
+ * Copyright (c) 2023-2024 Chevron
+ * Copyright (c) 2019-     GEOS/GEOSX Contributors
+ * All rights reserved
+ *
+ * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
+ * ------------------------------------------------------------------------------------------------------------
+ */
+
+/**
+ * @file CompositionalMultiphaseHybridFVMKernels_bcprops.cpp
+ *
+ * Auto-generated by generateKernels (kernelSpecs.json key: CompositionalMultiphaseHybridFVMKernels_bcprops).
+ * Do NOT edit manually.
+ *
+ * evaluateBCFaceProperties explicit instantiations for all NC in {1..5}, NP in {2,3}.
+ * NF-independent: one single translation unit.
+ */
+
+#include "physicsSolvers/fluidFlow/kernels/compositional/CompositionalMultiphaseHybridFVMKernels.hpp"
+#include "physicsSolvers/fluidFlow/kernels/compositional/CompositionalMultiphaseHybridFVMKernels_impl.hpp"
+
+namespace geos
+{
+namespace compositionalMultiphaseHybridFVMKernels
+{
+
+template void evaluateBCFaceProperties< 1, 2 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 1, 3 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 2, 2 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 2, 3 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 3, 2 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 3, 3 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 4, 2 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 4, 3 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 5, 2 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+template void evaluateBCFaceProperties< 5, 3 >(
+  integer const numPhases,
+  SortedArrayView< localIndex const > const & boundaryFaceSet,
+  arrayView1d< real64 const > const & facePres,
+  arrayView1d< real64 const > const & faceTemp,
+  arrayView2d< real64 const, compflow::USD_COMP > const & faceCompFrac,
+  arrayView2d< localIndex const > const & elemRegionList,
+  arrayView2d< localIndex const > const & elemSubRegionList,
+  arrayView2d< localIndex const > const & elemList,
+  localIndex const er,
+  localIndex const esr,
+  constitutive::MultiFluidBase & fluid,
+  constitutive::RelativePermeabilityBase & relperm,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMob,
+  arrayView2d< real64, compflow::USD_PHASE > const & facePhaseMassDens,
+  arrayView3d< real64, compflow::USD_PHASE_COMP > const & facePhaseCompFrac );
+
+} // namespace compositionalMultiphaseHybridFVMKernels
+} // namespace geos