From abf22cd4abb61fab0e2baf5c66a2bc085ec0a295 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Tue, 8 Oct 2019 15:23:55 +0200
Subject: [PATCH] [PL] CT; Remove unneeded getFlux implementation.
There is another one, which is used.
---
.../ComponentTransportFEM.h | 28 ++-----------------
1 file changed, 3 insertions(+), 25 deletions(-)
diff --git a/ProcessLib/ComponentTransport/ComponentTransportFEM.h b/ProcessLib/ComponentTransport/ComponentTransportFEM.h
index 4e19d84b500..59d6f9df50e 100644
--- a/ProcessLib/ComponentTransport/ComponentTransportFEM.h
+++ b/ProcessLib/ComponentTransport/ComponentTransportFEM.h
@@ -826,27 +826,6 @@ public:
auto const local_C = Eigen::Map<const NodalVectorType>(
&local_x[first_concentration_index], concentration_size);
- return calculateFlux(pnt_local_coords, t, local_p, local_C);
- }
-
- Eigen::Vector3d getFlux(
- MathLib::Point3d const& pnt_local_coords,
- double const t,
- std::vector<std::vector<double>> const& local_xs) const override
- {
- auto const local_p = MathLib::toVector(local_xs[hydraulic_process_id]);
- auto const local_C =
- MathLib::toVector(local_xs[first_transport_process_id]);
-
- return calculateFlux(pnt_local_coords, t, local_p, local_C);
- }
-
- Eigen::Vector3d calculateFlux(
- MathLib::Point3d const& pnt_local_coords,
- double const t,
- Eigen::Ref<const NodalVectorType> const& p_nodal_values,
- Eigen::Ref<const NodalVectorType> const& C_nodal_values) const
- {
// eval shape matrices at given point
auto const shape_matrices = [&]() {
auto const fe = NumLib::createIsoparametricFiniteElement<
@@ -874,11 +853,11 @@ public:
// local_x contains the local concentration and pressure values
NumLib::shapeFunctionInterpolate(
- C_nodal_values, shape_matrices.N,
+ local_C, shape_matrices.N,
std::get<double>(vars[static_cast<int>(
MaterialPropertyLib::Variable::concentration)]));
NumLib::shapeFunctionInterpolate(
- p_nodal_values, shape_matrices.N,
+ local_p, shape_matrices.N,
std::get<double>(vars[static_cast<int>(
MaterialPropertyLib::Variable::phase_pressure)]));
@@ -891,8 +870,7 @@ public:
.template value<double>(vars, pos, t);
GlobalDimMatrixType const K_over_mu = K / mu;
- GlobalDimVectorType q =
- -K_over_mu * shape_matrices.dNdx * p_nodal_values;
+ GlobalDimVectorType q = -K_over_mu * shape_matrices.dNdx * local_p;
auto const rho_w =
phase.property(MaterialPropertyLib::PropertyType::density)
.template value<double>(vars, pos, t);
--
GitLab