text: /// @ref core
  /// @file glm/detail/func_integer.inl
  
  #include "type_vec2.hpp"
  #include "type_vec3.hpp"
  #include "type_vec4.hpp"
  #include "type_int.hpp"
  #include "_vectorize.hpp"
  #if(GLM_ARCH & GLM_ARCH_X86 && GLM_COMPILER & GLM_COMPILER_VC)
  #	include <intrin.h>
  #	pragma intrinsic(_BitScanReverse)
  #endif//(GLM_ARCH & GLM_ARCH_X86 && GLM_COMPILER & GLM_COMPILER_VC)
  #include <limits>
  
  #if !GLM_HAS_EXTENDED_INTEGER_TYPE
  #	if GLM_COMPILER & GLM_COMPILER_GCC
  #		pragma GCC diagnostic ignored "-Wlong-long"
  #	endif
  #	if (GLM_COMPILER & GLM_COMPILER_CLANG)
  #		pragma clang diagnostic ignored "-Wc++11-long-long"
  #	endif
  #endif
  
  namespace glm{
  namespace detail
  {
  	template<typename T>
  	GLM_FUNC_QUALIFIER T mask(T Bits)
  	{
  		return Bits >= sizeof(T) * 8 ? ~static_cast<T>(0) : (static_cast<T>(1) << Bits) - static_cast<T>(1);
  	}
  
  	template<length_t L, typename T, glm::precision P, template<length_t, typename, precision> class vecType, bool Aligned, bool EXEC>
  	struct compute_bitfieldReverseStep
  	{
  		GLM_FUNC_QUALIFIER static vecType<L, T, P> call(vecType<L, T, P> const& v, T, T)
  		{
  			return v;
  		}
  	};
  
  	template<length_t L, typename T, glm::precision P, template<length_t, typename, precision> class vecType, bool Aligned>
  	struct compute_bitfieldReverseStep<L, T, P, vecType, Aligned, true>
  	{
  		GLM_FUNC_QUALIFIER static vecType<L, T, P> call(vecType<L, T, P> const& v, T Mask, T Shift)
  		{
  			return (v & Mask) << Shift | (v & (~Mask)) >> Shift;
  		}
  	};
  
  	template<length_t L, typename T, glm::precision P, template<length_t, typename, precision> class vecType, bool Aligned, bool EXEC>
  	struct compute_bitfieldBitCountStep
  	{
  		GLM_FUNC_QUALIFIER static vecType<L, T, P> call(vecType<L, T, P> const& v, T, T)
  		{
  			return v;
  		}
  	};
  
  	template<length_t L, typename T, glm::precision P, template<length_t, typename, precision> class vecType, bool Aligned>
  	struct compute_bitfieldBitCountStep<L, T, P, vecType, Aligned, true>
  	{
  		GLM_FUNC_QUALIFIER static vecType<L, T, P> call(vecType<L, T, P> const& v, T Mask, T Shift)
  		{
  			return (v & Mask) + ((v >> Shift) & Mask);
  		}
  	};
  
  	template<typename genIUType, size_t Bits>
  	struct compute_findLSB
  	{
  		GLM_FUNC_QUALIFIER static int call(genIUType Value)
  		{
  			if(Value == 0)
  				return -1;
  
  			return glm::bitCount(~Value & (Value - static_cast<genIUType>(1)));
  		}
  	};
  
  #	if GLM_HAS_BITSCAN_WINDOWS
  		template<typename genIUType>
  		struct compute_findLSB<genIUType, 32>
  		{
  			GLM_FUNC_QUALIFIER static int call(genIUType Value)
  			{
  				unsigned long Result(0);
  				unsigned char IsNotNull = _BitScanForward(&Result, *reinterpret_cast<unsigned long*>(&Value));
  				return IsNotNull ? int(Result) : -1;
  			}
  		};
  
  #		if !((GLM_COMPILER & GLM_COMPILER_VC) && (GLM_MODEL == GLM_MODEL_32))
  		template<typename genIUType>
  		struct compute_findLSB<genIUType, 64>
  		{
  			GLM_FUNC_QUALIFIER static int call(genIUType Value)
  			{
  				unsigned long Result(0);
  				unsigned char IsNotNull = _BitScanForward64(&Result, *reinterpret_cast<unsigned __int64*>(&Value));
  				return IsNotNull ? int(Result) : -1;
  			}
  		};
  #		endif
  #	endif//GLM_HAS_BITSCAN_WINDOWS
  
  	template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType, bool EXEC = true>
  	struct compute_findMSB_step_vec
  	{
  		GLM_FUNC_QUALIFIER static vecType<L, T, P> call(vecType<L, T, P> const & x, T Shift)
  		{
  			return x | (x >> Shift);
  		}
  	};
  
  	template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType>
  	struct compute_findMSB_step_vec<L, T, P, vecType, false>
  	{
  		GLM_FUNC_QUALIFIER static vecType<L, T, P> call(vecType<L, T, P> const& x, T)
  		{
  			return x;
  		}
  	};
  
  	template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType, int>
  	struct compute_findMSB_vec
  	{
  		GLM_FUNC_QUALIFIER static vecType<L, int, P> call(vecType<L, T, P> const& v)
  		{
  			vecType<L, T, P> x(v);
  			x = compute_findMSB_step_vec<L, T, P, vecType, sizeof(T) * 8 >=  8>::call(x, static_cast<T>( 1));
  			x = compute_findMSB_step_vec<L, T, P, vecType, sizeof(T) * 8 >=  8>::call(x, static_cast<T>( 2));
  			x = compute_findMSB_step_vec<L, T, P, vecType, sizeof(T) * 8 >=  8>::call(x, static_cast<T>( 4));
  			x = compute_findMSB_step_vec<L, T, P, vecType, sizeof(T) * 8 >= 16>::call(x, static_cast<T>( 8));
  			x = compute_findMSB_step_vec<L, T, P, vecType, sizeof(T) * 8 >= 32>::call(x, static_cast<T>(16));
  			x = compute_findMSB_step_vec<L, T, P, vecType, sizeof(T) * 8 >= 64>::call(x, static_cast<T>(32));
  			return vecType<L, int, P>(sizeof(T) * 8 - 1) - glm::bitCount(~x);
  		}
  	};
  
  #	if GLM_HAS_BITSCAN_WINDOWS
  		template<typename genIUType>
  		GLM_FUNC_QUALIFIER int compute_findMSB_32(genIUType Value)
  		{
  			unsigned long Result(0);
  			unsigned char IsNotNull = _BitScanReverse(&Result, *reinterpret_cast<unsigned long*>(&Value));
  			return IsNotNull ? int(Result) : -1;
  		}
  
  		template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType>
  		struct compute_findMSB_vec<L, T, P, vecType, 32>
  		{
  			GLM_FUNC_QUALIFIER static vecType<L, int, P> call(vecType<L, T, P> const& x)
  			{
  				return detail::functor1<L, int, T, P>::call(compute_findMSB_32, x);
  			}
  		};
  
  #		if !((GLM_COMPILER & GLM_COMPILER_VC) && (GLM_MODEL == GLM_MODEL_32))
  		template<typename genIUType>
  		GLM_FUNC_QUALIFIER int compute_findMSB_64(genIUType Value)
  		{
  			unsigned long Result(0);
  			unsigned char IsNotNull = _BitScanReverse64(&Result, *reinterpret_cast<unsigned __int64*>(&Value));
  			return IsNotNull ? int(Result) : -1;
  		}
  
  		template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType>
  		struct compute_findMSB_vec<L, T, P, vecType, 64>
  		{
  			GLM_FUNC_QUALIFIER static vecType<L, int, P> call(vecType<L, T, P> const& x)
  			{
  				return detail::functor1<L, int, T, P>::call(compute_findMSB_64, x);
  			}
  		};
  #		endif
  #	endif//GLM_HAS_BITSCAN_WINDOWS
  }//namespace detail
  
  	// uaddCarry
  	GLM_FUNC_QUALIFIER uint uaddCarry(uint const & x, uint const & y, uint & Carry)
  	{
  		uint64 const Value64(static_cast<uint64>(x) + static_cast<uint64>(y));
  		uint64 const Max32((static_cast<uint64>(1) << static_cast<uint64>(32)) - static_cast<uint64>(1));
  		Carry = Value64 > Max32 ? 1u : 0u;
  		return static_cast<uint32>(Value64 % (Max32 + static_cast<uint64>(1)));
  	}
  
  	template<length_t L, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, uint, P> uaddCarry(vecType<L, uint, P> const& x, vecType<L, uint, P> const& y, vecType<L, uint, P>& Carry)
  	{
  		vecType<L, uint64, P> Value64(vecType<L, uint64, P>(x) + vecType<L, uint64, P>(y));
  		vecType<L, uint64, P> Max32((static_cast<uint64>(1) << static_cast<uint64>(32)) - static_cast<uint64>(1));
  		Carry = mix(vecType<L, uint32, P>(0), vecType<L, uint32, P>(1), greaterThan(Value64, Max32));
  		return vecType<L, uint32,P>(Value64 % (Max32 + static_cast<uint64>(1)));
  	}
  
  	// usubBorrow
  	GLM_FUNC_QUALIFIER uint usubBorrow(uint const & x, uint const & y, uint & Borrow)
  	{
  		GLM_STATIC_ASSERT(sizeof(uint) == sizeof(uint32), "uint and uint32 size mismatch");
  
  		Borrow = x >= y ? static_cast<uint32>(0) : static_cast<uint32>(1);
  		if(y >= x)
  			return y - x;
  		else
  			return static_cast<uint32>((static_cast<int64>(1) << static_cast<int64>(32)) + (static_cast<int64>(y) - static_cast<int64>(x)));
  	}
  
  	template<length_t L, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, uint, P> usubBorrow(vecType<L, uint, P> const& x, vecType<L, uint, P> const& y, vecType<L, uint, P>& Borrow)
  	{
  		Borrow = mix(vecType<L, uint, P>(1), vecType<L, uint, P>(0), greaterThanEqual(x, y));
  		vecType<L, uint, P> const YgeX(y - x);
  		vecType<L, uint, P> const XgeY(vecType<L, uint32, P>((static_cast<int64>(1) << static_cast<int64>(32)) + (vecType<L, int64, P>(y) - vecType<L, int64, P>(x))));
  		return mix(XgeY, YgeX, greaterThanEqual(y, x));
  	}
  
  	// umulExtended
  	GLM_FUNC_QUALIFIER void umulExtended(uint const & x, uint const & y, uint & msb, uint & lsb)
  	{
  		GLM_STATIC_ASSERT(sizeof(uint) == sizeof(uint32), "uint and uint32 size mismatch");
  
  		uint64 Value64 = static_cast<uint64>(x) * static_cast<uint64>(y);
  		msb = static_cast<uint>(Value64 >> static_cast<uint64>(32));
  		lsb = static_cast<uint>(Value64);
  	}
  
  	template<length_t L, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER void umulExtended(vecType<L, uint, P> const& x, vecType<L, uint, P> const& y, vecType<L, uint, P>& msb, vecType<L, uint, P>& lsb)
  	{
  		GLM_STATIC_ASSERT(sizeof(uint) == sizeof(uint32), "uint and uint32 size mismatch");
  
  		vecType<L, uint64, P> Value64(vecType<L, uint64, P>(x) * vecType<L, uint64, P>(y));
  		msb = vecType<L, uint32, P>(Value64 >> static_cast<uint64>(32));
  		lsb = vecType<L, uint32, P>(Value64);
  	}
  
  	// imulExtended
  	GLM_FUNC_QUALIFIER void imulExtended(int x, int y, int& msb, int& lsb)
  	{
  		GLM_STATIC_ASSERT(sizeof(int) == sizeof(int32), "int and int32 size mismatch");
  
  		int64 Value64 = static_cast<int64>(x) * static_cast<int64>(y);
  		msb = static_cast<int>(Value64 >> static_cast<int64>(32));
  		lsb = static_cast<int>(Value64);
  	}
  
  	template<length_t L, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER void imulExtended(vecType<L, int, P> const& x, vecType<L, int, P> const& y, vecType<L, int, P>& msb, vecType<L, int, P>& lsb)
  	{
  		GLM_STATIC_ASSERT(sizeof(int) == sizeof(int32), "int and int32 size mismatch");
  
  		vecType<L, int64, P> Value64(vecType<L, int64, P>(x) * vecType<L, int64, P>(y));
  		lsb = vecType<L, int32, P>(Value64 & static_cast<int64>(0xFFFFFFFF));
  		msb = vecType<L, int32, P>((Value64 >> static_cast<int64>(32)) & static_cast<int64>(0xFFFFFFFF));
  	}
  
  	// bitfieldExtract
  	template<typename genIUType>
  	GLM_FUNC_QUALIFIER genIUType bitfieldExtract(genIUType Value, int Offset, int Bits)
  	{
  		return bitfieldExtract(vec<1, genIUType>(Value), Offset, Bits).x;
  	}
  
  	template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, T, P> bitfieldExtract(vecType<L, T, P> const& Value, int Offset, int Bits)
  	{
  		GLM_STATIC_ASSERT(std::numeric_limits<T>::is_integer, "'bitfieldExtract' only accept integer inputs");
  
  		return (Value >> static_cast<T>(Offset)) & static_cast<T>(detail::mask(Bits));
  	}
  
  	// bitfieldInsert
  	template<typename genIUType>
  	GLM_FUNC_QUALIFIER genIUType bitfieldInsert(genIUType const & Base, genIUType const & Insert, int Offset, int Bits)
  	{
  		return bitfieldInsert(vec<1, genIUType>(Base), vec<1, genIUType>(Insert), Offset, Bits).x;
  	}
  
  	template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, T, P> bitfieldInsert(vecType<L, T, P> const& Base, vecType<L, T, P> const& Insert, int Offset, int Bits)
  	{
  		GLM_STATIC_ASSERT(std::numeric_limits<T>::is_integer, "'bitfieldInsert' only accept integer values");
  
  		T const Mask = static_cast<T>(detail::mask(Bits) << Offset);
  		return (Base & ~Mask) | (Insert & Mask);
  	}
  
  	// bitfieldReverse
  	template<typename genType>
  	GLM_FUNC_QUALIFIER genType bitfieldReverse(genType x)
  	{
  		return bitfieldReverse(glm::vec<1, genType, glm::defaultp>(x)).x;
  	}
  
  	template<length_t L, typename T, glm::precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, T, P> bitfieldReverse(vecType<L, T, P> const& v)
  	{
  		vecType<L, T, P> x(v);
  		x = detail::compute_bitfieldReverseStep<L, T, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>=  2>::call(x, T(0x5555555555555555ull), static_cast<T>( 1));
  		x = detail::compute_bitfieldReverseStep<L, T, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>=  4>::call(x, T(0x3333333333333333ull), static_cast<T>( 2));
  		x = detail::compute_bitfieldReverseStep<L, T, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>=  8>::call(x, T(0x0F0F0F0F0F0F0F0Full), static_cast<T>( 4));
  		x = detail::compute_bitfieldReverseStep<L, T, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>= 16>::call(x, T(0x00FF00FF00FF00FFull), static_cast<T>( 8));
  		x = detail::compute_bitfieldReverseStep<L, T, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>= 32>::call(x, T(0x0000FFFF0000FFFFull), static_cast<T>(16));
  		x = detail::compute_bitfieldReverseStep<L, T, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>= 64>::call(x, T(0x00000000FFFFFFFFull), static_cast<T>(32));
  		return x;
  	}
  
  	// bitCount
  	template<typename genType>
  	GLM_FUNC_QUALIFIER int bitCount(genType x)
  	{
  		return bitCount(glm::vec<1, genType, glm::defaultp>(x)).x;
  	}
  
  	template<length_t L, typename T, glm::precision P, template<length_t, typename, glm::precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, int, P> bitCount(vecType<L, T, P> const& v)
  	{
  	#if GLM_COMPILER & GLM_COMPILER_VC
  	#pragma warning(push)
  	#pragma warning(disable : 4310) //cast truncates constant value
  	#endif
  		vecType<L, typename detail::make_unsigned<T>::type, P> x(*reinterpret_cast<vecType<L, typename detail::make_unsigned<T>::type, P> const *>(&v));
  		x = detail::compute_bitfieldBitCountStep<L, typename detail::make_unsigned<T>::type, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>=  2>::call(x, typename detail::make_unsigned<T>::type(0x5555555555555555ull), typename detail::make_unsigned<T>::type( 1));
  		x = detail::compute_bitfieldBitCountStep<L, typename detail::make_unsigned<T>::type, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>=  4>::call(x, typename detail::make_unsigned<T>::type(0x3333333333333333ull), typename detail::make_unsigned<T>::type( 2));
  		x = detail::compute_bitfieldBitCountStep<L, typename detail::make_unsigned<T>::type, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>=  8>::call(x, typename detail::make_unsigned<T>::type(0x0F0F0F0F0F0F0F0Full), typename detail::make_unsigned<T>::type( 4));
  		x = detail::compute_bitfieldBitCountStep<L, typename detail::make_unsigned<T>::type, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>= 16>::call(x, typename detail::make_unsigned<T>::type(0x00FF00FF00FF00FFull), typename detail::make_unsigned<T>::type( 8));
  		x = detail::compute_bitfieldBitCountStep<L, typename detail::make_unsigned<T>::type, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>= 32>::call(x, typename detail::make_unsigned<T>::type(0x0000FFFF0000FFFFull), typename detail::make_unsigned<T>::type(16));
  		x = detail::compute_bitfieldBitCountStep<L, typename detail::make_unsigned<T>::type, P, vecType, detail::is_aligned<P>::value, sizeof(T) * 8>= 64>::call(x, typename detail::make_unsigned<T>::type(0x00000000FFFFFFFFull), typename detail::make_unsigned<T>::type(32));
  		return vecType<L, int, P>(x);
  	#if GLM_COMPILER & GLM_COMPILER_VC
  	#pragma warning(pop)
  	#endif
  	}
  
  	// findLSB
  	template<typename genIUType>
  	GLM_FUNC_QUALIFIER int findLSB(genIUType Value)
  	{
  		GLM_STATIC_ASSERT(std::numeric_limits<genIUType>::is_integer, "'findLSB' only accept integer values");
  
  		return detail::compute_findLSB<genIUType, sizeof(genIUType) * 8>::call(Value);
  	}
  
  	template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, int, P> findLSB(vecType<L, T, P> const& x)
  	{
  		GLM_STATIC_ASSERT(std::numeric_limits<T>::is_integer, "'findLSB' only accept integer values");
  
  		return detail::functor1<L, int, T, P>::call(findLSB, x);
  	}
  
  	// findMSB
  	template<typename genIUType>
  	GLM_FUNC_QUALIFIER int findMSB(genIUType v)
  	{
  		GLM_STATIC_ASSERT(std::numeric_limits<genIUType>::is_integer, "'findMSB' only accept integer values");
  
  		return findMSB(vec<1, genIUType>(v)).x;
  	}
  
  	template<length_t L, typename T, precision P, template<length_t, typename, precision> class vecType>
  	GLM_FUNC_QUALIFIER vecType<L, int, P> findMSB(vecType<L, T, P> const& v)
  	{
  		GLM_STATIC_ASSERT(std::numeric_limits<T>::is_integer, "'findMSB' only accept integer values");
  
  		return detail::compute_findMSB_vec<L, T, P, vecType, sizeof(T) * 8>::call(v);
  	}
  }//namespace glm
  
  #if GLM_ARCH != GLM_ARCH_PURE && GLM_HAS_UNRESTRICTED_UNIONS
  #	include "func_integer_simd.inl"
  #endif
  
  