Newer
Older
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
/*
* TemplatePoint.h
*
* Created on: Jan 28, 2010
* Author: fischeth
*/
#ifndef TEMPLATEPOINT_H_
#define TEMPLATEPOINT_H_
#include <cassert>
#include <iostream>
#include <string>
#include <sstream>
#include "GeoObject.h"
namespace GEOLIB {
/**
* \ingroup GEOLIB
*
* \brief class-template for points can be instantiated by a numeric type.
* \param T the coordinate type
*/
template <class T> class TemplatePoint : public GeoObject
{
public:
/** default constructor */
TemplatePoint ();
/** constructor - constructs a TemplatePoint object
\param x1 value for the first coordinate
\param x2 value for the second coordinate
\param x3 value for the third coordinate
*/
TemplatePoint (T x1, T x2, T x3);
/** constructor - constructs a TemplatePoint object
\param x values for three coordinates
*/
TemplatePoint (T const* x);
/** virtual destructor */
virtual ~TemplatePoint() {};
/** \brief const access operator
* The access to the point coordinates is like the access to a field. Code example:
* \code
* Point<double> point (1.0, 2.0, 3.0);
* double sqrNrm2 = point[0] * point[0] + point[1] * point[1] + point[2] + point[2];
* \endcode
*/
const T& operator[] (size_t idx) const {
assert (idx <= 2);
return _x[idx];
}
/** \brief access operator (see book Effektiv C++ programmieren - subsection 1.3.2 ).
* \sa const T& operator[] (size_t idx) const
*/
T& operator[] (size_t idx) {
return const_cast<T&> (static_cast<const TemplatePoint&> (*this)[idx]);
}
/** returns an array containing the coordinates of the point */
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
/** write point coordinates into stream (used from operator<<)
* \param os a standard output stream
*/
virtual void write (std::ostream &os) const {
os << _x[0] << " " << _x[1] << " " << _x[2] << std::flush;
}
/**
* write point coordinates into string
*/
virtual std::string write () const {
std::ostringstream strStream;
strStream << _x[0] << " " << _x[1] << " " << _x[2];
return strStream.str();
}
/** read point coordinates into stream (used from operator>>) */
virtual void read (std::istream &is) {
is >> _x[0] >> _x[1] >> _x[2];
}
protected:
T _x[3];
};
template <class T> TemplatePoint<T>::TemplatePoint() :
GeoObject()
{
_x[0] = static_cast<T>(0);
_x[1] = static_cast<T>(0);
_x[2] = static_cast<T>(0);
}
template <class T> TemplatePoint<T>::TemplatePoint(T x1, T x2, T x3) :
GeoObject()
{
_x[0] = x1;
_x[1] = x2;
_x[2] = x3;
}
template <class T> TemplatePoint<T>::TemplatePoint (T const* x) :
GeoObject()
{
for (size_t k(0); k<3; k++) _x[k] = x[k];
}
/** overload the output operator for class Point */
template <class T>
std::ostream& operator<< (std::ostream &os, const TemplatePoint<T> &p)
{
p.write (os);
return os;
}
/** overload the input operator for class Point */
template <class T>
std::istream& operator>> (std::istream &is, TemplatePoint<T> &p)
{
p.read (is);
return is;
}
} // end namespace GEO
#endif /* TEMPLATEPOINT_H_ */