Commit 994580ec authored by wenqing's avatar wenqing
Browse files

[HC] Initialize specific_body_force

parent 3619bbe4
......@@ -193,12 +193,12 @@ std::unique_ptr<Process> createComponentTransportProcess(
}
/// \section parameterscomponenttransport Process Parameters
// Specific body force parameter.
Eigen::VectorXd specific_body_force;
std::vector<double> const b =
//! \ogs_file_param{prj__processes__process__ComponentTransport__specific_body_force}
config.getConfigParameter<std::vector<double>>("specific_body_force");
assert(!b.empty() && b.size() < 4);
// Specific body force parameter.
Eigen::VectorXd specific_body_force(b.size());
int const mesh_space_dimension =
MeshLib::getSpaceDimension(mesh.getNodes());
if (static_cast<int>(b.size()) != mesh_space_dimension)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment