diff --git a/Applications/ApplicationsLib/ProjectData.cpp b/Applications/ApplicationsLib/ProjectData.cpp index 2b42a54ecb2..cd6cefcd4f6 100644 --- a/Applications/ApplicationsLib/ProjectData.cpp +++ b/Applications/ApplicationsLib/ProjectData.cpp @@ -281,14 +281,13 @@ std::vector readRasters( BaseLib::ConfigTree const& config, std::string const& raster_directory, GeoLib::MinMaxPoints const& min_max_points) { - INFO("readRasters ... "); + INFO("readRasters ..."); std::vector 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 28e87338fee..fe7f4ac6b6c 100644 --- a/GeoLib/IO/NetCDFRasterReader.cpp +++ b/GeoLib/IO/NetCDFRasterReader.cpp @@ -78,7 +78,8 @@ std::vector readDataFromNetCDF( std::vector 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 readDataFromNetCDF( { continue; } - DBUG("variable '{}' has {} dimensions", variable_name, n_dims); + // DBUG("variable '{}' has {} dimensions", variable_name, n_dims); std::vector 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 start(n_dims, 0); if (dimension_number > counts[0]) { @@ -109,9 +110,9 @@ std::vector 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( @@ -129,10 +130,10 @@ std::vector readDataFromNetCDF( start[1] = static_cast(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 cfcd3f89f6e..30b7628f92a 100644 --- a/ProcessLib/BoundaryConditionAndSourceTerm/NeumannBoundaryConditionLocalAssembler.h +++ b/ProcessLib/BoundaryConditionAndSourceTerm/NeumannBoundaryConditionLocalAssembler.h @@ -102,8 +102,8 @@ class NeumannBoundaryConditionLocalAssembler final } 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 754f84a3941..b08af978e47 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 6057dbac705..a01d8bee9cb 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 10a9b6efa5d..3c081fc93e9 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 b5bede0610e..2889e61e1c4 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 dff4f59be19..eefa35b0911 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 3e49c1ebbb8..9412195da81 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