boundaryconstraints.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_CONSTRAINTS_BOUNDARYCONSTRAINTS_HH
#define DUNE_FUFEM_CONSTRAINTS_BOUNDARYCONSTRAINTS_HH
 
#include <variant>
 
#include <dune/functions/functionspacebases/interpolate.hh>
 
#include <dune/fufem/backends/traversal.hh>
#include <dune/fufem/functiontools/boundarydofs.hh>
#include <dune/fufem/constraints/affineconstraints.hh>
#include <dune/fufem/forms/domains.hh>
 
 
 
namespace Dune::Fufem {
 
template<class BV, class V, class MI, class C, class Basis, class Function, class IntersectionSet>
requires (Dune::Fufem::Forms::IsBoundaryDomain<IntersectionSet>::value or Dune::Fufem::Forms::IsSkeletonDomain<IntersectionSet>::value)
void computeBoundaryConstraints(AffineConstraints<BV, V, MI, C>& constraints, const Basis& basis, Function&& f, const IntersectionSet& intersectionSet)
{
  auto&& isConstrained = Dune::Functions::istlVectorBackend(constraints.isConstrained());
  auto isNewConstrained = constraints.isConstrained();
  Dune::Fufem::markIntersectionDofs(intersectionSet, basis, isNewConstrained);
  Dune::Fufem::recursiveForEachVectorEntry(isNewConstrained, [&](auto&& isNewConstrained_i, const auto& i)
  {
    isNewConstrained_i = isNewConstrained_i and not isConstrained[i];
    isConstrained[i] =  isNewConstrained_i or isConstrained[i];
  });
  Dune::Functions::interpolate(basis, constraints.constantTerm(), f, isNewConstrained);
  constraints.resolveDependencies();
}
 
 
template<class BV, class V, class MI, class C, class Basis, class Function>
void computeBoundaryConstraints(AffineConstraints<BV, V, MI, C>& constraints, const Basis& basis, Function&& f)
{
  computeBoundaryConstraints(constraints, basis, f, Dune::Fufem::Forms::boundary);
}
 
 
template<class BitVector=std::monostate, class Vector=std::monostate, class Basis, class Function, class IntersectionSet>
requires (Dune::Fufem::Forms::IsBoundaryDomain<IntersectionSet>::value or Dune::Fufem::Forms::IsSkeletonDomain<IntersectionSet>::value)
auto makeBoundaryConstraints(const Basis& basis, Function&& f, const IntersectionSet& intersectionSet)
{
  auto constraints = makeAffineConstraints<BitVector, Vector>(basis);
  computeBoundaryConstraints(constraints, basis, f, intersectionSet);
  return constraints;
}
 
 
template<class BitVector=std::monostate, class Vector=std::monostate, class Basis, class Function>
auto makeBoundaryConstraints(const Basis& basis, Function&& f)
{
  auto constraints = makeAffineConstraints<BitVector, Vector>(basis);
  computeBoundaryConstraints(constraints, basis, f);
  return constraints;
}
 
} // namespace Dune::Fufem
 
 
 
#endif // DUNE_FUFEM_CONSTRAINTS_BOUNDARYCONSTRAINTS_HH