From 0b508d57c0cf9184efff73a18153d49c4fe3253d Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Wed, 8 May 2019 10:59:13 +0200
Subject: [PATCH] [pointsource] Fix integration point source helper for box
---
.../multidomain/embedded/integrationpointsource.hh | 14 +++++++++++++-
1 file changed, 13 insertions(+), 1 deletion(-)
diff --git a/dumux/multidomain/embedded/integrationpointsource.hh b/dumux/multidomain/embedded/integrationpointsource.hh
index b4b3889b9c..f92c8730a7 100644
--- a/dumux/multidomain/embedded/integrationpointsource.hh
+++ b/dumux/multidomain/embedded/integrationpointsource.hh
@@ -74,6 +74,16 @@ public:
return integrationElement_;
}
+ void setQuadratureWeight(const Scalar w)
+ {
+ qpweight_ = w;
+ }
+
+ void setIntegrationElement(const Scalar ie)
+ {
+ integrationElement_ = ie;
+ }
+
const std::vector<std::size_t>& elementIndices() const
{
return elementIndices_;
@@ -170,8 +180,10 @@ public:
pointSourceMap.at(key).push_back(source);
else
pointSourceMap.insert({key, {source}});
+
+ // adjust the integration element
auto& s = pointSourceMap.at(key).back();
- s *= shapeValues[scv.indexInElement()];
+ s.setIntegrationElement(shapeValues[scv.indexInElement()]*s.integrationElement());
}
}
}
--
GitLab