GEOS
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
geos::CoupledSolver< SOLVERS > Class Template Reference
Inheritance diagram for geos::CoupledSolver< SOLVERS >:
Inheritance graph
[legend]

Public Member Functions

 CoupledSolver (const string &name, Group *const parent)
 main constructor for CoupledSolver Objects More...
 
 CoupledSolver (CoupledSolver const &)=delete
 deleted copy constructor
 
 CoupledSolver (CoupledSolver &&)=default
 default move constructor
 
CoupledSolveroperator= (CoupledSolver const &)=delete
 deleted assignment operator
 
CoupledSolveroperator= (CoupledSolver &&)=delete
 deleted move operator
 
void setSubSolvers ()
 Utility function to set the subsolvers pointers using the names provided by the user.
 
virtual void setupCoupling (DomainPartition const &domain, DofManager &dofManager) const
 Utility function to set the coupling between degrees of freedom. More...
 
virtual void assembleCouplingTerms (real64 const time_n, real64 const dt, DomainPartition const &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs)
 Utility function to compute coupling terms. More...
 
void setupDofs (DomainPartition const &domain, DofManager &dofManager) const override
 Populate degree-of-freedom manager with fields relevant to this solver. More...
 
virtual void implicitStepSetup (real64 const &time_n, real64 const &dt, DomainPartition &domain) override
 function to perform setup for implicit timestep More...
 
virtual void implicitStepComplete (real64 const &time_n, real64 const &dt, DomainPartition &domain) override
 perform cleanup for implicit timestep More...
 
virtual void assembleSystem (real64 const time_n, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) override
 function to assemble the linear system matrix and rhs More...
 
virtual void applySystemSolution (DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor, real64 const dt, DomainPartition &domain) override
 Function to apply the solution vector to the state. More...
 
virtual void updateState (DomainPartition &domain) override
 Recompute all dependent quantities from primary variables (including constitutive models) More...
 
virtual void resetStateToBeginningOfStep (DomainPartition &domain) override
 reset state of physics back to the beginning of the step. More...
 
real64 solverStep (real64 const &time_n, real64 const &dt, int const cycleNumber, DomainPartition &domain) override final
 
virtual real64 calculateResidualNorm (real64 const &time_n, real64 const &dt, DomainPartition const &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localRhs) override
 calculate the norm of the global system residual More...
 
virtual void applyBoundaryConditions (real64 const time_n, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) override
 apply boundary condition to system More...
 
virtual bool checkSystemSolution (DomainPartition &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor) override
 Function to check system solution for physical consistency and constraint violation. More...
 
virtual real64 scalingForSystemSolution (DomainPartition &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localSolution) override
 Function to determine if the solution vector should be scaled back in order to maintain a known constraint. More...
 
virtual real64 setNextDtBasedOnStateChange (real64 const &currentDt, DomainPartition &domain) override
 function to set the next dt based on state change More...
 
virtual real64 setNextDtBasedOnNewtonIter (real64 const &currentDt) override
 function to set the next time step size based on Newton convergence More...
 
virtual void cleanup (real64 const time_n, integer const cycleNumber, integer const eventCounter, real64 const eventProgress, DomainPartition &domain) override
 Called as the code exits the main run loop. More...
 
virtual bool checkSequentialSolutionIncrements (DomainPartition &domain) const override
 Check if the solution increments are ok to use. More...
 
virtual bool updateConfiguration (DomainPartition &domain) override
 updates the configuration (if needed) based on the state after a converged Newton loop. More...
 
virtual void outputConfigurationStatistics (DomainPartition const &domain) const override
 
virtual void resetConfigurationToBeginningOfStep (DomainPartition &domain) override
 resets the configuration to the beginning of the time-step. More...
 
virtual bool resetConfigurationToDefault (DomainPartition &domain) const override
 resets the configuration to the default value. More...
 
- Public Member Functions inherited from geos::PhysicsSolverBase
 PhysicsSolverBase (string const &name, Group *const parent)
 Constructor for PhysicsSolverBase. More...
 
 PhysicsSolverBase (PhysicsSolverBase &&)=default
 Move constructor for PhysicsSolverBase.
 
virtual ~PhysicsSolverBase () override
 Destructor for PhysicsSolverBase.
 
 PhysicsSolverBase ()=delete
 Deleted constructor.
 
 PhysicsSolverBase (PhysicsSolverBase const &)=delete
 Deleted copy constructor.
 
PhysicsSolverBaseoperator= (PhysicsSolverBase const &)=delete
 Deleted copy assignment operator.
 
PhysicsSolverBaseoperator= (PhysicsSolverBase &&)=delete
 Deleted move assignment operator.
 
virtual string getCatalogName () const =0
 
virtual void registerDataOnMesh (Group &MeshBodies) override
 Register wrappers that contain data on the mesh objects. More...
 
virtual void initialize_postMeshGeneration () override
 Initialization tasks after mesh generation is completed.
 
void generateMeshTargetsFromTargetRegions (Group const &meshBodies)
 Generate mesh targets from target regions. More...
 
virtual bool execute (real64 const time_n, real64 const dt, integer const cycleNumber, integer const eventCounter, real64 const eventProgress, DomainPartition &domain) override
 Main extension point of executable targets. More...
 
ParallelMatrixgetSystemMatrix ()
 Getter for system matrix. More...
 
ParallelMatrix const & getSystemMatrix () const
 Getter for system rhs vector. More...
 
ParallelVectorgetSystemRhs ()
 Getter for system rhs vector. More...
 
ParallelVector const & getSystemRhs () const
 Getter for system rhs vector. More...
 
ParallelVectorgetSystemSolution ()
 Getter for system solution vector. More...
 
ParallelVector const & getSystemSolution () const
 Getter for system solution vector. More...
 
DofManagergetDofManager ()
 Getter for degree-of-freedom manager. More...
 
DofManager const & getDofManager () const
 Getter for degree-of-freedom manager. More...
 
CRSMatrix< real64, globalIndex > & getLocalMatrix ()
 Getter for local matrix. More...
 
CRSMatrixView< real64 const, globalIndex const > getLocalMatrix () const
 Getter for local matrix. More...
 
virtual real64 solverStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
 entry function to perform a solver step More...
 
virtual real64 setNextDt (real64 const &currentDt, DomainPartition &domain)
 function to set the next time step size More...
 
virtual real64 setNextDtBasedOnCFL (real64 const &currentDt, DomainPartition &domain)
 function to set the next dt based on state change More...
 
virtual real64 explicitStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
 Entry function for an explicit time integration step. More...
 
virtual real64 nonlinearImplicitStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
 Function for a nonlinear implicit integration step. More...
 
virtual bool lineSearch (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, ParallelVector &rhs, ParallelVector &solution, real64 const scaleFactor, real64 &lastResidual)
 Function to perform line search. More...
 
virtual bool lineSearchWithParabolicInterpolation (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, ParallelVector &rhs, ParallelVector &solution, real64 const scaleFactor, real64 &lastResidual, real64 &residualNormT)
 Function to perform line search using a parabolic interpolation to find the scaling factor. More...
 
virtual real64 linearImplicitStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
 Function for a linear implicit integration step. More...
 
virtual void setupSystem (DomainPartition &domain, DofManager &dofManager, CRSMatrix< real64, globalIndex > &localMatrix, ParallelVector &rhs, ParallelVector &solution, bool const setSparsity=true)
 Set up the linear system (DOF indices and sparsity patterns) More...
 
void debugOutputSystem (real64 const &time, integer const cycleNumber, integer const nonlinearIteration, ParallelMatrix const &matrix, ParallelVector const &rhs) const
 Output the assembled linear system for debug purposes. More...
 
void debugOutputSolution (real64 const &time, integer const cycleNumber, integer const nonlinearIteration, ParallelVector const &solution) const
 Output the linear system solution for debug purposes. More...
 
virtual void solveLinearSystem (DofManager const &dofManager, ParallelMatrix &matrix, ParallelVector &rhs, ParallelVector &solution)
 function to apply a linear system solver to the assembled system. More...
 
virtual real64 getTimestepRequest (real64 const) override
 getter for the next timestep size More...
 
real64 getTimestepRequest ()
 getter for the next timestep size More...
 
virtual GroupcreateChild (string const &childKey, string const &childName) override
 creates a child group of of this PhysicsSolverBase instantiation More...
 
Timestamp getSystemSetupTimestamp () const
 getter for the timestamp of the system setup More...
 
Timestamp getMeshModificationTimestamp (DomainPartition &domain) const
 getter for the timestamp of the mesh modification on the mesh levels More...
 
void setSystemSetupTimestamp (Timestamp timestamp)
 set the timestamp of the system setup More...
 
R1Tensor const gravityVector () const
 return the value of the gravity vector specified in PhysicsSolverManager More...
 
virtual void saveSequentialIterationState (DomainPartition &domain)
 Save the state of the solver for sequential iteration. More...
 
LinearSolverParametersgetLinearSolverParameters ()
 accessor for the linear solver parameters. More...
 
LinearSolverParameters const & getLinearSolverParameters () const
 const accessor for the linear solver parameters. More...
 
NonlinearSolverParametersgetNonlinearSolverParameters ()
 accessor for the nonlinear solver parameters. More...
 
NonlinearSolverParameters const & getNonlinearSolverParameters () const
 const accessor for the nonlinear solver parameters. More...
 
localIndex targetRegionIndex (string const &regionName) const
 Get position of a given region within solver's target region list. More...
 
template<typename LAMBDA >
void forDiscretizationOnMeshTargets (Group const &meshBodies, LAMBDA &&lambda) const
 Loop over the target discretization on all mesh targets and apply callback. More...
 
template<typename LAMBDA >
void forDiscretizationOnMeshTargets (Group &meshBodies, LAMBDA &&lambda) const
 Loop over the target discretization on all mesh targets and apply callback. More...
 
string getDiscretizationName () const
 return the name of the discretization object More...
 
virtual bool registerCallback (void *func, const std::type_info &funcType) final override
 function to set the value of m_assemblyCallback More...
 
SolverStatisticsgetSolverStatistics ()
 accessor for the solver statistics. More...
 
SolverStatistics const & getSolverStatistics () const
 const accessor for the solver statistics. More...
 
map< std::pair< string, string >, array1d< string > > const & getMeshTargets () const
 accessor for m_meshTargets More...
 
- Public Member Functions inherited from geos::ExecutableGroup
virtual void signalToPrepareForExecution (real64 const time_n, real64 const dt, integer const cycle, DomainPartition &domain)
 Inform the object that it expects to execute during the next timestep. More...
 
void setTimesteppingBehavior (TimesteppingBehavior const timesteppingBehavior)
 Set the timestep behavior for a target. More...
 
TimesteppingBehavior getTimesteppingBehavior () const
 Get the target's time step behavior. More...
 
- Public Member Functions inherited from geos::dataRepository::Group
template<typename LOG_LEVEL_INFO >
std::enable_if_t< geos::is_log_level_info< LOG_LEVEL_INFO >, void > addLogLevel ()
 Append a levelCondition and a log description to the description of the wrapped object given a log info struct. Must be called in constructor. More...
 
virtual void reinit ()
 Performs re-initialization of certain variable depending on the solver being used.
 
 Group (string const &name, Group *const parent)
 Constructor. More...
 
 Group (string const &name, conduit::Node &rootNode)
 Constructor. More...
 
 Group (Group &&source)=default
 Move constructor. More...
 
virtual ~Group ()
 Destructor, deletes all Groups and Wrappers owned by this Group.
 
 Group ()=delete
 Deleted default constructor.
 
 Group (Group const &)=delete
 Deleted copy constructor.
 
Groupoperator= (Group const &)=delete
 Deleted copy assignment operator. More...
 
Groupoperator= (Group &&)=delete
 Deleted move assignment operator. More...
 
void printDataHierarchy (integer indent=0) const
 Prints the data hierarchy recursively. More...
 
string dumpInputOptions () const
 
string dumpSubGroupsNames () const
 
string dumpWrappersNames () const
 
template<typename T = Group>
T & registerGroup (string const &name, std::unique_ptr< T > newObject)
 Register a new Group as a sub-group of current Group. More...
 
template<typename T = Group>
T & registerGroup (string const &name, T *newObject)
 Register a new Group as a sub-group of current Group. More...
 
template<typename T = Group>
T & registerGroup (string const &name)
 Register a new Group as a sub-group of current Group. More...
 
template<typename T = Group>
T & registerGroup (subGroupMap::KeyIndex const &keyIndex)
 Register a new Group as a sub-group of current Group. More...
 
template<typename T = Group, typename TBASE = Group>
T & registerGroup (string const &name, string const &catalogName)
 Register a new Group as a sub-group of current Group. More...
 
void deregisterGroup (string const &name)
 Removes a child group from this group. More...
 
template<typename T = Group, typename KEY = void>
T * getGroupPointer (KEY const &key)
 Return a pointer to a sub-group of the current Group. More...
 
template<typename T = Group, typename KEY = void>
T const * getGroupPointer (KEY const &key) const
 Return a pointer to a sub-group of the current Group. More...
 
template<typename T = Group, typename KEY = void>
T & getGroup (KEY const &key)
 Return a reference to a sub-group of the current Group. More...
 
template<typename T = Group, typename KEY = void>
T const & getGroup (KEY const &key) const
 Return a reference to a sub-group of the current Group. More...
 
template<typename T = Group>
T & getGroupByPath (string const &path)
 Retrieve a group from the hierarchy using a path. More...
 
template<typename T = Group>
T const & getGroupByPath (string const &path) const
 Retrieve a group from the hierarchy using a path. More...
 
subGroupMapgetSubGroups ()
 Get the subgroups object. More...
 
subGroupMap const & getSubGroups () const
 Get the subgroups object. More...
 
localIndex numSubGroups () const
 return the number of sub groups in this Group More...
 
std::vector< stringgetSubGroupsNames () const
 
template<typename T = Group>
bool hasGroup (string const &name) const
 Check whether a sub-group exists. More...
 
template<typename T >
bool hasSubGroupOfType () const
 Check whether a sub-group exists by type. More...
 
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void forSubGroups (LAMBDA &&lambda)
 Apply the given functor to subgroups that can be casted to one of specified types. More...
 
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void forSubGroups (LAMBDA &&lambda) const
 Apply the given functor to subgroups that can be casted to one of specified types. More...
 
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void forSubGroupsIndex (LAMBDA &&lambda)
 Apply the given functor to subgroups that can be casted to one of specified types. More...
 
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void forSubGroupsIndex (LAMBDA &&lambda) const
 Apply the given functor to subgroups that can be casted to one of specified types. More...
 
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LOOKUP_CONTAINER , typename LAMBDA >
void forSubGroups (LOOKUP_CONTAINER const &subGroupKeys, LAMBDA &&lambda)
 Apply the given functor to subgroups that can be casted to one of specified types. More...
 
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LOOKUP_CONTAINER , typename LAMBDA >
void forSubGroups (LOOKUP_CONTAINER const &subGroupKeys, LAMBDA &&lambda) const
 Apply the given functor to subgroups that can be casted to one of specified types. More...
 
template<typename LAMBDA >
void forWrappers (LAMBDA &&lambda)
 Apply the given functor to wrappers. More...
 
template<typename LAMBDA >
void forWrappers (LAMBDA &&lambda) const
 Apply the given functor to wrappers. More...
 
template<typename TYPE , typename ... TYPES, typename LAMBDA >
void forWrappers (LAMBDA &&lambda)
 Apply the given functor to wrappers that can be cast to one of specified types. More...
 
template<typename TYPE , typename ... TYPES, typename LAMBDA >
void forWrappers (LAMBDA &&lambda) const
 Apply the given functor to wrappers that can be cast to one of specified types. More...
 
void initialize ()
 Run initialization functions on this and all subgroups. More...
 
virtual void initializationOrder (string_array &order)
 Sets the initialization order for sub-Groups. More...
 
void initializePostInitialConditions ()
 Initialization routine to be called after calling ApplyInitialConditions(). More...
 
void postRestartInitializationRecursive ()
 Initialization routine to be called after calling reading a restart file. More...
 
void processInputFileRecursive (xmlWrapper::xmlDocument &xmlDocument, xmlWrapper::xmlNode &targetNode)
 Recursively read values using ProcessInputFile() from the input file and put them into the wrapped values for this group. Also add the includes content to the xmlDocument when Include nodes are encountered. More...
 
void processInputFileRecursive (xmlWrapper::xmlDocument &xmlDocument, xmlWrapper::xmlNode &targetNode, xmlWrapper::xmlNodePos const &nodePos)
 Same as processInputFileRecursive(xmlWrapper::xmlDocument &, xmlWrapper::xmlNode &) but allow to reuse an existing xmlNodePos. More...
 
void postInputInitializationRecursive ()
 Recursively call postInputInitialization() to apply post processing after reading input values.
 
template<typename T , typename TBASE = T>
Wrapper< TBASE > & registerWrapper (string const &name, wrapperMap::KeyIndex::index_type *const rkey=nullptr)
 Create and register a Wrapper around a new object. More...
 
template<typename T , typename TBASE = T>
Wrapper< TBASE > & registerWrapper (Group::wrapperMap::KeyIndex const &viewKey)
 Create and register a Wrapper around a new object. More...
 
template<typename T >
Wrapper< T > & registerWrapper (string const &name, std::unique_ptr< T > newObject)
 Register a Wrapper around a given object and take ownership. More...
 
template<typename T >
Wrapper< T > & registerWrapper (string const &name, T *newObject)
 Register a Wrapper around an existing object, does not take ownership of the object. More...
 
WrapperBaseregisterWrapper (std::unique_ptr< WrapperBase > wrapper)
 Register and take ownership of an existing Wrapper. More...
 
void deregisterWrapper (string const &name)
 Removes a Wrapper from this group. More...
 
void generateDataStructureSkeleton (integer const level)
 Build a complete datastructure for schema generation. More...
 
virtual void expandObjectCatalogs ()
 Expand any catalogs in the data structure.
 
virtual void setSchemaDeviations (xmlWrapper::xmlNode schemaRoot, xmlWrapper::xmlNode schemaParent, integer documentationType)
 Inform the schema generator of any deviations between the xml and GEOS data structures. More...
 
virtual void registerDataOnMeshRecursive (Group &meshBodies)
 Calls RegisterDataOnMesh() recursively. More...
 
virtual localIndex packSize (string_array const &wrapperNames, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
 Get the size required to pack a list of wrappers. More...
 
virtual localIndex packSize (string_array const &wrapperNames, arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
 Get the size required to pack a list of indices within a list of wrappers. More...
 
localIndex packSize (arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
 Get the size required to pack a list of indices for all registered wrappers. More...
 
virtual localIndex pack (buffer_unit_type *&buffer, string_array const &wrapperNames, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
 Pack a list of wrappers to a buffer. More...
 
virtual localIndex pack (buffer_unit_type *&buffer, string_array const &wrapperNames, arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
 Pack a list of indices within a list of wrappers. More...
 
localIndex pack (buffer_unit_type *&buffer, arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
 Pack a list of indices for all registered wrappers. More...
 
virtual localIndex unpack (buffer_unit_type const *&buffer, arrayView1d< localIndex > &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events, MPI_Op op=MPI_REPLACE)
 Unpack a buffer. More...
 
template<typename KEY >
WrapperBase const & getWrapperBase (KEY const &key) const
 Return a reference to a WrapperBase stored in this group. More...
 
template<typename KEY >
WrapperBasegetWrapperBase (KEY const &key)
 Return a reference to a WrapperBase stored in this group. More...
 
indexType getWrapperIndex (string const &name) const
 
wrapperMap const & wrappers () const
 Get access to the internal wrapper storage. More...
 
wrapperMapwrappers ()
 Get access to the internal wrapper storage. More...
 
indexType numWrappers () const
 Return the number of wrappers. More...
 
std::vector< stringgetWrappersNames () const
 
template<typename LOOKUP_TYPE >
bool hasWrapper (LOOKUP_TYPE const &lookup) const
 Check if a wrapper exists. More...
 
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > const & getWrapper (LOOKUP_TYPE const &index) const
 Retrieve a Wrapper stored in this group. More...
 
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > & getWrapper (LOOKUP_TYPE const &index)
 Retrieve a Wrapper stored in this group. More...
 
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > const * getWrapperPointer (LOOKUP_TYPE const &index) const
 Retrieve a Wrapper stored in this group. More...
 
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > * getWrapperPointer (LOOKUP_TYPE const &index)
 Retrieve a Wrapper stored in this group. More...
 
template<typename T , typename LOOKUP_TYPE >
GEOS_DECLTYPE_AUTO_RETURN getReference (LOOKUP_TYPE const &lookup) const
 Look up a wrapper and get reference to wrapped object. More...
 
template<typename T , typename LOOKUP_TYPE >
T & getReference (LOOKUP_TYPE const &lookup)
 Look up a wrapper and get reference to wrapped object. More...
 
virtual void resize (localIndex const newSize)
 Resize the group and all contained wrappers that resize with parent. More...
 
virtual void reserve (indexType const newsize)
 Set the new capacity and reserve it in all wrappers that resize with parent. More...
 
localIndex capacity () const
 Get the "capacity" of the group, which determines the capacity of resizable wrappers. More...
 
localIndex size () const
 Get the "size" of the group, which determines the number of elements in resizable wrappers. More...
 
string const & getName () const
 Get group name. More...
 
string getPath () const
 Return the path of this Group in the data repository. Starts with '/' followed by the hierarchy of the children of the "Problem" in which the Group is. More...
 
DataContext const & getDataContext () const
 
template<typename KEY >
DataContext const & getWrapperDataContext (KEY key) const
 
GroupgetParent ()
 Access the group's parent. More...
 
Group const & getParent () const
 Access the group's parent. More...
 
bool hasParent () const
 
localIndex getIndexInParent () const
 Get the group's index within its parent group. More...
 
localIndex getSubGroupIndex (keyType const &key) const
 Get the index of a sub-Group within this group. More...
 
int sizedFromParent () const
 Check whether this Group is resized when its parent is resized. More...
 
GroupsetSizedFromParent (int val)
 Set whether this wrapper is resized when its parent is resized. More...
 
RestartFlags getRestartFlags () const
 Get flags that control restart output of this group. More...
 
void setRestartFlags (RestartFlags flags)
 Set flags that control restart output of this group. More...
 
InputFlags getInputFlags () const
 Get input flags for schema generation. More...
 
void setInputFlags (InputFlags flags)
 Set input flags for schema generation. More...
 
conduit::Node & getConduitNode ()
 Return the Conduit node object associated with this group. More...
 
conduit::Node const & getConduitNode () const
 Return the Conduit node object associated with this group. More...
 
void prepareToWrite ()
 Register the group and its wrappers with Conduit.
 
void finishWriting ()
 Write the group and its wrappers into Conduit.
 
void loadFromConduit ()
 Read the group and its wrappers from Conduit.
 
void enableLogLevelInput ()
 
void setLogLevel (integer const logLevel)
 Set verbosity level. More...
 
integer getLogLevel () const
 

Protected Member Functions

virtual real64 fullyCoupledSolverStep (real64 const &time_n, real64 const &dt, int const cycleNumber, DomainPartition &domain)
 Fully coupled solution approach solution step. More...
 
virtual real64 sequentiallyCoupledSolverStep (real64 const &time_n, real64 const &dt, int const cycleNumber, DomainPartition &domain)
 Sequentially coupled solver step. It solves a nonlinear system of equations using a sequential approach. More...
 
virtual void mapSolutionBetweenSolvers (DomainPartition &domain, integer const solverType)
 Maps the solution obtained from one solver to the fields used by the other solver(s) More...
 
virtual bool checkSequentialConvergence (int const &iter, real64 const &time_n, real64 const &dt, DomainPartition &domain)
 
virtual void postInputInitialization () override
 
virtual void validateNonlinearAcceleration ()
 
virtual void synchronizeNonlinearSolverParameters () override
 syncronize the nonlinear solver parameters.
 
virtual void startSequentialIteration (integer const &iter, DomainPartition &domain)
 
virtual void finishSequentialIteration (integer const &iter, DomainPartition &domain)
 
- Protected Member Functions inherited from geos::PhysicsSolverBase
real64 eisenstatWalker (real64 const newNewtonNorm, real64 const oldNewtonNorm, LinearSolverParameters::Krylov const &krylovParams)
 Eisenstat-Walker adaptive tolerance. More...
 
virtual void setConstitutiveNamesCallSuper (ElementSubRegionBase &subRegion) const
 This function sets constitutive name fields on an ElementSubRegionBase, and calls the base function it overrides. More...
 
- Protected Member Functions inherited from geos::dataRepository::Group
virtual void initializePreSubGroups ()
 Called by Initialize() prior to initializing sub-Groups.
 
virtual void initializePostSubGroups ()
 Called by Initialize() after to initializing sub-Groups.
 
virtual void initializePostInitialConditionsPreSubGroups ()
 Called by InitializePostInitialConditions() prior to initializing sub-Groups.
 
virtual void initializePostInitialConditionsPostSubGroups ()
 Called by InitializePostInitialConditions() after to initializing sub-Groups.
 
virtual void postRestartInitialization ()
 Performs initialization required after reading from a restart file.
 

Protected Attributes

std::tuple< SOLVERS *... > m_solvers
 Pointers of the single-physics solvers.
 
std::array< string, sizeof...(SOLVERS) > m_names
 Names of the single-physics solvers.
 
- Protected Attributes inherited from geos::PhysicsSolverBase
real64 m_cflFactor
 Courant–Friedrichs–Lewy factor for the timestep.
 
real64 m_maxStableDt
 maximum stable time step
 
real64 m_nextDt
 timestep of the next cycle
 
integer m_numTimestepsSinceLastDtCut
 Number of cycles since last timestep cut.
 
string m_discretizationName
 name of the FV discretization object in the data repository
 
DofManager m_dofManager
 Data structure to handle degrees of freedom.
 
ParallelMatrix m_matrix
 System matrix.
 
ParallelVector m_rhs
 System right-hand side vector.
 
ParallelVector m_solution
 System solution vector.
 
CRSMatrix< real64, globalIndexm_localMatrix
 Local system matrix and rhs.
 
std::unique_ptr< PreconditionerBase< LAInterface > > m_precond
 Custom preconditioner for the "native" iterative solver.
 
integer m_writeLinearSystem
 flag for debug output of matrix, rhs, and solution
 
LinearSolverParametersInput m_linearSolverParameters
 Linear solver parameters.
 
LinearSolverResult m_linearSolverResult
 Result of the last linear solve.
 
NonlinearSolverParameters m_nonlinearSolverParameters
 Nonlinear solver parameters.
 
SolverStatistics m_solverStatistics
 Solver statistics.
 
Timestamp m_systemSetupTimestamp
 Timestamp of the last call to setup system.
 
std::function< void(CRSMatrix< real64, globalIndex >, array1d< real64 >) > m_assemblyCallback
 Callback function for assembly step.
 
std::map< std::string, std::chrono::system_clock::duration > m_timers
 Timers for the aggregate profiling of the solver.
 

Additional Inherited Members

- Public Types inherited from geos::PhysicsSolverBase
using CatalogInterface = dataRepository::CatalogInterface< PhysicsSolverBase, string const &, Group *const >
 Type alias for catalog interface used by this class. See CatalogInterface.
 
- Public Types inherited from geos::ExecutableGroup
enum class  TimesteppingBehavior : integer { DeterminesTimeStepSize , DoesNotDetermineTimeStepSize }
 Timestepping type. More...
 
- Public Types inherited from geos::dataRepository::Group
using subGroupMap = MappedVector< Group, Group *, keyType, indexType >
 The template specialization of MappedVector to use for the collection of sub-Group objects.
 
using wrapperMap = MappedVector< WrapperBase, WrapperBase *, keyType, indexType >
 The template specialization of MappedVector to use for the collection wrappers objects.
 
using CatalogInterface = dataRepository::CatalogInterface< Group, string const &, Group *const >
 Type alias for catalog interface used by this class. See CatalogInterface.
 
- Static Public Member Functions inherited from geos::PhysicsSolverBase
static CatalogInterface::CatalogTypegetCatalog ()
 Get the singleton catalog for PhysicsSolverBase. More...
 
- Static Public Member Functions inherited from geos::dataRepository::Group
template<typename T0 , typename T1 , typename ... CASTTYPES, typename CONTAINERTYPE , typename LAMBDA >
static bool applyLambdaToContainer (CONTAINERTYPE &container, LAMBDA &&lambda)
 Apply a given functor to a container if the container can be cast to one of the specified types. More...
 
static CatalogInterface::CatalogTypegetCatalog ()
 Get the singleton catalog for this Group. More...
 
- Static Protected Member Functions inherited from geos::PhysicsSolverBase
template<typename CONSTITUTIVE_BASE_TYPE >
static string getConstitutiveName (ElementSubRegionBase const &subRegion)
 Get the Constitutive Name object. More...
 
template<typename CONSTITUTIVE_BASE_TYPE >
static string getConstitutiveName (ParticleSubRegionBase const &subRegion)
 Get the Constitutive Name object. More...
 
template<typename BASETYPE = constitutive::ConstitutiveBase, typename LOOKUP_TYPE >
static BASETYPE const & getConstitutiveModel (dataRepository::Group const &dataGroup, LOOKUP_TYPE const &key)
 Get the Constitutive Model object. More...
 
template<typename BASETYPE = constitutive::ConstitutiveBase, typename LOOKUP_TYPE >
static BASETYPE & getConstitutiveModel (dataRepository::Group &dataGroup, LOOKUP_TYPE const &key)
 Get the Constitutive Model object. More...
 

Detailed Description

template<typename ... SOLVERS>
class geos::CoupledSolver< SOLVERS >

Definition at line 33 of file CoupledSolver.hpp.

Constructor & Destructor Documentation

◆ CoupledSolver()

template<typename ... SOLVERS>
geos::CoupledSolver< SOLVERS >::CoupledSolver ( const string name,
Group *const  parent 
)
inline

main constructor for CoupledSolver Objects

Parameters
namethe name of this instantiation of CoupledSolver in the repository
parentthe parent group of this instantiation of CoupledSolver

Definition at line 43 of file CoupledSolver.hpp.

Member Function Documentation

◆ assembleCouplingTerms()

template<typename ... SOLVERS>
virtual void geos::CoupledSolver< SOLVERS >::assembleCouplingTerms ( real64 const  time_n,
real64 const  dt,
DomainPartition const &  domain,
DofManager const &  dofManager,
CRSMatrixView< real64, globalIndex const > const &  localMatrix,
arrayView1d< real64 > const &  localRhs 
)
inlinevirtual

Utility function to compute coupling terms.

Parameters
[in]time_nthe time at the beginning of the time step
[in]dtthe time step size
[in]domainthe domain partition @paran[in] dofManager the degree of freedom manager
[in]localMatrixthe local matrix
[in]localRhsthe local rhs

Reimplemented in geos::SinglePhaseReservoirAndWells< RESERVOIR_SOLVER >, and geos::CompositionalMultiphaseReservoirAndWells< RESERVOIR_SOLVER >.

Definition at line 119 of file CoupledSolver.hpp.

◆ checkSequentialSolutionIncrements()

template<typename ... SOLVERS>
virtual bool geos::CoupledSolver< SOLVERS >::checkSequentialSolutionIncrements ( DomainPartition domain) const
inlineoverridevirtual

Check if the solution increments are ok to use.

Parameters
domainthe domain partition
Returns
true if the solution increments are ok to use, false otherwise

Reimplemented from geos::PhysicsSolverBase.

Definition at line 352 of file CoupledSolver.hpp.

◆ fullyCoupledSolverStep()

template<typename ... SOLVERS>
virtual real64 geos::CoupledSolver< SOLVERS >::fullyCoupledSolverStep ( real64 const &  time_n,
real64 const &  dt,
int const  cycleNumber,
DomainPartition domain 
)
inlineprotectedvirtual

Fully coupled solution approach solution step.

Parameters
time_nthe current time
dttimestep size
cycleNumber
domainthe domain partition
Returns
real64 size of the accepted timestep

Definition at line 409 of file CoupledSolver.hpp.

◆ mapSolutionBetweenSolvers()

template<typename ... SOLVERS>
virtual void geos::CoupledSolver< SOLVERS >::mapSolutionBetweenSolvers ( DomainPartition domain,
integer const  solverType 
)
inlineprotectedvirtual

◆ outputConfigurationStatistics()

template<typename ... SOLVERS>
virtual void geos::CoupledSolver< SOLVERS >::outputConfigurationStatistics ( DomainPartition const &  domain) const
inlineoverridevirtual
Parameters
domainthe domain containing the mesh and fields

Reimplemented from geos::PhysicsSolverBase.

Definition at line 372 of file CoupledSolver.hpp.

◆ postInputInitialization()

template<typename ... SOLVERS>
virtual void geos::CoupledSolver< SOLVERS >::postInputInitialization ( )
inlineoverrideprotectedvirtual

◆ resetConfigurationToBeginningOfStep()

template<typename ... SOLVERS>
virtual void geos::CoupledSolver< SOLVERS >::resetConfigurationToBeginningOfStep ( DomainPartition domain)
inlineoverridevirtual

resets the configuration to the beginning of the time-step.

Parameters
domainthe domain containing the mesh and fields

Reimplemented from geos::PhysicsSolverBase.

Definition at line 380 of file CoupledSolver.hpp.

◆ resetConfigurationToDefault()

template<typename ... SOLVERS>
virtual bool geos::CoupledSolver< SOLVERS >::resetConfigurationToDefault ( DomainPartition domain) const
inlineoverridevirtual

resets the configuration to the default value.

Parameters
domainthe domain containing the mesh and fields
Returns
a bool that states whether the configuration was reset or not.

Reimplemented from geos::PhysicsSolverBase.

Definition at line 388 of file CoupledSolver.hpp.

◆ sequentiallyCoupledSolverStep()

template<typename ... SOLVERS>
virtual real64 geos::CoupledSolver< SOLVERS >::sequentiallyCoupledSolverStep ( real64 const &  time_n,
real64 const &  dt,
int const  cycleNumber,
DomainPartition domain 
)
inlineprotectedvirtual

Sequentially coupled solver step. It solves a nonlinear system of equations using a sequential approach.

Parameters
time_nthe current time
dttimestep size
cycleNumber
domainthe domain partition
Returns
real64 size of the accepted timestep

Sequential coupling loop

Definition at line 427 of file CoupledSolver.hpp.

◆ setupCoupling()

template<typename ... SOLVERS>
virtual void geos::CoupledSolver< SOLVERS >::setupCoupling ( DomainPartition const &  domain,
DofManager dofManager 
) const
inlinevirtual

◆ updateConfiguration()

template<typename ... SOLVERS>
virtual bool geos::CoupledSolver< SOLVERS >::updateConfiguration ( DomainPartition domain)
inlineoverridevirtual

updates the configuration (if needed) based on the state after a converged Newton loop.

Parameters
domainthe domain containing the mesh and fields
Returns
a bool that states whether the configuration used to solve the nonlinear loop is still valid or not.

Reimplemented from geos::PhysicsSolverBase.

Definition at line 362 of file CoupledSolver.hpp.


The documentation for this class was generated from the following file: