lsst.sphgeom gac51dbfff8+052faf71bd
NormalizedAngle.h
Go to the documentation of this file.
1/*
2 * LSST Data Management System
3 * Copyright 2014-2015 AURA/LSST.
4 *
5 * This product includes software developed by the
6 * LSST Project (http://www.lsst.org/).
7 *
8 * This program is free software: you can redistribute it and/or modify
9 * it under the terms of the GNU General Public License as published by
10 * the Free Software Foundation, either version 3 of the License, or
11 * (at your option) any later version.
12 *
13 * This program is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the LSST License Statement and
19 * the GNU General Public License along with this program. If not,
20 * see <https://www.lsstcorp.org/LegalNotices/>.
21 */
22
23#ifndef LSST_SPHGEOM_NORMALIZEDANGLE_H_
24#define LSST_SPHGEOM_NORMALIZEDANGLE_H_
25
28
29#include "Angle.h"
30
31
32namespace lsst {
33namespace sphgeom {
34
35// Forward declarations
36class LonLat;
37class Vector3d;
38
42public:
43 static NormalizedAngle nan() {
44 return NormalizedAngle(std::numeric_limits<double>::quiet_NaN());
45 }
46
47 static NormalizedAngle fromDegrees(double a) {
48 return NormalizedAngle(a * RAD_PER_DEG);
49 }
50
51 static NormalizedAngle fromRadians(double a) {
52 return NormalizedAngle(a);
53 }
54
63 NormalizedAngle const & b);
64
67 static NormalizedAngle center(NormalizedAngle const & a,
68 NormalizedAngle const & b);
69
72
74 NormalizedAngle(NormalizedAngle const & a) : _a(a._a) {}
75
77 explicit NormalizedAngle(Angle const & a) {
78 *this = NormalizedAngle(a.asRadians());
79 }
80
83 explicit NormalizedAngle(double a) {
84 // For really large |a|, the error in this reduction can exceed 2*PI
85 // (because PI is only an approximation to π).
86 if (a < 0.0) {
87 _a = Angle(std::fmod(a, 2.0 * PI) + 2.0 * PI);
88 } else if (a > 2 * PI) {
89 _a = Angle(std::fmod(a, 2.0 * PI));
90 } else {
91 _a = Angle(a);
92 }
93 }
94
97 NormalizedAngle(LonLat const &, LonLat const &);
98
101 NormalizedAngle(Vector3d const &, Vector3d const &);
102
106 operator Angle const & () const { return _a; }
107
108 // Comparison operators
109 bool operator==(Angle const & a) const { return _a == a; }
110 bool operator!=(Angle const & a) const { return _a != a; }
111 bool operator<(Angle const & a) const { return _a < a; }
112 bool operator>(Angle const & a) const { return _a > a; }
113 bool operator<=(Angle const & a) const { return _a <= a; }
114 bool operator>=(Angle const & a) const { return _a >= a; }
115
116 // Arithmetic operators
117 Angle operator-() const { return Angle(-_a); }
118 Angle operator+(Angle const & a) const { return _a + a; }
119 Angle operator-(Angle const & a) const { return _a - a; }
120 Angle operator*(double a) const { return _a * a; }
121 Angle operator/(double a) const { return _a / a; }
122 double operator/(Angle const & a) const { return _a / a; }
123
125 double asDegrees() const { return _a.asDegrees(); }
126
128 double asRadians() const { return _a.asRadians(); }
129
131 bool isNan() const { return _a.isNan(); }
132
141 double d = a.asRadians() - asRadians();
142 x._a = Angle((d < 0.0) ? 2.0 * PI + d : d);
143 return x;
144 }
145
146private:
147 Angle _a;
148};
149
150inline NormalizedAngle const & abs(NormalizedAngle const & a) { return a; }
151
152}} // namespace lsst::sphgeom
153
154#endif // LSST_SPHGEOM_NORMALIZEDANGLE_H_
This file declares a class for representing angles.
Definition: Angle.h:43
bool isNan() const
isNan returns true if the angle value is NaN.
Definition: Angle.h:91
double asDegrees() const
asDegrees returns the value of this angle in units of degrees.
Definition: Angle.h:82
double asRadians() const
asRadians returns the value of this angle in units of radians.
Definition: Angle.h:85
Definition: LonLat.h:48
Definition: NormalizedAngle.h:41
bool isNan() const
isNan returns true if the angle value is NaN.
Definition: NormalizedAngle.h:131
NormalizedAngle(double a)
Definition: NormalizedAngle.h:83
NormalizedAngle(NormalizedAngle const &a)
This constructor creates a copy of a.
Definition: NormalizedAngle.h:74
double asRadians() const
asRadians returns the value of this angle in units of radians.
Definition: NormalizedAngle.h:128
NormalizedAngle getAngleTo(NormalizedAngle const &a) const
Definition: NormalizedAngle.h:139
NormalizedAngle(Angle const &a)
This constructor creates a normalized copy of a.
Definition: NormalizedAngle.h:77
double asDegrees() const
asDegrees returns the value of this angle in units of degrees.
Definition: NormalizedAngle.h:125
NormalizedAngle()
This constructor creates a NormalizedAngle with a value of zero.
Definition: NormalizedAngle.h:71
static NormalizedAngle between(NormalizedAngle const &a, NormalizedAngle const &b)
Definition: NormalizedAngle.cc:35
static NormalizedAngle center(NormalizedAngle const &a, NormalizedAngle const &b)
Definition: NormalizedAngle.cc:45
Vector3d is a vector in ℝ³ with components stored in double precision.
Definition: Vector3d.h:44