boundarydofs.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set ts=4 sw=2 et 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_FUNCTIONTOOLS_BOUNDARYDOFS_HH
#define DUNE_FUFEM_FUNCTIONTOOLS_BOUNDARYDOFS_HH
 
#include <dune/common/concept.hh>
 
#include <dune/fufem/boundarypatch.hh>
#include <dune/fufem/forms/domains.hh>
 
#include <dune/functions/backends/concepts.hh>
#include <dune/functions/backends/istlvectorbackend.hh>
#include <dune/functions/functionspacebases/subentitydofs.hh>
 
 
 
namespace Dune::Fufem {
 
  template <class IntersectionSet, class Basis, class Vector>
  requires (Dune::Fufem::Forms::IsBoundaryDomain<IntersectionSet>::value or Dune::Fufem::Forms::IsSkeletonDomain<IntersectionSet>::value)
  void markIntersectionDofs(const IntersectionSet& intersectionSet,
                             const Basis& basis,
                             Vector& vector)
  {
    auto&& toVectorBackend = [&](auto& v) -> decltype(auto) {
      if constexpr (Dune::models<Dune::Functions::Concept::VectorBackend<Basis>, decltype(v)>()) {
        return v;
      } else {
        return Dune::Functions::istlVectorBackend(v);
      }
    };
    auto&& vectorBackend = toVectorBackend(vector);
    vectorBackend.resize(basis);
    vectorBackend = false;
 
    auto localView = basis.localView();
    auto seDOFs = Dune::Functions::subEntityDOFs(basis);
    const auto& gridView = basis.gridView();
    if constexpr (requires{ intersectionSet.begin(); intersectionSet.end(); })
      for(const auto& intersection: intersectionSet)
      {
        localView.bind(intersection.inside());
        for(auto localIndex: seDOFs.bind(localView,intersection))
          vectorBackend[localView.index(localIndex)] = true;
      }
    else
      for(auto&& element : elements(gridView))
      {
        if constexpr (requires{intersectionSet.containsFaceOf(element);})
          if (not intersectionSet.containsFaceOf(element))
            continue;
        localView.bind(element);
        for(const auto& intersection: intersections(gridView, element))
        {
          if constexpr (Dune::Fufem::Forms::IsBoundaryDomain<IntersectionSet>())
            if (not intersection.boundary())
              continue;
          if (intersectionSet.contains(intersection))
            for(auto localIndex: seDOFs.bind(localView,intersection))
              vectorBackend[localView.index(localIndex)] = true;
        }
      }
  }
 
 
 
  template <class GridView, class Basis, class Vector>
  void markBoundaryPatchDofs(const BoundaryPatch<GridView>& boundaryPatch,
                             const Basis& basis,
                             Vector& vector)
  {
    markIntersectionDofs(boundaryPatch, basis, vector);
  }
 
 
 
  template <class Basis, class Vector>
  void markBoundaryDofs(const Basis& basis, Vector& vector)
  {
    markIntersectionDofs(Dune::Fufem::Forms::boundary, basis, vector);
  }
 
 
} // namespace Dune::Fufem
 
 
 
#endif   // DUNE_FUFEM_FUNCTIONTOOLS_BOUNDARYDOFS_HH