dunefunctionsboundaryfunctionalassembler.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
// SPDX-FileCopyrightText: Copyright © DUNE-FUFEM Project contributors, see file AUTHORS.md
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception OR LGPL-3.0-or-later
 
#ifndef DUNE_FUFEM_ASSEMBERS_DUNEFUNCTIONSBOUNDARYFUNCTIONALASSEMBLER_HH
#define DUNE_FUFEM_ASSEMBERS_DUNEFUNCTIONSBOUNDARYFUNCTIONALASSEMBLER_HH
 
#ifndef DUNE_FUFEM_DISABLE_HEADER_DEPRECATION
#warning This header is deprecated and will be removed after 2.11. Use the global assemblers from dune/assembler/defaultglobalassembler.hh instead.
#endif
 
#include <dune/common/deprecated.hh>
 
#include <dune/istl/matrix.hh>
 
#include "dune/fufem/boundarypatch.hh"
 
 
namespace Dune {
namespace Fufem {
 
 
template <class TestBasis>
class
[[deprecated("This class will be removed after 2.11. Use Dune::Assembler::Assembler together with Dune::Fufem::Forms based local assemblers instead.")]]
DuneFunctionsBoundaryFunctionalAssembler
{
  using GridView = typename TestBasis::GridView;
 
public:
  [[deprecated("This class will be removed after 2.11. Use Dune::Assembler::Assembler together with Dune::Fufem::Forms based local assemblers instead.")]]
  DuneFunctionsBoundaryFunctionalAssembler(const TestBasis& tBasis, const BoundaryPatch<GridView>& boundaryPatch) :
      testBasis_(tBasis),
      boundaryPatch_(boundaryPatch)
  {}
 
 
 
  template <class VectorBackend, class LocalAssembler>
  void assembleBulkEntries(VectorBackend&& vectorBackend, LocalAssembler&& localAssembler) const
  {
    auto testLocalView     = testBasis_.localView();
 
    using Field = std::decay_t<decltype(vectorBackend[testLocalView.index(0)])>;
    using LocalVector = Dune::BlockVector<Dune::FieldVector<Field,1>>;
 
    auto localVector = LocalVector(testLocalView.maxSize());
 
    // iterate over boundary
    for (auto it =boundaryPatch_.begin(); it!= boundaryPatch_.end(); ++it)
    {
      const auto& element = it->inside();
      testLocalView.bind(element);
 
      for (size_t i=0; i<testLocalView.tree().size(); ++i)
      {
        auto localRow = testLocalView.tree().localIndex(i);
        localVector[localRow] = 0;
      }
 
      localAssembler(it, localVector, testLocalView);
 
      for (size_t i=0; i<testLocalView.tree().size(); ++i)
      {
        auto localRow = testLocalView.tree().localIndex(i);
        auto row = testLocalView.index(localRow);
        vectorBackend[row] += localVector[localRow];
      }
    }
  }
 
  template <class VectorBackend, class LocalAssembler>
  void assembleBulk(VectorBackend&& vectorBackend, LocalAssembler&& localAssembler) const
  {
    vectorBackend.resize(testBasis_);
    vectorBackend.vector() = 0.0;
    assembleBulkEntries(vectorBackend, std::forward<LocalAssembler>(localAssembler));
  }
 
 
protected:
  const TestBasis& testBasis_;
  const BoundaryPatch<GridView>& boundaryPatch_;
};
 
 
template <class TestBasis, class BP>
[[deprecated("This function will be removed after 2.11. Use Dune::Assembler::Assembler together with Dune::Fufem::Forms based local assemblers instead.")]]
auto duneFunctionsBoundaryFunctionalAssembler(const TestBasis& testBasis, const BP& bp)
{
DUNE_NO_DEPRECATED_BEGIN
  return DuneFunctionsBoundaryFunctionalAssembler<TestBasis>(testBasis, bp);
DUNE_NO_DEPRECATED_END
}
 
 
 
} // namespace Fufem
} // namespace Dune
 
 
#endif