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
|
/*************************************************************************
** Pair.hpp **
** **
** This file is part of dvisvgm -- a fast DVI to SVG converter **
** Copyright (C) 2005-2023 Martin Gieseking <martin.gieseking@uos.de> **
** **
** This program is free software; you can redistribute it and/or **
** modify it under the terms of the GNU General Public License as **
** published by the Free Software Foundation; either version 3 of **
** the License, or (at your option) any later version. **
** **
** This program is distributed in the hope that it will be useful, but **
** WITHOUT ANY WARRANTY; without even the implied warranty of **
** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the **
** GNU General Public License for more details. **
** **
** You should have received a copy of the GNU General Public License **
** along with this program; if not, see <http://www.gnu.org/licenses/>. **
*************************************************************************/
#ifndef PAIR_HPP
#define PAIR_HPP
#include <cmath>
#include <ostream>
#include "macros.hpp"
template <typename T>
class Pair {
public:
explicit Pair (T x=0, T y=0) : _x(x), _y(y) {}
template <typename U> Pair (const Pair<U> &p) : _x(U(p.x())), _y(U(p.y())) {}
Pair (const Pair &p) =default;
Pair (Pair &&p) =default;
Pair& operator = (const Pair &p) =default;
Pair& operator = (Pair &&p) =default;
Pair operator += (const Pair &p) {_x += p._x; _y += p._y; return *this;}
Pair operator -= (const Pair &p) {_x -= p._x; _y -= p._y; return *this;}
Pair operator *= (T c) {_x *= c; _y *= c; return *this;}
Pair operator /= (T c) {_x /= c; _y /= c; return *this;}
Pair operator - () const {return Pair(-_x, -_y);}
Pair ortho () const {return Pair(-_y, _x);}
double length () const {return std::hypot(_x, _y);}
bool operator == (const Pair &p) const {return _x == p._x && _y == p._y;}
bool operator != (const Pair &p) const {return _x != p._x || _y != p._y;}
T x () const {return _x;}
T y () const {return _y;}
void x (const T &xx) {_x = xx;}
void y (const T &yy) {_y = yy;}
std::ostream& write (std::ostream &os) const {return os << '(' << _x << ',' << _y << ')';}
private:
T _x, _y;
};
template <typename T>
inline Pair<T> abs (const Pair<T> &p) {
return Pair<T>(std::abs(p.x()), std::abs(p.y()));
}
/** Returns the dot product of two 2D vectors. */
template <typename T>
inline T dot (const Pair<T> &p1, const Pair<T> &p2) {
return p1.x()*p1.y() + p1.y()*p2.y();
}
/** Returns the determinant of two 2D vectors. */
template <typename T>
inline T det (const Pair<T> &p1, const Pair<T> &p2) {
return p1.x()*p2.y() - p1.y()*p2.x();
}
struct Pair32 : public Pair<int32_t> {
explicit Pair32 (int32_t x=0, int32_t y=0) : Pair<int32_t>(x, y) {}
explicit Pair32 (double x, double y) : Pair<int32_t>(lround(x), lround(y)) {}
explicit Pair32 (const Pair<double> &p) : Pair32(p.x(), p.y()) {}
Pair32 (const Pair<int32_t> &p) : Pair<int32_t>(p) {}
};
typedef Pair<double> DPair;
inline DPair round (const DPair &p) {
return DPair(std::lround(p.x()), std::lround(p.y()));
}
template<> template<>
inline Pair<int>::Pair (const Pair<double> &p) : _x(int(round(p.x()))), _y(int(round(p.y()))) {}
template <typename T>
IMPLEMENT_ARITHMETIC_OPERATOR(Pair<T>, +)
template <typename T>
IMPLEMENT_ARITHMETIC_OPERATOR(Pair<T>, -)
template <typename T>
IMPLEMENT_ARITHMETIC_OPERATOR2(Pair<T>, T, *)
template <typename T>
IMPLEMENT_ARITHMETIC_OPERATOR2(Pair<T>, T, /)
template <typename T>
IMPLEMENT_OUTPUT_OPERATOR(Pair<T>)
IMPLEMENT_ARITHMETIC_OPERATOR2(Pair32, int32_t, *)
#endif
|
