From 1cb993021573cdef88f729d922423d0d06f1b226 Mon Sep 17 00:00:00 2001
From: Lasse Hinrichsen <lh1887@mi.fu-berlin.de>
Date: Tue, 17 Jan 2017 11:17:31 +0100
Subject: [PATCH] Add BoundaryAssembler that works with dune-functions
Analogously to the DuneFunctionsFunctionalAssembler this adds the
possibility to use dune-functions with the assemblers for boundary
patches.
---
dune/fufem/assemblers/CMakeLists.txt | 1 +
...unefunctionsboundaryfunctionalassembler.hh | 93 +++++++++++++++++++
2 files changed, 94 insertions(+)
create mode 100644 dune/fufem/assemblers/dunefunctionsboundaryfunctionalassembler.hh
diff --git a/dune/fufem/assemblers/CMakeLists.txt b/dune/fufem/assemblers/CMakeLists.txt
index 24a2f374..f588564a 100644
--- a/dune/fufem/assemblers/CMakeLists.txt
+++ b/dune/fufem/assemblers/CMakeLists.txt
@@ -5,6 +5,7 @@ install(FILES
boundaryfunctionalassembler.hh
boundaryoperatorassembler.hh
dunefunctionsoperatorassembler.hh
+ dunefunctionsboundaryfunctionalassembler.hh
dunefunctionsfunctionalassembler.hh
functionalassembler.hh
integraloperatorassembler.hh
diff --git a/dune/fufem/assemblers/dunefunctionsboundaryfunctionalassembler.hh b/dune/fufem/assemblers/dunefunctionsboundaryfunctionalassembler.hh
new file mode 100644
index 00000000..af6cc24b
--- /dev/null
+++ b/dune/fufem/assemblers/dunefunctionsboundaryfunctionalassembler.hh
@@ -0,0 +1,93 @@
+// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
+// vi: set et ts=4 sw=2 sts=2:
+#ifndef DUNE_FUFEM_ASSEMBERS_DUNEFUNCTIONSBOUNDARYFUNCTIONALASSEMBLER_HH
+#define DUNE_FUFEM_ASSEMBERS_DUNEFUNCTIONSBOUNDARYFUNCTIONALASSEMBLER_HH
+
+#include <dune/istl/matrix.hh>
+#include <dune/istl/matrixindexset.hh>
+
+#include <dune/matrix-vector/axpy.hh>
+
+#include "dune/fufem/functionspacebases/functionspacebasis.hh"
+#include "dune/fufem/boundarypatch.hh"
+
+
+namespace Dune {
+namespace Fufem {
+
+
+//! Generic global assembler for functionals on a boundary
+template <class TrialBasis>
+class DuneFunctionsBoundaryFunctionalAssembler
+{
+ using GridView = typename TrialBasis::GridView;
+
+public:
+ //! create assembler for grid
+ DuneFunctionsBoundaryFunctionalAssembler(const TrialBasis& tBasis, const BoundaryPatch<GridView>& boundaryPatch) :
+ trialBasis_(tBasis),
+ boundaryPatch_(boundaryPatch)
+ {}
+
+
+
+ template <class VectorBackend, class LocalAssembler>
+ void assembleBulkEntries(VectorBackend&& vectorBackend, LocalAssembler&& localAssembler) const
+ {
+ auto trialLocalView = trialBasis_.localView();
+ auto trialLocalIndexSet = trialBasis_.localIndexSet();
+
+ using Field = std::decay_t<decltype(vectorBackend(trialLocalIndexSet.index(0)))>;
+ using LocalVector = Dune::BlockVector<Dune::FieldVector<Field,1>>;
+
+ auto localVector = LocalVector(trialLocalView.maxSize());
+
+ // iterate over boundary
+ for (auto it =boundaryPatch_.begin(); it!= boundaryPatch_.end(); ++it)
+ {
+ const auto& element = it->inside();
+ trialLocalView.bind(element);
+ trialLocalIndexSet.bind(trialLocalView);
+
+ localAssembler(it, localVector, trialLocalView);
+
+ for (size_t localRow=0; localRow<trialLocalIndexSet.size(); ++localRow)
+ {
+ auto row = trialLocalIndexSet.index(localRow);
+ vectorBackend(row) += localVector[localRow];
+ }
+ }
+ }
+
+ template <class VectorBackend, class LocalAssembler>
+ void assembleBulk(VectorBackend&& vectorBackend, LocalAssembler&& localAssembler) const
+ {
+ vectorBackend.resize(trialBasis_);
+ vectorBackend.vector() = 0.0;
+ assembleBulkEntries(vectorBackend, std::forward<LocalAssembler>(localAssembler));
+ }
+
+
+protected:
+ const TrialBasis& trialBasis_;
+ const BoundaryPatch<GridView>& boundaryPatch_;
+};
+
+
+/**
+ * \brief Create DuneFunctionsBoundaryFunctionalAssembler
+ */
+template <class TrialBasis, class BP>
+auto duneFunctionsBoundaryFunctionalAssembler(const TrialBasis& trialBasis, const BP& bp)
+{
+ return DuneFunctionsBoundaryFunctionalAssembler<TrialBasis>(trialBasis, bp);
+}
+
+
+
+} // namespace Fufem
+} // namespace Dune
+
+
+#endif
+
--
GitLab