11#ifndef MEB_FILTRATION_H_
12#define MEB_FILTRATION_H_
14namespace Gudhi::cech_complex {
35template<
typename Kernel,
typename SimplicialComplexForMEB,
typename Po
intRange>
37 using Point_d =
typename Kernel::Point_d;
38 using FT =
typename Kernel::FT;
39 using Sphere = std::pair<Point_d, FT>;
45 std::vector<Sphere> cache_;
46 std::vector<Point_d> pts;
47 CGAL::NT_converter<FT, Filtration_value> cvt;
49 auto fun = [&](Simplex_handle sh,
int dim){
55 auto vert_it = verts.begin();
58 auto&& pu = points[u];
59 Point_d m = k.midpoint_d_object()(pu, points[v]);
60 FT r = k.squared_distance_d_object()(m, pu);
61 if (exact) CGAL::exact(r);
64 cache_.emplace_back(std::move(m), std::move(r));
70 maxf = max(maxf, complex.
filtration(face_opposite_vertex.first));
72 auto key = complex.
key(face_opposite_vertex.first);
73 Sphere
const& sph = cache_[key];
74 if (k.squared_distance_d_object()(sph.first, points[face_opposite_vertex.second]) > sph.second)
continue;
86 pts.push_back(points[vertex]);
87 Point_d c = k.construct_circumcenter_d_object()(pts.begin(), pts.end());
88 FT r = k.squared_distance_d_object()(c, pts.front());
89 if (exact) CGAL::exact(r);
90 maxf = max(maxf, cvt(r));
93 cache_.emplace_back(std::move(c), std::move(r));
void assign_MEB_filtration(Kernel &&k, SimplicialComplexForMEB &complex, PointRange points, bool exact=false)
Given a simplicial complex and an embedding of its vertices, this assigns to each simplex a filtratio...
Definition MEB_filtration.h:36
Value type for a filtration function on a cell complex.
Definition FiltrationValue.h:20
Definition SimplicialComplexForMEB.h:22
unspecified Simplex_handle
Handle for a simplex.
Definition SimplicialComplexForMEB.h:24
Simplex_key key(Simplex_handle simplex)
Returns the key assigned to the 'simplex' with assign_key().
Filtration_value filtration(Simplex_handle simplex)
Returns the filtration value to the 'simplex'.
Boundary_opposite_vertex_simplex_range boundary_opposite_vertex_simplex_range(Simplex_handle simplex)
Returns a range of the pairs (simplex, opposite vertex) of the boundary of the 'simplex'.
unspecified Filtration_value
Type of filtration values.
Definition SimplicialComplexForMEB.h:29
void assign_key(Simplex_handle simplex, Simplex_key key)
Assigns this 'key' to the 'simplex'.
Simplex_vertex_range simplex_vertex_range(Simplex_handle simplex)
Returns a range over vertices (as Vertex_handle) of a given simplex.
unspecified Vertex_handle
Handle for a vertex. Must be a non-negative integer, it is also used as an index into the input list ...
Definition SimplicialComplexForMEB.h:27
int assign_filtration(Simplex_handle simplex, Filtration_value filtration)
Assigns this 'filtration' value to the 'simplex'.
void for_each_simplex(auto callback)
Calls callback(simplex, dim) for every simplex of the complex, with the guarantee that faces are visi...
Handle type for the vertices of a cell complex.
Definition VertexHandle.h:15