localfunctionadaptor.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_LOCALFUNCTIONADAPTOR_HH
#define DUNE_FUFEM_FORMS_LOCALFUNCTIONADAPTOR_HH
 
#include <utility>
 
#include <dune/fufem/forms/shapefunctioncache.hh>
 
 
 
namespace Dune::Fufem::Forms {
 
 
  template<class LocalOperator>
  class LocalFunctionAdaptor
  {
  public:
 
    using Element = typename LocalOperator::Element;
    using Range = typename LocalOperator::Range;
 
  private:
 
    static constexpr int dimension = Element::dimension;
    using CacheManager = typename Dune::Fufem::Forms::CacheManager<double, dimension>;
    using QuadratureRule = typename CacheManager::QuadratureRule;
    using QuadraturePoint = typename QuadratureRule::value_type;
    using Domain = typename Element::Geometry::LocalCoordinate;
 
  public:
 
    LocalFunctionAdaptor(LocalOperator&& localOperator) :
      localOperator_(std::move(localOperator)),
      internalRule_(),
      internalCacheManager_()
    {
      localOperator_.registerCaches(internalCacheManager_);
      internalRule_.emplace_back(typename QuadraturePoint::Vector(),0.0);
    }
 
    LocalFunctionAdaptor(const LocalOperator& localOperator) :
      localOperator_(localOperator),
      internalRule_(),
      internalCacheManager_()
    {
      localOperator_.registerCaches(internalCacheManager_);
      internalRule_.emplace_back(typename QuadraturePoint::Vector(),0.0);
    }
 
    LocalFunctionAdaptor(const LocalFunctionAdaptor& other) :
      LocalFunctionAdaptor(other.localOperator_)
    {}
 
    LocalFunctionAdaptor(LocalFunctionAdaptor&& other) :
      LocalFunctionAdaptor(std::move(other.localOperator_))
    {}
 
    void bind(const Element& element)
    {
      localOperator_.bind(element);
    }
 
    void unbind()
    {
      localOperator_.unbind();
    }
 
    bool bound() const
    {
      return localOperator_.bound();
    }
 
    auto operator()(const Domain& x) const
    {
      // The cache must always be fully invalidated,
      // because the quadrature rule was modified.
      // Since the internal cache was marked as non affine,
      // calling invalidate() will indeed remove all cached
      // entries.
      internalRule_[0] = QuadraturePoint{x, 0};
      internalCacheManager_.setRule(internalRule_);
      localOperator_.bindToCaches(internalCacheManager_);
      return localOperator_(0)();
    }
 
    const Element& localContext() const
    {
      return localOperator_.localContext();
    }
 
  private:
    mutable LocalOperator localOperator_;
    mutable QuadratureRule internalRule_;
    mutable CacheManager internalCacheManager_;
  };
 
 
} // namespace Dune::Fufem::Forms
 
 
#endif // DUNE_FUFEM_FORMS_LOCALFUNCTIONADAPTOR_HH