22 template <
int GlobalDim>
26 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
27 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
28 unsigned const integration_order,
29 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>&&
33 bool const use_monolithic_scheme)
34 :
Process(std::move(
name), mesh, std::move(jacobian_assembler), parameters,
35 integration_order, std::move(process_variables),
36 std::move(secondary_variables), use_monolithic_scheme),
37 _process_data(std::move(process_data))
41 template <
int GlobalDim>
45 _mesh_subset_all_nodes =
46 std::make_unique<MeshLib::MeshSubset>(_mesh, _mesh.getNodes());
49 _mesh_subset_base_nodes =
50 std::make_unique<MeshLib::MeshSubset>(_mesh, _base_nodes);
54 std::vector<MeshLib::MeshSubset> all_mesh_subsets_single_component{
55 *_mesh_subset_all_nodes};
56 _local_to_global_index_map_single_component =
57 std::make_unique<NumLib::LocalToGlobalIndexMap>(
58 std::move(all_mesh_subsets_single_component),
62 assert(_use_monolithic_scheme);
65 int const process_id = 0;
66 std::vector<MeshLib::MeshSubset> all_mesh_subsets(
67 getProcessVariables(process_id)[0]
69 .getNumberOfGlobalComponents(),
70 *_mesh_subset_all_nodes);
73 all_mesh_subsets.push_back(*_mesh_subset_base_nodes);
75 std::vector<int>
const vec_n_components{GlobalDim, 1};
77 _local_to_global_index_map =
78 std::make_unique<NumLib::LocalToGlobalIndexMap>(
79 std::move(all_mesh_subsets), vec_n_components,
81 assert(_local_to_global_index_map);
85 template <
int GlobalDim>
89 unsigned const integration_order)
91 int const process_id = 0;
102 _process_data.pressure_interpolated =
103 MeshLib::getOrCreateMeshProperty<double>(
108 template <
int GlobalDim>
111 assert(_use_monolithic_scheme);
113 int const process_id = 0;
114 initializeProcessBoundaryConditionsAndSourceTerms(
115 *_local_to_global_index_map, process_id);
119 template <
int GlobalDim>
124 assert(_use_monolithic_scheme);
126 auto const& l = *_local_to_global_index_map;
128 return {l.dofSizeWithoutGhosts(), l.dofSizeWithoutGhosts(),
129 &l.getGhostIndices(), &this->_sparsity_pattern};
133 template <
int GlobalDim>
135 const double t,
double const dt, std::vector<GlobalVector*>
const& x,
136 std::vector<GlobalVector*>
const& xdot,
int const process_id,
139 DBUG(
"Assemble StokesFlowProcess.");
141 std::vector<std::reference_wrapper<NumLib::LocalToGlobalIndexMap>>
143 assert(_use_monolithic_scheme);
145 dof_tables.push_back(std::ref(*_local_to_global_index_map));
156 template <
int GlobalDim>
158 const double ,
double const ,
159 std::vector<GlobalVector*>
const& ,
160 std::vector<GlobalVector*>
const& ,
const double ,
165 "Assembly of Jacobian matrix has not yet been implemented for "
166 "StokesFlowProcess.");
169 template <
int GlobalDim>
173 std::vector<GlobalVector*>
const& x,
175 int const process_id)
182 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
183 dof_tables.reserve(x.size());
184 assert(_use_monolithic_scheme);
186 dof_tables.push_back(_local_to_global_index_map.get());
196 template <
int GlobalDim>
198 std::vector<GlobalVector*>
const& x,
201 int const process_id)
208 std::vector<NumLib::LocalToGlobalIndexMap const*> dof_tables;
209 dof_tables.reserve(x.size());
210 assert(_use_monolithic_scheme);
212 dof_tables.push_back(_local_to_global_index_map.get());
221 template <
int GlobalDim>
225 assert(_use_monolithic_scheme);
227 return *_local_to_global_index_map;
void DBUG(char const *fmt, Args const &... args)
Global vector based on Eigen vector.
bool isAxiallySymmetric() const
unsigned getDimension() const
Returns the dimension of the mesh (determined by the maximum dimension over all elements).
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
virtual void postTimestep(std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, std::vector< GlobalVector * > const &x, double const t, double const dt)
virtual void computeSecondaryVariable(std::size_t const mesh_item_id, std::vector< NumLib::LocalToGlobalIndexMap const * > const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, GlobalVector const &x_dot, int const process_id)
unsigned getShapeFunctionOrder() const
std::vector< std::size_t > const & getActiveElementIDs() const
Handles configuration of several secondary variables from the project file.
StokesFlowProcess(std::string name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::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, StokesFlowProcessData &&process_data, SecondaryVariableCollection &&secondary_variables, bool const use_monolithic_scheme)
void initializeBoundaryConditions() override
void initializeConcreteProcess(NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, unsigned const integration_order) override
Process specific initialization called by initialize().
void computeSecondaryVariableConcrete(double const, double const, std::vector< GlobalVector * > const &x, GlobalVector const &, int const) override
void constructDofTable() override
void assembleConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b) override
void assembleWithJacobianConcreteProcess(const double t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, const double dxdot_dx, const double dx_dx, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b, GlobalMatrix &Jac) override
NumLib::LocalToGlobalIndexMap const & getDOFTable(const int) const override
void postTimestepConcreteProcess(std::vector< GlobalVector * > const &x, const double t, const double dt, int const process_id) override
MathLib::MatrixSpecifications getMatrixSpecifications(const int) const override
void assemble(std::size_t const mesh_item_id, LocalAssemblerInterface &local_assembler, std::vector< std::reference_wrapper< NumLib::LocalToGlobalIndexMap >> const &dof_tables, double const t, double const dt, std::vector< GlobalVector * > const &x, std::vector< GlobalVector * > const &xdot, int const process_id, GlobalMatrix &M, GlobalMatrix &K, GlobalVector &b)
std::vector< Node * > getBaseNodes(std::vector< Element * > const &elements)
@ BY_LOCATION
Ordering data by spatial location.
void createLocalAssemblers(const unsigned, std::vector< MeshLib::Element * > const &mesh_elements, NumLib::LocalToGlobalIndexMap const &dof_table, const unsigned shapefunction_order, std::vector< std::unique_ptr< LocalAssemblerInterface >> &local_assemblers, ExtraCtorArgs &&... extra_ctor_args)
static void executeSelectedMemberOnDereferenced(Method method, Container const &container, std::vector< std::size_t > const &active_container_ids, Args &&... args)
static void executeSelectedMemberDereferenced(Object &object, Method method, Container const &container, std::vector< std::size_t > const &active_container_ids, Args &&... args)