LSSTApplications  17.0+124,17.0+14,17.0+73,18.0.0+37,18.0.0+80,18.0.0-4-g68ffd23+4,18.1.0-1-g0001055+12,18.1.0-1-g03d53ef+5,18.1.0-1-g1349e88+55,18.1.0-1-g2505f39+44,18.1.0-1-g5315e5e+4,18.1.0-1-g5e4b7ea+14,18.1.0-1-g7e8fceb+4,18.1.0-1-g85f8cd4+48,18.1.0-1-g8ff0b9f+4,18.1.0-1-ga2c679d+1,18.1.0-1-gd55f500+35,18.1.0-10-gb58edde+2,18.1.0-11-g0997b02+4,18.1.0-13-gfe4edf0b+12,18.1.0-14-g259bd21+21,18.1.0-19-gdb69f3f+2,18.1.0-2-g5f9922c+24,18.1.0-2-gd3b74e5+11,18.1.0-2-gfbf3545+32,18.1.0-26-g728bddb4+5,18.1.0-27-g6ff7ca9+2,18.1.0-3-g52aa583+25,18.1.0-3-g8ea57af+9,18.1.0-3-gb69f684+42,18.1.0-3-gfcaddf3+6,18.1.0-32-gd8786685a,18.1.0-4-gf3f9b77+6,18.1.0-5-g1dd662b+2,18.1.0-5-g6dbcb01+41,18.1.0-6-gae77429+3,18.1.0-7-g9d75d83+9,18.1.0-7-gae09a6d+30,18.1.0-9-gc381ef5+4,w.2019.45
LSSTDataManagementBasePackage
UnitVector3d.cc
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 
25 
27 
28 #include <ostream>
29 
30 
31 namespace lsst {
32 namespace sphgeom {
33 
35  if (std::fabs(v.x()) > std::fabs(v.y())) {
36  return UnitVector3d(-v.z(), 0.0, v.x());
37  }
38  return UnitVector3d(0.0, v.z(), -v.y());
39 }
40 
42  Vector3d const & v2)
43 {
44  Vector3d n = (v2 + v1).cross(v2 - v1);
45  if (n.isZero()) {
46  return orthogonalTo(v1);
47  }
48  return UnitVector3d(n);
49 }
50 
52  Vector3d n(-v.x() * v.z(),
53  -v.y() * v.z(),
54  v.x() * v.x() + v.y() * v.y());
55  if (n.isZero()) {
56  UnitVector3d u;
57  u._v = Vector3d(-::copysign(1.0, v.z()), 0.0, 0.0);
58  return u;
59  }
60  return UnitVector3d(n);
61 }
62 
64  double sinLon = sin(lon);
65  double cosLon = cos(lon);
66  double sinLat = sin(lat);
67  double cosLat = cos(lat);
68  _v = Vector3d(cosLon * cosLat,
69  sinLon * cosLat,
70  sinLat);
71 }
72 
74  return os << static_cast<Vector3d const &>(v);
75 }
76 
77 }} // namespace lsst::sphgeom
static UnitVector3d orthogonalTo(Vector3d const &v)
orthogonalTo returns an arbitrary unit vector that is orthogonal to v.
Definition: UnitVector3d.cc:34
std::ostream & operator<<(std::ostream &, Angle const &)
Definition: Angle.cc:34
static UnitVector3d northFrom(Vector3d const &v)
northFrom returns the unit vector orthogonal to v that points "north" from v.
Definition: UnitVector3d.cc:51
Vector3d is a vector in ℝ³ with components stored in double precision.
Definition: Vector3d.h:44
double sin(Angle const &a)
Definition: Angle.h:102
double y() const
Definition: Vector3d.h:68
double cos(Angle const &a)
Definition: Angle.h:103
This file declares a class for representing unit vectors in ℝ³.
Vector3d cross(Vector3d const &v) const
cross returns the cross product of this unit vector and v.
Definition: UnitVector3d.h:155
double x() const
Definition: Vector3d.h:66
A base class for image defects.
bool isZero() const
isZero returns true if all the components of this vector are zero.
Definition: Vector3d.h:86
T fabs(T... args)
Angle represents an angle in radians.
Definition: Angle.h:43
double z() const
Definition: Vector3d.h:70
UnitVector3d is a unit vector in ℝ³ with components stored in double precision.
Definition: UnitVector3d.h:55
std::ostream * os
Definition: Schema.cc:746
STL class.
UnitVector3d()
The default constructor creates a unit vector equal to (1, 0, 0).
Definition: UnitVector3d.h:106