localinterpolation.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file LICENSE.md in module root
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
#ifndef DUNE_LOCALFUNCTIONS_ENRICHED_SIMPLEXP1BUBBLE_LOCALINTERPOLATION_HH
#define DUNE_LOCALFUNCTIONS_ENRICHED_SIMPLEXP1BUBBLE_LOCALINTERPOLATION_HH
 
#include <type_traits>
#include <vector>
 
namespace Dune
{
  template<class LB>
  class SimplexP1BubbleLocalInterpolation
  {
    static const int dim = LB::dimension;
    static const int numVertices = dim+1;
 
    using DomainType = typename LB::Traits::DomainType;
    using RangeType = typename LB::Traits::RangeType;
 
  public:
    template<class F, class C,
      class R = std::invoke_result_t<F, DomainType>,
      std::enable_if_t<std::is_convertible_v<R, C>, int> = 0>
    static constexpr void interpolate (const F& f, std::vector<C>& out)
    {
      out.resize(numVertices+1);
 
      // vertices
      DomainType x(0);
      out[0] = f(x);
 
      for (int i = 0; i < dim; ++i) {
        x = 0;
        x[i] = 1;
        out[i+1] = f(x);
      }
 
      // element bubble
      x = 1.0/(dim+1);
      R y = f(x);
 
      // evaluate the other shape functions in x and subtract this value
      std::vector<RangeType> sfValues;
      LB::evaluateFunction(x, sfValues);
 
      out[numVertices] = y;
      for (int i = 0; i < numVertices; ++i)
        out[numVertices] -= out[i]*sfValues[i];
    }
  };
 
} // end namespace Dune
 
#endif // DUNE_LOCALFUNCTIONS_ENRICHED_SIMPLEXP1BUBBLE_LOCALINTERPOLATION_HH