grid_tools_cache.cc

// ---------------------------------------------------------------------
//
// Copyright (C) 2017 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------

#include <deal.II/grid/grid_tools_cache.h>
#include <deal.II/grid/grid_tools.h>

DEAL_II_NAMESPACE_OPEN

namespace GridTools
{

  template<int dim, int spacedim>
  Cache<dim,spacedim>::Cache(
    const Triangulation<dim, spacedim> &tria,
    const Mapping<dim, spacedim> &mapping) :
    update_flags(update_all),
    tria(&tria),
    mapping(&mapping)
  {
    tria_signal = tria.signals.any_change.connect([&]()
    {
      mark_for_update(update_all);
    });
  }

  template<int dim, int spacedim>
  Cache<dim,spacedim>::~Cache()
  {
    // Make sure that the signals that was attached to the triangulation
    // is removed here.
    if (tria_signal.connected())
      tria_signal.disconnect();
  }



  template<int dim, int spacedim>
  void Cache<dim,spacedim>::mark_for_update(const CacheUpdateFlags &flags)
  {
    update_flags |= flags;
  }



  template<int dim, int spacedim>
  const std::vector<std::set<typename Triangulation<dim,spacedim>::active_cell_iterator> > &
  Cache<dim,spacedim>::get_vertex_to_cell_map() const
  {
    if (update_flags & update_vertex_to_cell_map)
      {
        vertex_to_cells = GridTools::vertex_to_cell_map(*tria);
        update_flags = update_flags & ~update_vertex_to_cell_map;
      }
    return vertex_to_cells;
  }



  template<int dim, int spacedim>
  const std::vector<std::vector<Tensor<1,spacedim>>> &
  Cache<dim,spacedim>::get_vertex_to_cell_centers_directions() const
  {
    if (update_flags & update_vertex_to_cell_centers_directions)
      {
        vertex_to_cell_centers =
          GridTools::vertex_to_cell_centers_directions(*tria, get_vertex_to_cell_map());
        update_flags = update_flags & ~update_vertex_to_cell_centers_directions;
      }
    return vertex_to_cell_centers;
  }



  template<int dim, int spacedim>
  const std::map<unsigned int, Point<spacedim> > &
  Cache<dim,spacedim>::get_used_vertices() const
  {
    if (update_flags & update_used_vertices)
      {
        used_vertices = GridTools::extract_used_vertices(*tria, *mapping);
        update_flags = update_flags & ~update_used_vertices;
      }
    return used_vertices;
  }



#ifdef DEAL_II_WITH_NANOFLANN
  template<int dim, int spacedim>
  const KDTree<spacedim> &Cache<dim,spacedim>::get_vertex_kdtree() const
  {
    if (update_flags & update_vertex_kdtree)
      {
        vertex_kdtree.set_points(tria->get_vertices());
        update_flags  = update_flags & ~update_vertex_kdtree;
      }
    return vertex_kdtree;
  }
#endif

#include "grid_tools_cache.inst"

} // GridTools

DEAL_II_NAMESPACE_CLOSE