vec4d.h

Go to the documentation of this file.

//
// Copyright 2016 Pixar
//
// Licensed under the Apache License, Version 2.0 (the "Apache License")
// with the following modification; you may not use this file except in
// compliance with the Apache License and the following modification to it:
// Section 6. Trademarks. is deleted and replaced with:
//
// 6. Trademarks. This License does not grant permission to use the trade
//    names, trademarks, service marks, or product names of the Licensor
//    and its affiliates, except as required to comply with Section 4(c) of
//    the License and to reproduce the content of the NOTICE file.
//
// You may obtain a copy of the Apache License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the Apache License with the above modification is
// distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
// KIND, either express or implied. See the Apache License for the specific
// language governing permissions and limitations under the Apache 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/tf/diagnostic.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/hash.h"

#include <cstddef>
#include <cmath>

#include <iosfwd>

PXR_NAMESPACE_OPEN_SCOPE

class GfVec4d;

template <>
struct GfIsGfVec<class GfVec4d> { static const bool value = true; };

class GfVec4d
{
public:
    typedef double ScalarType;
    static const size_t dimension = 4;

    GfVec4d() = default;

    constexpr explicit GfVec4d(double value)
        : _data{ value, value, value, value }
    {
    }

    constexpr GfVec4d(double s0, double s1, double s2, double s3)
        : _data{ s0, s1, s2, s3 }
    {
    }

    template <class Scl>
    constexpr explicit GfVec4d(Scl const *p)
        : _data{ p[0], p[1], p[2], 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) {
        _data[0] = s0;
        _data[1] = s1;
        _data[2] = s2;
        _data[3] = s3;
        return *this;
    }

    GfVec4d &Set(double const *a) {
        return Set(a[0], a[1], a[2], a[3]);
    }

    double const *data() const { return _data; }
    double *data() { return _data; }
    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[0], vec[1], vec[2], vec[3]);
    }

    bool operator==(GfVec4d const &other) const {
        return _data[0] == other[0] &&
               _data[1] == other[1] &&
               _data[2] == other[2] &&
               _data[3] == other[3];
    }
    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(-_data[0], -_data[1], -_data[2], -_data[3]);
    }

    GfVec4d &operator+=(GfVec4d const &other) {
        _data[0] += other[0];
        _data[1] += other[1];
        _data[2] += other[2];
        _data[3] += other[3];
        return *this;
    }
    friend GfVec4d operator+(GfVec4d const &l, GfVec4d const &r) {
        return GfVec4d(l) += r;
    }

    GfVec4d &operator-=(GfVec4d const &other) {
        _data[0] -= other[0];
        _data[1] -= other[1];
        _data[2] -= other[2];
        _data[3] -= other[3];
        return *this;
    }
    friend GfVec4d operator-(GfVec4d const &l, GfVec4d const &r) {
        return GfVec4d(l) -= r;
    }

    GfVec4d &operator*=(double s) {
        _data[0] *= s;
        _data[1] *= s;
        _data[2] *= s;
        _data[3] *= 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 _data[0] * v[0] + _data[1] * v[1] + _data[2] * v[2] + _data[3] * v[3];
    }

    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;
    }

  
private:
    double _data[4];
};

GF_API std::ostream& operator<<(std::ostream &, GfVec4d const &);


PXR_NAMESPACE_CLOSE_SCOPE

#include "pxr/base/gf/vec4f.h"
#include "pxr/base/gf/vec4h.h"
#include "pxr/base/gf/vec4i.h"

PXR_NAMESPACE_OPEN_SCOPE

inline
GfVec4d::GfVec4d(class GfVec4f const &other)
{
    _data[0] = other[0];
    _data[1] = other[1];
    _data[2] = other[2];
    _data[3] = other[3];
}
inline
GfVec4d::GfVec4d(class GfVec4h const &other)
{
    _data[0] = other[0];
    _data[1] = other[1];
    _data[2] = other[2];
    _data[3] = other[3];
}
inline
GfVec4d::GfVec4d(class GfVec4i const &other)
{
    _data[0] = other[0];
    _data[1] = other[1];
    _data[2] = other[2];
    _data[3] = other[3];
}

inline GfVec4d
GfCompMult(GfVec4d const &v1, GfVec4d const &v2) {
    return GfVec4d(
        v1[0] * v2[0],
        v1[1] * v2[1],
        v1[2] * v2[2],
        v1[3] * v2[3]
        );
}

inline GfVec4d
GfCompDiv(GfVec4d const &v1, GfVec4d const &v2) {
    return GfVec4d(
        v1[0] / v2[0],
        v1[1] / v2[1],
        v1[2] / v2[2],
        v1[3] / v2[3]
        );
}

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