cellLocator.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) 2017-2020 Zuse Institute Berlin                            */
/*                                                                           */
/*  KASKADE 7 is distributed under the terms of the ZIB Academic License.    */
/*    see $KASKADE/academic.txt                                              */
/*                                                                           */
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
#ifndef CELLLOCATOR_HH
#define CELLLOCATOR_HH
 
#include <memory>
#include <optional>
#include <tuple>
 
#include "dune/common/fvector.hh"
#include "dune/geometry/type.hh"
 
#include "utilities/timing.hh"
 
namespace Kaskade
{
  // forward declaration for PIMPL idiome
  namespace CellLocatorDetail
  {
    template <class GridView, int dimw>
    class SpatialIndex;    
  }
  
  // ----------------------------------------------------------------------------------------------
 
 
 
  template <class GridView, int dimw=GridView::Grid::dimension>
  class CellLocator
  {
    using Index = CellLocatorDetail::SpatialIndex<GridView,dimw>;
    using Grid = typename GridView::Grid;
 
  public:
    static int const dimension = GridView::Grid::dimension;
    static int const dimension_world = dimw;
    using ctype    = typename GridView::ctype;
    using Position = Dune::FieldVector<ctype,dimension_world>;
    using Cell = typename GridView::template Codim<0>::Entity;
 
    
    CellLocator(GridView const& gridView);
 
    // we need this because the default destructor doesn't know about the size of SpatialIndex when deallocating via unique_ptr
    ~CellLocator();
 
    ctype diameter() const 
    {
      return diam;
    }
    std::vector<Cell> cells() const;
    
    std::pair<Cell,double> closestCell(Position const& pos) const;
 
  private:
    std::unique_ptr<CellLocatorDetail::SpatialIndex<GridView,dimension_world>> spatialIndex;
    ctype diam;
 
    void init(GridView const& gridView);
    
  };
}
 
#endif