diff --git a/Applications/ApplicationsLib/ProjectData.cpp b/Applications/ApplicationsLib/ProjectData.cpp
index 147c4db819f19dbb6cd516cec5e472360b009e8f..492ef56ef856c2feff34f4c8ac21acf4b0ee4ab7 100644
--- a/Applications/ApplicationsLib/ProjectData.cpp
+++ b/Applications/ApplicationsLib/ProjectData.cpp
@@ -265,16 +265,7 @@ void ProjectData::parseProcessVariables(
BaseLib::ConfigTree const& process_variables_config)
{
DBUG("Parse process variables:")
- if (_geoObjects == nullptr)
- {
- ERR("Geometric objects are required to define process variables.");
- ERR("No geometric objects present.");
- return;
- }
- // TODO at the moment we have only one mesh, later there
- // can be several meshes. Then we have to check for correct mesh here and
- // assign the referenced mesh below.
if (_mesh_vec.empty() || _mesh_vec[0] == nullptr)
{
ERR("A mesh is required to define process variables.");
@@ -287,9 +278,8 @@ void ProjectData::parseProcessVariables(
//! \ogs_file_param{prj__process_variables__process_variable}
: process_variables_config.getConfigSubtreeList("process_variable"))
{
- // TODO Extend to referenced meshes.
- auto pv = ProcessLib::ProcessVariable{var_config, *_mesh_vec[0],
- *_geoObjects, _parameters};
+ auto pv =
+ ProcessLib::ProcessVariable{var_config, _mesh_vec, _parameters};
if (!names.insert(pv.getName()).second)
OGS_FATAL("A process variable with name `%s' already exists.",
pv.getName().c_str());
diff --git a/ProcessLib/ProcessVariable.cpp b/ProcessLib/ProcessVariable.cpp
index 48d45cd23714ee5b1fe49743e6610847ae032155..5407654b45109677bf05e2af33b1d729885ef5fd 100644
--- a/ProcessLib/ProcessVariable.cpp
+++ b/ProcessLib/ProcessVariable.cpp
@@ -12,19 +12,19 @@
#include <utility>
#include <logog/include/logog.hpp>
-#include "GeoLib/GEOObjects.h"
#include "MeshLib/Mesh.h"
#include "ProcessLib/Utils/ProcessUtils.h"
namespace ProcessLib
{
ProcessVariable::ProcessVariable(
- BaseLib::ConfigTree const& config, MeshLib::Mesh& mesh,
- GeoLib::GEOObjects const& geometries,
+ BaseLib::ConfigTree const& config,
+ std::vector<MeshLib::Mesh*> const& meshes,
std::vector<std::unique_ptr<ParameterBase>> const& parameters)
: //! \ogs_file_param{prj__process_variables__process_variable__name}
_name(config.getConfigParameter<std::string>("name")),
- _mesh(mesh),
+ _mesh(*meshes[0]), // Using the first mesh as the main mesh.
+ // TODO (naumov) potentially extend to named meshes.
//! \ogs_file_param{prj__process_variables__process_variable__components}
_n_components(config.getConfigParameter<int>("components")),
//! \ogs_file_param{prj__process_variables__process_variable__order}
@@ -56,18 +56,23 @@ ProcessVariable::ProcessVariable(
//! \ogs_file_param{prj__process_variables__process_variable__boundary_conditions__boundary_condition__geometry}
bc_config.getConfigParameter<std::string>("geometry");
- GeoLib::GeoObject const* const geometry =
- geometries.getGeoObject(geometrical_set_name, geometry_name);
-
- if (! geometry)
- OGS_FATAL(
- "No geometry with name `%s' has been found in the "
- "geometrical set `%s'.",
- geometry_name.c_str(), geometrical_set_name.c_str());
-
- DBUG(
- "Found geometry type \"%s\"",
- GeoLib::convertGeoTypeToString(geometry->getGeoType()).c_str());
+ auto const full_geometry_name =
+ geometrical_set_name + "_" + geometry_name;
+ auto const mesh_it =
+ std::find_if(begin(meshes), end(meshes),
+ [&full_geometry_name](MeshLib::Mesh* const mesh) {
+ assert(mesh != nullptr);
+ return mesh->getName() == full_geometry_name;
+ });
+ if (mesh_it == end(meshes))
+ {
+ OGS_FATAL("Required mesh with name '%s' not found.",
+ full_geometry_name.c_str());
+ }
+ MeshLib::Mesh& bc_mesh = **mesh_it;
+
+ DBUG("Found mesh '%s' with id %d.", bc_mesh.getName().c_str(),
+ bc_mesh.getID());
auto component_id =
//! \ogs_file_param{prj__process_variables__process_variable__boundary_conditions__boundary_condition__component}
@@ -77,7 +82,7 @@ ProcessVariable::ProcessVariable(
// default value for single component vars.
component_id = 0;
- _bc_configs.emplace_back(std::move(bc_config), *geometry,
+ _bc_configs.emplace_back(std::move(bc_config), bc_mesh,
component_id);
}
} else {
@@ -92,6 +97,7 @@ ProcessVariable::ProcessVariable(
//! \ogs_file_param{prj__process_variables__process_variable__source_terms__source_term}
sts_config->getConfigSubtreeList("source_term"))
{
+ // TODO (naumov) Remove code duplication with the bc_config parsing.
auto const geometrical_set_name =
//! \ogs_file_param{prj__process_variables__process_variable__source_terms__source_term__geometrical_set}
st_config.getConfigParameter<std::string>("geometrical_set");
@@ -99,18 +105,23 @@ ProcessVariable::ProcessVariable(
//! \ogs_file_param{prj__process_variables__process_variable__source_terms__source_term__geometry}
st_config.getConfigParameter<std::string>("geometry");
- GeoLib::GeoObject const* const geometry =
- geometries.getGeoObject(geometrical_set_name, geometry_name);
-
- if (! geometry)
- OGS_FATAL(
- "No geometry with name `%s' has been found in the "
- "geometrical set `%s'.",
- geometry_name.c_str(), geometrical_set_name.c_str());
-
- DBUG(
- "Found geometry type \"%s\"",
- GeoLib::convertGeoTypeToString(geometry->getGeoType()).c_str());
+ auto const full_geometry_name =
+ geometrical_set_name + "_" + geometry_name;
+ auto const mesh_it =
+ std::find_if(begin(meshes), end(meshes),
+ [&full_geometry_name](MeshLib::Mesh* const mesh) {
+ assert(mesh != nullptr);
+ return mesh->getName() == full_geometry_name;
+ });
+ if (mesh_it == end(meshes))
+ {
+ OGS_FATAL("Required mesh with name '%s' not found.",
+ full_geometry_name.c_str());
+ }
+ MeshLib::Mesh& bc_mesh = **mesh_it;
+
+ DBUG("Found mesh '%s' with id %d.", bc_mesh.getName().c_str(),
+ bc_mesh.getID());
auto component_id =
//! \ogs_file_param{prj__process_variables__process_variable__source_terms__source_term__component}
@@ -120,7 +131,7 @@ ProcessVariable::ProcessVariable(
// default value for single component vars.
component_id = 0;
- _source_term_configs.emplace_back(std::move(st_config), *geometry,
+ _source_term_configs.emplace_back(std::move(st_config), bc_mesh,
component_id);
}
} else {
diff --git a/ProcessLib/ProcessVariable.h b/ProcessLib/ProcessVariable.h
index 4a33f59bb18346915eea1cf62369e8ce6a4541eb..11343c425d81274966bab713b81b6c3c7adfaa75 100644
--- a/ProcessLib/ProcessVariable.h
+++ b/ProcessLib/ProcessVariable.h
@@ -31,8 +31,8 @@ class ProcessVariable
{
public:
ProcessVariable(
- BaseLib::ConfigTree const& config, MeshLib::Mesh& mesh,
- GeoLib::GEOObjects const& geometries,
+ BaseLib::ConfigTree const& config,
+ std::vector<MeshLib::Mesh*> const& meshes,
std::vector<std::unique_ptr<ParameterBase>> const& parameters);
ProcessVariable(ProcessVariable&&);