diff --git a/Applications/ApplicationsLib/ProjectData.cpp b/Applications/ApplicationsLib/ProjectData.cpp
index 2b42a54ecb256e08b1e2a97bc0a639fa345b7e95..cd6cefcd4f6f25e2d91cad1598a9000fdeed195a 100644
--- a/Applications/ApplicationsLib/ProjectData.cpp
+++ b/Applications/ApplicationsLib/ProjectData.cpp
@@ -281,14 +281,13 @@ std::vector<GeoLib::NamedRaster> readRasters(
BaseLib::ConfigTree const& config, std::string const& raster_directory,
GeoLib::MinMaxPoints const& min_max_points)
{
- INFO("readRasters ... ");
+ INFO("readRasters ...");
std::vector<GeoLib::NamedRaster> named_rasters;
//! \ogs_file_param{prj__rasters}
auto optional_rasters = config.getConfigSubtreeOptional("rasters");
if (optional_rasters)
{
- INFO("Reading rasters.");
//! \ogs_file_param{prj__rasters__raster}
auto const configs = optional_rasters->getConfigSubtreeList("raster");
std::transform(
diff --git a/GeoLib/IO/NetCDFRasterReader.cpp b/GeoLib/IO/NetCDFRasterReader.cpp
index 28e87338fee4ccf363b3932c9e65383fdfea23fe..fe7f4ac6b6cb99161f0ec06366576a72b41b717d 100644
--- a/GeoLib/IO/NetCDFRasterReader.cpp
+++ b/GeoLib/IO/NetCDFRasterReader.cpp
@@ -78,7 +78,8 @@ std::vector<double> readDataFromNetCDF(
std::vector<double> raster_data;
for (auto [variable_name, variable] : variables)
{
- DBUG("checking variable '{}'", variable_name);
+ // uncomment for debugging
+ // DBUG("checking variable '{}'", variable_name);
if (variable_name == var_name)
{
auto const n_dims = variable.getDimCount();
@@ -86,14 +87,14 @@ std::vector<double> readDataFromNetCDF(
{
continue;
}
- DBUG("variable '{}' has {} dimensions", variable_name, n_dims);
+ // DBUG("variable '{}' has {} dimensions", variable_name, n_dims);
std::vector<std::size_t> counts(n_dims, 1);
for (int i = 0; i < n_dims; ++i)
{
counts[i] = variable.getDim(i).getSize();
}
- DBUG("counts {}", fmt::join(counts, ", "));
+ // DBUG("counts {}", fmt::join(counts, ", "));
std::vector<std::size_t> start(n_dims, 0);
if (dimension_number > counts[0])
{
@@ -109,9 +110,9 @@ std::vector<double> readDataFromNetCDF(
// y dimension - netcdf dimension number 1
// With the counts[1] and cell_size the 'usual' lower y-coordinate
// of the origin is computed.
- DBUG("reset header.origin[1]: original y0: {}, new y0: {} ",
- header.origin[1],
- header.origin[1] - counts[1] * header.cell_size);
+ // DBUG("reset header.origin[1]: original y0: {}, new y0: {} ",
+ // header.origin[1],
+ // header.origin[1] - counts[1] * header.cell_size);
header.origin[1] -= counts[1] * header.cell_size;
counts[2] = static_cast<int>(
@@ -129,10 +130,10 @@ std::vector<double> readDataFromNetCDF(
start[1] = static_cast<int>(std::floor(
(min_max_points.min[1] - header.origin[1]) / header.cell_size));
- DBUG("cut-out raster: pixel in x dimension: {} ", counts[2]);
- DBUG("cut-out raster: pixel in y dimension: {} ", counts[1]);
- DBUG("cut-out raster: start index x dimension: {}", start[2]);
- DBUG("cut-out raster: start index y dimension: {}", start[1]);
+ // DBUG("cut-out raster: pixel in x dimension: {} ", counts[2]);
+ // DBUG("cut-out raster: pixel in y dimension: {} ", counts[1]);
+ // DBUG("cut-out raster: start index x dimension: {}", start[2]);
+ // DBUG("cut-out raster: start index y dimension: {}", start[1]);
header.n_cols = counts[2];
header.n_rows = counts[1];
diff --git a/ProcessLib/BoundaryConditionAndSourceTerm/NeumannBoundaryConditionLocalAssembler.h b/ProcessLib/BoundaryConditionAndSourceTerm/NeumannBoundaryConditionLocalAssembler.h
index cfcd3f89f6e51c1e79395c20da790e7a677201cf..30b7628f92a6365645274377a1fd11b108098ac1 100644
--- a/ProcessLib/BoundaryConditionAndSourceTerm/NeumannBoundaryConditionLocalAssembler.h
+++ b/ProcessLib/BoundaryConditionAndSourceTerm/NeumannBoundaryConditionLocalAssembler.h
@@ -102,8 +102,8 @@ public:
}
else
{
- auto const value = _data.neumann_bc_parameter(t, position);
- _local_rhs.noalias() += N * value[0] * w * integral_measure;
+ auto const value = _data.neumann_bc_parameter(t, position)[0];
+ _local_rhs.noalias() += N * value * w * integral_measure;
}
}
diff --git a/ProcessLib/LargeDeformation/LargeDeformationFEM.h b/ProcessLib/LargeDeformation/LargeDeformationFEM.h
index 754f84a3941c12036c923252bbc3c2d015452f8f..b08af978e470707b5f91cd1cd44c787599daa00f 100644
--- a/ProcessLib/LargeDeformation/LargeDeformationFEM.h
+++ b/ProcessLib/LargeDeformation/LargeDeformationFEM.h
@@ -70,7 +70,6 @@ computeSigmaGeom(Eigen::Matrix3d const& sigma_tensor)
if constexpr (DisplacementDim == 3)
{
- // Assuming sigma_geom_op is defined for 3D
return sigma_geom_op(sigma_tensor);
}
}
diff --git a/ProcessLib/LargeDeformation/LocalAssemblerInterface.h b/ProcessLib/LargeDeformation/LocalAssemblerInterface.h
index 6057dbac7054acc1b3d373bdfc1ac9270a14b9ed..a01d8bee9cb84d5d24fcc0b911c6b77855a57022 100644
--- a/ProcessLib/LargeDeformation/LocalAssemblerInterface.h
+++ b/ProcessLib/LargeDeformation/LocalAssemblerInterface.h
@@ -78,7 +78,6 @@ struct LargeDeformationLocalAssemblerInterface
if (name.starts_with("material_state_variable_"))
{
std::string const variable_name = name.substr(24, name.size() - 24);
- DBUG("Setting material state variable '{:s}'", variable_name);
auto const& internal_variables =
solid_material_.getInternalVariables();
@@ -88,6 +87,7 @@ struct LargeDeformationLocalAssemblerInterface
{ return iv.name == variable_name; });
iv != end(internal_variables))
{
+ DBUG("Setting material state variable '{:s}'", variable_name);
return ProcessLib::
setIntegrationPointDataMaterialStateVariables(
values, material_states_,
diff --git a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
index 10a9b6efa5dce70b91fe52eb6376b83d77df7631..3c081fc93e9fdda7b56e1a7a34accff7fa5d14ad 100644
--- a/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
+++ b/ProcessLib/RichardsMechanics/RichardsMechanicsFEM-impl.h
@@ -251,7 +251,6 @@ std::size_t RichardsMechanicsLocalAssembler<
if (name.starts_with("material_state_variable_") && name.ends_with("_ip"))
{
std::string const variable_name = name.substr(24, name.size() - 24 - 3);
- DBUG("Setting material state variable '{:s}'", variable_name);
// Using first ip data for solid material. TODO (naumov) move solid
// material into element, store only material state in IPs.
@@ -263,6 +262,7 @@ std::size_t RichardsMechanicsLocalAssembler<
{ return iv.name == variable_name; });
iv != end(internal_variables))
{
+ DBUG("Setting material state variable '{:s}'", variable_name);
return ProcessLib::setIntegrationPointDataMaterialStateVariables(
values, _ip_data, &IpData::material_state_variables,
iv->reference);
diff --git a/ProcessLib/SmallDeformation/LocalAssemblerInterface.h b/ProcessLib/SmallDeformation/LocalAssemblerInterface.h
index b5bede0610ea78eeede1b706d0b9e3156cc5d6fc..2889e61e1c48ecf25e33d279dcc5ed6c434608c5 100644
--- a/ProcessLib/SmallDeformation/LocalAssemblerInterface.h
+++ b/ProcessLib/SmallDeformation/LocalAssemblerInterface.h
@@ -81,7 +81,6 @@ struct SmallDeformationLocalAssemblerInterface
if (name.starts_with("material_state_variable_"))
{
std::string const variable_name = name.substr(24, name.size() - 24);
- DBUG("Setting material state variable '{:s}'", variable_name);
auto const& internal_variables =
solid_material_.getInternalVariables();
@@ -91,6 +90,7 @@ struct SmallDeformationLocalAssemblerInterface
{ return iv.name == variable_name; });
iv != end(internal_variables))
{
+ DBUG("Setting material state variable '{:s}'", variable_name);
return ProcessLib::
setIntegrationPointDataMaterialStateVariables(
values, material_states_,
diff --git a/ProcessLib/ThermoRichardsMechanics/LocalAssemblerInterface.h b/ProcessLib/ThermoRichardsMechanics/LocalAssemblerInterface.h
index dff4f59be197aab84534f4d67a262f88bcfd1605..eefa35b0911e397a941ec94564242ad1533d59be 100644
--- a/ProcessLib/ThermoRichardsMechanics/LocalAssemblerInterface.h
+++ b/ProcessLib/ThermoRichardsMechanics/LocalAssemblerInterface.h
@@ -84,7 +84,6 @@ struct LocalAssemblerInterface : public ProcessLib::LocalAssemblerInterface,
if (name.starts_with("material_state_variable_"))
{
std::string const variable_name = name.substr(24, name.size() - 24);
- DBUG("Setting material state variable '{:s}'", variable_name);
auto const& internal_variables =
solid_material_.getInternalVariables();
@@ -94,6 +93,7 @@ struct LocalAssemblerInterface : public ProcessLib::LocalAssemblerInterface,
{ return iv.name == variable_name; });
iv != end(internal_variables))
{
+ DBUG("Setting material state variable '{:s}'", variable_name);
return ProcessLib::
setIntegrationPointDataMaterialStateVariables(
values, material_states_,
diff --git a/web/content/docs/benchmarks/liquid-flow/time-dependent-heterogeneous-source-term-and-boundary-conditions/index.md b/web/content/docs/benchmarks/liquid-flow/time-dependent-heterogeneous-source-term-and-boundary-conditions/index.md
index 3e49c1ebbb83d35c60193bdcf7b63e647d19aa34..9412195da81f4ce212ad59cf70412c03a2dbf70b 100644
--- a/web/content/docs/benchmarks/liquid-flow/time-dependent-heterogeneous-source-term-and-boundary-conditions/index.md
+++ b/web/content/docs/benchmarks/liquid-flow/time-dependent-heterogeneous-source-term-and-boundary-conditions/index.md
@@ -51,7 +51,7 @@ This simple example should demonstrate the use of the time dependent
heterogeneous parameter. We start with homogeneous parabolic problem:
$$
\begin{equation}
-s\\;\frac{\partial p}{\partial t} + k\; \Delta p = q(t,x) \quad \text{in }\Omega
+s\\;\frac{\partial p}{\partial t} + k\\; \Delta p = q(t,x) \quad \text{in }\Omega
\end{equation}
$$
w.r.t boundary conditions