From d55bf37cdbd430264850c0aefe937a657cf990fb Mon Sep 17 00:00:00 2001
From: Kilian Weishaupt <kilian.weishaupt@iws.uni-stuttgart.de>
Date: Tue, 14 Aug 2018 18:02:04 +0200
Subject: [PATCH] [md][couplingmanager] Store tuple of problem weak ptrs in the
base class
---
dumux/multidomain/couplingmanager.hh | 46 ++++++++++++++++++++++++++--
1 file changed, 44 insertions(+), 2 deletions(-)
diff --git a/dumux/multidomain/couplingmanager.hh b/dumux/multidomain/couplingmanager.hh
index bb58de44a6..1008a16bb2 100644
--- a/dumux/multidomain/couplingmanager.hh
+++ b/dumux/multidomain/couplingmanager.hh
@@ -25,11 +25,14 @@
#ifndef DUMUX_MULTIDOMAIN_COUPLING_MANAGER_HH
#define DUMUX_MULTIDOMAIN_COUPLING_MANAGER_HH
+#include <memory>
+#include <tuple>
#include <vector>
-#include <dune/common/indices.hh>
#include <dune/common/exceptions.hh>
-#include <dumux/common/typetraits/typetraits.hh>
+#include <dune/common/indices.hh>
#include <dumux/assembly/numericepsilon.hh>
+#include <dumux/common/properties.hh>
+#include <dumux/common/typetraits/typetraits.hh>
namespace Dumux {
@@ -45,6 +48,9 @@ class CouplingManager
template<std::size_t id> using PrimaryVariables = typename GET_PROP_TYPE(SubDomainTypeTag<id>, PrimaryVariables);
template<std::size_t id> using GridView = typename GET_PROP_TYPE(SubDomainTypeTag<id>, FVGridGeometry)::GridView;
template<std::size_t id> using Element = typename GridView<id>::template Codim<0>::Entity;
+ template<std::size_t id> using Problem = typename GET_PROP_TYPE(SubDomainTypeTag<id>, Problem);
+ template<std::size_t id> using ProblemWeakPtr = std::weak_ptr<const Problem<id>>;
+ using Problems = typename Traits::template MakeTuple<ProblemWeakPtr>;
public:
//! default type used for coupling element stencils
@@ -231,6 +237,36 @@ public:
return NumericEpsilon<typename Traits::Scalar, numEq>(paramGroup);
}
+ /*!
+ * \brief set the pointers to the sub problems
+ * \param problems A tuple of shared pointers to the sub problems
+ */
+ template<typename... SubProblems>
+ void setSubProblems(const std::tuple<std::shared_ptr<SubProblems>...>& problems)
+ { problems_ = problems; }
+
+ /*!
+ * \brief set a pointer to one of the sub problems
+ * \param problem a pointer to the sub problem
+ * \param domainIdx the domain index of the sub problem
+ */
+ template<class SubProblem, std::size_t i>
+ void setSubProblem(std::shared_ptr<SubProblem> problem, Dune::index_constant<i> domainIdx)
+ { std::get<i>(problems_) = problem; }
+
+ /*!
+ * \brief Return a reference to the sub problem
+ * \param domainIdx The domain index
+ */
+ template<std::size_t i>
+ const Problem<i>& problem(Dune::index_constant<i> domainIdx) const
+ {
+ if (!std::get<i>(problems_).expired())
+ return *std::get<i>(problems_).lock();
+ else
+ DUNE_THROW(Dune::InvalidStateException, "The problem pointer was not set or has already expired. Use setSubProblems() before calling this function");
+ }
+
protected:
/*!
@@ -254,6 +290,12 @@ private:
*/
SolutionVector curSol_;
+ /*!
+ * \brief A tuple of std::weak_ptrs to the sub problems
+ * \note these are weak pointers and not shared pointers to break the cyclic dependency between coupling manager and problems
+ */
+ Problems problems_;
+
};
} //end namespace Dumux
--
GitLab