LSSTApplications  17.0+11,17.0+34,17.0+56,17.0+57,17.0+59,17.0+7,17.0-1-g377950a+33,17.0.1-1-g114240f+2,17.0.1-1-g4d4fbc4+28,17.0.1-1-g55520dc+49,17.0.1-1-g5f4ed7e+52,17.0.1-1-g6dd7d69+17,17.0.1-1-g8de6c91+11,17.0.1-1-gb9095d2+7,17.0.1-1-ge9fec5e+5,17.0.1-1-gf4e0155+55,17.0.1-1-gfc65f5f+50,17.0.1-1-gfc6fb1f+20,17.0.1-10-g87f9f3f+1,17.0.1-11-ge9de802+16,17.0.1-16-ga14f7d5c+4,17.0.1-17-gc79d625+1,17.0.1-17-gdae4c4a+8,17.0.1-2-g26618f5+29,17.0.1-2-g54f2ebc+9,17.0.1-2-gf403422+1,17.0.1-20-g2ca2f74+6,17.0.1-23-gf3eadeb7+1,17.0.1-3-g7e86b59+39,17.0.1-3-gb5ca14a,17.0.1-3-gd08d533+40,17.0.1-30-g596af8797,17.0.1-4-g59d126d+4,17.0.1-4-gc69c472+5,17.0.1-6-g5afd9b9+4,17.0.1-7-g35889ee+1,17.0.1-7-gc7c8782+18,17.0.1-9-gc4bbfb2+3,w.2019.22
LSSTDataManagementBasePackage
Vector3d.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_VECTOR3D_H_
24 #define LSST_SPHGEOM_VECTOR3D_H_
25 
28 
29 #include <cmath>
30 #include <iosfwd>
31 #include <limits>
32 #include <stdexcept>
33 
34 
35 namespace lsst {
36 namespace sphgeom {
37 
38 // Forward declarations
39 class Angle;
40 class UnitVector3d;
41 
42 
44 class Vector3d {
45 public:
47  Vector3d() { _v[0] = 0.0; _v[1] = 0.0; _v[2] = 0.0; }
48 
50  Vector3d(double x, double y, double z) { _v[0] = x; _v[1] = y; _v[2] = z; }
51 
52  bool operator==(Vector3d const & v) const {
53  return _v[0] == v._v[0] && _v[1] == v._v[1] && _v[2] == v._v[2];
54  }
55 
56  bool operator!=(Vector3d const & v) const {
57  return _v[0] != v._v[0] || _v[1] != v._v[1] || _v[2] != v._v[2];
58  }
59 
61  double const * getData() const { return _v; }
62 
64  double operator()(int i) const { return _v[i]; }
65 
66  double x() const { return _v[0]; }
67 
68  double y() const { return _v[1]; }
69 
70  double z() const { return _v[2]; }
71 
73  double dot(Vector3d const & v) const {
74  return _v[0] * v._v[0] + _v[1] * v._v[1] + _v[2] * v._v[2];
75  }
76 
78  double getSquaredNorm() const { return dot(*this); }
79 
81  double getNorm() const {
82  return std::sqrt(getSquaredNorm());
83  }
84 
86  bool isZero() const { return *this == Vector3d(); }
87 
93  double normalize();
94 
96  bool isNormalized() const {
97  return std::fabs(1.0 - getSquaredNorm()) <= 1e-15;
98  }
99 
101  Vector3d cross(Vector3d const & v) const {
102  return Vector3d(_v[1] * v._v[2] - _v[2] * v._v[1],
103  _v[2] * v._v[0] - _v[0] * v._v[2],
104  _v[0] * v._v[1] - _v[1] * v._v[0]);
105  }
106 
108  Vector3d operator-() const {
109  return Vector3d(-_v[0],
110  -_v[1],
111  -_v[2]);
112  }
113 
116  Vector3d operator*(double s) const {
117  return Vector3d(_v[0] * s,
118  _v[1] * s,
119  _v[2] * s);
120  }
121 
124  Vector3d operator/(double s) const {
125  return Vector3d(_v[0] / s,
126  _v[1] / s,
127  _v[2] / s);
128  }
129 
131  Vector3d operator+(Vector3d const & v) const {
132  return Vector3d(_v[0] + v._v[0],
133  _v[1] + v._v[1],
134  _v[2] + v._v[2]);
135  }
136 
138  Vector3d operator-(Vector3d const & v) const {
139  return Vector3d(_v[0] - v._v[0],
140  _v[1] - v._v[1],
141  _v[2] - v._v[2]);
142  }
143 
144  Vector3d & operator*=(double s) { *this = *this * s; return *this; }
145  Vector3d & operator/=(double s) { *this = *this / s; return *this; }
146  Vector3d & operator+=(Vector3d const & v) { *this = *this + v; return *this; }
147  Vector3d & operator-=(Vector3d const & v) { *this = *this - v; return *this; }
148 
150  Vector3d cwiseProduct(Vector3d const & v) const {
151  return Vector3d(_v[0] * v._v[0],
152  _v[1] * v._v[1],
153  _v[2] * v._v[2]);
154  }
155 
158  Vector3d rotatedAround(UnitVector3d const & k, Angle a) const;
159 
160 private:
161  double _v[3];
162 };
163 
164 
165 inline Vector3d operator*(double s, Vector3d const & v) { return v * s; }
166 
168 
169 }} // namespace lsst::sphgeom
170 
171 #endif // LSST_SPHGEOM_VECTOR3D_H_
double normalize()
normalize scales this vector to have unit norm and returns its norm prior to scaling.
Definition: Vector3d.cc:41
Vector3d cwiseProduct(Vector3d const &v) const
cwiseProduct returns the component-wise product of this vector and v.
Definition: Vector3d.h:150
bool isNormalized() const
isNormalized returns true if this vectors norm is very close to 1.
Definition: Vector3d.h:96
Vector3d()
The default constructor creates a zero vector.
Definition: Vector3d.h:47
bool operator!=(Vector3d const &v) const
Definition: Vector3d.h:56
Vector3d & operator+=(Vector3d const &v)
Definition: Vector3d.h:146
table::Key< int > a
std::ostream & operator<<(std::ostream &, Angle const &)
Definition: Angle.cc:34
Vector3d operator*(double s) const
The multiplication operator returns the component-wise product of this vector with scalar s...
Definition: Vector3d.h:116
Vector3d rotatedAround(UnitVector3d const &k, Angle a) const
rotatedAround returns a copy of this vector, rotated around the unit vector k by angle a according to...
Definition: Vector3d.cc:125
double operator()(int i) const
The function call operator returns the i-th component of this vector.
Definition: Vector3d.h:64
Vector3d is a vector in ℝ³ with components stored in double precision.
Definition: Vector3d.h:44
double y() const
Definition: Vector3d.h:68
bool operator==(Vector3d const &v) const
Definition: Vector3d.h:52
double x() const
Definition: Vector3d.h:66
A base class for image defects.
Vector3d & operator/=(double s)
Definition: Vector3d.h:145
bool isZero() const
isZero returns true if all the components of this vector are zero.
Definition: Vector3d.h:86
T fabs(T... args)
solver_t * s
Vector3d operator-() const
The unary minus operator negates every component of this vector.
Definition: Vector3d.h:108
Angle represents an angle in radians.
Definition: Angle.h:43
Vector3d operator+(Vector3d const &v) const
The addition operator returns the sum of this vector and v.
Definition: Vector3d.h:131
Vector3d cross(Vector3d const &v) const
cross returns the cross product of this vector and v.
Definition: Vector3d.h:101
Vector3d & operator*=(double s)
Definition: Vector3d.h:144
double z() const
Definition: Vector3d.h:70
Vector3d operator-(Vector3d const &v) const
The subtraction operator returns the difference between this vector and v.
Definition: Vector3d.h:138
Vector3d(double x, double y, double z)
This constructor creates a vector with the given components.
Definition: Vector3d.h:50
lsst::geom::Angle Angle
Definition: misc.h:33
Vector3d & operator-=(Vector3d const &v)
Definition: Vector3d.h:147
UnitVector3d is a unit vector in ℝ³ with components stored in double precision.
Definition: UnitVector3d.h:55
T sqrt(T... args)
double getSquaredNorm() const
getSquaredNorm returns the inner product of this vector with itself.
Definition: Vector3d.h:78
STL class.
Vector3d operator/(double s) const
The division operator returns the component-wise quotient of this vector with scalar s...
Definition: Vector3d.h:124
double const * getData() const
data returns a pointer to the 3 components of this vector.
Definition: Vector3d.h:61
double dot(Vector3d const &v) const
dot returns the inner product of this vector and v.
Definition: Vector3d.h:73
double getNorm() const
getNorm returns the L2 norm of this vector.
Definition: Vector3d.h:81