![]() |
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.
|
#include <trilinos_elimination_builder_and_solver.h>
Public Member Functions | |
Life Cycle | |
TrilinosResidualBasedEliminationBuilderAndSolver (Epetra_MpiComm &Comm, int guess_row_size, typename TLinearSolver::Pointer pNewLinearSystemSolver) | |
virtual | ~TrilinosResidualBasedEliminationBuilderAndSolver () |
Operators | |
void | Build (typename TSchemeType::Pointer pScheme, ModelPart &r_model_part, TSystemMatrixType &A, TSystemVectorType &b) override |
equivalent (but generally faster) then performing BuildLHS and BuildRHS More... | |
void | BuildLHS (typename TSchemeType::Pointer pScheme, ModelPart &r_model_part, TSystemMatrixType &A) override |
Function to perform the building of the LHS, depending on the implementation chosen the size of the matrix could be equal to the total number of Dofs or to the number unrestrained dofs. More... | |
void | SystemSolveWithPhysics (TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b, ModelPart &rModelPart) |
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 &r_model_part, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) 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 &r_model_part, TSystemVectorType &b) 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 | SetUpDofSet (typename TSchemeType::Pointer pScheme, ModelPart &r_model_part) override |
Builds the list of the DofSets involved in the problem by "asking" to each element and condition its Dofs. More... | |
void | SetUpSystem (ModelPart &r_model_part) override |
It organises the dofset in order to speed up the building phase. More... | |
void | UpdateGhostDofs (ModelPart &rThisModelPart) |
void | ResizeAndInitializeVectors (typename TSchemeType::Pointer pScheme, TSystemMatrixPointerType &pA, TSystemVectorPointerType &pDx, TSystemVectorPointerType &pb, ModelPart &rModelPart) override |
This method initializes and resizes the system of equations. More... | |
void | InitializeSolutionStep (ModelPart &r_model_part, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) override |
It applies certain operations at the system of equations at the beginning of the solution step. More... | |
void | FinalizeSolutionStep (ModelPart &r_model_part, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) override |
It applies certain operations at the system of equations at the end of the solution step. More... | |
void | CalculateReactions (typename TSchemeType::Pointer pScheme, ModelPart &r_model_part, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) override |
It computes the reactions of the system. More... | |
void | BuildLHS_CompleteOnFreeRows (typename TSchemeType::Pointer pScheme, ModelPart &r_model_part, TSystemMatrixType &A) override |
It builds 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 | ApplyDirichletConditions (typename TSchemeType::Pointer pScheme, ModelPart &r_model_part, TSystemMatrixType &A, TSystemVectorType &Dx, TSystemVectorType &b) override |
It 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... | |
![]() | |
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 () |
virtual ClassType::Pointer | Create (typename TLinearSolver::Pointer pNewLinearSystemSolver, Parameters ThisParameters) const |
Create method. More... | |
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... | |
virtual std::string | Info () const |
Turn back information as a string. More... | |
virtual void | PrintInfo (std::ostream &rOStream) const |
Print information about this object. More... | |
virtual void | PrintData (std::ostream &rOStream) const |
Print object's data. 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 | SystemSolve (TSystemMatrixType &rA, TSystemVectorType &rDx, TSystemVectorType &rb) |
This is a call to the linear system solver. 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 | Clear () |
This function is intended to be called at the end of the solution step to clean up memory storage not needed. More... | |
virtual int | Check (ModelPart &rModelPart) |
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... | |
virtual Parameters | GetDefaultParameters () const |
This method provides the defaults parameters to avoid conflicts between the different constructors. More... | |
Protected Attributes | |
Protected member Variables | |
Epetra_MpiComm & | mrComm |
int | mguess_row_size |
int | mFirstMyId |
int | mLastMyId |
![]() | |
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 |
Type Definitions | |
typedef BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > | BaseType |
typedef TSparseSpace | SparseSpaceType |
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 BaseType::NodesArrayType | NodesArrayType |
typedef BaseType::ElementsArrayType | ElementsArrayType |
typedef BaseType::ConditionsArrayType | ConditionsArrayType |
typedef BaseType::ElementsContainerType | ElementsContainerType |
KRATOS_CLASS_POINTER_DEFINITION (TrilinosResidualBasedEliminationBuilderAndSolver) | |
Additional Inherited Members | |
![]() | |
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... | |
![]() | |
virtual Parameters | ValidateAndAssignParameters (Parameters ThisParameters, const Parameters DefaultParameters) const |
This method validate and assign default parameters. More... | |
virtual void | AssignSettings (const Parameters ThisParameters) |
This method assigns settings to member variables. More... | |
Short class definition.
Detail class definition.
Current class provides an implementation for standard builder and solving operations.
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.
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
\URL[Example of use html]{ extended_documentation/no_ex_of_use.html}
\URL[Example of use pdf]{ extended_documentation/no_ex_of_use.pdf}
\URL[Example of use doc]{ extended_documentation/no_ex_of_use.doc}
\URL[Example of use ps]{ extended_documentation/no_ex_of_use.ps}
\URL[Extended documentation html]{ extended_documentation/no_ext_doc.html}
\URL[Extended documentation pdf]{ extended_documentation/no_ext_doc.pdf}
\URL[Extended documentation doc]{ extended_documentation/no_ext_doc.doc}
\URL[Extended documentation ps]{ extended_documentation/no_ext_doc.ps}
typedef BuilderAndSolver<TSparseSpace,TDenseSpace, TLinearSolver > Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::BaseType |
typedef BaseType::ConditionsArrayType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::ConditionsArrayType |
typedef BaseType::DofsArrayType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::DofsArrayType |
typedef BaseType::ElementsArrayType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::ElementsArrayType |
typedef BaseType::ElementsContainerType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::ElementsContainerType |
typedef BaseType::LocalSystemMatrixType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::LocalSystemMatrixType |
typedef BaseType::LocalSystemVectorType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::LocalSystemVectorType |
typedef BaseType::NodesArrayType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::NodesArrayType |
typedef TSparseSpace Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::SparseSpaceType |
typedef BaseType::TDataType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TDataType |
typedef BaseType::TSchemeType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSchemeType |
typedef BaseType::TSystemMatrixPointerType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSystemMatrixPointerType |
typedef BaseType::TSystemMatrixType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSystemMatrixType |
typedef BaseType::TSystemVectorPointerType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSystemVectorPointerType |
typedef BaseType::TSystemVectorType Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::TSystemVectorType |
|
inline |
Constructor.
|
inlinevirtual |
Destructor.
|
inlineoverridevirtual |
It 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.
For explanation of how it works for a particular implementation the user should refer to the particular Builder And Solver chosen
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
rDx | The vector of unknowns |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
equivalent (but generally faster) then performing BuildLHS and BuildRHS
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
Function to perform the building and solving phase at the same time.
It is ideally the fastest and safer function to use when it is possible to solve just after building
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
rDx | The vector of unknowns |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
Function to perform the building of the LHS, depending on the implementation chosen the size of the matrix could be equal to the total number of Dofs or to the number unrestrained dofs.
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
It builds 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.
This matrix is obtained by building the total matrix without the lines corresponding to the fixed degrees of freedom (but keeping the columns!!)
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
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.
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
Corresponds to the previews, but the System's matrix is considered already built and only the RHS is built again.
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
rDx | The vector of unknowns |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
It computes the reactions of the system.
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
rDx | The vector of unknowns |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
It applies certain operations at the system of equations at the end of the solution step.
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
rDx | The vector of unknowns |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
It applies certain operations at the system of equations at the beginning of the solution step.
rModelPart | The model part to compute |
rA | The LHS matrix of the system of equations |
rDx | The vector of unknowns |
rb | The RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
Kratos::TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >::KRATOS_CLASS_POINTER_DEFINITION | ( | TrilinosResidualBasedEliminationBuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver > | ) |
|
inlineoverridevirtual |
This method initializes and resizes the system of equations.
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
pA | The pointer to LHS matrix of the system of equations |
pDx | The pointer to vector of unknowns |
pb | The pointer to RHS vector of the system of equations |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
Builds the list of the DofSets involved in the problem by "asking" to each element and condition its Dofs.
The list of dofs is stores inside the BuilderAndSolver as it is closely connected to the way the matrix and RHS are built
pScheme | The pointer to the integration scheme |
rModelPart | The model part to compute |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inlineoverridevirtual |
It organises the dofset in order to speed up the building phase.
rModelPart | The model part to compute |
Reimplemented from Kratos::BuilderAndSolver< TSparseSpace, TDenseSpace, TLinearSolver >.
|
inline |
|
inline |
|
protected |
|
protected |
|
protected |
|
protected |