diff --git a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
index 58b3ccd1dd171c6856e0a30518bccad445c16cf8..a4bede7634e26cdbdb8b6fe7d4dad821de08d341 100644
--- a/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
+++ b/ProcessLib/HydroMechanics/CreateHydroMechanicsProcess.cpp
@@ -55,8 +55,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
DBUG("Associate displacement with process variable \'%s\'.",
process_variables[1].get().getName().c_str());
- if (process_variables[1].get().getNumberOfComponents() !=
- DisplacementDim)
+ if (process_variables[1].get().getNumberOfComponents() != DisplacementDim)
{
OGS_FATAL(
"Number of components of the process variable '%s' is different "
@@ -77,7 +76,6 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
process_variables[0].get().getNumberOfComponents());
}
-
// Constitutive relation.
// read type;
auto const constitutive_relation_config =
@@ -107,8 +105,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
auto& intrinsic_permeability = findParameter<double>(
config,
//! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__intrinsic_permeability}
- "intrinsic_permeability",
- parameters, 1);
+ "intrinsic_permeability", parameters, 1);
DBUG("Use \'%s\' as intrinsic conductivity parameter.",
intrinsic_permeability.name.c_str());
@@ -126,8 +123,7 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
auto& fluid_viscosity = findParameter<double>(
config,
//! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__fluid_viscosity}
- "fluid_viscosity",
- parameters, 1);
+ "fluid_viscosity", parameters, 1);
DBUG("Use \'%s\' as fluid viscosity parameter.",
fluid_viscosity.name.c_str());
@@ -135,17 +131,14 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
auto& fluid_density = findParameter<double>(
config,
//! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__fluid_density}
- "fluid_density",
- parameters, 1);
- DBUG("Use \'%s\' as fluid density parameter.",
- fluid_density.name.c_str());
+ "fluid_density", parameters, 1);
+ DBUG("Use \'%s\' as fluid density parameter.", fluid_density.name.c_str());
// Biot coefficient
auto& biot_coefficient = findParameter<double>(
config,
//! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__biot_coefficient}
- "biot_coefficient",
- parameters, 1);
+ "biot_coefficient", parameters, 1);
DBUG("Use \'%s\' as Biot coefficient parameter.",
biot_coefficient.name.c_str());
@@ -153,19 +146,15 @@ std::unique_ptr<Process> createHydroMechanicsProcess(
auto& porosity = findParameter<double>(
config,
//! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__porosity}
- "porosity",
- parameters, 1);
- DBUG("Use \'%s\' as porosity parameter.",
- porosity.name.c_str());
+ "porosity", parameters, 1);
+ DBUG("Use \'%s\' as porosity parameter.", porosity.name.c_str());
// Solid density
auto& solid_density = findParameter<double>(
config,
//! \ogs_file_param_special{prj__processes__process__HYDRO_MECHANICS__solid_density}
- "solid_density",
- parameters, 1);
- DBUG("Use \'%s\' as solid density parameter.",
- solid_density.name.c_str());
+ "solid_density", parameters, 1);
+ DBUG("Use \'%s\' as solid density parameter.", solid_density.name.c_str());
// Specific body force
Eigen::Matrix<double, DisplacementDim, 1> specific_body_force;
diff --git a/ProcessLib/HydroMechanics/CreateLocalAssemblers.h b/ProcessLib/HydroMechanics/CreateLocalAssemblers.h
index 305e574426937de343d87067dae3393120309eeb..3201b94494c067948a253cce9485d1517eca9e0a 100644
--- a/ProcessLib/HydroMechanics/CreateLocalAssemblers.h
+++ b/ProcessLib/HydroMechanics/CreateLocalAssemblers.h
@@ -64,8 +64,9 @@ void createLocalAssemblers(
* The first two template parameters cannot be deduced from the arguments.
* Therefore they always have to be provided manually.
*/
-template <int DisplacementDim, template <typename, typename, typename, unsigned, int>
- class LocalAssemblerImplementation,
+template <int DisplacementDim,
+ template <typename, typename, typename, unsigned, int>
+ class LocalAssemblerImplementation,
typename LocalAssemblerInterface, typename... ExtraCtorArgs>
void createLocalAssemblers(
const unsigned dimension,
diff --git a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
index 0cb5c7bcb033b84fdd248bc1f17ead10369d8bd6..d49106882831eef9a65d19d5c94e80c0f6caefe5 100644
--- a/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
+++ b/ProcessLib/HydroMechanics/HydroMechanicsFEM.h
@@ -61,7 +61,7 @@ struct IntegrationPointData final
typename ShapeMatrixTypeDisplacement::template MatrixType<
DisplacementDim, NPoints * DisplacementDim>
N_u;
- //typename ShapeMatrixTypeDisplacement::NodalRowVectorType N_u;
+ // typename ShapeMatrixTypeDisplacement::NodalRowVectorType N_u;
typename BMatricesType::BMatrixType b_matrices;
typename BMatricesType::KelvinVectorType sigma_eff, sigma_eff_prev;
typename BMatricesType::KelvinVectorType eps, eps_prev;
@@ -207,8 +207,8 @@ public:
auto& ip_data = _ip_data[ip];
auto const& sm = shape_matrices_u[ip];
_ip_data[ip].integration_weight =
- _integration_method.getWeightedPoint(ip).getWeight() * sm.integralMeasure *
- shape_matrices_u[ip].detJ;
+ _integration_method.getWeightedPoint(ip).getWeight() *
+ sm.integralMeasure * shape_matrices_u[ip].detJ;
ip_data.b_matrices.resize(
kelvin_vector_size,
ShapeFunctionDisplacement::NPOINTS * DisplacementDim);
@@ -327,8 +327,7 @@ public:
auto const& C = _ip_data[ip].C;
- double const S =
- _process_data.specific_storage(t, x_position)[0];
+ double const S = _process_data.specific_storage(t, x_position)[0];
double const K_over_mu =
_process_data.intrinsic_permeability(t, x_position)[0] /
_process_data.fluid_viscosity(t, x_position)[0];
@@ -394,8 +393,7 @@ public:
.noalias() += Kup.transpose() / dt;
// pressure equation
- local_rhs.template segment<pressure_size>(pressure_index)
- .noalias() -=
+ local_rhs.template segment<pressure_size>(pressure_index).noalias() -=
laplace_p * p + storage_p * p_dot + Kup.transpose() * u_dot;
// displacement equation
@@ -604,7 +602,8 @@ private:
IntegrationMethod _integration_method;
MeshLib::Element const& _element;
- SecondaryData<typename ShapeMatricesTypeDisplacement::ShapeMatrices::ShapeType>
+ SecondaryData<
+ typename ShapeMatricesTypeDisplacement::ShapeMatrices::ShapeType>
_secondary_data;
std::vector<std::vector<double>> _darcy_velocities =
@@ -612,7 +611,6 @@ private:
DisplacementDim,
std::vector<double>(_integration_method.getNumberOfPoints()));
-
static const int pressure_index = 0;
static const int pressure_size = ShapeFunctionPressure::NPOINTS;
static const int displacement_index = ShapeFunctionPressure::NPOINTS;
diff --git a/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp b/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp
index 5882b72e349f2a0168d6638d121d67c3c89146e0..12614d40370c3eb923b5988244f79f9cd2aa7306 100644
--- a/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp
+++ b/ProcessLib/HydroMechanics/HydroMechanicsProcess.cpp
@@ -7,16 +7,15 @@
*
*/
-#include "HydroMechanicsProcess-fwd.h"
#include "HydroMechanicsProcess.h"
+#include "HydroMechanicsProcess-fwd.h"
namespace ProcessLib
{
namespace HydroMechanics
{
-
template class HydroMechanicsProcess<2>;
template class HydroMechanicsProcess<3>;
-} // namespace HydroMechanics
-} // namespace ProcessLib
+} // namespace HydroMechanics
+} // namespace ProcessLib
diff --git a/ProcessLib/HydroMechanics/HydroMechanicsProcess.h b/ProcessLib/HydroMechanics/HydroMechanicsProcess.h
index 45e3db37db7e6b4f39344c968b05c061afd97d15..4496c7fa41c5e4ef114593806b0d74cab0ca0eb2 100644
--- a/ProcessLib/HydroMechanics/HydroMechanicsProcess.h
+++ b/ProcessLib/HydroMechanics/HydroMechanicsProcess.h
@@ -139,7 +139,8 @@ private:
getExtrapolator(), _local_assemblers,
&HydroMechanicsLocalAssemblerInterface::getIntPtSigmaXY));
- if (DisplacementDim == 3) {
+ if (DisplacementDim == 3)
+ {
Base::_secondary_variables.addSecondaryVariable(
"sigma_xz", 1,
makeExtrapolator(
@@ -196,11 +197,9 @@ private:
getIntPtDarcyVelocityZ));
}
- void assembleConcreteProcess(const double t, GlobalVector const& x,
- GlobalMatrix& M, GlobalMatrix& K,
- GlobalVector& b,
- StaggeredCouplingTerm const&
- coupling_term) override
+ void assembleConcreteProcess(
+ const double t, GlobalVector const& x, GlobalMatrix& M, GlobalMatrix& K,
+ GlobalVector& b, StaggeredCouplingTerm const& coupling_term) override
{
DBUG("Assemble HydroMechanicsProcess.");
diff --git a/ProcessLib/HydroMechanics/LocalDataInitializer.h b/ProcessLib/HydroMechanics/LocalDataInitializer.h
index 2869c4cd017b5ddeaf264ee926d8fa0662f80781..d2b68e03be5c71678234a36c58979cb49286f3b9 100644
--- a/ProcessLib/HydroMechanics/LocalDataInitializer.h
+++ b/ProcessLib/HydroMechanics/LocalDataInitializer.h
@@ -11,9 +11,9 @@
#include <functional>
#include <memory>
+#include <type_traits>
#include <typeindex>
#include <typeinfo>
-#include <type_traits>
#include <unordered_map>
#include "MeshLib/Elements/Elements.h"