integrate.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_FORMS_INTEGRATE_HH
#define DUNE_FUFEM_FORMS_INTEGRATE_HH
 
#include <type_traits>
#include <utility>
 
#include <dune/common/typetraits.hh>
 
#include <dune/fufem/forms/baseclass.hh>
 
#include <dune/fufem/forms/domains.hh>
#include <dune/fufem/forms/integratedlinearform.hh>
#include <dune/fufem/forms/integratedbilinearform.hh>
#include <dune/fufem/forms/integratedboundarylinearform.hh>
#include <dune/fufem/forms/integratedboundarybilinearform.hh>
#include <dune/fufem/forms/integratedskeletonbilinearform.hh>
#include <dune/fufem/forms/integratedskeletonlinearform.hh>
#include <dune/fufem/forms/quadraturehints.hh>
 
 
 
namespace Dune::Fufem::Forms {
 
 
  template<class MultilinearOperator, class Domain, class QuadratureHint>
  requires isOperatorOrSumOperator_v<MultilinearOperator> && IsDomain<Domain>::value &&  QuadratureHints::IsQuadratureHint<QuadratureHint>::value
  auto integrate(MultilinearOperator op, const Domain& domain, QuadratureHint hint)
  {
    auto sumOperator = SumOperator(op);
    if constexpr(MultilinearOperator::arity==1)
    {
      if constexpr (IsBulkDomain<Domain>())
        return IntegratedLinearForm(sumOperator, domain, hint);
      if constexpr (IsBoundaryDomain<Domain>())
        return IntegratedBoundaryLinearForm(sumOperator, domain, hint);
      if constexpr (IsSkeletonDomain<Domain>())
        return IntegratedSkeletonLinearForm(sumOperator, domain, hint);
    }
    else if constexpr(MultilinearOperator::arity==2)
    {
      if constexpr (IsBulkDomain<Domain>())
        return IntegratedBilinearForm(sumOperator, domain, hint);
      if constexpr (IsBoundaryDomain<Domain>())
        return IntegratedBoundaryBilinearForm(sumOperator, domain, hint);
      if constexpr (IsSkeletonDomain<Domain>())
        return IntegratedSkeletonBilinearForm(sumOperator, domain, hint);
    }
  }
 
 
 
  template<class MultilinearOperator, class Domain>
  requires (isOperatorOrSumOperator_v<MultilinearOperator> and IsDomain<Domain>::value)
  auto integrate(MultilinearOperator op, const Domain& domain)
  {
    auto sumOperator = SumOperator(op);
    if constexpr(MultilinearOperator::arity==1)
    {
      if constexpr (IsBulkDomain<Domain>())
        return IntegratedLinearForm(sumOperator, domain);
      if constexpr (IsBoundaryDomain<Domain>())
        return IntegratedBoundaryLinearForm(sumOperator, domain);
      if constexpr (IsSkeletonDomain<Domain>())
        return IntegratedSkeletonLinearForm(sumOperator, domain);
    }
    else if constexpr(MultilinearOperator::arity==2)
    {
      if constexpr (IsBulkDomain<Domain>())
        return IntegratedBilinearForm(sumOperator, domain);
      if constexpr (IsBoundaryDomain<Domain>())
        return IntegratedBoundaryBilinearForm(sumOperator, domain);
      if constexpr (IsSkeletonDomain<Domain>())
        return IntegratedSkeletonBilinearForm(sumOperator, domain);
    }
  }
 
 
 
  template<class MultilinearOperator, class QuadratureHint>
  requires isOperatorOrSumOperator_v<MultilinearOperator> && QuadratureHints::IsQuadratureHint<QuadratureHint>::value
  auto integrate(MultilinearOperator op, QuadratureHint hint)
  {
    auto sumOperator = SumOperator(op);
    if constexpr(MultilinearOperator::arity==1)
      return IntegratedLinearForm(sumOperator, Dune::Fufem::Forms::bulk, hint);
    else if constexpr(MultilinearOperator::arity==2)
      return IntegratedBilinearForm(sumOperator, Dune::Fufem::Forms::bulk, hint);
  }
 
 
 
  template<class MultilinearOperator>
  requires isOperatorOrSumOperator_v<MultilinearOperator>
  auto integrate(MultilinearOperator op)
  {
    auto sumOperator = SumOperator(op);
    if constexpr(MultilinearOperator::arity==1)
      return IntegratedLinearForm(sumOperator, Dune::Fufem::Forms::bulk);
    else if constexpr(MultilinearOperator::arity==2)
      return IntegratedBilinearForm(sumOperator, Dune::Fufem::Forms::bulk);
  }
 
 
 
} // namespace Dune::Fufem::Forms
 
 
#endif // DUNE_FUFEM_FORMS_INTEGRATE_HH