24 namespace RichardsFlow
30 std::array
const required_medium_properties = {
38 std::array
const required_liquid_properties = {
42 std::array<MaterialPropertyLib::PropertyType, 0>
const
43 required_solid_properties{};
46 mesh, media_map, required_medium_properties, required_solid_properties,
47 required_liquid_properties);
53 std::unique_ptr<ProcessLib::AbstractJacobianAssembler>&& jacobian_assembler,
54 std::vector<ProcessVariable>
const& variables,
55 std::vector<std::unique_ptr<ParameterLib::ParameterBase>>
const& parameters,
56 unsigned const integration_order,
59 std::unique_ptr<MathLib::PiecewiseLinearInterpolation>>
const&
61 std::map<
int, std::shared_ptr<MaterialPropertyLib::Medium>>
const& media)
66 DBUG(
"Create RichardsFlowProcess.");
77 std::vector<std::vector<std::reference_wrapper<ProcessVariable>>>
79 process_variables.push_back(std::move(per_process_variables));
86 std::vector<double>
const b =
89 assert(!b.empty() && b.size() < 4);
90 Eigen::VectorXd specific_body_force(b.size());
94 std::copy_n(b.data(), b.size(), specific_body_force.data());
104 DBUG(
"Check the media properties of RichardsFlow process ...");
106 DBUG(
"Media properties verified.");
109 std::move(media_map), specific_body_force, has_gravity, mass_lumping};
111 return std::make_unique<RichardsFlowProcess>(
112 std::move(
name), mesh, std::move(jacobian_assembler), parameters,
113 integration_order, std::move(process_variables),
114 std::move(process_data), std::move(secondary_variables));
void DBUG(char const *fmt, Args const &... args)
void checkConfigParameter(std::string const ¶m, T const &value) const
T getConfigParameter(std::string const ¶m) const
ConfigTree getConfigSubtree(std::string const &root) const
Handles configuration of several secondary variables from the project file.
std::unique_ptr< MaterialSpatialDistributionMap > createMaterialSpatialDistributionMap(std::map< int, std::shared_ptr< Medium >> const &media, MeshLib::Mesh const &mesh)
void checkMaterialSpatialDistributionMap(MeshLib::Mesh const &mesh, MaterialPropertyLib::MaterialSpatialDistributionMap const &media_map, ContainerMedium const &required_properties_medium, ContainerSolid const &required_properties_solid_phase, ContainerLiquid const &required_properties_liquid_phase)
Eigen::Map< const Vector > toVector(std::vector< double > const &data, Eigen::VectorXd::Index size)
Creates an Eigen mapped vector from the given data vector.
std::unique_ptr< Process > createRichardsFlowProcess(std::string name, MeshLib::Mesh &mesh, std::unique_ptr< ProcessLib::AbstractJacobianAssembler > &&jacobian_assembler, std::vector< ProcessVariable > const &variables, std::vector< std::unique_ptr< ParameterLib::ParameterBase >> const ¶meters, unsigned const integration_order, BaseLib::ConfigTree const &config, std::map< std::string, std::unique_ptr< MathLib::PiecewiseLinearInterpolation >> const &, std::map< int, std::shared_ptr< MaterialPropertyLib::Medium >> const &media)
void checkMPLProperties(MeshLib::Mesh const &mesh, MaterialPropertyLib::MaterialSpatialDistributionMap const &media_map)
std::vector< std::reference_wrapper< ProcessVariable > > findProcessVariables(std::vector< ProcessVariable > const &variables, BaseLib::ConfigTree const &pv_config, std::initializer_list< std::string > tags)
void createSecondaryVariables(BaseLib::ConfigTree const &config, SecondaryVariableCollection &secondary_variables)