vec4d.h

Go to the documentation of this file.

//
// Copyright 2016 Pixar
//
// Licensed under the terms set forth in the LICENSE.txt file available at
// https://openusd.org/license.
//
// This file is generated by a script.  Do not edit directly.  Edit the
// vec.template.h file to make changes.
 
#ifndef PXR_BASE_GF_VEC4D_H
#define PXR_BASE_GF_VEC4D_H
 
 
#include "pxr/pxr.h"
#include "pxr/base/gf/api.h"
#include "pxr/base/gf/limits.h"
#include "pxr/base/gf/traits.h"
#include "pxr/base/gf/math.h"
#include "pxr/base/tf/diagnostic.h"
#include "pxr/base/tf/hash.h"
 
#include <cstddef>
#include <type_traits>
#include <cmath>
 
#include <iosfwd>
 
PXR_NAMESPACE_OPEN_SCOPE
 
class GfVec2d;
class GfVec3d;
 
template <>
struct GfIsGfVec<class GfVec4d> { static const bool value = true; };
 
class GfVec4d
{
public:
    typedef double ScalarType;
    static const size_t dimension = 4;
 
    double x, y, z, w;
 
    GfVec4d() = default;
 
    constexpr explicit GfVec4d(double value)
        : x(value), y(value), z(value), w(value)
    {
    }
 
    constexpr GfVec4d(double s0, double s1, double s2, double s3)
        : x(s0), y(s1), z(s2), w(s3)
    {
    }
 
    template <class Scl>
    constexpr explicit GfVec4d(Scl const *p)
        : x(p[0]), y(p[1]), z(p[2]), w(p[3])
    {
    }
 
    GfVec4d(class GfVec4f const &other);
 
    GfVec4d(class GfVec4h const &other);
 
    GfVec4d(class GfVec4i const &other);
    static GfVec4d XAxis() {
        GfVec4d result(0);
        result[0] = 1;
        return result;
    }
 
    static GfVec4d YAxis() {
        GfVec4d result(0);
        result[1] = 1;
        return result;
    }
 
    static GfVec4d ZAxis() {
        GfVec4d result(0);
        result[2] = 1;
        return result;
    }
 
    static GfVec4d WAxis() {
        GfVec4d result(0);
        result[3] = 1;
        return result;
    }
 
    static GfVec4d Axis(size_t i) {
        GfVec4d result(0);
        if (i < 4)
            result[i] = 1;
        return result;
    }
 
    GfVec4d &Set(double s0, double s1, double s2, double s3) {
        x = s0;
        y = s1;
        z = s2;
        w = s3;
        return *this;
    }
 
    GfVec4d &Set(double const *a) {
        return Set(a[0], a[1], a[2], a[3]);
    }
 
    double const *data() const { return &x; }
    double *data() { return &x; }
    double const *GetArray() const { return data(); }
 
    double const &operator[](size_t i) const { return data()[i]; }
    double &operator[](size_t i) { return data()[i]; }
 
    friend inline size_t hash_value(GfVec4d const &vec) {
        return TfHash::Combine(vec.x, vec.y, vec.z, vec.w);
    }
 
    bool operator==(GfVec4d const &other) const {
        return x == other.x &&
               y == other.y &&
               z == other.z &&
               w == other.w;
    }
    bool operator!=(GfVec4d const &other) const {
        return !(*this == other);
    }
 
    // TODO Add inequality for other vec types...
    GF_API
    bool operator==(class GfVec4f const &other) const;
    GF_API
    bool operator==(class GfVec4h const &other) const;
    GF_API
    bool operator==(class GfVec4i const &other) const;
    
    GfVec4d operator-() const {
        return GfVec4d(-x, -y, -z, -w);
    }
 
    GfVec4d &operator+=(GfVec4d const &other) {
        x += other.x;
        y += other.y;
        z += other.z;
        w += other.w;
        return *this;
    }
    friend GfVec4d operator+(GfVec4d const &l, GfVec4d const &r) {
        return GfVec4d(l) += r;
    }
 
    GfVec4d &operator-=(GfVec4d const &other) {
        x -= other.x;
        y -= other.y;
        z -= other.z;
        w -= other.w;
        return *this;
    }
    friend GfVec4d operator-(GfVec4d const &l, GfVec4d const &r) {
        return GfVec4d(l) -= r;
    }
 
    GfVec4d &operator*=(double s) {
        x *= s;
        y *= s;
        z *= s;
        w *= s;
        return *this;
    }
    GfVec4d operator*(double s) const {
        return GfVec4d(*this) *= s;
    }
    friend GfVec4d operator*(double s, GfVec4d const &v) {
        return v * s;
    }
 
    // TODO should divide by the scalar type.
    GfVec4d &operator/=(double s) {
        // TODO This should not multiply by 1/s, it should do the division.
        // Doing the division is more numerically stable when s is close to
        // zero.
        return *this *= (1.0 / s);
    }
    GfVec4d operator/(double s) const {
        return *this * (1.0 / s);
    }
    
    double operator*(GfVec4d const &v) const {
        return x * v.x + y * v.y + z * v.z + w * v.w;
    }
 
    GfVec4d GetProjection(GfVec4d const &v) const {
        return v * (*this * v);
    }
 
    GfVec4d GetComplement(GfVec4d const &b) const {
        return *this - this->GetProjection(b);
    }
 
    double GetLengthSq() const {
        return *this * *this;
    }
 
    double GetLength() const {
        return GfSqrt(GetLengthSq());
    }
 
    double Normalize(double eps = GF_MIN_VECTOR_LENGTH) {
        // TODO this seems suspect...  suggest dividing by length so long as
        // length is not zero.
        double length = GetLength();
        *this /= (length > eps) ? length : eps;
        return length;
    }
 
    GfVec4d GetNormalized(double eps = GF_MIN_VECTOR_LENGTH) const {
        GfVec4d normalized(*this);
        normalized.Normalize(eps);
        return normalized;
    }
 
  
#if defined (doxygen)
    constexpr GfVec2d xy() const;
    constexpr GfVec2d xz() const;
    constexpr GfVec2d xw() const;
    constexpr GfVec2d yx() const;
    constexpr GfVec2d yz() const;
    constexpr GfVec2d yw() const;
    constexpr GfVec2d zx() const;
    constexpr GfVec2d zy() const;
    constexpr GfVec2d zw() const;
    constexpr GfVec2d wx() const;
    constexpr GfVec2d wy() const;
    constexpr GfVec2d wz() const;
    constexpr GfVec3d xyz() const;
    constexpr GfVec3d xyw() const;
    constexpr GfVec3d xzy() const;
    constexpr GfVec3d xzw() const;
    constexpr GfVec3d xwy() const;
    constexpr GfVec3d xwz() const;
    constexpr GfVec3d yxz() const;
    constexpr GfVec3d yxw() const;
    constexpr GfVec3d yzx() const;
    constexpr GfVec3d yzw() const;
    constexpr GfVec3d ywx() const;
    constexpr GfVec3d ywz() const;
    constexpr GfVec3d zxy() const;
    constexpr GfVec3d zxw() const;
    constexpr GfVec3d zyx() const;
    constexpr GfVec3d zyw() const;
    constexpr GfVec3d zwx() const;
    constexpr GfVec3d zwy() const;
    constexpr GfVec3d wxy() const;
    constexpr GfVec3d wxz() const;
    constexpr GfVec3d wyx() const;
    constexpr GfVec3d wyz() const;
    constexpr GfVec3d wzx() const;
    constexpr GfVec3d wzy() const;
    constexpr GfVec4d xyzw() const;
    constexpr GfVec4d xywz() const;
    constexpr GfVec4d xzyw() const;
    constexpr GfVec4d xzwy() const;
    constexpr GfVec4d xwyz() const;
    constexpr GfVec4d xwzy() const;
    constexpr GfVec4d yxzw() const;
    constexpr GfVec4d yxwz() const;
    constexpr GfVec4d yzxw() const;
    constexpr GfVec4d yzwx() const;
    constexpr GfVec4d ywxz() const;
    constexpr GfVec4d ywzx() const;
    constexpr GfVec4d zxyw() const;
    constexpr GfVec4d zxwy() const;
    constexpr GfVec4d zyxw() const;
    constexpr GfVec4d zywx() const;
    constexpr GfVec4d zwxy() const;
    constexpr GfVec4d zwyx() const;
    constexpr GfVec4d wxyz() const;
    constexpr GfVec4d wxzy() const;
    constexpr GfVec4d wyxz() const;
    constexpr GfVec4d wyzx() const;
    constexpr GfVec4d wzxy() const;
    constexpr GfVec4d wzyx() const;
#else // !doxygen
    // Note that these are templated on their own return type so that the
    // declarations cost nothing in compile time & debug info unless they are
    // actually called in a TU.  Naming a different type is an error and
    // static_assert()ed against.
    template <class Vec = GfVec2d> constexpr Vec xy() const;
    template <class Vec = GfVec2d> constexpr Vec xz() const;
    template <class Vec = GfVec2d> constexpr Vec xw() const;
    template <class Vec = GfVec2d> constexpr Vec yx() const;
    template <class Vec = GfVec2d> constexpr Vec yz() const;
    template <class Vec = GfVec2d> constexpr Vec yw() const;
    template <class Vec = GfVec2d> constexpr Vec zx() const;
    template <class Vec = GfVec2d> constexpr Vec zy() const;
    template <class Vec = GfVec2d> constexpr Vec zw() const;
    template <class Vec = GfVec2d> constexpr Vec wx() const;
    template <class Vec = GfVec2d> constexpr Vec wy() const;
    template <class Vec = GfVec2d> constexpr Vec wz() const;
    template <class Vec = GfVec3d> constexpr Vec xyz() const;
    template <class Vec = GfVec3d> constexpr Vec xyw() const;
    template <class Vec = GfVec3d> constexpr Vec xzy() const;
    template <class Vec = GfVec3d> constexpr Vec xzw() const;
    template <class Vec = GfVec3d> constexpr Vec xwy() const;
    template <class Vec = GfVec3d> constexpr Vec xwz() const;
    template <class Vec = GfVec3d> constexpr Vec yxz() const;
    template <class Vec = GfVec3d> constexpr Vec yxw() const;
    template <class Vec = GfVec3d> constexpr Vec yzx() const;
    template <class Vec = GfVec3d> constexpr Vec yzw() const;
    template <class Vec = GfVec3d> constexpr Vec ywx() const;
    template <class Vec = GfVec3d> constexpr Vec ywz() const;
    template <class Vec = GfVec3d> constexpr Vec zxy() const;
    template <class Vec = GfVec3d> constexpr Vec zxw() const;
    template <class Vec = GfVec3d> constexpr Vec zyx() const;
    template <class Vec = GfVec3d> constexpr Vec zyw() const;
    template <class Vec = GfVec3d> constexpr Vec zwx() const;
    template <class Vec = GfVec3d> constexpr Vec zwy() const;
    template <class Vec = GfVec3d> constexpr Vec wxy() const;
    template <class Vec = GfVec3d> constexpr Vec wxz() const;
    template <class Vec = GfVec3d> constexpr Vec wyx() const;
    template <class Vec = GfVec3d> constexpr Vec wyz() const;
    template <class Vec = GfVec3d> constexpr Vec wzx() const;
    template <class Vec = GfVec3d> constexpr Vec wzy() const;
    template <class Vec = GfVec4d> constexpr Vec xyzw() const;
    template <class Vec = GfVec4d> constexpr Vec xywz() const;
    template <class Vec = GfVec4d> constexpr Vec xzyw() const;
    template <class Vec = GfVec4d> constexpr Vec xzwy() const;
    template <class Vec = GfVec4d> constexpr Vec xwyz() const;
    template <class Vec = GfVec4d> constexpr Vec xwzy() const;
    template <class Vec = GfVec4d> constexpr Vec yxzw() const;
    template <class Vec = GfVec4d> constexpr Vec yxwz() const;
    template <class Vec = GfVec4d> constexpr Vec yzxw() const;
    template <class Vec = GfVec4d> constexpr Vec yzwx() const;
    template <class Vec = GfVec4d> constexpr Vec ywxz() const;
    template <class Vec = GfVec4d> constexpr Vec ywzx() const;
    template <class Vec = GfVec4d> constexpr Vec zxyw() const;
    template <class Vec = GfVec4d> constexpr Vec zxwy() const;
    template <class Vec = GfVec4d> constexpr Vec zyxw() const;
    template <class Vec = GfVec4d> constexpr Vec zywx() const;
    template <class Vec = GfVec4d> constexpr Vec zwxy() const;
    template <class Vec = GfVec4d> constexpr Vec zwyx() const;
    template <class Vec = GfVec4d> constexpr Vec wxyz() const;
    template <class Vec = GfVec4d> constexpr Vec wxzy() const;
    template <class Vec = GfVec4d> constexpr Vec wyxz() const;
    template <class Vec = GfVec4d> constexpr Vec wyzx() const;
    template <class Vec = GfVec4d> constexpr Vec wzxy() const;
    template <class Vec = GfVec4d> constexpr Vec wzyx() const;
#endif // doxygen
 };
 
GF_API std::ostream& operator<<(std::ostream &, GfVec4d const &);
 
PXR_NAMESPACE_CLOSE_SCOPE
 
// Other scalar types of the same dimension, for the conversions above.
#include "pxr/base/gf/vec4f.h"
#include "pxr/base/gf/vec4h.h"
#include "pxr/base/gf/vec4i.h"
// Lower dimensions of the same scalar type, for the swizzles above.  Note that
// lower dimensions never include higher ones, so this introduces no cycle.
#include "pxr/base/gf/vec2d.h"
#include "pxr/base/gf/vec3d.h"
 
PXR_NAMESPACE_OPEN_SCOPE
 
inline
GfVec4d::GfVec4d(class GfVec4f const &other)
{
    x = other.x;
    y = other.y;
    z = other.z;
    w = other.w;
}
inline
GfVec4d::GfVec4d(class GfVec4h const &other)
{
    x = other.x;
    y = other.y;
    z = other.z;
    w = other.w;
}
inline
GfVec4d::GfVec4d(class GfVec4i const &other)
{
    x = other.x;
    y = other.y;
    z = other.z;
    w = other.w;
}
 
template <class Vec>
constexpr Vec
GfVec4d::xy() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(x, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xz() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(x, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xw() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(x, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yx() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(y, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yz() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(y, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yw() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(y, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zx() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(z, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zy() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(z, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zw() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(z, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wx() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(w, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wy() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(w, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wz() const
{
    static_assert(std::is_same<Vec, GfVec2d>::value);
    return Vec(w, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xyz() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(x, y, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xyw() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(x, y, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xzy() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(x, z, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xzw() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(x, z, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xwy() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(x, w, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xwz() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(x, w, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yxz() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(y, x, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yxw() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(y, x, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yzx() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(y, z, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yzw() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(y, z, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::ywx() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(y, w, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::ywz() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(y, w, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zxy() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(z, x, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zxw() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(z, x, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zyx() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(z, y, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zyw() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(z, y, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zwx() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(z, w, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zwy() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(z, w, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wxy() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(w, x, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wxz() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(w, x, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wyx() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(w, y, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wyz() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(w, y, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wzx() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(w, z, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wzy() const
{
    static_assert(std::is_same<Vec, GfVec3d>::value);
    return Vec(w, z, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xyzw() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(x, y, z, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xywz() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(x, y, w, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xzyw() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(x, z, y, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xzwy() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(x, z, w, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xwyz() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(x, w, y, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::xwzy() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(x, w, z, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yxzw() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(y, x, z, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yxwz() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(y, x, w, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yzxw() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(y, z, x, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::yzwx() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(y, z, w, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::ywxz() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(y, w, x, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::ywzx() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(y, w, z, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zxyw() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(z, x, y, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zxwy() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(z, x, w, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zyxw() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(z, y, x, w);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zywx() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(z, y, w, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zwxy() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(z, w, x, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::zwyx() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(z, w, y, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wxyz() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(w, x, y, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wxzy() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(w, x, z, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wyxz() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(w, y, x, z);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wyzx() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(w, y, z, x);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wzxy() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(w, z, x, y);
}
 
template <class Vec>
constexpr Vec
GfVec4d::wzyx() const
{
    static_assert(std::is_same<Vec, GfVec4d>::value);
    return Vec(w, z, y, x);
}
 
 
inline GfVec4d
GfCompMult(GfVec4d const &v1, GfVec4d const &v2) {
    return GfVec4d(
        v1.x * v2.x,
        v1.y * v2.y,
        v1.z * v2.z,
        v1.w * v2.w
        );
}
 
inline GfVec4d
GfCompDiv(GfVec4d const &v1, GfVec4d const &v2) {
    return GfVec4d(
        v1.x / v2.x,
        v1.y / v2.y,
        v1.z / v2.z,
        v1.w / v2.w
        );
}
 
inline double
GfDot(GfVec4d const &v1, GfVec4d const &v2) {
    return v1 * v2;
}
 
 
inline double
GfGetLength(GfVec4d const &v)
{
    return v.GetLength();
}
 
inline double
GfNormalize(GfVec4d *v, double eps = GF_MIN_VECTOR_LENGTH)
{
    return v->Normalize(eps);
}
 
inline GfVec4d
GfGetNormalized(GfVec4d const &v, double eps = GF_MIN_VECTOR_LENGTH)
{
    return v.GetNormalized(eps);
}
 
inline GfVec4d
GfGetProjection(GfVec4d const &a, GfVec4d const &b)
{
    return a.GetProjection(b);
}
 
inline GfVec4d
GfGetComplement(GfVec4d const &a, GfVec4d const &b)
{
    return a.GetComplement(b);
}
 
inline bool
GfIsClose(GfVec4d const &v1, GfVec4d const &v2, double tolerance)
{
    GfVec4d delta = v1 - v2;
    return delta.GetLengthSq() <= tolerance * tolerance;
}
 
 
 
PXR_NAMESPACE_CLOSE_SCOPE
 
#endif // PXR_BASE_GF_VEC4D_H