lsst::sphgeom::UnitVector3d Class Reference

UnitVector3d is a unit vector in ℝ³ with components stored in double precision. More...

#include <UnitVector3d.h>

Public Member Functions
	UnitVector3d ()
	The default constructor creates a unit vector equal to (1, 0, 0).
	UnitVector3d (UnitVector3d const &v)=default
	UnitVector3d (LonLat const &p)
	This constructor creates the unit vector corresponding to the point p on the unit sphere.
	UnitVector3d (Angle lon, Angle lat)
	This constructor creates a unit vector corresponding to the given spherical coordinates.
	operator Vector3d const & () const
	This conversion operator returns a const reference to the underlying Vector3d.
bool	operator== (Vector3d const &v) const
bool	operator!= (Vector3d const &v) const
double const *	getData () const
	getData returns a pointer to the 3 components of this unit vector.
double	operator() (int i) const
	The function call operator returns the i-th component of this vector.
double	x () const
double	y () const
double	z () const
double	dot (Vector3d const &v) const
	dot returns the inner product of this unit vector and v.
Vector3d	cross (Vector3d const &v) const
	cross returns the cross product of this unit vector and v.
Vector3d	robustCross (UnitVector3d const &v) const
	a.robustCross(b) is (b + a).cross(b - a) - twice the cross product of a and b.
UnitVector3d	operator- () const
	The unary minus operator negates every component of this unit vector.
Vector3d	operator* (double s) const
	The multiplication operator returns the component-wise product of this unit vector with scalar s.
Vector3d	operator/ (double s) const
	The division operator returns the component-wise quotient of this unit vector with scalar s.
Vector3d	operator+ (Vector3d const &v) const
	The addition operator returns the sum of this unit vector and v.
Vector3d	operator- (Vector3d const &v) const
	The subtraction operator returns the difference between this unit vector and v.
Vector3d	cwiseProduct (Vector3d const &v) const
	cwiseProduct returns the component-wise product of this unit vector and v.
UnitVector3d	rotatedAround (UnitVector3d const &k, Angle a) const
	rotatedAround returns a copy of this unit vector, rotated around the unit vector k by angle a according to the right hand rule.
	UnitVector3d (Vector3d const &v)
	UnitVector3d (double x, double y, double z)

Static Public Member Functions
static UnitVector3d	orthogonalTo (Vector3d const &v)
	orthogonalTo returns an arbitrary unit vector that is orthogonal to v.
static UnitVector3d	orthogonalTo (Vector3d const &v1, Vector3d const &v2)
	orthogonalTo returns the unit vector n orthogonal to both v1 and v2 and with the same direction as the cross product of v1 and v2.
static UnitVector3d	northFrom (Vector3d const &v)
	northFrom returns the unit vector orthogonal to v that points "north" from v.
static UnitVector3d	orthogonalTo (NormalizedAngle const &a)
	orthogonalTo returns the unit vector orthogonal to the meridian with the given longitude.
static UnitVector3d	fromNormalized (Vector3d const &v)
	fromNormalized returns the unit vector equal to v, which is assumed to be normalized.
static UnitVector3d	fromNormalized (double x, double y, double z)
	fromNormalized returns the unit vector with the given components, which are assumed to correspond to those of a normalized vector.
static UnitVector3d	X ()
static UnitVector3d	Y ()
static UnitVector3d	Z ()

Detailed Description

UnitVector3d is a unit vector in ℝ³ with components stored in double precision.

Conceptually, a UnitVector3d is a Vector3d. However, it does not inherit from Vector3d because all its mutators must correspond to isometries. Propagating this restriction to the Vector3d base class would make the base class unduly difficult to use. Furthermore, the base class assignment operator would have to be made virtual - otherwise it would be possible to cast a UnitVector3d reference to a Vector3d reference and then assign components yielding non-unit norm. For a class this compact and performance critical, the addition of a vtable pointer per instance and the potential for virtual call overhead on assignment is deemed prohibitive.

This class supports conversion from longitude and latitude angles. All methods that accept a LonLat or equivalent angles shall convert from longitude/latitude to a unit vector according to the following conventions:

• (0°, 0°) → (1, 0, 0)
• (90°, 0°) → (0, 1, 0)
• (0°, +90°) → (0, 0, 1)

Definition at line 62 of file UnitVector3d.h.

Constructor & Destructor Documentation

UnitVector3d() [1/6]

lsst::sphgeom::UnitVector3d::UnitVector3d ( )

inline

The default constructor creates a unit vector equal to (1, 0, 0).

Definition at line 113 of file UnitVector3d.h.

: _v(1.0, 0.0, 0.0) {}

UnitVector3d() [2/6]

lsst::sphgeom::UnitVector3d::UnitVector3d ( UnitVector3d const & v )

default

UnitVector3d() [3/6]

lsst::sphgeom::UnitVector3d::UnitVector3d ( Vector3d const & v )

inlineexplicit

This constructor creates a unit vector with the given direction.

Definition at line 119 of file UnitVector3d.h.

                                              : _v(v) {
        _v.normalize();
    }

UnitVector3d() [4/6]

lsst::sphgeom::UnitVector3d::UnitVector3d ( double x, double y, double z )

inline

This constructor creates a unit vector with the given direction.

Definition at line 123 of file UnitVector3d.h.

                                               : _v(x, y, z) {
        _v.normalize();
    }

UnitVector3d() [5/6]

lsst::sphgeom::UnitVector3d::UnitVector3d ( LonLat const & p )

inlineexplicit

This constructor creates the unit vector corresponding to the point p on the unit sphere.

Definition at line 130 of file UnitVector3d.h.

                                            {
        *this = UnitVector3d(p.getLon(), p.getLat());
    }

UnitVector3d() [6/6]

lsst::sphgeom::UnitVector3d::UnitVector3d	(	Angle	lon,
		Angle	lat )

This constructor creates a unit vector corresponding to the given spherical coordinates.

Definition at line 70 of file UnitVector3d.cc.

                                               {
    double sinLon = sin(lon);
    double cosLon = cos(lon);
    double sinLat = sin(lat);
    double cosLat = cos(lat);
    _v = Vector3d(cosLon * cosLat,
                  sinLon * cosLat,
                  sinLat);
}

Member Function Documentation

cross()

Vector3d lsst::sphgeom::UnitVector3d::cross ( Vector3d const & v ) const

inline

cross returns the cross product of this unit vector and v.

Definition at line 162 of file UnitVector3d.h.

{ return _v.cross(v); }

cwiseProduct()

Vector3d lsst::sphgeom::UnitVector3d::cwiseProduct ( Vector3d const & v ) const

inline

cwiseProduct returns the component-wise product of this unit vector and v.

Definition at line 194 of file UnitVector3d.h.

                                                    {
        return _v.cwiseProduct(v);
    }

dot()

double lsst::sphgeom::UnitVector3d::dot ( Vector3d const & v ) const

inline

dot returns the inner product of this unit vector and v.

Definition at line 159 of file UnitVector3d.h.

{ return _v.dot(v); }

fromNormalized() [1/2]

static UnitVector3d lsst::sphgeom::UnitVector3d::fromNormalized ( double x, double y, double z )

inlinestatic

fromNormalized returns the unit vector with the given components, which are assumed to correspond to those of a normalized vector.

Use with caution - this assumption is not verified!

Definition at line 96 of file UnitVector3d.h.

                                                                     {
        UnitVector3d u; u._v = Vector3d(x, y, z); return u;
    }

fromNormalized() [2/2]

static UnitVector3d lsst::sphgeom::UnitVector3d::fromNormalized ( Vector3d const & v )

inlinestatic

fromNormalized returns the unit vector equal to v, which is assumed to be normalized.

Use with caution - this assumption is not verified!

Definition at line 89 of file UnitVector3d.h.

                                                           {
        UnitVector3d u; u._v = v; return u;
    }

getData()

double const * lsst::sphgeom::UnitVector3d::getData ( ) const

inline

getData returns a pointer to the 3 components of this unit vector.

Definition at line 147 of file UnitVector3d.h.

{ return _v.getData(); }

northFrom()

UnitVector3d lsst::sphgeom::UnitVector3d::northFrom ( Vector3d const & v )

static

northFrom returns the unit vector orthogonal to v that points "north" from v.

More precisely, it returns a unit vector orthogonal to v that lies in the plane defined by v, the north pole (0, 0, 1), and the origin. If v is colinear with the z axis, then (-1, 0, 0) is returned if v points north, and (1, 0, 0) is returned if v points south.

Definition at line 58 of file UnitVector3d.cc.

                                                       {
    Vector3d n(-v.x() * v.z(),
               -v.y() * v.z(),
               v.x() * v.x() + v.y() * v.y());
    if (n.isZero()) {
        UnitVector3d u;
        u._v = Vector3d(-::copysign(1.0, v.z()), 0.0, 0.0);
        return u;
    }
    return UnitVector3d(n);
}

operator Vector3d const &()

lsst::sphgeom::UnitVector3d::operator Vector3d const & ( ) const

inline

This conversion operator returns a const reference to the underlying Vector3d.

It allows a UnitVector3d to transparently replace a Vector3d as an argument in most function calls.

Definition at line 141 of file UnitVector3d.h.

{ return _v; }

operator!=()

bool lsst::sphgeom::UnitVector3d::operator!= ( Vector3d const & v ) const

inline

Definition at line 144 of file UnitVector3d.h.

{ return _v != v; }

operator()()

double lsst::sphgeom::UnitVector3d::operator() ( int i ) const

inline

The function call operator returns the i-th component of this vector.

Definition at line 150 of file UnitVector3d.h.

{ return _v(i); }

operator*()

Vector3d lsst::sphgeom::UnitVector3d::operator* ( double s ) const

inline

The multiplication operator returns the component-wise product of this unit vector with scalar s.

Definition at line 179 of file UnitVector3d.h.

{ return _v * s; }

operator+()

Vector3d lsst::sphgeom::UnitVector3d::operator+ ( Vector3d const & v ) const

inline

The addition operator returns the sum of this unit vector and v.

Definition at line 186 of file UnitVector3d.h.

{ return _v + v; }

operator-() [1/2]

UnitVector3d lsst::sphgeom::UnitVector3d::operator- ( ) const

inline

The unary minus operator negates every component of this unit vector.

Definition at line 173 of file UnitVector3d.h.

                                   {
        UnitVector3d u; u._v = -_v; return u;
    }

operator-() [2/2]

Vector3d lsst::sphgeom::UnitVector3d::operator- ( Vector3d const & v ) const

inline

The subtraction operator returns the difference between this unit vector and v.

Definition at line 190 of file UnitVector3d.h.

{ return _v - v; }

operator/()

Vector3d lsst::sphgeom::UnitVector3d::operator/ ( double s ) const

inline

The division operator returns the component-wise quotient of this unit vector with scalar s.

Definition at line 183 of file UnitVector3d.h.

{ return _v / s; }

operator==()

bool lsst::sphgeom::UnitVector3d::operator== ( Vector3d const & v ) const

inline

Definition at line 143 of file UnitVector3d.h.

{ return _v == v; }

orthogonalTo() [1/3]

static UnitVector3d lsst::sphgeom::UnitVector3d::orthogonalTo ( NormalizedAngle const & a )

inlinestatic

orthogonalTo returns the unit vector orthogonal to the meridian with the given longitude.

Definition at line 83 of file UnitVector3d.h.

                                                                {
        UnitVector3d u; u._v = Vector3d(-sin(a), cos(a), 0.0); return u;
    }

orthogonalTo() [2/3]

UnitVector3d lsst::sphgeom::UnitVector3d::orthogonalTo ( Vector3d const & v )

static

orthogonalTo returns an arbitrary unit vector that is orthogonal to v.

Definition at line 41 of file UnitVector3d.cc.

                                                          {
    if (std::fabs(v.x()) > std::fabs(v.y())) {
        return UnitVector3d(-v.z(), 0.0, v.x());
    }
    return UnitVector3d(0.0, v.z(), -v.y());
}

orthogonalTo() [3/3]

UnitVector3d lsst::sphgeom::UnitVector3d::orthogonalTo ( Vector3d const & v1, Vector3d const & v2 )

static

orthogonalTo returns the unit vector n orthogonal to both v1 and v2 and with the same direction as the cross product of v1 and v2.

If v1 and v2 are nearly (anti)parallel, then an arbitrary unit vector orthogonal to v1 is returned.

Definition at line 48 of file UnitVector3d.cc.

{
    Vector3d n = (v2 + v1).cross(v2 - v1);
    if (n.isZero()) {
        return orthogonalTo(v1);
    }
    return UnitVector3d(n);
}

robustCross()

Vector3d lsst::sphgeom::UnitVector3d::robustCross ( UnitVector3d const & v ) const

inline

a.robustCross(b) is (b + a).cross(b - a) - twice the cross product of a and b.

The result is almost orthogonal to a and b even for nearly (anti-)parallel unit vectors. The idea comes from the Google S2 library.

Definition at line 168 of file UnitVector3d.h.

                                                       {
        return (v + *this).cross(v - *this);
    }

rotatedAround()

UnitVector3d lsst::sphgeom::UnitVector3d::rotatedAround ( UnitVector3d const & k, Angle a ) const

inline

rotatedAround returns a copy of this unit vector, rotated around the unit vector k by angle a according to the right hand rule.

Definition at line 200 of file UnitVector3d.h.

                                                                      {
        return UnitVector3d(_v.rotatedAround(k, a));
    }

X()

static UnitVector3d lsst::sphgeom::UnitVector3d::X ( )

inlinestatic

Definition at line 100 of file UnitVector3d.h.

                            {
        return UnitVector3d();
    }

x()

double lsst::sphgeom::UnitVector3d::x ( ) const

inline

Definition at line 152 of file UnitVector3d.h.

{ return _v.x(); }

Y()

static UnitVector3d lsst::sphgeom::UnitVector3d::Y ( )

inlinestatic

Definition at line 104 of file UnitVector3d.h.

                            {
        UnitVector3d u; u._v = Vector3d(0.0, 1.0, 0.0); return u;
    }

y()

double lsst::sphgeom::UnitVector3d::y ( ) const

inline

Definition at line 154 of file UnitVector3d.h.

{ return _v.y(); }

Z()

static UnitVector3d lsst::sphgeom::UnitVector3d::Z ( )

inlinestatic

Definition at line 108 of file UnitVector3d.h.

                            {
        UnitVector3d u; u._v = Vector3d(0.0, 0.0, 1.0); return u;
    }

z()

double lsst::sphgeom::UnitVector3d::z ( ) const

inline

Definition at line 156 of file UnitVector3d.h.

{ return _v.z(); }

---

The documentation for this class was generated from the following files:

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-9.0.0/Linux64/sphgeom/g781aacb6e4+a8ce1bb630/include/lsst/sphgeom/UnitVector3d.h
• /j/snowflake/release/lsstsw/stack/lsst-scipipe-9.0.0/Linux64/sphgeom/g781aacb6e4+a8ce1bb630/src/UnitVector3d.cc