Feature/byer3/well debug

inputFiles/compositionalMultiphaseWell/compositional_multiphase_wells_1d.xml

@@ -26,6 +26,7 @@
     <CompositionalMultiphaseWell
       name="compositionalMultiphaseWell"
       logLevel="1"
+      writeSegDebug="1"
       targetRegions="{ wellRegion1, wellRegion2 }">
       <WellControls
         name="wellControls1"

src/cmake/GeosxOptions.cmake

@@ -90,7 +90,7 @@ option( GEOS_INSTALL_SCHEMA "Enables schema generation and installation" ON )
 
 option( GEOS_BUILD_OBJ_LIBS "Builds coreComponent modules as object libraries" OFF )
 
-option( GEOS_BUILD_SHARED_LIBS "Builds geosx_core as a shared library " ON )
+option( GEOS_BUILD_SHARED_LIBS "Builds geosx_core as a shared library " OFF )
 
 set( GEOS_PARALLEL_COMPILE_JOBS "" CACHE STRING "Maximum number of concurrent compilation jobs" )
 if( GEOS_PARALLEL_COMPILE_JOBS )

src/coreComponents/physicsSolvers/SolverStatistics.cpp

@@ -38,6 +38,18 @@ SolverStatistics::SolverStatistics( string const & name, Group * const parent )
     setApplyDefaultValue( 0 ).
     setDescription( "Number of time step cuts" );
 
+  registerWrapper( viewKeyStruct::numCurrentOuterLoopIterationsString(), &m_currentNumOuterLoopIterations ).
+    setApplyDefaultValue( 0 ).
+    setDescription( "Cumulative number of current outer loop iterations" );
+
+  registerWrapper( viewKeyStruct::numCurrentNonlinearIterationsString(), &m_currentNumNonlinearIterations ).
+    setApplyDefaultValue( 0 ).
+    setDescription( "Cumulative number of current nonlinear iterations" );
+
+  registerWrapper( viewKeyStruct::numCurrentLinearIterationsString(), &m_currentNumLinearIterations ).
+    setApplyDefaultValue( 0 ).
+    setDescription( "Cumulative number of current linear iterations" );
+
   registerWrapper( viewKeyStruct::numSuccessfulOuterLoopIterationsString(), &m_numSuccessfulOuterLoopIterations ).
     setApplyDefaultValue( 0 ).
     setDescription( "Cumulative number of successful outer loop iterations" );

src/coreComponents/physicsSolvers/SolverStatistics.hpp

@@ -140,6 +140,14 @@ class SolverStatistics : public dataRepository::Group
     /// String key for the number of time step cuts
     static constexpr char const * numTimeStepCutsString() { return "numTimeStepCuts"; }
 
+
+    /// String key for the current number of outer loop iterations
+    static constexpr char const * numCurrentOuterLoopIterationsString() { return "numCurrentOuterLoopIterations"; }
+    /// String key for the current number of nonlinear iterations
+    static constexpr char const * numCurrentNonlinearIterationsString() { return "numCurrentNonlinearIterations"; }
+    /// String key for the current number of linear iterations
+    static constexpr char const * numCurrentLinearIterationsString() { return "numCurrentLinearIterations"; }
+
     /// String key for the successful number of outer loop iterations
     static constexpr char const * numSuccessfulOuterLoopIterationsString() { return "numSuccessfulOuterLoopIterations"; }
     /// String key for the successful number of nonlinear iterations

src/coreComponents/physicsSolvers/fluidFlow/wells/CompositionalMultiphaseWell.hpp

@@ -241,6 +241,26 @@ class CompositionalMultiphaseWell : public WellSolverBase
                                           CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                           arrayView1d< real64 > const & localRhs ) override;
 
+  /**
+   * @brief apply a special treatment to the wells that are shut
+   * @param time_n the time at the previous converged time step
+   * @param dt the time step size
+   * @param domain the physical domain object
+   * @param dofManager degree-of-freedom manager associated with the linear system
+   * @param matrix the system matrix
+   * @param rhs the system right-hand side vector
+   */
+
+
+  virtual void outputWellDebug( real64 const time,
+                                real64 const dt,
+                                integer num_timesteps,
+                                integer current_newton_iteration,
+                                integer num_timestep_cuts,
+                                DomainPartition & domain,
+                                DofManager const & dofManager,
+                                CRSMatrixView< real64, globalIndex const > const & localMatrix,
+                                arrayView1d< real64 > const & localRhs ) override;
   /**
    * @brief Sets all the negative component densities (if any) to zero.
    * @param domain the physical domain object
@@ -351,6 +371,13 @@ class CompositionalMultiphaseWell : public WellSolverBase
                    real64 const & dt,
                    DomainPartition & domain ) override;
 
+  void printSegRates( real64 const & time_n,
+                      real64 const & dt,
+                      integer num_timesteps,
+                      integer num_timestep_cuts,
+                      integer current_newton_iteration,
+                      DomainPartition & domain );
+
 private:
 
   /**

src/coreComponents/physicsSolvers/fluidFlow/wells/SinglePhaseWell.cpp

@@ -614,6 +614,27 @@ void SinglePhaseWell::assembleVolumeBalanceTerms( DomainPartition const & GEOS_U
   // not implemented for single phase flow
 }
 
+
+void SinglePhaseWell::outputWellDebug( real64 const time,
+                                       real64 const dt,
+                                       integer num_timesteps,
+                                       integer current_newton_iteration,
+                                       integer num_timestep_cuts,
+                                       DomainPartition & domain,
+                                       DofManager const & dofManager,
+                                       CRSMatrixView< real64, globalIndex const > const & localMatrix,
+                                       arrayView1d< real64 > const & localRhs )
+{
+  GEOS_UNUSED_VAR( time );
+  GEOS_UNUSED_VAR( dt );
+  GEOS_UNUSED_VAR( num_timesteps );
+  GEOS_UNUSED_VAR( current_newton_iteration );
+  GEOS_UNUSED_VAR( num_timestep_cuts );
+  GEOS_UNUSED_VAR( domain );
+  GEOS_UNUSED_VAR( dofManager );
+  GEOS_UNUSED_VAR( localMatrix );
+  GEOS_UNUSED_VAR( localRhs );
+}
 void SinglePhaseWell::computePerforationRates( real64 const & time_n,
                                                real64 const & dt, DomainPartition & domain )
 {

src/coreComponents/physicsSolvers/fluidFlow/wells/SinglePhaseWell.hpp

@@ -224,6 +224,18 @@ class SinglePhaseWell : public WellSolverBase
                                           CRSMatrixView< real64, globalIndex const > const & localMatrix,
                                           arrayView1d< real64 > const & localRhs ) override;
 
+
+
+  virtual void outputWellDebug( real64 const time,
+                                real64 const dt,
+                                integer num_timesteps,
+                                integer current_newton_iteration,
+                                integer num_timestep_cuts,
+                                DomainPartition & domain,
+                                DofManager const & dofManager,
+                                CRSMatrixView< real64, globalIndex const > const & localMatrix,
+                                arrayView1d< real64 > const & localRhs )override;
+
   struct viewKeyStruct : WellSolverBase::viewKeyStruct
   {
     static constexpr char const * dofFieldString() { return "singlePhaseWellVars"; }