1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
// ---------------------------------------------------------------------
//
// 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
|
