GEOS
Classes | Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Protected Member Functions | Protected Attributes | List of all members
geos::SolidMechanicsMPM Class Reference

#include <SolidMechanicsMPM.hpp>

Inheritance diagram for geos::SolidMechanicsMPM:
Inheritance graph
[legend]

Classes

struct  viewKeyStruct
 

Public Types

enum class  TimeIntegrationOption : integer { QuasiStatic , ImplicitDynamic , ExplicitDynamic }
 
enum class  BoundaryConditionOption : integer { OUTFLOW , SYMMETRY }
 
- 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.
 

Public Member Functions

 SolidMechanicsMPM (const string &name, Group *const parent)
 
 SolidMechanicsMPM (SolidMechanicsMPM const &)=delete
 
 SolidMechanicsMPM (SolidMechanicsMPM &&)=default
 
SolidMechanicsMPMoperator= (SolidMechanicsMPM const &)=delete
 
SolidMechanicsMPMoperator= (SolidMechanicsMPM &&)=delete
 
virtual ~SolidMechanicsMPM () override
 
string getCatalogName () const override
 
virtual void initializePreSubGroups () override
 Called by Initialize() prior to initializing sub-Groups.
 
virtual void registerDataOnMesh (Group &meshBodies) override final
 Register wrappers that contain data on the mesh objects. More...
 
virtual real64 solverStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain) override
 entry function to perform a solver step More...
 
virtual real64 explicitStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain) override
 Entry function for an explicit time integration step. More...
 
virtual void updateState (DomainPartition &domain) override final
 Recompute all dependent quantities from primary variables (including constitutive models) More...
 
virtual bool execute (real64 const time_n, real64 const dt, integer const cycleNumber, integer const eventCounter, real64 const eventProgress, DomainPartition &domain) override
 
template<typename CONSTITUTIVE_BASE , typename KERNEL_WRAPPER , typename ... PARAMS>
void assemblyLaunch (DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs, PARAMS &&... params)
 
template<typename ... PARAMS>
real64 explicitKernelDispatch (MeshLevel &mesh, arrayView1d< string const > const &targetRegions, string const &finiteElementName, real64 const dt, std::string const &elementListName)
 
void initialize (NodeManager &nodeManager, ParticleManager &particleManager, SpatialPartition &partition)
 
void resizeGrid (SpatialPartition &partition, NodeManager &nodeManager, real64 const dt)
 
void syncGridFields (std::vector< std::string > const &fieldNames, DomainPartition &domain, NodeManager &nodeManager, MeshLevel &mesh, MPI_Op op)
 
void singleFaceVectorFieldSymmetryBC (const int face, arrayView3d< real64 > const &vectorMultiField, arrayView2d< real64 const, nodes::REFERENCE_POSITION_USD > const gridPosition, Group &nodeSets)
 
void enforceGridVectorFieldSymmetryBC (arrayView3d< real64 > const &vectorMultiField, arrayView2d< real64 const, nodes::REFERENCE_POSITION_USD > const gridPosition, Group &nodeSets)
 
void applyEssentialBCs (const real64 dt, const real64 time_n, NodeManager &nodeManager)
 
void computeGridSurfaceNormals (ParticleManager &particleManager, NodeManager &nodeManager)
 
void normalizeGridSurfaceNormals (arrayView2d< real64 const > const &gridMass, arrayView3d< real64 > const &gridSurfaceNormal)
 
void computeContactForces (real64 const dt, arrayView2d< real64 const > const &gridMass, arrayView2d< real64 const > const &gridDamage, arrayView2d< real64 const > const &gridMaxDamage, arrayView3d< real64 const > const &gridVelocity, arrayView3d< real64 const > const &gridMomentum, arrayView3d< real64 const > const &gridSurfaceNormal, arrayView3d< real64 const > const &gridMaterialPosition, arrayView3d< real64 > const &gridContactForce)
 
void computePairwiseNodalContactForce (int const &separable, real64 const &dt, real64 const &mA, real64 const &mB, arraySlice1d< real64 const > const vA, arraySlice1d< real64 const > const GEOS_UNUSED_PARAM(vB), arraySlice1d< real64 const > const qA, arraySlice1d< real64 const > const qB, arraySlice1d< real64 const > const nA, arraySlice1d< real64 const > const nB, arraySlice1d< real64 const > const xA, arraySlice1d< real64 const > const xB, arraySlice1d< real64 > const fA, arraySlice1d< real64 > const fB)
 
void computeOrthonormalBasis (const real64 *e1, real64 *e2, real64 *e3)
 
void setGridFieldLabels (NodeManager &nodeManager)
 
void solverProfiling (std::string label)
 
void solverProfilingIf (std::string label, bool condition)
 
real64 computeNeighborList (ParticleManager &particleManager)
 
void optimizeBinSort (ParticleManager &particleManager)
 
real64 kernel (real64 const &r)
 
void kernelGradient (arraySlice1d< real64 const > const x, std::vector< real64 > &xp, real64 const &r, real64 *result)
 
real64 computeKernelField (arraySlice1d< real64 const > const x, arrayView2d< real64 const > const xp, arrayView1d< real64 const > const Vp, arrayView1d< real64 const > const fp)
 
void computeKernelFieldGradient (arraySlice1d< real64 const > const x, std::vector< std::vector< real64 > > &xp, std::vector< real64 > &Vp, std::vector< real64 > &fp, arraySlice1d< real64 > const result)
 
void computeKernelVectorGradient (arraySlice1d< real64 const > const x, std::vector< std::vector< real64 > > &xp, std::vector< real64 > &Vp, std::vector< std::vector< real64 > > &fp, arraySlice2d< real64 > const result)
 
void computeDamageFieldGradient (ParticleManager &particleManager)
 
void updateSurfaceFlagOverload (ParticleManager &particleManager)
 
void projectDamageFieldGradientToGrid (ParticleManager &particleManager, NodeManager &nodeManager)
 
void updateDeformationGradient (real64 dt, ParticleManager &particleManager)
 
void updateConstitutiveModelDependencies (ParticleManager &particleManager)
 
void updateStress (real64 dt, ParticleManager &particleManager)
 
void particleKinematicUpdate (ParticleManager &particleManager)
 
void computeAndWriteBoxAverage (const real64 dt, const real64 time_n, ParticleManager &particleManager)
 
void initializeGridFields (NodeManager &nodeManager)
 
void boundaryConditionUpdate (real64 dt, real64 time_n)
 
void particleToGrid (ParticleManager &particleManager, NodeManager &nodeManager)
 
void gridTrialUpdate (real64 dt, NodeManager &nodeManager)
 
void enforceContact (real64 dt, DomainPartition &domain, ParticleManager &particleManager, NodeManager &nodeManager, MeshLevel &mesh)
 
void interpolateFTable (real64 dt, real64 time_n)
 
void gridToParticle (real64 dt, ParticleManager &particleManager, NodeManager &nodeManager)
 
void updateSolverDependencies (ParticleManager &particleManager)
 
real64 getStableTimeStep (ParticleManager &particleManager)
 
void deleteBadParticles (ParticleManager &particleManager)
 
void printProfilingResults ()
 
void computeSurfaceFlags (ParticleManager &particleManager)
 
void computeSphF (ParticleManager &particleManager)
 
int evaluateSeparabilityCriterion (localIndex const &A, localIndex const &B, real64 const &damageA, real64 const &damageB, real64 const &maxDamageA, real64 const &maxDamageB)
 
void flagOutOfRangeParticles (ParticleManager &particleManager)
 
void computeRVectors (ParticleManager &particleManager)
 
void cpdiDomainScaling (ParticleManager &particleManager)
 
void resizeMappingArrays (ParticleManager &particleManager)
 
void populateMappingArrays (ParticleManager &particleManager, NodeManager &nodeManager)
 
- 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 void initialize_postMeshGeneration () override
 Initialization tasks after mesh generation is completed.
 
void generateMeshTargetsFromTargetRegions (Group const &meshBodies)
 Generate mesh targets from target regions. 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...
 
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 setNextDt (real64 const &currentDt, DomainPartition &domain)
 function to set the next time step size More...
 
virtual real64 setNextDtBasedOnNewtonIter (real64 const &currentDt)
 function to set the next time step size based on Newton convergence More...
 
virtual real64 setNextDtBasedOnStateChange (real64 const &currentDt, DomainPartition &domain)
 function to set the next dt based on state change More...
 
virtual real64 setNextDtBasedOnCFL (real64 const &currentDt, DomainPartition &domain)
 function to set the next dt based on state change 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 implicitStepSetup (real64 const &time_n, real64 const &dt, DomainPartition &domain)
 function to perform setup for implicit timestep More...
 
virtual void setupDofs (DomainPartition const &domain, DofManager &dofManager) const
 Populate degree-of-freedom manager with fields relevant to this solver. 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 assembleSystem (real64 const time, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs)
 function to assemble the linear system matrix and rhs More...
 
virtual void applyBoundaryConditions (real64 const time, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs)
 apply boundary condition to 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 real64 calculateResidualNorm (real64 const &time, real64 const &dt, DomainPartition const &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localRhs)
 calculate the norm of the global system residual 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 bool checkSystemSolution (DomainPartition &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor)
 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)
 Function to determine if the solution vector should be scaled back in order to maintain a known constraint. More...
 
virtual void applySystemSolution (DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor, real64 const dt, DomainPartition &domain)
 Function to apply the solution vector to the state. More...
 
virtual bool updateConfiguration (DomainPartition &domain)
 updates the configuration (if needed) based on the state after a converged Newton loop. More...
 
virtual void outputConfigurationStatistics (DomainPartition const &domain) const
 
virtual void resetConfigurationToBeginningOfStep (DomainPartition &domain)
 resets the configuration to the beginning of the time-step. More...
 
virtual bool resetConfigurationToDefault (DomainPartition &domain) const
 resets the configuration to the default value. More...
 
virtual void resetStateToBeginningOfStep (DomainPartition &domain)
 reset state of physics back to the beginning of the step. More...
 
virtual void implicitStepComplete (real64 const &time, real64 const &dt, DomainPartition &domain)
 perform cleanup for implicit timestep 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 bool checkSequentialSolutionIncrements (DomainPartition &domain) const
 Check if the solution increments are ok to use. 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...
 
virtual void synchronizeNonlinearSolverParameters ()
 syncronize the nonlinear solver parameters.
 
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
 

Static Public Member Functions

static string catalogName ()
 
- 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...
 

Public Attributes

geos::SolidMechanicsMPM::viewKeyStruct solidMechanicsViewKeys
 

Protected Member Functions

virtual void postInputInitialization () override final
 
- 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 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::vector< array2d< localIndex > > m_mappedNodes
 
std::vector< array2d< real64 > > m_shapeFunctionValues
 
std::vector< array3d< real64 > > m_shapeFunctionGradientValues
 
int m_solverProfiling
 
std::vector< real64m_profilingTimes
 
std::vector< std::string > m_profilingLabels
 
TimeIntegrationOption m_timeIntegrationOption
 
int m_prescribedBcTable
 
array1d< int > m_boundaryConditionTypes
 
array2d< real64m_bcTable
 
int m_prescribedBoundaryFTable
 
Path m_fTablePath
 
int m_fTableInterpType
 
array2d< real64m_fTable
 
array1d< real64m_domainF
 
array1d< real64m_domainL
 
int m_boxAverageHistory
 
int m_reactionHistory
 
int m_needsNeighborList
 
real64 m_neighborRadius
 
int m_binSizeMultiplier
 
int m_useDamageAsSurfaceFlag
 
int m_cpdiDomainScaling
 
real64 m_smallMass
 
int m_numContactGroups
 
int m_numContactFlags
 
int m_numVelocityFields
 
real64 m_separabilityMinDamage
 
int m_treatFullyDamagedAsSingleField
 
int m_surfaceDetection
 
int m_damageFieldPartitioning
 
int m_contactGapCorrection
 
real64 m_frictionCoefficient
 
int m_planeStrain
 
int m_numDims
 
array1d< real64m_hEl
 
array1d< real64m_xLocalMin
 
array1d< real64m_xLocalMax
 
array1d< real64m_xLocalMinNoGhost
 
array1d< real64m_xLocalMaxNoGhost
 
array1d< real64m_xGlobalMin
 
array1d< real64m_xGlobalMax
 
array1d< real64m_partitionExtent
 
array1d< real64m_domainExtent
 
array1d< int > m_nEl
 
array3d< int > m_ijkMap
 
- 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

- 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

This class implements a material point method solution to the equations of motion.

Definition at line 45 of file SolidMechanicsMPM.hpp.

Member Enumeration Documentation

◆ BoundaryConditionOption

The options for essential boundary conditions

Enumerator
OUTFLOW 

Outflow.

SYMMETRY 

Symmetry.

Definition at line 66 of file SolidMechanicsMPM.hpp.

◆ TimeIntegrationOption

The options for time integration

Enumerator
QuasiStatic 

QuasiStatic.

ImplicitDynamic 

ImplicitDynamic.

ExplicitDynamic 

ExplicitDynamic.

Definition at line 54 of file SolidMechanicsMPM.hpp.

Constructor & Destructor Documentation

◆ SolidMechanicsMPM()

geos::SolidMechanicsMPM::SolidMechanicsMPM ( const string name,
Group *const  parent 
)

Constructor

Parameters
nameThe name of the solver instance
parentthe parent group of the solver

◆ ~SolidMechanicsMPM()

virtual geos::SolidMechanicsMPM::~SolidMechanicsMPM ( )
overridevirtual

destructor

Member Function Documentation

◆ catalogName()

static string geos::SolidMechanicsMPM::catalogName ( )
inlinestatic
Returns
The string that may be used to generate a new instance from the PhysicsSolverBase::CatalogInterface::CatalogType

Definition at line 95 of file SolidMechanicsMPM.hpp.

◆ execute()

virtual bool geos::SolidMechanicsMPM::execute ( real64 const  time_n,
real64 const  dt,
integer const  cycleNumber,
integer const  eventCounter,
real64 const  eventProgress,
DomainPartition domain 
)
overridevirtual

This method is called when its host event is triggered

Reimplemented from geos::PhysicsSolverBase.

◆ getCatalogName()

string geos::SolidMechanicsMPM::getCatalogName ( ) const
inlineoverridevirtual

Returns
Get the final class Catalog name

Implements geos::PhysicsSolverBase.

Definition at line 99 of file SolidMechanicsMPM.hpp.

◆ postInputInitialization()

virtual void geos::SolidMechanicsMPM::postInputInitialization ( )
finaloverrideprotectedvirtual

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

Reimplemented from geos::dataRepository::Group.

◆ registerDataOnMesh()

virtual void geos::SolidMechanicsMPM::registerDataOnMesh ( Group MeshBodies)
finaloverridevirtual

Register wrappers that contain data on the mesh objects.

Parameters
MeshBodiesthe group of mesh bodies

Reimplemented from geos::PhysicsSolverBase.


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