41 NumLib::ODESystemTag::FirstOrderImplicitQuasilinear,
42 NumLib::NonlinearSolverTag::Newton>,
50 std::unique_ptr<AbstractJacobianAssembler>&& jacobian_assembler,
51 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const&
53 unsigned const integration_order,
54 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
57 const bool use_monolithic_scheme =
true);
60 void preTimestep(std::vector<GlobalVector*>
const& x,
const double t,
61 const double delta_t,
const int process_id);
65 std::vector<GlobalVector*>
const& x_prev,
const double t,
66 const double delta_t,
int const process_id);
71 std::vector<GlobalVector*>
const& x_prev,
72 const double t,
double const dt,
73 int const process_id);
80 std::vector<GlobalVector*>
const& x,
82 int const process_id);
87 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const&
91 std::vector<GlobalVector*>& process_solutions,
92 std::vector<GlobalVector*>
const& process_solutions_prev,
94 int const process_id);
97 const int process_id)
const override;
108 std::vector<GlobalVector*>
const& ,
109 std::vector<GlobalVector*>& )
116 void assemble(
const double t,
double const dt,
117 std::vector<GlobalVector*>
const& x,
118 std::vector<GlobalVector*>
const& x_prev,
123 std::vector<GlobalVector*>
const& x,
124 std::vector<GlobalVector*>
const& x_prev,
135 void preOutput(
const double t,
double const dt,
136 std::vector<GlobalVector*>
const& x,
137 std::vector<GlobalVector*>
const& x_prev,
138 int const process_id);
140 std::vector<NumLib::IndexValueVector<GlobalIndexType>>
const*
142 int const process_id)
const final
154 std::vector<std::reference_wrapper<ProcessVariable>>
const&
170 std::vector<std::unique_ptr<MeshLib::IntegrationPointWriter>>
const&
182 std::vector<GlobalVector*>
const& )
const
184 return Eigen::Vector3d{};
188 std::vector<GlobalVector*>& ,
189 std::vector<GlobalVector*>
const& ,
double const ,
213 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const&
221 unsigned const integration_order) = 0;
226 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const&
230 std::vector<GlobalVector*>& ,
243 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
244 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
248 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
249 std::vector<GlobalVector*>
const& x_prev,
int const process_id,
254 std::vector<GlobalVector*>
const& ,
262 std::vector<GlobalVector*>
const& ,
263 std::vector<GlobalVector*>
const& ,
271 std::vector<GlobalVector*>
const& ,
272 std::vector<GlobalVector*>
const& ,
const double ,
273 double const ,
int const )
285 std::vector<GlobalVector*>
const& ,
298 const double ,
double const ,
299 std::vector<GlobalVector*>
const& ,
300 std::vector<GlobalVector*>
const& ,
int const )
325 const int specified_process_id);
335 virtual std::tuple<NumLib::LocalToGlobalIndexMap*, bool>
340 std::vector<GlobalIndexType>
376 std::vector<std::unique_ptr<MeshLib::IntegrationPointWriter>>
386 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>
MathLib::SparsityPattern< GlobalIndexType > GlobalSparsityPattern
Global vector based on Eigen vector.
std::vector< NumLib::IndexValueVector< GlobalIndexType > > const * getKnownSolutions(double const t, GlobalVector const &x, int const process_id) const final
std::vector< std::unique_ptr< MeshLib::IntegrationPointWriter > > const & getIntegrationPointWriters() const
MeshLib::Mesh & getMesh() const
virtual NumLib::IterationResult postIterationConcreteProcess(GlobalVector const &)
virtual void assembleWithJacobianConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac)=0
virtual Eigen::Vector3d getFlux(std::size_t, MathLib::Point3d const &, double const, std::vector< GlobalVector * > const &) const
virtual void preTimestepConcreteProcess(std::vector< GlobalVector * > const &, const double, const double, const int)
std::vector< BoundaryConditionCollection > _boundary_conditions
std::vector< std::unique_ptr< MeshLib::IntegrationPointWriter > > _integration_point_writer
SecondaryVariableCollection const & getSecondaryVariables() const
virtual void preIterationConcreteProcess(const unsigned, GlobalVector const &)
void constructDofTableOfSpecifiedProcessStaggeredScheme(const int specified_process_id)
static PROCESSLIB_EXPORT const std::string constant_one_parameter_name
void preAssemble(const double t, double const dt, GlobalVector const &x) final
virtual void initializeBoundaryConditions(std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
void initialize(std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
void initializeProcessBoundaryConditionsAndSourceTerms(const NumLib::LocalToGlobalIndexMap &dof_table, const int process_id, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium > > const &media)
virtual void solveReactionEquation(std::vector< GlobalVector * > &, std::vector< GlobalVector * > const &, double const, double const, NumLib::EquationSystem &, int const)
virtual void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order)=0
Process specific initialization called by initialize().
std::vector< std::size_t > _ids_of_active_elements
Union of active element ids per process variable.
virtual bool isMonolithicSchemeUsed() const
void assembleWithJacobian(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) final
void assemble(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) final
virtual void postTimestepConcreteProcess(std::vector< GlobalVector * > const &, std::vector< GlobalVector * > const &, const double, const double, int const)
Process(std::string name_, MeshLib::Mesh &mesh, std::unique_ptr< AbstractJacobianAssembler > &&jacobian_assembler, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const ¶meters, unsigned const integration_order, std::vector< std::vector< std::reference_wrapper< ProcessVariable > > > &&process_variables, SecondaryVariableCollection &&secondary_variables, const bool use_monolithic_scheme=true)
virtual void computeSecondaryVariableConcrete(double const, double const, std::vector< GlobalVector * > const &, GlobalVector const &, int const)
std::unique_ptr< MeshLib::MeshSubset const > _mesh_subset_all_nodes
void preOutput(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id)
virtual void assembleConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b)=0
std::vector< std::vector< std::reference_wrapper< ProcessVariable > > > _process_variables
std::vector< std::size_t > const & getActiveElementIDs() const
NumLib::LocalToGlobalIndexMap const & getSingleComponentDOFTable() const
virtual NumLib::LocalToGlobalIndexMap const & getDOFTable(const int) const
virtual void setInitialConditionsConcreteProcess(std::vector< GlobalVector * > &, double const, int const)
NumLib::IterationResult postIteration(GlobalVector const &x) final
void postTimestep(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, const double t, const double delta_t, int const process_id)
Postprocessing after a complete timestep.
SecondaryVariableCollection _secondary_variables
void preTimestep(std::vector< GlobalVector * > const &x, const double t, const double delta_t, const int process_id)
Preprocessing before starting assembly for new timestep.
void computeSecondaryVariable(double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_prev, int const process_id)
compute secondary variables for the coupled equations or for output.
void initializeExtrapolator()
virtual void postNonLinearSolverConcreteProcess(std::vector< GlobalVector * > const &, std::vector< GlobalVector * > const &, const double, double const, int const)
std::vector< std::reference_wrapper< ProcessVariable > > const & getProcessVariables(const int process_id) const
virtual void constructDofTable()
void constructMonolithicProcessDofTable()
VectorMatrixAssembler _global_assembler
virtual void preAssembleConcreteProcess(const double, double const, GlobalVector const &)
void postNonLinearSolver(std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &x_prev, const double t, double const dt, int const process_id)
std::vector< GlobalIndexType > getIndicesOfResiduumWithoutInitialCompensation() const override
unsigned const _integration_order
std::unique_ptr< NumLib::LocalToGlobalIndexMap > _local_to_global_index_map
std::vector< SourceTermCollection > _source_term_collections
virtual void extrapolateIntegrationPointValuesToNodes(const double, std::vector< GlobalVector * > const &, std::vector< GlobalVector * > &)
std::unique_ptr< ProcessLib::AbstractJacobianAssembler > _jacobian_assembler
void updateDeactivatedSubdomains(double const time, const int process_id)
void computeSparsityPattern()
MathLib::MatrixSpecifications getMatrixSpecifications(const int process_id) const override
void setInitialConditions(std::vector< GlobalVector * > &process_solutions, std::vector< GlobalVector * > const &process_solutions_prev, double const t, int const process_id)
virtual void preOutputConcreteProcess(const double, double const, std::vector< GlobalVector * > const &, std::vector< GlobalVector * > const &, int const)
ExtrapolatorData _extrapolator_data
NumLib::Extrapolator & getExtrapolator() const
GlobalSparsityPattern _sparsity_pattern
void preIteration(const unsigned iter, GlobalVector const &x) final
const bool _use_monolithic_scheme
virtual std::tuple< NumLib::LocalToGlobalIndexMap *, bool > getDOFTableForExtrapolatorData() const
Handles configuration of several secondary variables from the project file.
IterationResult
Status flags telling the NonlinearSolver if an iteration succeeded.
#define PROCESSLIB_EXPORT