blob: 7005f02c0a0742777e56cd0b0070573dc4485de3 (
plain)
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
|
/*
* Edge_profile.h
*
* Created on: Feb 13, 2014
* Author: dsalinas
*/
#ifndef GUDHI_EDGE_PROFILE_H_
#define GUDHI_EDGE_PROFILE_H_
//#include "combinatorics/Skeleton_blocker/Simplex.h"
namespace Gudhi{
namespace contraction {
template<typename GeometricSimplifiableComplex> class Edge_profile{
public:
typedef GeometricSimplifiableComplex Complex;
typedef typename Complex::GT GT;
typedef typename GeometricSimplifiableComplex::Vertex_handle Vertex_handle;
typedef typename GeometricSimplifiableComplex::Root_vertex_handle Root_vertex_handle;
typedef typename GeometricSimplifiableComplex::Edge_handle Edge_handle;
typedef typename GeometricSimplifiableComplex::Graph_vertex Graph_vertex;
typedef typename GeometricSimplifiableComplex::Graph_edge Graph_edge;
typedef typename GeometricSimplifiableComplex::Point Point;
Edge_profile( GeometricSimplifiableComplex& complex,Edge_handle edge):complex_(complex),edge_handle_(edge),
v0_(complex_.first_vertex(edge_handle_)),v1_(complex_.second_vertex(edge_handle_))
{
assert(complex_.get_address(complex_[edge_handle_].first()));
assert(complex_.get_address(complex_[edge_handle_].second()));
assert(complex_.contains_edge(v0_handle(),v1_handle()));
assert(v0_handle() != v1_handle());
}
virtual ~Edge_profile(){ }
GeometricSimplifiableComplex& complex() const {
return complex_;
}
Edge_handle edge_handle() const{
return edge_handle_;
}
Graph_edge& edge() const{
return complex_[edge_handle_];
}
Graph_vertex& v0() const{return complex_[v0_handle()];}
Graph_vertex& v1() const{return complex_[v1_handle()];}
Vertex_handle v0_handle() const{
return v0_;
// Root_vertex_handle root = complex_[edge_handle_].first();
// assert(complex_.get_address(root));
// return *complex_.get_address(root);
}
Vertex_handle v1_handle() const{
return v1_;
// Root_vertex_handle root = complex_[edge_handle_].second();
// assert(complex_.get_address(root));
// return *complex_.get_address(root);
}
const Point& p0() const {return complex_.point(v0_handle());}
const Point& p1() const {return complex_.point(v1_handle());}
friend std::ostream& operator << (std::ostream& o, const Edge_profile & v){
o << "v0:"<<v.v0_handle() << " v1:"<<v.v1_handle();
return o;
}
private:
GeometricSimplifiableComplex& complex_;
Edge_handle edge_handle_;
Vertex_handle v0_;
Vertex_handle v1_;
};
template<typename EdgeProfile> class Edge_profile_factory{
public:
typedef typename EdgeProfile::Edge_handle Edge_handle_;
typedef typename EdgeProfile::Complex Complex_;
virtual EdgeProfile make_profile(
Complex_& complex,
Edge_handle_ edge) const{
return EdgeProfile(complex,edge);
}
virtual ~Edge_profile_factory(){};
};
} // namespace contraction
} // namespace GUDHI
#endif /* GUDHI_EDGE_PROFILE_H_ */
|