MRVector4.h

Go to the documentation of this file.

#pragma once
 
#include <cmath>
#include <iosfwd>
#include <type_traits>
#include "MRPch/MRBindingMacros.h"
#include "MRMesh/MRMacros.h"
#include "MRVector3.h"
 
namespace MR
{
 
 
#ifdef _MSC_VER
#pragma warning(push)
#pragma warning(disable: 4804) 
#pragma warning(disable: 4146) 
#endif

template <typename T>
struct Vector4
{
    using ValueType = T;
    using MatrixType = Matrix4<T>;
    using SymMatrixType = SymMatrix4<T>;
    static constexpr int elements = 4;
 
    T x, y, z, w;
 
    constexpr Vector4() noexcept : x( 0 ), y( 0 ), z( 0 ), w( 0 )
    {
        static_assert( sizeof( Vector4<ValueType> ) == elements * sizeof( ValueType ), "Struct size invalid" );
        static_assert( elements == 4, "Invalid number of elements" );
    }
    explicit Vector4( NoInit ) noexcept { }
    constexpr Vector4( T x, T y, T z, T w ) noexcept : x( x ), y( y ), z( z ), w( w ) { }
    static constexpr Vector4 diagonal( T a ) noexcept
    {
        return Vector4( a, a, a, a );
    }

    template <typename U> MR_REQUIRES_IF_SUPPORTED( !std::is_same_v<T, U> )
    constexpr explicit Vector4( const Vector4<U> & v ) noexcept : x( T( v.x ) ), y( T( v.y ) ), z( T( v.z ) ), w( T( v.w ) )
    {
    }
 
    constexpr const T & operator []( int e ) const noexcept { return *( ( ValueType *)this + e ); }
    constexpr       T & operator []( int e )       noexcept { return *( ( ValueType* )this + e ); }
 
    T lengthSq() const
    {
        return x * x + y * y + z * z + w * w;
    }
    auto length() const
    {
        using std::sqrt;
        return sqrt( lengthSq() );
    }
 
    Vector4 normalized() const MR_REQUIRES_IF_SUPPORTED( !std::is_integral_v<T> )
    {
        auto len = length();
        if ( len <= 0 )
            return {};
        return ( 1 / len ) * ( *this );
    }

    Vector3<T> proj3d() const MR_REQUIRES_IF_SUPPORTED( !std::is_integral_v<T> )
    {
        return { x / w, y / w, z / w };
    }
 
    [[nodiscard]] bool isFinite() const MR_REQUIRES_IF_SUPPORTED( std::is_floating_point_v<T> )
    {
        return std::isfinite( x ) && std::isfinite( y ) && std::isfinite( z ) && std::isfinite( w );
    }
 
    [[nodiscard]] friend constexpr bool operator ==( const Vector4<T> & a, const Vector4<T> & b ) { return a.x == b.x && a.y == b.y && a.z == b.z && a.w == b.w; }
    [[nodiscard]] friend constexpr bool operator !=( const Vector4<T> & a, const Vector4<T> & b ) { return !( a == b ); }

 
    [[nodiscard]] friend constexpr const Vector4<T> & operator +( const Vector4<T> & a ) { return a; }
    [[nodiscard]] friend constexpr auto operator -( const Vector4<T> & a ) -> Vector4<decltype( -std::declval<T>() )> { return { -a.x, -a.y, -a.z, -a.w }; }
 
    [[nodiscard]] friend constexpr auto operator +( const Vector4<T> & a, const Vector4<T> & b ) -> Vector4<decltype( std::declval<T>() + std::declval<T>() )> { return { a.x + b.x, a.y + b.y, a.z + b.z, a.w + b.w }; }
    [[nodiscard]] friend constexpr auto operator -( const Vector4<T> & a, const Vector4<T> & b ) -> Vector4<decltype( std::declval<T>() - std::declval<T>() )> { return { a.x - b.x, a.y - b.y, a.z - b.z, a.w - b.w }; }
    [[nodiscard]] friend constexpr auto operator *(               T    a, const Vector4<T> & b ) -> Vector4<decltype( std::declval<T>() * std::declval<T>() )> { return { a * b.x, a * b.y, a * b.z, a * b.w }; }
    [[nodiscard]] friend constexpr auto operator *( const Vector4<T> & b,               T    a ) -> Vector4<decltype( std::declval<T>() * std::declval<T>() )> { return { a * b.x, a * b.y, a * b.z, a * b.w }; }
    [[nodiscard]] friend constexpr auto operator /(       Vector4<T>   b,               T    a ) -> Vector4<decltype( std::declval<T>() / std::declval<T>() )>
    {
        if constexpr ( std::is_integral_v<T> )
            return { b.x / a, b.y / a, b.z / a, b.w / a };
        else
            return b * ( 1 / a );
    }
 
    friend constexpr Vector4<T> & operator +=( Vector4<T> & a, const Vector4<T> & b ) { a.x += b.x; a.y += b.y; a.z += b.z; a.w += b.w; return a; }
    friend constexpr Vector4<T> & operator -=( Vector4<T> & a, const Vector4<T> & b ) { a.x -= b.x; a.y -= b.y; a.z -= b.z; a.w -= b.w; return a; }
    friend constexpr Vector4<T> & operator *=( Vector4<T> & a,               T    b ) { a.x *= b; a.y *= b; a.z *= b; a.w *= b; return a; }
    friend constexpr Vector4<T> & operator /=( Vector4<T> & a,               T    b )
    {
        if constexpr ( std::is_integral_v<T> )
            { a.x /= b; a.y /= b; a.z /= b; a.w /= b; return a; }
        else
            return a *= ( 1 / b );
    }
 
    friend std::ostream& operator<<( std::ostream& s, const Vector4& vec )
    {
        return s << vec.x << ' ' << vec.y << ' ' << vec.z << ' ' << vec.w;
    }
 
    friend std::istream& operator>>( std::istream& s, Vector4& vec )
    {
        return s >> vec.x >> vec.y >> vec.z >> vec.w;
    }
 

 
    MR_BIND_IGNORE friend auto begin( const Vector4& v ) { return &v[0]; }
    MR_BIND_IGNORE friend auto begin( Vector4& v ) { return &v[0]; }
    MR_BIND_IGNORE friend auto end( const Vector4& v ) { return &v[4]; }
    MR_BIND_IGNORE friend auto end( Vector4& v ) { return &v[4]; }
};


template <typename T>
inline T distanceSq( const Vector4<T> & a, const Vector4<T> & b )
{
    return ( a - b ).lengthSq();
}

template <typename T>
inline T distance( const Vector4<T> & a, const Vector4<T> & b )
{
    return ( a - b ).length();
}

template <typename T>
inline auto dot( const Vector4<T> & a, const Vector4<T> & b ) -> decltype( a.x * b.x )
{
    return a.x * b.x + a.y * b.y + a.z * b.z + a.w * b.w;
}

template <typename T>
inline T sqr( const Vector4<T> & a )
{
    return a.lengthSq();
}

template <typename T>
inline Vector4<T> mult( const Vector4<T>& a, const Vector4<T>& b )
{
    return { a.x * b.x, a.y * b.y, a.z * b.z, a.w * b.w };
}

template <typename T>
inline Vector4<T> div( const Vector4<T>& a, const Vector4<T>& b )
{
    return { a.x / b.x, a.y / b.y, a.z / b.z, a.w / b.w };
}

 
#ifdef _MSC_VER
#pragma warning(pop)
#endif
 
}