File: floatpoint.hpp

package info (click to toggle)
gamera 1:3.4.2+git20160808.1725654-1
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 22,312 kB
  • ctags: 24,991
  • sloc: xml: 122,324; ansic: 52,869; cpp: 50,664; python: 35,034; makefile: 118; sh: 101
file content (124 lines) | stat: -rw-r--r-- 3,163 bytes parent folder | download | duplicates (4)
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
/*
 *
 * Copyright (C) 2005 Michael Droettboom
 *
 * 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 2
 * 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, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
 */

#ifndef mgd_floatpoint
#define mgd_floatpoint

#include "gamera.hpp"
#include <math.h>

template<class T>
class PointBase {
private:
  T m_x, m_y;
public:
  typedef PointBase<T> self;
  
  // Default constructor
  inline PointBase() : m_x(), m_y() {}

  // Basic constructor
  inline PointBase(T x, T y) : m_x(x), m_y(y) {}

  // From standard Gamera point
  inline PointBase(const Gamera::Point& p) : m_x(T(p.x())), m_y(T(p.y())) {}

  // Copy constructor
  inline PointBase(const self& p) : m_x(p.x()), m_y(p.y()) {}

  // Return the x coordinate
  inline T x() const {
    return m_x;
  }

  // Return the y coordinate
  inline T y() const {
    return m_y;
  }

  // Equality
  template<class OtherPoint>
  inline bool operator==(const OtherPoint& p) const {
    T e = std::numeric_limits<T>::epsilon();
    return (std::abs(m_x - T(p.x())) < e &&
	    std::abs(m_y - T(p.y())) < e);
  }

  // Inequality
  template<class OtherPoint>
  inline bool operator!=(const OtherPoint& p) const {
    return (m_x != T(p.x()) || m_y != T(p.y()));
  }

  // Addition
  template<class OtherPoint>
  inline self operator+(const OtherPoint& p) const {
    return self(m_x + T(p.x()), m_y + T(p.y()));
  } 

  // Subtraction
  template<class OtherPoint>
  inline self operator-(const OtherPoint& p) const {
    return self(m_x - T(p.x()), m_y - T(p.y()));
  } 

  // Multiplication
  template<class OtherPoint>
  inline self operator*(const OtherPoint& p) const {
    return self(m_x * T(p.x()), m_y * T(p.y()));
  } 

  // Division
  template<class OtherPoint>
  inline self operator/(const OtherPoint& p) const {
    return self(m_x / T(p.x()), m_y / T(p.y()));
  } 

  // Unary minus
  inline self operator-() const {
    return self(-m_x, -m_y);
  }

  // Unary plus
  inline self operator+() const {
    return self(+m_x, +m_y);
  }
  
  // Euclidean distance
  template<class OtherPoint>
  inline T distance(const OtherPoint& p) const {
    T x_dist = m_x - T(p.x());
    T y_dist = m_y - T(p.y());
    return std::sqrt(x_dist*x_dist + y_dist*y_dist);
  }
};

template<class T>
inline static PointBase<T> abs(PointBase<T> p) {
  return PointBase<T>(std::abs(p.x()), std::abs(p.y()));
}

template<class T>
std::ostream& operator<<(std::ostream& s, const PointBase<T>& point) {
  return s << "FloatPoint(" << point.x() << ", " << point.y() << ")";
}

typedef PointBase<double> FloatPoint;

#endif