Inheritance diagram for geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >:

[legend]

Classes
struct	viewKeyStruct

Public Types
enum class	SolverType : integer { Flow = 0 , SolidMechanics = 1 }

using	Base = CoupledSolver< FLOW_SOLVER, MECHANICS_SOLVER >

Public Types inherited from geos::PhysicsSolverBase
enum class	StatsOutputType : integer { none , iteration , convergence , all }
	Type of the stat output.

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.

Public Member Functions
	PoromechanicsSolver (const string &name, dataRepository::Group *const parent)
	main constructor for CoupledSolver Objects More...

virtual void	postInputInitialization () override

virtual void	setConstitutiveNamesCallSuper (ElementSubRegionBase &subRegion) const override final
	This function sets constitutive name fields on an ElementSubRegionBase, and calls the base function it overrides. More...

virtual void	initializePreSubGroups () override
	Called by Initialize() prior to initializing sub-Groups.

virtual void	registerDataOnMesh (dataRepository::Group &meshBodies) override
	Register wrappers that contain data on the mesh objects. More...

virtual void	implicitStepSetup (real64 const &time_n, real64 const &dt, DomainPartition &domain) override
	function to perform setup for implicit timestep More...

virtual void	setupDofs (DomainPartition const &domain, DofManager &dofManager) const override
	Populate degree-of-freedom manager with fields relevant to this solver. More...

virtual void	setupCoupling (DomainPartition const &GEOS_UNUSED_PARAM(domain), DofManager &dofManager) const override

virtual void	assembleSystem (real64 const time, 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	assembleElementBasedTerms (real64 const time_n, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs)=0

virtual bool	checkSequentialConvergence (integer const cycleNumber, integer const iter, real64 const &time_n, real64 const &dt, DomainPartition &domain) override

MECHANICS_SOLVER *	solidMechanicsSolver () const
	accessor for the pointer to the solid mechanics solver More...

FLOW_SOLVER *	flowSolver () const
	accessor for the pointer to the flow solver More...

void	setStressInitialization (bool const performStressInitialization)

void	updateStabilizationParameters (DomainPartition &domain) const

void	updateBulkDensity (DomainPartition &domain)

Public Member Functions inherited from geos::CoupledSolver< FLOW_SOLVER, SolidMechanicsLagrangianFEM >
	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

CoupledSolver &	operator= (CoupledSolver const &)=delete
	deleted assignment operator

CoupledSolver &	operator= (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	implicitStepComplete (real64 const &time_n, real64 const &dt, DomainPartition &domain) override
	perform cleanup for implicit timestep 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 void	updateAndWriteConvergenceStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, integer const iteration) override
	Update the convergence information and write then into a CSV file. More...

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	setNextDt (real64 const &currentTime, real64 const &currentDt, DomainPartition &domain) override
	function to set the next time step size More...

virtual void	cleanup (real64 const time_n, integer const cycleNumber, integer const eventCounter, real64 const eventProgress, DomainPartition &domain) override

virtual bool	checkSequentialSolutionIncrements (DomainPartition &domain) const override
	Check if the solution increments are ok to use. More...

virtual bool	updateConfiguration (DomainPartition &domain, integer const configurationLoopIter) 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...

virtual void	synchronizeNonlinearSolverParameters () override
	synchronize the nonlinear solver parameters.

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.

PhysicsSolverBase &	operator= (PhysicsSolverBase const &)=delete
	Deleted copy assignment operator.

PhysicsSolverBase &	operator= (PhysicsSolverBase &&)=delete
	Deleted move assignment operator.

virtual string	getCatalogName () const =0

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...

ParallelMatrix &	getSystemMatrix ()
	Getter for system matrix. More...

ParallelMatrix const &	getSystemMatrix () const
	Getter for system rhs vector. More...

ParallelVector &	getSystemRhs ()
	Getter for system rhs vector. More...

ParallelVector const &	getSystemRhs () const
	Getter for system rhs vector. More...

ParallelVector &	getSystemSolution ()
	Getter for system solution vector. More...

ParallelVector const &	getSystemSolution () const
	Getter for system solution vector. More...

DofManager &	getDofManager ()
	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	setNextDtBasedOnIterNumber (real64 const &currentDt)
	function to set the next time step size based on convergence More...

virtual real64	setNextDtBasedOnStateChange (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, integer const newtonIter, 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, integer const newtonIter, 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...

virtual void	setSparsityPattern (DomainPartition &domain, DofManager &dofManager, CRSMatrix< real64, globalIndex > &localMatrix, SparsityPattern< globalIndex > &pattern)
	Set the sparsity pattern of the linear system matrix. More...

virtual std::unique_ptr< PreconditionerBase< LAInterface > >	createPreconditioner (DomainPartition &domain) const
	Create a preconditioner for this solver's linear system. 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 Group *	createChild (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...

LinearSolverParameters &	getLinearSolverParameters ()
	accessor for the linear solver parameters. More...

LinearSolverParameters const &	getLinearSolverParameters () const
	const accessor for the linear solver parameters. More...

NonlinearSolverParameters &	getNonlinearSolverParameters ()
	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...

string_array const &	getTargetRegionNames () const
	return the list of target regions 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...

IterationsStatistics &	getIterationStats ()

IterationsStatistics const &	getIterationStats () const

ConvergenceStatistics &	getConvergenceStats ()

ConvergenceStatistics const &	getConvergenceStats () const

SolverStatistics &	getSolverStatistics ()
	accessor for the solver statistics. More...

SolverStatistics const &	getSolverStatistics () const
	const accessor for the solver statistics. More...

map< std::pair< string, string >, string_array > const &	getMeshTargets () const
	accessor for m_meshTargets More...

bool	detectOscillations () const
	Detect oscillations in the solution. More...

Public Member Functions inherited from geos::ExecutableGroup
	ExecutableGroup ()=delete
	deleted default constructor

	ExecutableGroup (ExecutableGroup const &)=delete
	deleted copy constructor

	ExecutableGroup (ExecutableGroup &&)=default
	default move constructor

ExecutableGroup &	operator= (ExecutableGroup const &)=delete
	deleted assignment operator

ExecutableGroup &	operator= (ExecutableGroup &&)=delete
	deleted move operator

virtual	~ExecutableGroup () override=default
	Default destructor for the 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.

Group &	operator= (Group const &)=delete
	Deleted copy assignment operator. More...

Group &	operator= (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...

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...

subGroupMap &	getSubGroups ()
	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...

stdVector< string >	getSubGroupsNames () 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...

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...

WrapperBase &	registerWrapper (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 >
WrapperBase &	getWrapperBase (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...

wrapperMap &	wrappers ()
	Get access to the internal wrapper storage. More...

indexType	numWrappers () const
	Return the number of wrappers. More...

stdVector< string >	getWrappersNames () 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

Group &	getParent ()
	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...

Group &	setSizedFromParent (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	setLogLevel (integer const logLevel)
	Set verbosity level. More...

integer	getLogLevel () const

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 &targetNodePos)
	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.

Static Public Member Functions
static string	coupledSolverAttributePrefix ()
	String used to form the solverName used to register solvers in CoupledSolver.

Static Public Member Functions inherited from geos::PhysicsSolverBase
static CatalogInterface::CatalogType &	getCatalog ()
	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::CatalogType &	getCatalog ()
	Get the singleton catalog for this Group. More...

static string	processInputName (xmlWrapper::xmlNode const &targetNode, xmlWrapper::xmlNodePos const &targetNodePos, string_view parentNodeName, xmlWrapper::xmlNodePos const &parentNodePos, std::set< string > &siblingNames)

Public Attributes
std::tuple< SOLVERS *... >	m_solvers
	Pointers of the single-physics solvers.

Protected Member Functions
virtual void	initializePostInitialConditionsPreSubGroups () override
	Called by InitializePostInitialConditions() prior to initializing sub-Groups.

template<typename TYPE_LIST , typename KERNEL_WRAPPER , typename ... PARAMS>
real64	assemblyLaunch (MeshLevel &mesh, DofManager const &dofManager, string_array const &regionNames, string const &materialNamesString, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs, real64 const dt, PARAMS &&... params)

void	recordAverageMeanTotalStressIncrement (DomainPartition &domain, array1d< real64 > &averageMeanTotalStressIncrement)

void	applyAcceleratedAverageMeanTotalStressIncrement (DomainPartition &domain, array1d< real64 > &averageMeanTotalStressIncrement)

real64	computeAitkenRelaxationFactor (array1d< real64 > const &s0, array1d< real64 > const &s1, array1d< real64 > const &s1_tilde, array1d< real64 > const &s2_tilde, real64 const omega0)

array1d< real64 >	computeUpdate (array1d< real64 > const &s1, array1d< real64 > const &s2_tilde, real64 const omega1)

void	startSequentialIteration (integer const &iter, DomainPartition &domain) override

void	finishSequentialIteration (integer const &iter, DomainPartition &domain) override

virtual void	mapSolutionBetweenSolvers (DomainPartition &domain, integer const solverType) override
	Maps the solution obtained from one solver to the fields used by the other solver(s) More...

void	averageMeanTotalStressIncrement (DomainPartition &domain)
	Helper function to average the mean total stress increment over quadrature points. More...

virtual void	setMGRStrategy ()=0

virtual void	updateBulkDensity (ElementSubRegionBase &subRegion)=0

virtual string	getFlowDofKey () const =0

virtual void	validateNonlinearAcceleration () override

Protected Member Functions inherited from geos::CoupledSolver< FLOW_SOLVER, SolidMechanicsLagrangianFEM >
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, integer const cycleNumber, DomainPartition &domain)
	Sequentially coupled solver step. It solves a nonlinear system of equations using a sequential approach. More...

Protected Member Functions inherited from geos::PhysicsSolverBase
real64	eisenstatWalker (real64 const newNewtonNorm, real64 const oldNewtonNorm, LinearSolverParameters::Krylov const &krylovParams)
	Eisenstat-Walker adaptive tolerance. More...

template<typename CONSTITUTIVE >
void	setConstitutiveName (ElementSubRegionBase &subRegion, string const &wrapperName, string const &constitutiveType) const
	Register wrapper with given name and store constitutive model name on the subregion. More...

Protected Member Functions inherited from geos::dataRepository::Group
virtual void	initializePostSubGroups ()
	Called by Initialize() after 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
integer	m_isThermal
	Flag to determine whether or not this is a thermal simulation.

bool	m_performStressInitialization
	Flag to indicate that the solver is going to perform stress initialization.

stabilization::StabilizationType	m_stabilizationType
	Type of stabilization used.

string_array	m_stabilizationRegionNames
	Names of regions where stabilization applied.

real64	m_stabilizationMultiplier
	Stabilization Multiplier.

array1d< real64 >	m_s0
	Member variables needed for Nonlinear Acceleration ( Aitken ). Naming convention follows ( Jiang & Tchelepi, 2019 )

array1d< real64 >	m_s1

array1d< real64 >	m_s1_tilde

array1d< real64 >	m_s2

array1d< real64 >	m_s2_tilde

real64	m_omega0

real64	m_omega1

Protected Attributes inherited from geos::CoupledSolver< FLOW_SOLVER, SolidMechanicsLagrangianFEM >
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_nextDt
	timestep of the next cycle

integer	m_allowNonConvergedLinearSolverSolution
	behavior in case of linear solver failure

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.

ParallelVector	m_scaling
	Diagonal scaling vector D (Ahat = D * A * D, bhat = D * b, x = D * xhat)

integer	m_usePhysicsScaling
	Flag to decide whether to apply physics-based scaling to the linear system.

CRSMatrix< real64, globalIndex >	m_localMatrix
	Local system matrix and rhs.

std::unique_ptr< LinearSolverBase< LAInterface > >	m_linearSolver
	Custom linear solver for the "native" solver type.

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

StatsOutputType	m_writeStatisticsCSV

LinearSolverParametersInput	m_linearSolverParameters
	Linear solver parameters.

LinearSolverResult	m_linearSolverResult
	Result of the last linear solver.

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.

ArrayOfArrays< real64 >	m_solutionHistory
	History of the solution vector, used for oscillation detection.

Additional Inherited Members
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...

template<typename CONSTITUTIVE_TYPE >
static CONSTITUTIVE_TYPE &	getConstitutiveModel (ElementSubRegionBase &subRegion)
	Get the Constitutive Model object. More...

Detailed Description

template<typename FLOW_SOLVER, typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>
class geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >

Definition at line 55 of file PoromechanicsSolver.hpp.

Constructor & Destructor Documentation

◆ PoromechanicsSolver()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::PoromechanicsSolver	(	const string &	name,
		dataRepository::Group *const	parent
	)

inline

main constructor for CoupledSolver Objects

Parameters

name	the name of this instantiation of CoupledSolver in the repository
parent	the parent group of this instantiation of CoupledSolver

Definition at line 80 of file PoromechanicsSolver.hpp.

Member Function Documentation

◆ assembleSystem()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

virtual void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::assembleSystem	(	real64 const	time_n,
		real64 const	dt,
		DomainPartition &	domain,
		DofManager const &	dofManager,
		CRSMatrixView< real64, globalIndex const > const &	localMatrix,
		arrayView1d< real64 > const &	localRhs
	)

inlineoverridevirtual

function to assemble the linear system matrix and rhs

Parameters

time	the time at the beginning of the step
dt	the desired timestep
domain	the domain partition
dofManager	degree-of-freedom manager associated with the linear system
localMatrix	the system matrix
localRhs	the system right-hand side vector

This function assembles the residual and the jacobian of the residual wrt the primary variables. In a stand alone physics solver, this function will fill a single block in the block system. However the capability to query the block system structure for any coupled blocks may be implemented to fill in off diagonal blocks of the system to enable coupling between solvers.

Note: This function must be overridden in the derived physics solver in order to use an implict solution method such as LinearImplicitStep() or NonlinearImplicitStep().

Fully-coupled assembly.

Reimplemented from geos::CoupledSolver< FLOW_SOLVER, SolidMechanicsLagrangianFEM >.

Reimplemented in geos::SinglePhasePoromechanicsConformingFractures< FLOW_SOLVER >, geos::SinglePhasePoromechanicsEmbeddedFractures, geos::SinglePhasePoromechanics< FLOW_SOLVER, MECHANICS_SOLVER >, geos::SinglePhasePoromechanics< SinglePhaseBase, SolidMechanicsEmbeddedFractures >, geos::SinglePhasePoromechanics< SinglePhaseBase, SolidMechanicsLagrangeContact >, geos::MultiphasePoromechanics< FLOW_SOLVER, MECHANICS_SOLVER >, geos::MultiphasePoromechanics< CompositionalMultiphaseBase, SolidMechanicsLagrangeContact >, and geos::HydrofractureSolver< POROMECHANICS_SOLVER >.

Definition at line 267 of file PoromechanicsSolver.hpp.

◆ averageMeanTotalStressIncrement()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::averageMeanTotalStressIncrement ( DomainPartition & domain )

inlineprotected

Helper function to average the mean total stress increment over quadrature points.

Parameters

[in] domain the domain partition

Definition at line 698 of file PoromechanicsSolver.hpp.

◆ flowSolver()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

FLOW_SOLVER* geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::flowSolver ( ) const

inline

accessor for the pointer to the flow solver

Returns: a pointer to the flow solver

Definition at line 344 of file PoromechanicsSolver.hpp.

◆ implicitStepSetup()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

virtual void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::implicitStepSetup	(	real64 const &	time_n,
		real64 const &	dt,
		DomainPartition &	domain
	)

inlineoverridevirtual

function to perform setup for implicit timestep

Parameters

time_n	the time at the beginning of the step
dt	the desired timestep
domain	the domain partition

This function should contain any step level initialization required to perform an implicit step.

Note: This function must be overridden in the derived physics solver in order to use an implict solution method such as LinearImplicitStep() or NonlinearImplicitStep().

Reimplemented from geos::CoupledSolver< FLOW_SOLVER, SolidMechanicsLagrangianFEM >.

Reimplemented in geos::HydrofractureSolver< POROMECHANICS_SOLVER >.

Definition at line 234 of file PoromechanicsSolver.hpp.

◆ mapSolutionBetweenSolvers()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

virtual void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::mapSolutionBetweenSolvers	(	DomainPartition &	domain,
		integer const	solverType
	)

inlineoverrideprotectedvirtual

Maps the solution obtained from one solver to the fields used by the other solver(s)

Parameters

domain	the domain partition
solverType	the index of the solver withing this coupled solver.

After the flow solver

After the solid mechanics solver

Reimplemented from geos::CoupledSolver< FLOW_SOLVER, SolidMechanicsLagrangianFEM >.

Reimplemented in geos::PoromechanicsConformingFractures< MultiphasePoromechanics, CompositionalMultiphaseBase >, and geos::PoromechanicsConformingFractures< SinglePhasePoromechanics, SinglePhaseBase >.

Definition at line 652 of file PoromechanicsSolver.hpp.

◆ postInputInitialization()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

virtual void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::postInputInitialization ( )

inlineoverridevirtual

This function provides capability to post process input values prior to any other initialization operations.

Reimplemented from geos::CoupledSolver< FLOW_SOLVER, SolidMechanicsLagrangianFEM >.

Reimplemented in geos::SinglePhasePoromechanicsEmbeddedFractures, geos::HydrofractureSolver< POROMECHANICS_SOLVER >, and geos::MultiphasePoromechanicsConformingFractures< FLOW_SOLVER >.

Definition at line 113 of file PoromechanicsSolver.hpp.

◆ registerDataOnMesh()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

virtual void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::registerDataOnMesh ( dataRepository::Group & MeshBodies )

inlineoverridevirtual

Register wrappers that contain data on the mesh objects.

Parameters

MeshBodies the group of mesh bodies

Reimplemented from geos::PhysicsSolverBase.

Reimplemented in geos::HydrofractureSolver< POROMECHANICS_SOLVER >, and geos::SinglePhasePoromechanicsEmbeddedFractures.

Definition at line 183 of file PoromechanicsSolver.hpp.

◆ setConstitutiveNamesCallSuper()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

virtual void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::setConstitutiveNamesCallSuper ( ElementSubRegionBase & subRegion ) const

inlinefinaloverridevirtual

This function sets constitutive name fields on an ElementSubRegionBase, and calls the base function it overrides.

Parameters

subRegion The ElementSubRegionBase that will have constitutive names set.

Reimplemented from geos::PhysicsSolverBase.

Definition at line 133 of file PoromechanicsSolver.hpp.

◆ setupDofs()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

virtual void geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::setupDofs	(	DomainPartition const &	domain,
		DofManager &	dofManager
	)		const

inlineoverridevirtual

Populate degree-of-freedom manager with fields relevant to this solver.

Parameters

domain	the domain containing the mesh and fields
dofManager	degree-of-freedom manager associated with the linear system

Reimplemented from geos::PhysicsSolverBase.

Definition at line 248 of file PoromechanicsSolver.hpp.

◆ solidMechanicsSolver()

template<typename FLOW_SOLVER , typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>

MECHANICS_SOLVER* geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >::solidMechanicsSolver ( ) const

inline

accessor for the pointer to the solid mechanics solver

Returns: a pointer to the solid mechanics solver

Definition at line 335 of file PoromechanicsSolver.hpp.

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

coreComponents/physicsSolvers/multiphysics/PoromechanicsSolver.hpp

Classes

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

template<typename FLOW_SOLVER, typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM> class geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >

Constructor & Destructor Documentation

◆ PoromechanicsSolver()

Member Function Documentation

◆ assembleSystem()

◆ averageMeanTotalStressIncrement()

◆ flowSolver()

◆ implicitStepSetup()

◆ mapSolutionBetweenSolvers()

◆ postInputInitialization()

◆ registerDataOnMesh()

◆ setConstitutiveNamesCallSuper()

◆ setupDofs()

◆ solidMechanicsSolver()

template<typename FLOW_SOLVER, typename MECHANICS_SOLVER = SolidMechanicsLagrangianFEM>
class geos::PoromechanicsSolver< FLOW_SOLVER, MECHANICS_SOLVER >