From 73aa6821200ce22a19fec8580d2377939bcf01e6 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Mon, 27 Apr 2015 17:54:47 +0200
Subject: [PATCH] [PL] Remove wrongly inserted scaling value.
The Neumann BC is a*du/dn = g_N(x), and the user
provides only the g_N(x) function.
---
ProcessLib/GroundwaterFlowProcess.h | 1 -
ProcessLib/NeumannBc.h | 6 +-----
2 files changed, 1 insertion(+), 6 deletions(-)
diff --git a/ProcessLib/GroundwaterFlowProcess.h b/ProcessLib/GroundwaterFlowProcess.h
index be13d4103cf..5ef915b9db9 100644
--- a/ProcessLib/GroundwaterFlowProcess.h
+++ b/ProcessLib/GroundwaterFlowProcess.h
@@ -152,7 +152,6 @@ public:
hydraulic_head_mesh_element_searcher,
_global_setup,
_integration_order,
- _hydraulic_conductivity,
*_local_to_global_index_map,
*_mesh_subset_all_nodes);
}
diff --git a/ProcessLib/NeumannBc.h b/ProcessLib/NeumannBc.h
index 5797f3cb4f9..1a5d24f8649 100644
--- a/ProcessLib/NeumannBc.h
+++ b/ProcessLib/NeumannBc.h
@@ -36,12 +36,10 @@ public:
NeumannBc(
NeumannBcConfig* bc,
unsigned const integration_order,
- double const scaling_value,
AssemblerLib::LocalToGlobalIndexMap const& local_to_global_index_map,
MeshLib::MeshSubset const& mesh_subset_all_nodes
)
:
- _scaling_value(scaling_value),
_function(*bc->getFunction()),
_integration_order(integration_order)
{
@@ -97,7 +95,7 @@ public:
auto elementValueLookup = [this](MeshLib::Element const& e)
{
- return _function() * _scaling_value;
+ return _function();
};
DBUG("Calling local Neumann assembler builder for neumann boundary elements.");
@@ -122,8 +120,6 @@ public:
private:
- double const _scaling_value;
-
MathLib::ConstantFunction<double> const _function;
std::vector<MeshLib::Element const*> _elements;
--
GitLab