![]() |
KratosMultiphysics
KRATOS Multiphysics (Kratos) is a framework for building parallel, multi-disciplinary simulation software, aiming at modularity, extensibility, and high performance. Kratos is written in C++, and counts with an extensive Python interface.
|
Current class provides an implementation for contact builder and solving operations. (elimination) More...
#include <contact_residualbased_elimination_builder_and_solver.h>
Public Member Functions | |
Life Cycle | |
ContactResidualBasedEliminationBuilderAndSolver () | |
Default constructor. More... | |
ContactResidualBasedEliminationBuilderAndSolver (typename TLinearSolver::Pointer pNewLinearSystemSolver, Parameters ThisParameters) | |
Default constructor. (with parameters) More... | |
ContactResidualBasedEliminationBuilderAndSolver (typename TLinearSolver::Pointer pNewLinearSystemSolver) | |
~ContactResidualBasedEliminationBuilderAndSolver () override | |
![]() | |
ResidualBasedEliminationBuilderAndSolver () | |
Default constructor. More... | |
ResidualBasedEliminationBuilderAndSolver (typename TLinearSolver::Pointer pNewLinearSystemSolver, Parameters ThisParameters) | |
Default constructor. (with parameters) More... | |
ResidualBasedEliminationBuilderAndSolver (typename TLinearSolver::Pointer pNewLinearSystemSolver) | |
Constructor. More... | |
~ResidualBasedEliminationBuilderAndSolver () override | |
std::string | Info () const override |
Turn back information as a string. More... | |
void | PrintInfo (std::ostream &rOStream) const override |
Print information about this object. More... | |
void | PrintData (std::ostream &rOStream) const override |
Print object's data. More... | |
KRATOS_CLASS_POINTER_DEFINITION (ResidualBasedEliminationBuilderAndSolver) | |
Pointer definition of ResidualBasedEliminationBuilderAndSolverWithConstraints. More... | |
void | Build (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rb) override |
Function to perform the build of the RHS. The vector could be sized as the total number of dofs or as the number of unrestrained ones. More... | |
void | BuildLHS (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA) override |
Function to perform the building of the LHS. More... | |
void | BuildLHS_CompleteOnFreeRows (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA) override |
Build a rectangular matrix of size n*N where "n" is the number of unrestrained degrees of freedom and "N" is the total number of degrees of freedom involved. More... | |
void | SystemSolve (TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) override |
This is a call to the linear system solver. More... | |
void | SystemSolveWithPhysics (TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb, ModelPart &rModelPart) |
This is a call to the linear system solver (taking into account some physical particularities of the problem) More... | |
void | BuildAndSolve (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) override |
Function to perform the building and solving phase at the same time. More... | |
void | BuildRHSAndSolve (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) override |
Corresponds to the previews, but the System's matrix is considered already built and only the RHS is built again. More... | |
void | BuildRHS (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemVectorType &rb) override |
Function to perform the build of the RHS. More... | |
void | SetUpDofSet (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart) override |
Builds the list of the DofSets involved in the problem by "asking" to each element and condition its Dofs. More... | |
void | ResizeAndInitializeVectors (typename TSchemeType::Pointer pScheme, TSystemMatrixPointerType &pA, TSystemVectorPointerType &pDx, TSystemVectorPointerType &pb, ModelPart &rModelPart) override |
This method resize and initializes the system of euqations. More... | |
void | CalculateReactions (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) override |
This method computes the reactions. More... | |
void | ApplyDirichletConditions (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) override |
Applies the dirichlet conditions. This operation may be very heavy or completely unexpensive depending on the implementation chosen and on how the System Matrix is built. More... | |
void | Clear () override |
This function is intended to be called at the end of the solution step to clean up memory storage not needed. More... | |
int | Check (ModelPart &rModelPart) override |
This function is designed to be called once to perform all the checks needed on the input provided. Checks can be "expensive" as the function is designed to catch user's errors. More... | |
![]() | |
BuilderAndSolver () | |
Default constructor. More... | |
BuilderAndSolver (typename TLinearSolver::Pointer pNewLinearSystemSolver, Parameters ThisParameters) | |
Default constructor with Parameters. More... | |
BuilderAndSolver (typename TLinearSolver::Pointer pNewLinearSystemSolver) | |
Default constructor. More... | |
virtual | ~BuilderAndSolver () |
void | SetEchoLevel (int Level) |
It sets the level of echo for the solving strategy. More... | |
int | GetEchoLevel () const |
It returns the echo level. More... | |
virtual TSparseSpace::MatrixType & | GetConstraintRelationMatrix () |
This method returns constraint relation (T) matrix. More... | |
virtual TSparseSpace::VectorType & | GetConstraintConstantVector () |
This method returns constraint constant vector. More... | |
KRATOS_CLASS_POINTER_DEFINITION (BuilderAndSolver) | |
Pointer definition of BuilderAndSolver. More... | |
bool | GetCalculateReactionsFlag () const |
This method returns the flag mCalculateReactionsFlag. More... | |
void | SetCalculateReactionsFlag (bool flag) |
This method sets the flag mCalculateReactionsFlag. More... | |
bool | GetDofSetIsInitializedFlag () const |
This method returns the flag mDofSetIsInitialized. More... | |
void | SetDofSetIsInitializedFlag (bool DofSetIsInitialized) |
This method sets the flag mDofSetIsInitialized. More... | |
bool | GetReshapeMatrixFlag () const |
This method returns the flag mReshapeMatrixFlag. More... | |
void | SetReshapeMatrixFlag (bool ReshapeMatrixFlag) |
This method sets the flag mReshapeMatrixFlag. More... | |
unsigned int | GetEquationSystemSize () const |
This method returns the value mEquationSystemSize. More... | |
TLinearSolver::Pointer | GetLinearSystemSolver () const |
This method return the linear solver used. More... | |
void | SetLinearSystemSolver (typename TLinearSolver::Pointer pLinearSystemSolver) |
This method sets the linear solver to be used. More... | |
virtual void | BuildLHS_Complete (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA) |
It builds a matrix of size N*N where "N" is the total number of degrees of freedom involved. More... | |
virtual void | BuildAndSolveLinearizedOnPreviousIteration (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb, const bool MoveMesh) |
Function to perform the building and solving phase at the same time Linearizing with the database at the old iteration. More... | |
virtual void | ApplyDirichletConditions_LHS (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx) |
The same of the precedent but affecting only the LHS. More... | |
virtual void | ApplyDirichletConditions_RHS (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemVectorType &rDx, TSystemVectorType &rb) |
The same of the precedent but affecting only the RHS. More... | |
virtual void | ApplyRHSConstraints (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemVectorType &rb) |
Applies the constraints with master-slave relation matrix (RHS only) More... | |
virtual void | ApplyConstraints (typename TSchemeType::Pointer pScheme, ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rb) |
Applies the constraints with master-slave relation matrix. More... | |
virtual DofsArrayType & | GetDofSet () |
It allows to get the list of Dofs from the element. More... | |
virtual const DofsArrayType & | GetDofSet () const |
It allows to get the list of Dofs from the element. More... | |
virtual void | InitializeSolutionStep (ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) |
It applies certain operations at the system of equations at the beginning of the solution step. More... | |
virtual void | FinalizeSolutionStep (ModelPart &rModelPart, TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) |
It applies certain operations at the system of equations at the end of the solution step. More... | |
Protected Member Functions | |
Protected Operations | |
void | AssignSettings (const Parameters ThisParameters) override |
This method assigns settings to member variables. More... | |
![]() | |
void | Assemble (TSystemMatrixType &rA, TSystemVectorType &rb, const LocalSystemMatrixType &rLHSContribution, const LocalSystemVectorType &rRHSContribution, const Element::EquationIdVectorType &rEquationId) |
This method assembles the system. More... | |
virtual void | ConstructMatrixStructure (typename TSchemeType::Pointer pScheme, TSystemMatrixType &rA, ModelPart &rModelPart) |
This method constructs the relationship between the DoF. More... | |
void | AssembleLHS (TSystemMatrixType &rA, LocalSystemMatrixType &rLHSContribution, EquationIdVectorType &rEquationId) |
This method assembles the LHS of the system. More... | |
void | AssembleRowContributionFreeDofs (TSystemMatrixType &rA, const Matrix &rALocal, const IndexType i, const IndexType i_local, const Element::EquationIdVectorType &EquationId) |
This function is equivalent to the AssembleRowContribution of the block builder and solver. More... | |
IndexType | ForwardFind (const IndexType id_to_find, const IndexType start, const IndexType *index_vector) |
IndexType | BackwardFind (const IndexType id_to_find, const IndexType start, const IndexType *index_vector) |
![]() | |
virtual Parameters | ValidateAndAssignParameters (Parameters ThisParameters, const Parameters DefaultParameters) const |
This method validate and assign default parameters. More... | |
Type Definitions | |
using | BaseBuilderAndSolverType = BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > |
Builder and solver base class. More... | |
using | BaseType = ResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > |
Definitions dependent on the base class. More... | |
using | ClassType = ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > |
The definition of the current class. More... | |
using | TSchemeType = typename BaseType::TSchemeType |
Base types definitions. More... | |
using | TDataType = typename BaseType::TDataType |
using | DofsArrayType = typename BaseType::DofsArrayType |
using | TSystemMatrixType = typename BaseType::TSystemMatrixType |
using | TSystemVectorType = typename BaseType::TSystemVectorType |
using | DofType = typename ModelPart::DofType |
The definition of the dof type. More... | |
using | SizeType = std::size_t |
The size type. More... | |
using | IndexType = std::size_t |
The index type. More... | |
using | IndexSetType = std::unordered_set< IndexType > |
Index set definition. More... | |
KRATOS_CLASS_POINTER_DEFINITION (ContactResidualBasedEliminationBuilderAndSolver) | |
Pointer definition of ContactResidualBasedEliminationBuilderAndSolver. More... | |
Operations | |
BaseBuilderAndSolverType::Pointer | Create (typename TLinearSolver::Pointer pNewLinearSystemSolver, Parameters ThisParameters) const override |
Create method. More... | |
void | SetUpSystem (ModelPart &rModelPart) override |
It organises the dofset in order to speed up the building phase. More... | |
Parameters | GetDefaultParameters () const override |
This method provides the defaults parameters to avoid conflicts between the different constructors. More... | |
static std::string | Name () |
Returns the name of the class as used in the settings (snake_case format) More... | |
Additional Inherited Members | |
![]() | |
typedef BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > | BaseType |
Definition of the base class. More... | |
typedef ResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > | ClassType |
The definition of the current class. More... | |
typedef BaseType::SizeType | SizeType |
Definition of the classes from the base class. More... | |
typedef BaseType::IndexType | IndexType |
typedef BaseType::TSchemeType | TSchemeType |
typedef BaseType::TDataType | TDataType |
typedef BaseType::DofsArrayType | DofsArrayType |
typedef BaseType::TSystemMatrixType | TSystemMatrixType |
typedef BaseType::TSystemVectorType | TSystemVectorType |
typedef BaseType::LocalSystemVectorType | LocalSystemVectorType |
typedef BaseType::LocalSystemMatrixType | LocalSystemMatrixType |
typedef BaseType::TSystemMatrixPointerType | TSystemMatrixPointerType |
typedef BaseType::TSystemVectorPointerType | TSystemVectorPointerType |
typedef Element::EquationIdVectorType | EquationIdVectorType |
Definition of the equation id vector. More... | |
typedef Element::DofsVectorType | DofsVectorType |
typedef Node | NodeType |
Node definition. More... | |
typedef BaseType::NodesArrayType | NodesArrayType |
Containers definition. More... | |
typedef BaseType::ElementsArrayType | ElementsArrayType |
typedef BaseType::ConditionsArrayType | ConditionsArrayType |
typedef BaseType::ElementsContainerType | ElementsContainerType |
![]() | |
typedef BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > | ClassType |
The definition of the current class. More... | |
typedef std::size_t | SizeType |
Definition of the size type. More... | |
typedef std::size_t | IndexType |
Definition of the index type. More... | |
typedef TSparseSpace::DataType | TDataType |
Definition of the data type. More... | |
typedef TSparseSpace::MatrixType | TSystemMatrixType |
Definition of the sparse matrix. More... | |
typedef TSparseSpace::VectorType | TSystemVectorType |
Definition of the vector size. More... | |
typedef TSparseSpace::MatrixPointerType | TSystemMatrixPointerType |
Definition of the pointer to the sparse matrix. More... | |
typedef TSparseSpace::VectorPointerType | TSystemVectorPointerType |
Definition of the pointer to the vector. More... | |
typedef TDenseSpace::MatrixType | LocalSystemMatrixType |
The local matrix definition. More... | |
typedef TDenseSpace::VectorType | LocalSystemVectorType |
The local vector definition. More... | |
typedef Scheme< TSparseSpace, TDenseSpace > | TSchemeType |
Definition of the scheme type. More... | |
typedef ModelPart::DofType | TDofType |
Definition of the DoF class. More... | |
typedef ModelPart::DofsArrayType | DofsArrayType |
Definition of the DoF array type. More... | |
typedef ModelPart::NodesContainerType | NodesArrayType |
The containers of the entities. More... | |
typedef ModelPart::ElementsContainerType | ElementsArrayType |
typedef ModelPart::ConditionsContainerType | ConditionsArrayType |
typedef PointerVectorSet< Element, IndexedObject > | ElementsContainerType |
The definition of the element container type. More... | |
![]() | |
static std::string | Name () |
Returns the name of the class as used in the settings (snake_case format) More... | |
![]() | |
static std::string | Name () |
Returns the name of the class as used in the settings (snake_case format) More... | |
![]() | |
double | mScaleFactor = 1.0 |
SCALING_DIAGONAL | mScalingDiagonal = SCALING_DIAGONAL::NO_SCALING |
The manually set scale factor. More... | |
![]() | |
TLinearSolver::Pointer | mpLinearSystemSolver = nullptr |
DofsArrayType | mDofSet |
Pointer to the linear solver. More... | |
bool | mReshapeMatrixFlag = false |
The set containing the DoF of the system. More... | |
bool | mDofSetIsInitialized = false |
If the matrix is reshaped each step. More... | |
bool | mCalculateReactionsFlag = false |
Flag taking care if the dof set was initialized ot not. More... | |
unsigned int | mEquationSystemSize |
Flag taking in account if it is needed or not to calculate the reactions. More... | |
int | mEchoLevel = 0 |
Number of degrees of freedom of the problem to be solve. More... | |
TSystemVectorPointerType | mpReactionsVector |
Current class provides an implementation for contact builder and solving operations. (elimination)
The RHS is constituted by the unbalanced loads (residual). Degrees of freedom are reordered putting the restrained degrees of freedom at the end of the system ordered in reverse order with respect to the DofSet and not considered the inactive ones. Imposition of the dirichlet conditions is naturally dealt with as the residual already contains this information. Calculation of the reactions involves a cost very similiar to the calculation of the total residual
TSparseSpace | The sparse matrix system considered |
TDenseSpace | The dense matrix system |
TLinearSolver | The type of linear solver considered |
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::BaseBuilderAndSolverType = BuilderAndSolver<TSparseSpace, TDenseSpace, TLinearSolver> |
Builder and solver base class.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::BaseType = ResidualBasedEliminationBuilderAndSolver<TSparseSpace, TDenseSpace, TLinearSolver> |
Definitions dependent on the base class.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::ClassType = ContactResidualBasedEliminationBuilderAndSolver<TSparseSpace, TDenseSpace, TLinearSolver> |
The definition of the current class.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::DofsArrayType = typename BaseType::DofsArrayType |
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::DofType = typename ModelPart::DofType |
The definition of the dof type.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::IndexSetType = std::unordered_set<IndexType> |
Index set definition.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::IndexType = std::size_t |
The index type.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::SizeType = std::size_t |
The size type.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TDataType = typename BaseType::TDataType |
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSchemeType = typename BaseType::TSchemeType |
Base types definitions.
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSystemMatrixType = typename BaseType::TSystemMatrixType |
using Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSystemVectorType = typename BaseType::TSystemVectorType |
|
inlineexplicit |
Default constructor.
|
inlineexplicit |
Default constructor. (with parameters)
|
inline |
Constructor.
|
inlineoverride |
Destructor.
|
inlineoverrideprotectedvirtual |
This method assigns settings to member variables.
ThisParameters | Parameters that are assigned to the member variables |
Reimplemented from Kratos::ResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
Create method.
pNewLinearSystemSolver | The linear solver for the system of equations |
ThisParameters | The configuration parameters |
Reimplemented from Kratos::ResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
This method provides the defaults parameters to avoid conflicts between the different constructors.
Reimplemented from Kratos::ResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
Kratos::ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::KRATOS_CLASS_POINTER_DEFINITION | ( | ContactResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > | ) |
Pointer definition of ContactResidualBasedEliminationBuilderAndSolver.
|
inlinestatic |
Returns the name of the class as used in the settings (snake_case format)
|
inlineoverridevirtual |
It organises the dofset in order to speed up the building phase.
rModelPart | The model part to compute |
Idem to the not contact version, except that if we fix the displacement in one slave node we should fix the corresponding LM for consistency
Reimplemented from Kratos::ResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.