Vector4.hh

Go to the documentation of this file.

/*
 * Copyright (C) 2012 Open Source Robotics Foundation
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 *
*/
#ifndef GZ_MATH_VECTOR4_HH_
#define GZ_MATH_VECTOR4_HH_
 
#include <algorithm>
#include <cmath>
#include <limits>
 
#include <gz/math/Matrix4.hh>
#include <gz/math/Helpers.hh>
#include <gz/math/config.hh>
 
namespace gz::math
{
  // Inline bracket to help doxygen filtering.
  inline namespace GZ_MATH_VERSION_NAMESPACE {
  //
  template<typename T>
  class Vector4
  {
    public: static const Vector4<T> &Zero;
 
    public: static const Vector4<T> &One;
 
    public: static const Vector4 &NaN;
 
    public: constexpr Vector4()
    : data{0, 0, 0, 0}
    {
    }
 
    public: constexpr Vector4(const T &_x, const T &_y, const T &_z,
                              const T &_w)
    : data{_x, _y, _z, _w}
    {
    }
 
    public: T Distance(const Vector4<T> &_pt) const
    {
      return static_cast<T>(sqrt(
                  (this->data[0]-_pt[0])*(this->data[0]-_pt[0]) +
                  (this->data[1]-_pt[1])*(this->data[1]-_pt[1]) +
                  (this->data[2]-_pt[2])*(this->data[2]-_pt[2]) +
                  (this->data[3]-_pt[3])*(this->data[3]-_pt[3])));
    }
 
    public: T Distance(T _x, T _y, T _z, T _w) const
    {
      return this->Distance(Vector4(_x, _y, _z, _w));
    }
 
    public: T Length() const
    {
      return static_cast<T>(sqrt(this->SquaredLength()));
    }
 
    public: T SquaredLength() const
    {
      return
        this->data[0] * this->data[0] +
        this->data[1] * this->data[1] +
        this->data[2] * this->data[2] +
        this->data[3] * this->data[3];
    }
 
    public: void Round()
    {
      this->data[0] = static_cast<T>(std::nearbyint(this->data[0]));
      this->data[1] = static_cast<T>(std::nearbyint(this->data[1]));
      this->data[2] = static_cast<T>(std::nearbyint(this->data[2]));
      this->data[3] = static_cast<T>(std::nearbyint(this->data[3]));
    }
 
    public: Vector4 Rounded() const
    {
      Vector4<T> result = *this;
      result.Round();
      return result;
    }
 
    public: inline void Correct()
    {
      // std::isfinite works with floating point values,
      // need to explicit cast to avoid ambiguity in vc++.
      if (!std::isfinite(static_cast<double>(this->data[0])))
        this->data[0] = 0;
      if (!std::isfinite(static_cast<double>(this->data[1])))
        this->data[1] = 0;
      if (!std::isfinite(static_cast<double>(this->data[2])))
        this->data[2] = 0;
      if (!std::isfinite(static_cast<double>(this->data[3])))
        this->data[3] = 0;
    }
 
    public: void Normalize()
    {
      T d = this->Length();
 
      if (!equal<T>(d, static_cast<T>(0.0)))
      {
        this->data[0] /= d;
        this->data[1] /= d;
        this->data[2] /= d;
        this->data[3] /= d;
      }
    }
 
    public: Vector4 Normalized() const
    {
      Vector4<T> result = *this;
      result.Normalize();
      return result;
    }
 
    public: T Dot(const Vector4<T> &_v) const
    {
      return this->data[0] * _v[0] +
             this->data[1] * _v[1] +
             this->data[2] * _v[2] +
             this->data[3] * _v[3];
    }
 
    public: T AbsDot(const Vector4<T> &_v) const
    {
      return std::abs(this->data[0] * _v[0]) +
             std::abs(this->data[1] * _v[1]) +
             std::abs(this->data[2] * _v[2]) +
             std::abs(this->data[3] * _v[3]);
    }
 
    public: Vector4 Abs() const
    {
      return Vector4(std::abs(this->data[0]),
                     std::abs(this->data[1]),
                     std::abs(this->data[2]),
                     std::abs(this->data[3]));
    }
 
    public: void Set(T _x = 0, T _y = 0, T _z = 0, T _w = 0)
    {
      this->data[0] = _x;
      this->data[1] = _y;
      this->data[2] = _z;
      this->data[3] = _w;
    }
 
    public: void Max(const Vector4<T> &_v)
    {
      this->data[0] = std::max(_v[0], this->data[0]);
      this->data[1] = std::max(_v[1], this->data[1]);
      this->data[2] = std::max(_v[2], this->data[2]);
      this->data[3] = std::max(_v[3], this->data[3]);
    }
 
    public: void Min(const Vector4<T> &_v)
    {
      this->data[0] = std::min(_v[0], this->data[0]);
      this->data[1] = std::min(_v[1], this->data[1]);
      this->data[2] = std::min(_v[2], this->data[2]);
      this->data[3] = std::min(_v[3], this->data[3]);
    }
 
    public: T Max() const
    {
      return *std::max_element(this->data, this->data+4);
    }
 
    public: T Min() const
    {
      return *std::min_element(this->data, this->data+4);
    }
 
    public: T Sum() const
    {
      return this->data[0] + this->data[1] + this->data[2] + this->data[3];
    }
 
    public: Vector4<T> &operator=(T _value)
    {
      this->data[0] = _value;
      this->data[1] = _value;
      this->data[2] = _value;
      this->data[3] = _value;
 
      return *this;
    }
 
    public: Vector4<T> operator+(const Vector4<T> &_v) const
    {
      return Vector4<T>(this->data[0] + _v[0],
                        this->data[1] + _v[1],
                        this->data[2] + _v[2],
                        this->data[3] + _v[3]);
    }
 
    public: const Vector4<T> &operator+=(const Vector4<T> &_v)
    {
      this->data[0] += _v[0];
      this->data[1] += _v[1];
      this->data[2] += _v[2];
      this->data[3] += _v[3];
 
      return *this;
    }
 
    public: inline Vector4<T> operator+(const T _s) const
    {
      return Vector4<T>(this->data[0] + _s,
                        this->data[1] + _s,
                        this->data[2] + _s,
                        this->data[3] + _s);
    }
 
    public: friend inline Vector4<T> operator+(const T _s,
                                               const Vector4<T> &_v)
    {
      return _v + _s;
    }
 
    public: const Vector4<T> &operator+=(const T _s)
    {
      this->data[0] += _s;
      this->data[1] += _s;
      this->data[2] += _s;
      this->data[3] += _s;
 
      return *this;
    }
 
    public: inline Vector4 operator-() const
    {
      return Vector4(-this->data[0], -this->data[1],
                     -this->data[2], -this->data[3]);
    }
 
    public: Vector4<T> operator-(const Vector4<T> &_v) const
    {
      return Vector4<T>(this->data[0] - _v[0],
                        this->data[1] - _v[1],
                        this->data[2] - _v[2],
                        this->data[3] - _v[3]);
    }
 
    public: const Vector4<T> &operator-=(const Vector4<T> &_v)
    {
      this->data[0] -= _v[0];
      this->data[1] -= _v[1];
      this->data[2] -= _v[2];
      this->data[3] -= _v[3];
 
      return *this;
    }
 
    public: inline Vector4<T> operator-(const T _s) const
    {
      return Vector4<T>(this->data[0] - _s,
                        this->data[1] - _s,
                        this->data[2] - _s,
                        this->data[3] - _s);
    }
 
    public: friend inline Vector4<T> operator-(const T _s,
                                               const Vector4<T> &_v)
    {
      return {_s - _v.X(), _s - _v.Y(), _s - _v.Z(), _s - _v.W()};
    }
 
    public: const Vector4<T> &operator-=(const T _s)
    {
      this->data[0] -= _s;
      this->data[1] -= _s;
      this->data[2] -= _s;
      this->data[3] -= _s;
 
      return *this;
    }
 
    public: const Vector4<T> operator/(const Vector4<T> &_v) const
    {
      return Vector4<T>(this->data[0] / _v[0],
                        this->data[1] / _v[1],
                        this->data[2] / _v[2],
                        this->data[3] / _v[3]);
    }
 
    public: const Vector4<T> &operator/=(const Vector4<T> &_v)
    {
      this->data[0] /= _v[0];
      this->data[1] /= _v[1];
      this->data[2] /= _v[2];
      this->data[3] /= _v[3];
 
      return *this;
    }
 
    public: const Vector4<T> operator/(T _v) const
    {
      return Vector4<T>(this->data[0] / _v, this->data[1] / _v,
          this->data[2] / _v, this->data[3] / _v);
    }
 
    public: const Vector4<T> &operator/=(T _v)
    {
      this->data[0] /= _v;
      this->data[1] /= _v;
      this->data[2] /= _v;
      this->data[3] /= _v;
 
      return *this;
    }
 
    public: const Vector4<T> operator*(const Vector4<T> &_pt) const
    {
      return Vector4<T>(this->data[0] * _pt[0],
                        this->data[1] * _pt[1],
                        this->data[2] * _pt[2],
                        this->data[3] * _pt[3]);
    }
 
    public: const Vector4<T> operator*(const Matrix4<T> &_m) const
    {
      return Vector4<T>(
          this->data[0]*_m(0, 0) + this->data[1]*_m(1, 0) +
          this->data[2]*_m(2, 0) + this->data[3]*_m(3, 0),
          this->data[0]*_m(0, 1) + this->data[1]*_m(1, 1) +
          this->data[2]*_m(2, 1) + this->data[3]*_m(3, 1),
          this->data[0]*_m(0, 2) + this->data[1]*_m(1, 2) +
          this->data[2]*_m(2, 2) + this->data[3]*_m(3, 2),
          this->data[0]*_m(0, 3) + this->data[1]*_m(1, 3) +
          this->data[2]*_m(2, 3) + this->data[3]*_m(3, 3));
    }
 
    public: const Vector4<T> &operator*=(const Vector4<T> &_pt)
    {
      this->data[0] *= _pt[0];
      this->data[1] *= _pt[1];
      this->data[2] *= _pt[2];
      this->data[3] *= _pt[3];
 
      return *this;
    }
 
    public: const Vector4<T> operator*(T _v) const
    {
      return Vector4<T>(this->data[0] * _v, this->data[1] * _v,
          this->data[2] * _v, this->data[3] * _v);
    }
 
    public: friend inline const Vector4 operator*(const T _s,
                                                  const Vector4 &_v)
    {
      return Vector4(_v * _s);
    }
 
    public: const Vector4<T> &operator*=(T _v)
    {
      this->data[0] *= _v;
      this->data[1] *= _v;
      this->data[2] *= _v;
      this->data[3] *= _v;
 
      return *this;
    }
 
    public: bool Equal(const Vector4 &_v, const T &_tol) const
    {
      return equal<T>(this->data[0], _v[0], _tol)
          && equal<T>(this->data[1], _v[1], _tol)
          && equal<T>(this->data[2], _v[2], _tol)
          && equal<T>(this->data[3], _v[3], _tol);
    }
 
    public: bool operator==(const Vector4<T> &_v) const
    {
      return this->Equal(_v, static_cast<T>(1e-6));
    }
 
    public: bool operator!=(const Vector4<T> &_pt) const
    {
      return !(*this == _pt);
    }
 
    public: bool IsFinite() const
    {
      // std::isfinite works with floating point values,
      // need to explicit cast to avoid ambiguity in vc++.
      return std::isfinite(static_cast<double>(this->data[0])) &&
             std::isfinite(static_cast<double>(this->data[1])) &&
             std::isfinite(static_cast<double>(this->data[2])) &&
             std::isfinite(static_cast<double>(this->data[3]));
    }
 
    public: T &operator[](const std::size_t _index)
    {
      return this->data[clamp(_index, GZ_ZERO_SIZE_T, GZ_THREE_SIZE_T)];
    }
 
    public: T operator[](const std::size_t _index) const
    {
      return this->data[clamp(_index, GZ_ZERO_SIZE_T, GZ_THREE_SIZE_T)];
    }
 
    public: T &X()
    {
      return this->data[0];
    }
 
    public: T &Y()
    {
      return this->data[1];
    }
 
    public: T &Z()
    {
      return this->data[2];
    }
 
    public: T &W()
    {
      return this->data[3];
    }
 
    public: T X() const
    {
      return this->data[0];
    }
 
    public: T Y() const
    {
      return this->data[1];
    }
 
    public: T Z() const
    {
      return this->data[2];
    }
 
    public: T W() const
    {
      return this->data[3];
    }
 
    public: inline void X(const T &_v)
    {
      this->data[0] = _v;
    }
 
    public: inline void Y(const T &_v)
    {
      this->data[1] = _v;
    }
 
    public: inline void Z(const T &_v)
    {
      this->data[2] = _v;
    }
 
    public: inline void W(const T &_v)
    {
      this->data[3] = _v;
    }
 
    public: bool operator<(const Vector4<T> &_pt) const
    {
      return this->data[0] < _pt[0] || this->data[1] < _pt[1] ||
             this->data[2] < _pt[2] || this->data[3] < _pt[3];
    }
 
    public: friend std::ostream &operator<<(
                std::ostream &_out, const gz::math::Vector4<T> &_pt)
    {
      for (auto i : {0, 1, 2, 3})
      {
        if (i > 0)
          _out << " ";
 
        appendToStream(_out, _pt[i]);
      }
      return _out;
    }
 
    public: friend std::istream &operator>>(
                std::istream &_in, gz::math::Vector4<T> &_pt)
    {
      T x, y, z, w;
 
      // Skip white spaces
      _in.setf(std::ios_base::skipws);
      _in >> x >> y >> z >> w;
      if (!_in.fail())
      {
        _pt.Set(x, y, z, w);
      }
      return _in;
    }
 
    private: T data[4];
  };
 
  namespace detail {
 
    template<typename T>
    constexpr Vector4<T> gVector4Zero(0, 0, 0, 0);
 
    template<typename T>
    constexpr Vector4<T> gVector4One(1, 1, 1, 1);
 
    template<typename T>
    constexpr Vector4<T> gVector4NaN(
        std::numeric_limits<T>::quiet_NaN(),
        std::numeric_limits<T>::quiet_NaN(),
        std::numeric_limits<T>::quiet_NaN(),
        std::numeric_limits<T>::quiet_NaN());
 
  }  // namespace detail
 
  template<typename T>
  const Vector4<T> &Vector4<T>::Zero = detail::gVector4Zero<T>;
 
  template<typename T>
  const Vector4<T> &Vector4<T>::One = detail::gVector4One<T>;
 
  template<typename T>
  const Vector4<T> &Vector4<T>::NaN = detail::gVector4NaN<T>;
 
  typedef Vector4<int> Vector4i;
  typedef Vector4<double> Vector4d;
  typedef Vector4<float> Vector4f;
  }  // namespace GZ_MATH_VERSION_NAMESPACE
}  // namespace gz::math
#endif  // GZ_MATH_VERSION_NAMESPACE