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
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
|
// Copyright (c) 2014 INRIA Sophia-Antipolis (France), INRIA Lorraine LORIA.
// All rights reserved.
//
// This file is part of CGAL (www.cgal.org).
//
// $URL: https://github.com/CGAL/cgal/blob/v6.1.1/Optimal_transportation_reconstruction_2/include/CGAL/OTR_2/Sample.h $
// $Id: include/CGAL/OTR_2/Sample.h 08b27d3db14 $
// SPDX-License-Identifier: GPL-3.0-or-later OR LicenseRef-Commercial
//
// Author(s) : Fernando de Goes, Pierre Alliez, Ivo Vigan, Clément Jamin
#ifndef CGAL_OTR2_SAMPLE_H
#define CGAL_OTR2_SAMPLE_H
#include <CGAL/license/Optimal_transportation_reconstruction_2.h>
/// \cond SKIP_IN_MANUAL
namespace CGAL {
namespace OTR_2 {
template <class Traits_>
class Sample
{
public:
typedef typename Traits_::FT FT;
typedef typename Traits_::Point_2 Point;
private:
Point m_point;
FT m_mass;
FT m_dist2_to_edge;
FT m_coordinate;
FT m_backup_dist2;
FT m_backup_coord;
public:
Sample()
{}
Sample(const Point& point,
const FT mass = FT(1))
: m_point(point),
m_mass(mass),
m_dist2_to_edge(0),
m_coordinate(0),
m_backup_dist2(0),
m_backup_coord(0)
{
}
Sample(const Sample& sample)
: m_point(sample.point()),
m_mass(sample.mass()),
m_dist2_to_edge(0),
m_coordinate(0),
m_backup_dist2(0),
m_backup_coord(0)
{
}
~Sample() { }
const Point& point() const { return m_point; }
Point& point() { return m_point; }
const FT& mass() const { return m_mass; }
FT& mass() { return m_mass; }
const FT& distance2() const { return m_dist2_to_edge; }
FT& distance2() { return m_dist2_to_edge; }
const FT& coordinate() const { return m_coordinate; }
FT& coordinate() { return m_coordinate; }
void backup()
{
m_backup_dist2 = m_dist2_to_edge;
m_backup_coord = m_coordinate;
}
void restore()
{
m_dist2_to_edge = m_backup_dist2;
m_coordinate = m_backup_coord;
}
};
template <class Sample_>
class Sample_with_priority
{
public:
typedef typename Sample_::FT FT;
private:
int m_sample;
FT m_priority;
public:
Sample_with_priority(int sample, const FT priority = FT(0))
: m_sample(sample), m_priority(priority)
{}
Sample_with_priority(const Sample_with_priority& psample)
: m_sample(psample.sample()), m_priority(psample.priority())
{}
~Sample_with_priority()
{}
Sample_with_priority& operator = (const Sample_with_priority& psample)
{
m_sample = psample.sample();
m_priority = psample.priority();
return *this;
}
int sample() const { return m_sample; }
const FT priority() const { return m_priority; }
};
template <class T>
struct greater_priority
{
bool operator() (const T& a, const T& b) const
{
return ( a.priority() > b.priority() );
}
};
} } //end namespaces
/// \endcond
#endif // CGAL_OTR2_SAMPLE_H
|
