polynomialcollection.hh

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
/*                                                                           */
/*  This file is part of the library KASKADE 7                               */
/*  https://www.zib.de/research/projects/kaskade7-finite-element-toolbox     */
/*                                                                           */
/*  Copyright (C) 2002-2012 Zuse Institute Berlin                            */
/*                                                                           */
/*  KASKADE 7 is distributed under the terms of the ZIB Academic License.    */
/*    see $KASKADE/academic.txt                                              */
/*                                                                           */
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
#ifndef INTERPOLATION_POLYNOMIAL_COLLECTION_HH
#define INTERPOLATION_POLYNOMIAL_COLLECTION_HH
 
#include <vector>
 
#include "utilities/interpolation/tools.hh"
 
namespace Kaskade
{
  template<class InterpolationPolynomial, template <typename> class PolynomialPolicy = Policy::ShapeFunctionSetPolicy>
  class InterpolationPolynomialCollection
  : private PolynomialPolicy<InterpolationPolynomial>
  {
    typedef typename InterpolationPolynomial::GridView GridView;
    typedef typename GridView::template Codim<0>::Entity Entity;
    typedef typename GridView::template Codim<0>::Iterator CellIterator;
    typedef PolynomialPolicy<InterpolationPolynomial> Policy;
  public:
    typedef typename InterpolationPolynomial::range_type ValueType;
 
    // no copying
    InterpolationPolynomialCollection(InterpolationPolynomialCollection const&) = delete;
    InterpolationPolynomialCollection& operator=(InterpolationPolynomialCollection const&) = delete;
 
 
    template <typename... Parameters>
    InterpolationPolynomialCollection(GridView const& gridView_, Parameters const&... params) :
    gridView(gridView_),
    interpolationPolynomials(gridView.size(0))
    {
      CellIterator end = gridView.template end<0>();
      for(CellIterator ci = gridView.template begin<0>(); ci!=end; ++ci) interpolationPolynomials[gridView.indexSet().index(*ci)] = Policy::init(*ci, gridView, params...);
    }
 
 
    template <class Cell>
    ValueType value(Cell const& cell, Dune::FieldVector<typename Cell::Geometry::ctype, Cell::dimension> const& localCoordinate) const
    {
      return Policy::evaluate( interpolationPolynomials[ gridView.indexSet().index(cell) ], localCoordinate );
    }
 
  private:
    GridView const& gridView;
    std::vector<InterpolationPolynomial> interpolationPolynomials;
    using PolynomialPolicy<InterpolationPolynomial>::init;
  };
} /* end of namespace Kaskade */
#endif