pixo 0.1.0

data/ext/pixo/ext/glm/glm/gtc/quaternion_simd.inl

@@ -0,0 +1,198 @@
+/// @ref core
+/// @file glm/gtc/quaternion_simd.inl
+
+#if GLM_ARCH & GLM_ARCH_SSE2_BIT
+
+namespace glm{
+namespace detail
+{
+/*
+	template <precision P>
+	struct compute_quat_mul<float, P, true>
+	{
+		static tquat<float, P> call(tquat<float, P> const& q1, tquat<float, P> const& q2)
+		{
+			// SSE2 STATS: 11 shuffle, 8 mul, 8 add
+			// SSE4 STATS: 3 shuffle, 4 mul, 4 dpps
+
+			__m128 const mul0 = _mm_mul_ps(q1.Data, _mm_shuffle_ps(q2.Data, q2.Data, _MM_SHUFFLE(0, 1, 2, 3)));
+			__m128 const mul1 = _mm_mul_ps(q1.Data, _mm_shuffle_ps(q2.Data, q2.Data, _MM_SHUFFLE(1, 0, 3, 2)));
+			__m128 const mul2 = _mm_mul_ps(q1.Data, _mm_shuffle_ps(q2.Data, q2.Data, _MM_SHUFFLE(2, 3, 0, 1)));
+			__m128 const mul3 = _mm_mul_ps(q1.Data, q2.Data);
+
+#			if GLM_ARCH & GLM_ARCH_SSE41_BIT
+				__m128 const add0 = _mm_dp_ps(mul0, _mm_set_ps(1.0f, -1.0f,  1.0f,  1.0f), 0xff);
+				__m128 const add1 = _mm_dp_ps(mul1, _mm_set_ps(1.0f,  1.0f,  1.0f, -1.0f), 0xff);
+				__m128 const add2 = _mm_dp_ps(mul2, _mm_set_ps(1.0f,  1.0f, -1.0f,  1.0f), 0xff);
+				__m128 const add3 = _mm_dp_ps(mul3, _mm_set_ps(1.0f, -1.0f, -1.0f, -1.0f), 0xff);
+#			else
+				__m128 const mul4 = _mm_mul_ps(mul0, _mm_set_ps(1.0f, -1.0f,  1.0f,  1.0f));
+				__m128 const add0 = _mm_add_ps(mul0, _mm_movehl_ps(mul4, mul4));
+				__m128 const add4 = _mm_add_ss(add0, _mm_shuffle_ps(add0, add0, 1));
+
+				__m128 const mul5 = _mm_mul_ps(mul1, _mm_set_ps(1.0f,  1.0f,  1.0f, -1.0f));
+				__m128 const add1 = _mm_add_ps(mul1, _mm_movehl_ps(mul5, mul5));
+				__m128 const add5 = _mm_add_ss(add1, _mm_shuffle_ps(add1, add1, 1));
+
+				__m128 const mul6 = _mm_mul_ps(mul2, _mm_set_ps(1.0f,  1.0f, -1.0f,  1.0f));
+				__m128 const add2 = _mm_add_ps(mul6, _mm_movehl_ps(mul6, mul6));
+				__m128 const add6 = _mm_add_ss(add2, _mm_shuffle_ps(add2, add2, 1));
+
+				__m128 const mul7 = _mm_mul_ps(mul3, _mm_set_ps(1.0f, -1.0f, -1.0f, -1.0f));
+				__m128 const add3 = _mm_add_ps(mul3, _mm_movehl_ps(mul7, mul7));
+				__m128 const add7 = _mm_add_ss(add3, _mm_shuffle_ps(add3, add3, 1));
+		#endif
+
+			// This SIMD code is a politically correct way of doing this, but in every test I've tried it has been slower than
+			// the final code below. I'll keep this here for reference - maybe somebody else can do something better...
+			//
+			//__m128 xxyy = _mm_shuffle_ps(add4, add5, _MM_SHUFFLE(0, 0, 0, 0));
+			//__m128 zzww = _mm_shuffle_ps(add6, add7, _MM_SHUFFLE(0, 0, 0, 0));
+			//
+			//return _mm_shuffle_ps(xxyy, zzww, _MM_SHUFFLE(2, 0, 2, 0));
+
+			tquat<float, P> Result(uninitialize);
+			_mm_store_ss(&Result.x, add4);
+			_mm_store_ss(&Result.y, add5);
+			_mm_store_ss(&Result.z, add6);
+			_mm_store_ss(&Result.w, add7);
+			return Result;
+		}
+	};
+*/
+
+	template <precision P>
+	struct compute_dot<tquat, float, P, true>
+	{
+		static GLM_FUNC_QUALIFIER float call(tquat<float, P> const& x, tquat<float, P> const& y)
+		{
+			return _mm_cvtss_f32(glm_vec1_dot(x.data, y.data));
+		}
+	};
+
+	template <precision P>
+	struct compute_quat_add<float, P, true>
+	{
+		static tquat<float, P> call(tquat<float, P> const& q, tquat<float, P> const& p)
+		{
+			tquat<float, P> Result(uninitialize);
+			Result.data = _mm_add_ps(q.data, p.data);
+			return Result;
+		}
+	};
+
+#	if GLM_ARCH & GLM_ARCH_AVX_BIT
+	template <precision P>
+	struct compute_quat_add<double, P, true>
+	{
+		static tquat<double, P> call(tquat<double, P> const & a, tquat<double, P> const & b)
+		{
+			tquat<double, P> Result(uninitialize);
+			Result.data = _mm256_add_pd(a.data, b.data);
+			return Result;
+		}
+	};
+#	endif
+
+	template <precision P>
+	struct compute_quat_sub<float, P, true>
+	{
+		static tquat<float, P> call(tquat<float, P> const& q, tquat<float, P> const& p)
+		{
+			tvec4<float, P> Result(uninitialize);
+			Result.data = _mm_sub_ps(q.data, p.data);
+			return Result;
+		}
+	};
+
+#	if GLM_ARCH & GLM_ARCH_AVX_BIT
+	template <precision P>
+	struct compute_quat_sub<double, P, true>
+	{
+		static tquat<double, P> call(tquat<double, P> const & a, tquat<double, P> const & b)
+		{
+			tquat<double, P> Result(uninitialize);
+			Result.data = _mm256_sub_pd(a.data, b.data);
+			return Result;
+		}
+	};
+#	endif
+
+	template <precision P>
+	struct compute_quat_mul_scalar<float, P, true>
+	{
+		static tquat<float, P> call(tquat<float, P> const& q, float s)
+		{
+			tvec4<float, P> Result(uninitialize);
+			Result.data = _mm_mul_ps(q.data, _mm_set_ps1(s));
+			return Result;
+		}
+	};
+
+#	if GLM_ARCH & GLM_ARCH_AVX_BIT
+	template <precision P>
+	struct compute_quat_mul_scalar<double, P, true>
+	{
+		static tquat<double, P> call(tquat<double, P> const& q, double s)
+		{
+			tquat<double, P> Result(uninitialize);
+			Result.data = _mm256_mul_pd(q.data, _mm_set_ps1(s));
+			return Result;
+		}
+	};
+#	endif
+
+	template <precision P>
+	struct compute_quat_div_scalar<float, P, true>
+	{
+		static tquat<float, P> call(tquat<float, P> const& q, float s)
+		{
+			tvec4<float, P> Result(uninitialize);
+			Result.data = _mm_div_ps(q.data, _mm_set_ps1(s));
+			return Result;
+		}
+	};
+
+#	if GLM_ARCH & GLM_ARCH_AVX_BIT
+	template <precision P>
+	struct compute_quat_div_scalar<double, P, true>
+	{
+		static tquat<double, P> call(tquat<double, P> const& q, double s)
+		{
+			tquat<double, P> Result(uninitialize);
+			Result.data = _mm256_div_pd(q.data, _mm_set_ps1(s));
+			return Result;
+		}
+	};
+#	endif
+
+	template <precision P>
+	struct compute_quat_mul_vec4<float, P, true>
+	{
+		static tvec4<float, P> call(tquat<float, P> const& q, tvec4<float, P> const& v)
+		{
+			__m128 const q_wwww = _mm_shuffle_ps(q.data, q.data, _MM_SHUFFLE(3, 3, 3, 3));
+			__m128 const q_swp0 = _mm_shuffle_ps(q.data, q.data, _MM_SHUFFLE(3, 0, 2, 1));
+			__m128 const q_swp1 = _mm_shuffle_ps(q.data, q.data, _MM_SHUFFLE(3, 1, 0, 2));
+			__m128 const v_swp0 = _mm_shuffle_ps(v.data, v.data, _MM_SHUFFLE(3, 0, 2, 1));
+			__m128 const v_swp1 = _mm_shuffle_ps(v.data, v.data, _MM_SHUFFLE(3, 1, 0, 2));
+	
+			__m128 uv      = _mm_sub_ps(_mm_mul_ps(q_swp0, v_swp1), _mm_mul_ps(q_swp1, v_swp0));
+			__m128 uv_swp0 = _mm_shuffle_ps(uv, uv, _MM_SHUFFLE(3, 0, 2, 1));
+			__m128 uv_swp1 = _mm_shuffle_ps(uv, uv, _MM_SHUFFLE(3, 1, 0, 2));
+			__m128 uuv     = _mm_sub_ps(_mm_mul_ps(q_swp0, uv_swp1), _mm_mul_ps(q_swp1, uv_swp0));
+
+			__m128 const two = _mm_set1_ps(2.0f);
+			uv  = _mm_mul_ps(uv, _mm_mul_ps(q_wwww, two));
+			uuv = _mm_mul_ps(uuv, two);
+
+			tvec4<float, P> Result(uninitialize);
+			Result.data = _mm_add_ps(v.Data, _mm_add_ps(uv, uuv));
+			return Result;
+		}
+	};
+}//namespace detail
+}//namespace glm
+
+#endif//GLM_ARCH & GLM_ARCH_SSE2_BIT
+

data/ext/pixo/ext/glm/glm/gtc/random.hpp

@@ -0,0 +1,98 @@
+/// @ref gtc_random
+/// @file glm/gtc/random.hpp
+///
+/// @see core (dependence)
+/// @see gtc_half_float (dependence)
+/// @see gtx_random (extended)
+///
+/// @defgroup gtc_random GLM_GTC_random
+/// @ingroup gtc
+///
+/// @brief Generate random number from various distribution methods.
+///
+/// <glm/gtc/random.hpp> need to be included to use these functionalities.
+
+#pragma once
+
+// Dependency:
+#include "../vec2.hpp"
+#include "../vec3.hpp"
+
+#if GLM_MESSAGES == GLM_MESSAGES_ENABLED && !defined(GLM_EXT_INCLUDED)
+#	pragma message("GLM: GLM_GTC_random extension included")
+#endif
+
+namespace glm
+{
+	/// @addtogroup gtc_random
+	/// @{
+	
+	/// Generate random numbers in the interval [Min, Max], according a linear distribution 
+	/// 
+	/// @param Min 
+	/// @param Max 
+	/// @tparam genType Value type. Currently supported: float or double scalars.
+	/// @see gtc_random
+	template <typename genTYpe>
+	GLM_FUNC_DECL genTYpe linearRand(
+		genTYpe Min,
+		genTYpe Max);
+
+	/// Generate random numbers in the interval [Min, Max], according a linear distribution 
+	/// 
+	/// @param Min 
+	/// @param Max 
+	/// @tparam T Value type. Currently supported: float or double.
+	/// @tparam vecType A vertor type: tvec1, tvec2, tvec3, tvec4 or compatible
+	/// @see gtc_random
+	template <typename T, precision P, template <typename, precision> class vecType>
+	GLM_FUNC_DECL vecType<T, P> linearRand(
+		vecType<T, P> const & Min,
+		vecType<T, P> const & Max);
+
+	/// Generate random numbers in the interval [Min, Max], according a gaussian distribution 
+	/// 
+	/// @param Mean
+	/// @param Deviation
+	/// @see gtc_random
+	template <typename genType>
+	GLM_FUNC_DECL genType gaussRand(
+		genType Mean,
+		genType Deviation);
+	
+	/// Generate a random 2D vector which coordinates are regulary distributed on a circle of a given radius
+	/// 
+	/// @param Radius 
+	/// @see gtc_random
+	template <typename T>
+	GLM_FUNC_DECL tvec2<T, defaultp> circularRand(
+		T Radius);
+	
+	/// Generate a random 3D vector which coordinates are regulary distributed on a sphere of a given radius
+	/// 
+	/// @param Radius
+	/// @see gtc_random
+	template <typename T>
+	GLM_FUNC_DECL tvec3<T, defaultp> sphericalRand(
+		T Radius);
+	
+	/// Generate a random 2D vector which coordinates are regulary distributed within the area of a disk of a given radius
+	/// 
+	/// @param Radius
+	/// @see gtc_random
+	template <typename T>
+	GLM_FUNC_DECL tvec2<T, defaultp> diskRand(
+		T Radius);
+	
+	/// Generate a random 3D vector which coordinates are regulary distributed within the volume of a ball of a given radius
+	/// 
+	/// @param Radius
+	/// @see gtc_random
+	template <typename T>
+	GLM_FUNC_DECL tvec3<T, defaultp> ballRand(
+		T Radius);
+	
+	/// @}
+}//namespace glm
+
+#include "random.inl"

data/ext/pixo/ext/glm/glm/gtc/random.inl

@@ -0,0 +1,350 @@
+/// @ref gtc_random
+/// @file glm/gtc/random.inl
+
+#include "../geometric.hpp"
+#include "../exponential.hpp"
+#include <cstdlib>
+#include <ctime>
+#include <cassert>
+
+namespace glm{
+namespace detail
+{
+	template <typename T, precision P, template <class, precision> class vecType>
+	struct compute_rand
+	{
+		GLM_FUNC_QUALIFIER static vecType<T, P> call();
+	};
+
+	template <precision P>
+	struct compute_rand<uint8, P, tvec1>
+	{
+		GLM_FUNC_QUALIFIER static tvec1<uint8, P> call()
+		{
+			return tvec1<uint8, P>(
+				std::rand() % std::numeric_limits<uint8>::max());
+		}
+	};
+
+	template <precision P>
+	struct compute_rand<uint8, P, tvec2>
+	{
+		GLM_FUNC_QUALIFIER static tvec2<uint8, P> call()
+		{
+			return tvec2<uint8, P>(
+				std::rand() % std::numeric_limits<uint8>::max(),
+				std::rand() % std::numeric_limits<uint8>::max());
+		}
+	};
+
+	template <precision P>
+	struct compute_rand<uint8, P, tvec3>
+	{
+		GLM_FUNC_QUALIFIER static tvec3<uint8, P> call()
+		{
+			return tvec3<uint8, P>(
+				std::rand() % std::numeric_limits<uint8>::max(),
+				std::rand() % std::numeric_limits<uint8>::max(),
+				std::rand() % std::numeric_limits<uint8>::max());
+		}
+	};
+
+	template <precision P>
+	struct compute_rand<uint8, P, tvec4>
+	{
+		GLM_FUNC_QUALIFIER static tvec4<uint8, P> call()
+		{
+			return tvec4<uint8, P>(
+				std::rand() % std::numeric_limits<uint8>::max(),
+				std::rand() % std::numeric_limits<uint8>::max(),
+				std::rand() % std::numeric_limits<uint8>::max(),
+				std::rand() % std::numeric_limits<uint8>::max());
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_rand<uint16, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<uint16, P> call()
+		{
+			return
+				(vecType<uint16, P>(compute_rand<uint8, P, vecType>::call()) << static_cast<uint16>(8)) |
+				(vecType<uint16, P>(compute_rand<uint8, P, vecType>::call()) << static_cast<uint16>(0));
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_rand<uint32, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<uint32, P> call()
+		{
+			return
+				(vecType<uint32, P>(compute_rand<uint16, P, vecType>::call()) << static_cast<uint32>(16)) |
+				(vecType<uint32, P>(compute_rand<uint16, P, vecType>::call()) << static_cast<uint32>(0));
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_rand<uint64, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<uint64, P> call()
+		{
+			return
+				(vecType<uint64, P>(compute_rand<uint32, P, vecType>::call()) << static_cast<uint64>(32)) |
+				(vecType<uint64, P>(compute_rand<uint32, P, vecType>::call()) << static_cast<uint64>(0));
+		}
+	};
+
+	template <typename T, precision P, template <class, precision> class vecType>
+	struct compute_linearRand
+	{
+		GLM_FUNC_QUALIFIER static vecType<T, P> call(vecType<T, P> const & Min, vecType<T, P> const & Max);
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<int8, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<int8, P> call(vecType<int8, P> const & Min, vecType<int8, P> const & Max)
+		{
+			return (vecType<int8, P>(compute_rand<uint8, P, vecType>::call() % vecType<uint8, P>(Max + static_cast<int8>(1) - Min))) + Min;
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<uint8, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<uint8, P> call(vecType<uint8, P> const & Min, vecType<uint8, P> const & Max)
+		{
+			return (compute_rand<uint8, P, vecType>::call() % (Max + static_cast<uint8>(1) - Min)) + Min;
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<int16, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<int16, P> call(vecType<int16, P> const & Min, vecType<int16, P> const & Max)
+		{
+			return (vecType<int16, P>(compute_rand<uint16, P, vecType>::call() % vecType<uint16, P>(Max + static_cast<int16>(1) - Min))) + Min;
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<uint16, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<uint16, P> call(vecType<uint16, P> const & Min, vecType<uint16, P> const & Max)
+		{
+			return (compute_rand<uint16, P, vecType>::call() % (Max + static_cast<uint16>(1) - Min)) + Min;
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<int32, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<int32, P> call(vecType<int32, P> const & Min, vecType<int32, P> const & Max)
+		{
+			return (vecType<int32, P>(compute_rand<uint32, P, vecType>::call() % vecType<uint32, P>(Max + static_cast<int32>(1) - Min))) + Min;
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<uint32, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<uint32, P> call(vecType<uint32, P> const & Min, vecType<uint32, P> const & Max)
+		{
+			return (compute_rand<uint32, P, vecType>::call() % (Max + static_cast<uint32>(1) - Min)) + Min;
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<int64, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<int64, P> call(vecType<int64, P> const & Min, vecType<int64, P> const & Max)
+		{
+			return (vecType<int64, P>(compute_rand<uint64, P, vecType>::call() % vecType<uint64, P>(Max + static_cast<int64>(1) - Min))) + Min;
+		}
+	};
+
+	template <precision P, template <class, precision> class vecType>
+	struct compute_linearRand<uint64, P, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<uint64, P> call(vecType<uint64, P> const & Min, vecType<uint64, P> const & Max)
+		{
+			return (compute_rand<uint64, P, vecType>::call() % (Max + static_cast<uint64>(1) - Min)) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<float, lowp, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<float, lowp> call(vecType<float, lowp> const & Min, vecType<float, lowp> const & Max)
+		{
+			return vecType<float, lowp>(compute_rand<uint8, lowp, vecType>::call()) / static_cast<float>(std::numeric_limits<uint8>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<float, mediump, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<float, mediump> call(vecType<float, mediump> const & Min, vecType<float, mediump> const & Max)
+		{
+			return vecType<float, mediump>(compute_rand<uint16, mediump, vecType>::call()) / static_cast<float>(std::numeric_limits<uint16>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<float, highp, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<float, highp> call(vecType<float, highp> const & Min, vecType<float, highp> const & Max)
+		{
+			return vecType<float, highp>(compute_rand<uint32, highp, vecType>::call()) / static_cast<float>(std::numeric_limits<uint32>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<double, lowp, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<double, lowp> call(vecType<double, lowp> const & Min, vecType<double, lowp> const & Max)
+		{
+			return vecType<double, lowp>(compute_rand<uint16, lowp, vecType>::call()) / static_cast<double>(std::numeric_limits<uint16>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<double, mediump, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<double, mediump> call(vecType<double, mediump> const & Min, vecType<double, mediump> const & Max)
+		{
+			return vecType<double, mediump>(compute_rand<uint32, mediump, vecType>::call()) / static_cast<double>(std::numeric_limits<uint32>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<double, highp, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<double, highp> call(vecType<double, highp> const & Min, vecType<double, highp> const & Max)
+		{
+			return vecType<double, highp>(compute_rand<uint64, highp, vecType>::call()) / static_cast<double>(std::numeric_limits<uint64>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<long double, lowp, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<long double, lowp> call(vecType<long double, lowp> const & Min, vecType<long double, lowp> const & Max)
+		{
+			return vecType<long double, lowp>(compute_rand<uint32, lowp, vecType>::call()) / static_cast<long double>(std::numeric_limits<uint32>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<long double, mediump, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<long double, mediump> call(vecType<long double, mediump> const & Min, vecType<long double, mediump> const & Max)
+		{
+			return vecType<long double, mediump>(compute_rand<uint64, mediump, vecType>::call()) / static_cast<long double>(std::numeric_limits<uint64>::max()) * (Max - Min) + Min;
+		}
+	};
+
+	template <template <class, precision> class vecType>
+	struct compute_linearRand<long double, highp, vecType>
+	{
+		GLM_FUNC_QUALIFIER static vecType<long double, highp> call(vecType<long double, highp> const & Min, vecType<long double, highp> const & Max)
+		{
+			return vecType<long double, highp>(compute_rand<uint64, highp, vecType>::call()) / static_cast<long double>(std::numeric_limits<uint64>::max()) * (Max - Min) + Min;
+		}
+	};
+}//namespace detail
+
+	template <typename genType>
+	GLM_FUNC_QUALIFIER genType linearRand(genType Min, genType Max)
+	{
+		return detail::compute_linearRand<genType, highp, tvec1>::call(
+			tvec1<genType, highp>(Min),
+			tvec1<genType, highp>(Max)).x;
+	}
+
+	template <typename T, precision P, template <typename, precision> class vecType>
+	GLM_FUNC_QUALIFIER vecType<T, P> linearRand(vecType<T, P> const & Min, vecType<T, P> const & Max)
+	{
+		return detail::compute_linearRand<T, P, vecType>::call(Min, Max);
+	}
+
+	template <typename genType>
+	GLM_FUNC_QUALIFIER genType gaussRand(genType Mean, genType Deviation)
+	{
+		genType w, x1, x2;
+	
+		do
+		{
+			x1 = linearRand(genType(-1), genType(1));
+			x2 = linearRand(genType(-1), genType(1));
+		
+			w = x1 * x1 + x2 * x2;
+		} while(w > genType(1));
+	
+		return x2 * Deviation * Deviation * sqrt((genType(-2) * log(w)) / w) + Mean;
+	}
+
+	template <typename T, precision P, template <typename, precision> class vecType>
+	GLM_FUNC_QUALIFIER vecType<T, P> gaussRand(vecType<T, P> const & Mean, vecType<T, P> const & Deviation)
+	{
+		return detail::functor2<T, P, vecType>::call(gaussRand, Mean, Deviation);
+	}
+
+	template <typename T>
+	GLM_FUNC_QUALIFIER tvec2<T, defaultp> diskRand(T Radius)
+	{		
+		tvec2<T, defaultp> Result(T(0));
+		T LenRadius(T(0));
+		
+		do
+		{
+			Result = linearRand(
+				tvec2<T, defaultp>(-Radius),
+				tvec2<T, defaultp>(Radius));
+			LenRadius = length(Result);
+		}
+		while(LenRadius > Radius);
+		
+		return Result;
+	}
+	
+	template <typename T>
+	GLM_FUNC_QUALIFIER tvec3<T, defaultp> ballRand(T Radius)
+	{		
+		tvec3<T, defaultp> Result(T(0));
+		T LenRadius(T(0));
+		
+		do
+		{
+			Result = linearRand(
+				tvec3<T, defaultp>(-Radius),
+				tvec3<T, defaultp>(Radius));
+			LenRadius = length(Result);
+		}
+		while(LenRadius > Radius);
+		
+		return Result;
+	}
+	
+	template <typename T>
+	GLM_FUNC_QUALIFIER tvec2<T, defaultp> circularRand(T Radius)
+	{
+		T a = linearRand(T(0), T(6.283185307179586476925286766559f));
+		return tvec2<T, defaultp>(cos(a), sin(a)) * Radius;		
+	}
+	
+	template <typename T>
+	GLM_FUNC_QUALIFIER tvec3<T, defaultp> sphericalRand(T Radius)
+	{
+		T z = linearRand(T(-1), T(1));
+		T a = linearRand(T(0), T(6.283185307179586476925286766559f));
+	
+		T r = sqrt(T(1) - z * z);
+	
+		T x = r * cos(a);
+		T y = r * sin(a);
+	
+		return tvec3<T, defaultp>(x, y, z) * Radius;	
+	}
+}//namespace glm