LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
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
std::ostream * os
Definition: Schema.cc:746
This file declares a class for representing unit vectors in ℝ³.
Angle represents an angle in radians.
Definition: Angle.h:43
UnitVector3d is a unit vector in ℝ³ with components stored in double precision.
Definition: UnitVector3d.h:55
UnitVector3d()
The default constructor creates a unit vector equal to (1, 0, 0).
Definition: UnitVector3d.h:106
static UnitVector3d orthogonalTo(Vector3d const &v)
orthogonalTo returns an arbitrary unit vector that is orthogonal to v.
Definition: UnitVector3d.cc:34
Vector3d cross(Vector3d const &v) const
cross returns the cross product of this unit vector and v.
Definition: UnitVector3d.h:155
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 x() const
Definition: Vector3d.h:66
bool isZero() const
isZero returns true if all the components of this vector are zero.
Definition: Vector3d.h:86
double y() const
Definition: Vector3d.h:68
double z() const
Definition: Vector3d.h:70
T fabs(T... args)
std::ostream & operator<<(std::ostream &, Angle const &)
Definition: Angle.cc:34
double sin(Angle const &a)
Definition: Angle.h:102
double cos(Angle const &a)
Definition: Angle.h:103
A base class for image defects.