Christophe Riccio
ago%!(EXTRA string=14 years)
commit
fe41daa938
47 changed files with 3781 additions and 1034 deletions
@ -0,0 +1,81 @@ |
|||||||
|
///////////////////////////////////////////////////////////////////////////////////
|
||||||
|
/// OpenGL Mathematics (glm.g-truc.net)
|
||||||
|
///
|
||||||
|
/// Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net)
|
||||||
|
/// Permission is hereby granted, free of charge, to any person obtaining a copy
|
||||||
|
/// of this software and associated documentation files (the "Software"), to deal
|
||||||
|
/// in the Software without restriction, including without limitation the rights
|
||||||
|
/// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
||||||
|
/// copies of the Software, and to permit persons to whom the Software is
|
||||||
|
/// furnished to do so, subject to the following conditions:
|
||||||
|
///
|
||||||
|
/// The above copyright notice and this permission notice shall be included in
|
||||||
|
/// all copies or substantial portions of the Software.
|
||||||
|
///
|
||||||
|
/// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
||||||
|
/// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
||||||
|
/// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
||||||
|
/// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
||||||
|
/// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
||||||
|
/// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
|
||||||
|
/// THE SOFTWARE.
|
||||||
|
///
|
||||||
|
/// @ref gtc_noise
|
||||||
|
/// @file glm/gtc/noise.hpp
|
||||||
|
/// @date 2011-04-21 / 2011-09-27
|
||||||
|
/// @author Christophe Riccio
|
||||||
|
///
|
||||||
|
/// @see core (dependence)
|
||||||
|
///
|
||||||
|
/// @defgroup gtx_noise GLM_GTX_noise: Procedural noise functions
|
||||||
|
/// @ingroup gtx
|
||||||
|
///
|
||||||
|
/// Defines 2D, 3D and 4D procedural noise functions
|
||||||
|
/// Based on the work of Stefan Gustavson and Ashima Arts on "webgl-noise":
|
||||||
|
/// https://github.com/ashima/webgl-noise
|
||||||
|
/// Following Stefan Gustavson's paper "Simplex noise demystified":
|
||||||
|
/// http://www.itn.liu.se/~stegu/simplexnoise/simplexnoise.pdf
|
||||||
|
/// Defines the half-precision floating-point type, along with various typedefs for vectors and matrices.
|
||||||
|
/// <glm/gtx/noise.hpp> need to be included to use these functionalities.
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////
|
||||||
|
|
||||||
|
#ifndef GLM_GTC_noise |
||||||
|
#define GLM_GTC_noise GLM_VERSION |
||||||
|
|
||||||
|
// Dependency:
|
||||||
|
#include "../glm.hpp" |
||||||
|
|
||||||
|
#if(defined(GLM_MESSAGES) && !defined(glm_ext)) |
||||||
|
# pragma message("GLM: GLM_GTC_noise extension included") |
||||||
|
#endif |
||||||
|
|
||||||
|
namespace glm |
||||||
|
{ |
||||||
|
/// @addtogroup gtx_noise
|
||||||
|
/// @{
|
||||||
|
|
||||||
|
/// Classic perlin noise.
|
||||||
|
/// From GLM_GTC_noise extension.
|
||||||
|
template <typename T, template<typename> class vecType>
|
||||||
|
T perlin( |
||||||
|
vecType<T> const & p); |
||||||
|
|
||||||
|
/// Periodic perlin noise.
|
||||||
|
/// From GLM_GTC_noise extension.
|
||||||
|
template <typename T, template<typename> class vecType>
|
||||||
|
T perlin( |
||||||
|
vecType<T> const & p,
|
||||||
|
vecType<T> const & rep); |
||||||
|
|
||||||
|
/// Simplex noise.
|
||||||
|
/// From GLM_GTC_noise extension.
|
||||||
|
template <typename T, template<typename> class vecType>
|
||||||
|
T simplex( |
||||||
|
vecType<T> const & p); |
||||||
|
|
||||||
|
/// @}
|
||||||
|
}//namespace glm
|
||||||
|
|
||||||
|
#include "noise.inl" |
||||||
|
|
||||||
|
#endif//GLM_GTC_noise
|
||||||
@ -0,0 +1,852 @@ |
|||||||
|
/////////////////////////////////////////////////////////////////////////////////////////////////// |
||||||
|
// OpenGL Mathematics Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net) |
||||||
|
/////////////////////////////////////////////////////////////////////////////////////////////////// |
||||||
|
// Based on the work of Stefan Gustavson and Ashima Arts on "webgl-noise": |
||||||
|
// https://github.com/ashima/webgl-noise |
||||||
|
// Following Stefan Gustavson's paper "Simplex noise demystified": |
||||||
|
// http://www.itn.liu.se/~stegu/simplexnoise/simplexnoise.pdf |
||||||
|
/////////////////////////////////////////////////////////////////////////////////////////////////// |
||||||
|
// Created : 2011-04-21 |
||||||
|
// Updated : 2011-09-27 |
||||||
|
// Licence : This source is under MIT License |
||||||
|
// File : glm/gtc/noise.inl |
||||||
|
/////////////////////////////////////////////////////////////////////////////////////////////////// |
||||||
|
// Dependency: |
||||||
|
// - GLM core |
||||||
|
/////////////////////////////////////////////////////////////////////////////////////////////////// |
||||||
|
|
||||||
|
namespace glm{ |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T mod289(T const & x) |
||||||
|
{ |
||||||
|
return x - floor(x * T(1.0 / 289.0)) * T(289.0); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T permute(T const & x) |
||||||
|
{ |
||||||
|
return mod289(((x * T(34)) + T(1)) * x); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T, template<typename> class vecType> |
||||||
|
GLM_FUNC_QUALIFIER vecType<T> permute(vecType<T> const & x) |
||||||
|
{ |
||||||
|
return mod289(((x * T(34)) + T(1)) * x); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T taylorInvSqrt(T const & r) |
||||||
|
{ |
||||||
|
return T(1.79284291400159) - T(0.85373472095314) * r; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T, template<typename> class vecType> |
||||||
|
GLM_FUNC_QUALIFIER vecType<T> taylorInvSqrt(vecType<T> const & r) |
||||||
|
{ |
||||||
|
return T(1.79284291400159) - T(0.85373472095314) * r; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T, template <typename> class vecType> |
||||||
|
GLM_FUNC_QUALIFIER vecType<T> fade(vecType<T> const & t) |
||||||
|
{ |
||||||
|
return t * t * t * (t * (t * T(6) - T(15)) + T(10)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<T> grad4(T const & j, detail::tvec4<T> const & ip) |
||||||
|
{ |
||||||
|
detail::tvec3<T> pXYZ = floor(fract(detail::tvec3<T>(j) * detail::tvec3<T>(ip)) * T(7)) * ip[2] - T(1); |
||||||
|
T pW = T(1.5) - dot(abs(pXYZ), detail::tvec3<T>(1)); |
||||||
|
detail::tvec4<T> s = detail::tvec4<T>(lessThan(detail::tvec4<T>(pXYZ, pW), detail::tvec4<T>(0.0))); |
||||||
|
pXYZ = pXYZ + (detail::tvec3<T>(s) * T(2) - T(1)) * s.w; |
||||||
|
return detail::tvec4<T>(pXYZ, pW); |
||||||
|
} |
||||||
|
|
||||||
|
// Classic Perlin noise |
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T perlin(detail::tvec2<T> const & P) |
||||||
|
{ |
||||||
|
detail::tvec4<T> Pi = glm::floor(detail::tvec4<T>(P.x, P.y, P.x, P.y)) + detail::tvec4<T>(0.0, 0.0, 1.0, 1.0); |
||||||
|
detail::tvec4<T> Pf = glm::fract(detail::tvec4<T>(P.x, P.y, P.x, P.y)) - detail::tvec4<T>(0.0, 0.0, 1.0, 1.0); |
||||||
|
Pi = mod(Pi, T(289)); // To avoid truncation effects in permutation |
||||||
|
detail::tvec4<T> ix(Pi.x, Pi.z, Pi.x, Pi.z); |
||||||
|
detail::tvec4<T> iy(Pi.y, Pi.y, Pi.w, Pi.w); |
||||||
|
detail::tvec4<T> fx(Pf.x, Pf.z, Pf.x, Pf.z); |
||||||
|
detail::tvec4<T> fy(Pf.y, Pf.y, Pf.w, Pf.w); |
||||||
|
|
||||||
|
detail::tvec4<T> i = glm::permute(glm::permute(ix) + iy); |
||||||
|
|
||||||
|
detail::tvec4<T> gx = T(2) * glm::fract(i / T(41)) - T(1); |
||||||
|
detail::tvec4<T> gy = glm::abs(gx) - T(0.5); |
||||||
|
detail::tvec4<T> tx = glm::floor(gx + T(0.5)); |
||||||
|
gx = gx - tx; |
||||||
|
|
||||||
|
detail::tvec2<T> g00(gx.x, gy.x); |
||||||
|
detail::tvec2<T> g10(gx.y, gy.y); |
||||||
|
detail::tvec2<T> g01(gx.z, gy.z); |
||||||
|
detail::tvec2<T> g11(gx.w, gy.w); |
||||||
|
|
||||||
|
detail::tvec4<T> norm = taylorInvSqrt(detail::tvec4<T>(dot(g00, g00), dot(g01, g01), dot(g10, g10), dot(g11, g11))); |
||||||
|
g00 *= norm.x; |
||||||
|
g01 *= norm.y; |
||||||
|
g10 *= norm.z; |
||||||
|
g11 *= norm.w; |
||||||
|
|
||||||
|
T n00 = dot(g00, detail::tvec2<T>(fx.x, fy.x)); |
||||||
|
T n10 = dot(g10, detail::tvec2<T>(fx.y, fy.y)); |
||||||
|
T n01 = dot(g01, detail::tvec2<T>(fx.z, fy.z)); |
||||||
|
T n11 = dot(g11, detail::tvec2<T>(fx.w, fy.w)); |
||||||
|
|
||||||
|
detail::tvec2<T> fade_xy = fade(detail::tvec2<T>(Pf.x, Pf.y)); |
||||||
|
detail::tvec2<T> n_x = mix(detail::tvec2<T>(n00, n01), detail::tvec2<T>(n10, n11), fade_xy.x); |
||||||
|
T n_xy = mix(n_x.x, n_x.y, fade_xy.y); |
||||||
|
return T(2.3) * n_xy; |
||||||
|
} |
||||||
|
|
||||||
|
// Classic Perlin noise |
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T perlin(detail::tvec3<T> const & P) |
||||||
|
{ |
||||||
|
detail::tvec3<T> Pi0 = floor(P); // Integer part for indexing |
||||||
|
detail::tvec3<T> Pi1 = Pi0 + T(1); // Integer part + 1 |
||||||
|
Pi0 = mod289(Pi0); |
||||||
|
Pi1 = mod289(Pi1); |
||||||
|
detail::tvec3<T> Pf0 = fract(P); // Fractional part for interpolation |
||||||
|
detail::tvec3<T> Pf1 = Pf0 - T(1); // Fractional part - 1.0 |
||||||
|
detail::tvec4<T> ix(Pi0.x, Pi1.x, Pi0.x, Pi1.x); |
||||||
|
detail::tvec4<T> iy = detail::tvec4<T>(detail::tvec2<T>(Pi0.y), detail::tvec2<T>(Pi1.y)); |
||||||
|
detail::tvec4<T> iz0(Pi0.z); |
||||||
|
detail::tvec4<T> iz1(Pi1.z); |
||||||
|
|
||||||
|
detail::tvec4<T> ixy = permute(permute(ix) + iy); |
||||||
|
detail::tvec4<T> ixy0 = permute(ixy + iz0); |
||||||
|
detail::tvec4<T> ixy1 = permute(ixy + iz1); |
||||||
|
|
||||||
|
detail::tvec4<T> gx0 = ixy0 * T(1.0 / 7.0); |
||||||
|
detail::tvec4<T> gy0 = fract(floor(gx0) * T(1.0 / 7.0)) - T(0.5); |
||||||
|
gx0 = fract(gx0); |
||||||
|
detail::tvec4<T> gz0 = detail::tvec4<T>(0.5) - abs(gx0) - abs(gy0); |
||||||
|
detail::tvec4<T> sz0 = step(gz0, detail::tvec4<T>(0.0)); |
||||||
|
gx0 -= sz0 * (step(T(0), gx0) - T(0.5)); |
||||||
|
gy0 -= sz0 * (step(T(0), gy0) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx1 = ixy1 * T(1.0 / 7.0); |
||||||
|
detail::tvec4<T> gy1 = fract(floor(gx1) * T(1.0 / 7.0)) - T(0.5); |
||||||
|
gx1 = fract(gx1); |
||||||
|
detail::tvec4<T> gz1 = detail::tvec4<T>(0.5) - abs(gx1) - abs(gy1); |
||||||
|
detail::tvec4<T> sz1 = step(gz1, detail::tvec4<T>(0.0)); |
||||||
|
gx1 -= sz1 * (step(T(0), gx1) - T(0.5)); |
||||||
|
gy1 -= sz1 * (step(T(0), gy1) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec3<T> g000(gx0.x, gy0.x, gz0.x); |
||||||
|
detail::tvec3<T> g100(gx0.y, gy0.y, gz0.y); |
||||||
|
detail::tvec3<T> g010(gx0.z, gy0.z, gz0.z); |
||||||
|
detail::tvec3<T> g110(gx0.w, gy0.w, gz0.w); |
||||||
|
detail::tvec3<T> g001(gx1.x, gy1.x, gz1.x); |
||||||
|
detail::tvec3<T> g101(gx1.y, gy1.y, gz1.y); |
||||||
|
detail::tvec3<T> g011(gx1.z, gy1.z, gz1.z); |
||||||
|
detail::tvec3<T> g111(gx1.w, gy1.w, gz1.w); |
||||||
|
|
||||||
|
detail::tvec4<T> norm0 = taylorInvSqrt(detail::tvec4<T>(dot(g000, g000), dot(g010, g010), dot(g100, g100), dot(g110, g110))); |
||||||
|
g000 *= norm0.x; |
||||||
|
g010 *= norm0.y; |
||||||
|
g100 *= norm0.z; |
||||||
|
g110 *= norm0.w; |
||||||
|
detail::tvec4<T> norm1 = taylorInvSqrt(detail::tvec4<T>(dot(g001, g001), dot(g011, g011), dot(g101, g101), dot(g111, g111))); |
||||||
|
g001 *= norm1.x; |
||||||
|
g011 *= norm1.y; |
||||||
|
g101 *= norm1.z; |
||||||
|
g111 *= norm1.w; |
||||||
|
|
||||||
|
T n000 = dot(g000, Pf0); |
||||||
|
T n100 = dot(g100, detail::tvec3<T>(Pf1.x, Pf0.y, Pf0.z)); |
||||||
|
T n010 = dot(g010, detail::tvec3<T>(Pf0.x, Pf1.y, Pf0.z)); |
||||||
|
T n110 = dot(g110, detail::tvec3<T>(Pf1.x, Pf1.y, Pf0.z)); |
||||||
|
T n001 = dot(g001, detail::tvec3<T>(Pf0.x, Pf0.y, Pf1.z)); |
||||||
|
T n101 = dot(g101, detail::tvec3<T>(Pf1.x, Pf0.y, Pf1.z)); |
||||||
|
T n011 = dot(g011, detail::tvec3<T>(Pf0.x, Pf1.y, Pf1.z)); |
||||||
|
T n111 = dot(g111, Pf1); |
||||||
|
|
||||||
|
detail::tvec3<T> fade_xyz = fade(Pf0); |
||||||
|
detail::tvec4<T> n_z = mix(detail::tvec4<T>(n000, n100, n010, n110), detail::tvec4<T>(n001, n101, n011, n111), fade_xyz.z); |
||||||
|
detail::tvec2<T> n_yz = mix(detail::tvec2<T>(n_z.x, n_z.y), detail::tvec2<T>(n_z.z, n_z.w), fade_xyz.y); |
||||||
|
T n_xyz = mix(n_yz.x, n_yz.y, fade_xyz.x); |
||||||
|
return T(2.2) * n_xyz; |
||||||
|
} |
||||||
|
/* |
||||||
|
// Classic Perlin noise |
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T perlin(detail::tvec3<T> const & P) |
||||||
|
{ |
||||||
|
detail::tvec3<T> Pi0 = floor(P); // Integer part for indexing |
||||||
|
detail::tvec3<T> Pi1 = Pi0 + T(1); // Integer part + 1 |
||||||
|
Pi0 = mod(Pi0, T(289)); |
||||||
|
Pi1 = mod(Pi1, T(289)); |
||||||
|
detail::tvec3<T> Pf0 = fract(P); // Fractional part for interpolation |
||||||
|
detail::tvec3<T> Pf1 = Pf0 - T(1); // Fractional part - 1.0 |
||||||
|
detail::tvec4<T> ix(Pi0.x, Pi1.x, Pi0.x, Pi1.x); |
||||||
|
detail::tvec4<T> iy(Pi0.y, Pi0.y, Pi1.y, Pi1.y); |
||||||
|
detail::tvec4<T> iz0(Pi0.z); |
||||||
|
detail::tvec4<T> iz1(Pi1.z); |
||||||
|
|
||||||
|
detail::tvec4<T> ixy = permute(permute(ix) + iy); |
||||||
|
detail::tvec4<T> ixy0 = permute(ixy + iz0); |
||||||
|
detail::tvec4<T> ixy1 = permute(ixy + iz1); |
||||||
|
|
||||||
|
detail::tvec4<T> gx0 = ixy0 / T(7); |
||||||
|
detail::tvec4<T> gy0 = fract(floor(gx0) / T(7)) - T(0.5); |
||||||
|
gx0 = fract(gx0); |
||||||
|
detail::tvec4<T> gz0 = detail::tvec4<T>(0.5) - abs(gx0) - abs(gy0); |
||||||
|
detail::tvec4<T> sz0 = step(gz0, detail::tvec4<T>(0.0)); |
||||||
|
gx0 -= sz0 * (step(0.0, gx0) - T(0.5)); |
||||||
|
gy0 -= sz0 * (step(0.0, gy0) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx1 = ixy1 / T(7); |
||||||
|
detail::tvec4<T> gy1 = fract(floor(gx1) / T(7)) - T(0.5); |
||||||
|
gx1 = fract(gx1); |
||||||
|
detail::tvec4<T> gz1 = detail::tvec4<T>(0.5) - abs(gx1) - abs(gy1); |
||||||
|
detail::tvec4<T> sz1 = step(gz1, detail::tvec4<T>(0.0)); |
||||||
|
gx1 -= sz1 * (step(T(0), gx1) - T(0.5)); |
||||||
|
gy1 -= sz1 * (step(T(0), gy1) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec3<T> g000(gx0.x, gy0.x, gz0.x); |
||||||
|
detail::tvec3<T> g100(gx0.y, gy0.y, gz0.y); |
||||||
|
detail::tvec3<T> g010(gx0.z, gy0.z, gz0.z); |
||||||
|
detail::tvec3<T> g110(gx0.w, gy0.w, gz0.w); |
||||||
|
detail::tvec3<T> g001(gx1.x, gy1.x, gz1.x); |
||||||
|
detail::tvec3<T> g101(gx1.y, gy1.y, gz1.y); |
||||||
|
detail::tvec3<T> g011(gx1.z, gy1.z, gz1.z); |
||||||
|
detail::tvec3<T> g111(gx1.w, gy1.w, gz1.w); |
||||||
|
|
||||||
|
detail::tvec4<T> norm0 = taylorInvSqrt(detail::tvec4<T>(dot(g000, g000), dot(g010, g010), dot(g100, g100), dot(g110, g110))); |
||||||
|
g000 *= norm0.x; |
||||||
|
g010 *= norm0.y; |
||||||
|
g100 *= norm0.z; |
||||||
|
g110 *= norm0.w; |
||||||
|
detail::tvec4<T> norm1 = taylorInvSqrt(detail::tvec4<T>(dot(g001, g001), dot(g011, g011), dot(g101, g101), dot(g111, g111))); |
||||||
|
g001 *= norm1.x; |
||||||
|
g011 *= norm1.y; |
||||||
|
g101 *= norm1.z; |
||||||
|
g111 *= norm1.w; |
||||||
|
|
||||||
|
T n000 = dot(g000, Pf0); |
||||||
|
T n100 = dot(g100, detail::tvec3<T>(Pf1.x, Pf0.y, Pf0.z)); |
||||||
|
T n010 = dot(g010, detail::tvec3<T>(Pf0.x, Pf1.y, Pf0.z)); |
||||||
|
T n110 = dot(g110, detail::tvec3<T>(Pf1.x, Pf1.y, Pf0.z)); |
||||||
|
T n001 = dot(g001, detail::tvec3<T>(Pf0.x, Pf0.y, Pf1.z)); |
||||||
|
T n101 = dot(g101, detail::tvec3<T>(Pf1.x, Pf0.y, Pf1.z)); |
||||||
|
T n011 = dot(g011, detail::tvec3<T>(Pf0.x, Pf1.y, Pf1.z)); |
||||||
|
T n111 = dot(g111, Pf1); |
||||||
|
|
||||||
|
detail::tvec3<T> fade_xyz = fade(Pf0); |
||||||
|
detail::tvec4<T> n_z = mix(detail::tvec4<T>(n000, n100, n010, n110), detail::tvec4<T>(n001, n101, n011, n111), fade_xyz.z); |
||||||
|
detail::tvec2<T> n_yz = mix( |
||||||
|
detail::tvec2<T>(n_z.x, n_z.y), |
||||||
|
detail::tvec2<T>(n_z.z, n_z.w), fade_xyz.y); |
||||||
|
T n_xyz = mix(n_yz.x, n_yz.y, fade_xyz.x); |
||||||
|
return T(2.2) * n_xyz; |
||||||
|
} |
||||||
|
*/ |
||||||
|
// Classic Perlin noise |
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T perlin(detail::tvec4<T> const & P) |
||||||
|
{ |
||||||
|
detail::tvec4<T> Pi0 = floor(P); // Integer part for indexing |
||||||
|
detail::tvec4<T> Pi1 = Pi0 + T(1); // Integer part + 1 |
||||||
|
Pi0 = mod(Pi0, T(289)); |
||||||
|
Pi1 = mod(Pi1, T(289)); |
||||||
|
detail::tvec4<T> Pf0 = fract(P); // Fractional part for interpolation |
||||||
|
detail::tvec4<T> Pf1 = Pf0 - T(1); // Fractional part - 1.0 |
||||||
|
detail::tvec4<T> ix(Pi0.x, Pi1.x, Pi0.x, Pi1.x); |
||||||
|
detail::tvec4<T> iy(Pi0.y, Pi0.y, Pi1.y, Pi1.y); |
||||||
|
detail::tvec4<T> iz0(Pi0.z); |
||||||
|
detail::tvec4<T> iz1(Pi1.z); |
||||||
|
detail::tvec4<T> iw0(Pi0.w); |
||||||
|
detail::tvec4<T> iw1(Pi1.w); |
||||||
|
|
||||||
|
detail::tvec4<T> ixy = permute(permute(ix) + iy); |
||||||
|
detail::tvec4<T> ixy0 = permute(ixy + iz0); |
||||||
|
detail::tvec4<T> ixy1 = permute(ixy + iz1); |
||||||
|
detail::tvec4<T> ixy00 = permute(ixy0 + iw0); |
||||||
|
detail::tvec4<T> ixy01 = permute(ixy0 + iw1); |
||||||
|
detail::tvec4<T> ixy10 = permute(ixy1 + iw0); |
||||||
|
detail::tvec4<T> ixy11 = permute(ixy1 + iw1); |
||||||
|
|
||||||
|
detail::tvec4<T> gx00 = ixy00 / T(7); |
||||||
|
detail::tvec4<T> gy00 = floor(gx00) / T(7); |
||||||
|
detail::tvec4<T> gz00 = floor(gy00) / T(6); |
||||||
|
gx00 = fract(gx00) - T(0.5); |
||||||
|
gy00 = fract(gy00) - T(0.5); |
||||||
|
gz00 = fract(gz00) - T(0.5); |
||||||
|
detail::tvec4<T> gw00 = detail::tvec4<T>(0.75) - abs(gx00) - abs(gy00) - abs(gz00); |
||||||
|
detail::tvec4<T> sw00 = step(gw00, detail::tvec4<T>(0.0)); |
||||||
|
gx00 -= sw00 * (step(T(0), gx00) - T(0.5)); |
||||||
|
gy00 -= sw00 * (step(T(0), gy00) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx01 = ixy01 / T(7); |
||||||
|
detail::tvec4<T> gy01 = floor(gx01) / T(7); |
||||||
|
detail::tvec4<T> gz01 = floor(gy01) / T(6); |
||||||
|
gx01 = fract(gx01) - T(0.5); |
||||||
|
gy01 = fract(gy01) - T(0.5); |
||||||
|
gz01 = fract(gz01) - T(0.5); |
||||||
|
detail::tvec4<T> gw01 = detail::tvec4<T>(0.75) - abs(gx01) - abs(gy01) - abs(gz01); |
||||||
|
detail::tvec4<T> sw01 = step(gw01, detail::tvec4<T>(0.0)); |
||||||
|
gx01 -= sw01 * (step(T(0), gx01) - T(0.5)); |
||||||
|
gy01 -= sw01 * (step(T(0), gy01) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx10 = ixy10 / T(7); |
||||||
|
detail::tvec4<T> gy10 = floor(gx10) / T(7); |
||||||
|
detail::tvec4<T> gz10 = floor(gy10) / T(6); |
||||||
|
gx10 = fract(gx10) - T(0.5); |
||||||
|
gy10 = fract(gy10) - T(0.5); |
||||||
|
gz10 = fract(gz10) - T(0.5); |
||||||
|
detail::tvec4<T> gw10 = detail::tvec4<T>(0.75) - abs(gx10) - abs(gy10) - abs(gz10); |
||||||
|
detail::tvec4<T> sw10 = step(gw10, detail::tvec4<T>(0)); |
||||||
|
gx10 -= sw10 * (step(T(0), gx10) - T(0.5)); |
||||||
|
gy10 -= sw10 * (step(T(0), gy10) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx11 = ixy11 / T(7); |
||||||
|
detail::tvec4<T> gy11 = floor(gx11) / T(7); |
||||||
|
detail::tvec4<T> gz11 = floor(gy11) / T(6); |
||||||
|
gx11 = fract(gx11) - T(0.5); |
||||||
|
gy11 = fract(gy11) - T(0.5); |
||||||
|
gz11 = fract(gz11) - T(0.5); |
||||||
|
detail::tvec4<T> gw11 = detail::tvec4<T>(0.75) - abs(gx11) - abs(gy11) - abs(gz11); |
||||||
|
detail::tvec4<T> sw11 = step(gw11, detail::tvec4<T>(0.0)); |
||||||
|
gx11 -= sw11 * (step(T(0), gx11) - T(0.5)); |
||||||
|
gy11 -= sw11 * (step(T(0), gy11) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> g0000(gx00.x, gy00.x, gz00.x, gw00.x); |
||||||
|
detail::tvec4<T> g1000(gx00.y, gy00.y, gz00.y, gw00.y); |
||||||
|
detail::tvec4<T> g0100(gx00.z, gy00.z, gz00.z, gw00.z); |
||||||
|
detail::tvec4<T> g1100(gx00.w, gy00.w, gz00.w, gw00.w); |
||||||
|
detail::tvec4<T> g0010(gx10.x, gy10.x, gz10.x, gw10.x); |
||||||
|
detail::tvec4<T> g1010(gx10.y, gy10.y, gz10.y, gw10.y); |
||||||
|
detail::tvec4<T> g0110(gx10.z, gy10.z, gz10.z, gw10.z); |
||||||
|
detail::tvec4<T> g1110(gx10.w, gy10.w, gz10.w, gw10.w); |
||||||
|
detail::tvec4<T> g0001(gx01.x, gy01.x, gz01.x, gw01.x); |
||||||
|
detail::tvec4<T> g1001(gx01.y, gy01.y, gz01.y, gw01.y); |
||||||
|
detail::tvec4<T> g0101(gx01.z, gy01.z, gz01.z, gw01.z); |
||||||
|
detail::tvec4<T> g1101(gx01.w, gy01.w, gz01.w, gw01.w); |
||||||
|
detail::tvec4<T> g0011(gx11.x, gy11.x, gz11.x, gw11.x); |
||||||
|
detail::tvec4<T> g1011(gx11.y, gy11.y, gz11.y, gw11.y); |
||||||
|
detail::tvec4<T> g0111(gx11.z, gy11.z, gz11.z, gw11.z); |
||||||
|
detail::tvec4<T> g1111(gx11.w, gy11.w, gz11.w, gw11.w); |
||||||
|
|
||||||
|
detail::tvec4<T> norm00 = taylorInvSqrt(detail::tvec4<T>(dot(g0000, g0000), dot(g0100, g0100), dot(g1000, g1000), dot(g1100, g1100))); |
||||||
|
g0000 *= norm00.x; |
||||||
|
g0100 *= norm00.y; |
||||||
|
g1000 *= norm00.z; |
||||||
|
g1100 *= norm00.w; |
||||||
|
|
||||||
|
detail::tvec4<T> norm01 = taylorInvSqrt(detail::tvec4<T>(dot(g0001, g0001), dot(g0101, g0101), dot(g1001, g1001), dot(g1101, g1101))); |
||||||
|
g0001 *= norm01.x; |
||||||
|
g0101 *= norm01.y; |
||||||
|
g1001 *= norm01.z; |
||||||
|
g1101 *= norm01.w; |
||||||
|
|
||||||
|
detail::tvec4<T> norm10 = taylorInvSqrt(detail::tvec4<T>(dot(g0010, g0010), dot(g0110, g0110), dot(g1010, g1010), dot(g1110, g1110))); |
||||||
|
g0010 *= norm10.x; |
||||||
|
g0110 *= norm10.y; |
||||||
|
g1010 *= norm10.z; |
||||||
|
g1110 *= norm10.w; |
||||||
|
|
||||||
|
detail::tvec4<T> norm11 = taylorInvSqrt(detail::tvec4<T>(dot(g0011, g0011), dot(g0111, g0111), dot(g1011, g1011), dot(g1111, g1111))); |
||||||
|
g0011 *= norm11.x; |
||||||
|
g0111 *= norm11.y; |
||||||
|
g1011 *= norm11.z; |
||||||
|
g1111 *= norm11.w; |
||||||
|
|
||||||
|
T n0000 = dot(g0000, Pf0); |
||||||
|
T n1000 = dot(g1000, detail::tvec4<T>(Pf1.x, Pf0.y, Pf0.z, Pf0.w)); |
||||||
|
T n0100 = dot(g0100, detail::tvec4<T>(Pf0.x, Pf1.y, Pf0.z, Pf0.w)); |
||||||
|
T n1100 = dot(g1100, detail::tvec4<T>(Pf1.x, Pf1.y, Pf0.z, Pf0.w)); |
||||||
|
T n0010 = dot(g0010, detail::tvec4<T>(Pf0.x, Pf0.y, Pf1.z, Pf0.w)); |
||||||
|
T n1010 = dot(g1010, detail::tvec4<T>(Pf1.x, Pf0.y, Pf1.z, Pf0.w)); |
||||||
|
T n0110 = dot(g0110, detail::tvec4<T>(Pf0.x, Pf1.y, Pf1.z, Pf0.w)); |
||||||
|
T n1110 = dot(g1110, detail::tvec4<T>(Pf1.x, Pf1.y, Pf1.z, Pf0.w)); |
||||||
|
T n0001 = dot(g0001, detail::tvec4<T>(Pf0.x, Pf0.y, Pf0.z, Pf1.w)); |
||||||
|
T n1001 = dot(g1001, detail::tvec4<T>(Pf1.x, Pf0.y, Pf0.z, Pf1.w)); |
||||||
|
T n0101 = dot(g0101, detail::tvec4<T>(Pf0.x, Pf1.y, Pf0.z, Pf1.w)); |
||||||
|
T n1101 = dot(g1101, detail::tvec4<T>(Pf1.x, Pf1.y, Pf0.z, Pf1.w)); |
||||||
|
T n0011 = dot(g0011, detail::tvec4<T>(Pf0.x, Pf0.y, Pf1.z, Pf1.w)); |
||||||
|
T n1011 = dot(g1011, detail::tvec4<T>(Pf1.x, Pf0.y, Pf1.z, Pf1.w)); |
||||||
|
T n0111 = dot(g0111, detail::tvec4<T>(Pf0.x, Pf1.y, Pf1.z, Pf1.w)); |
||||||
|
T n1111 = dot(g1111, Pf1); |
||||||
|
|
||||||
|
detail::tvec4<T> fade_xyzw = fade(Pf0); |
||||||
|
detail::tvec4<T> n_0w = mix(detail::tvec4<T>(n0000, n1000, n0100, n1100), detail::tvec4<T>(n0001, n1001, n0101, n1101), fade_xyzw.w); |
||||||
|
detail::tvec4<T> n_1w = mix(detail::tvec4<T>(n0010, n1010, n0110, n1110), detail::tvec4<T>(n0011, n1011, n0111, n1111), fade_xyzw.w); |
||||||
|
detail::tvec4<T> n_zw = mix(n_0w, n_1w, fade_xyzw.z); |
||||||
|
detail::tvec2<T> n_yzw = mix(detail::tvec2<T>(n_zw.x, n_zw.y), detail::tvec2<T>(n_zw.z, n_zw.w), fade_xyzw.y); |
||||||
|
T n_xyzw = mix(n_yzw.x, n_yzw.y, fade_xyzw.x); |
||||||
|
return T(2.2) * n_xyzw; |
||||||
|
} |
||||||
|
|
||||||
|
// Classic Perlin noise, periodic variant |
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T perlin(detail::tvec2<T> const & P, detail::tvec2<T> const & rep) |
||||||
|
{ |
||||||
|
detail::tvec4<T> Pi = floor(detail::tvec4<T>(P.x, P.y, P.x, P.y)) + detail::tvec4<T>(0.0, 0.0, 1.0, 1.0); |
||||||
|
detail::tvec4<T> Pf = fract(detail::tvec4<T>(P.x, P.y, P.x, P.y)) - detail::tvec4<T>(0.0, 0.0, 1.0, 1.0); |
||||||
|
Pi = mod(Pi, detail::tvec4<T>(rep.x, rep.y, rep.x, rep.y)); // To create noise with explicit period |
||||||
|
Pi = mod(Pi, T(289)); // To avoid truncation effects in permutation |
||||||
|
detail::tvec4<T> ix(Pi.x, Pi.z, Pi.x, Pi.z); |
||||||
|
detail::tvec4<T> iy(Pi.y, Pi.y, Pi.w, Pi.w); |
||||||
|
detail::tvec4<T> fx(Pf.x, Pf.z, Pf.x, Pf.z); |
||||||
|
detail::tvec4<T> fy(Pf.y, Pf.y, Pf.w, Pf.w); |
||||||
|
|
||||||
|
detail::tvec4<T> i = permute(permute(ix) + iy); |
||||||
|
|
||||||
|
detail::tvec4<T> gx = T(2) * fract(i / T(41)) - T(1); |
||||||
|
detail::tvec4<T> gy = abs(gx) - T(0.5); |
||||||
|
detail::tvec4<T> tx = floor(gx + T(0.5)); |
||||||
|
gx = gx - tx; |
||||||
|
|
||||||
|
detail::tvec2<T> g00(gx.x, gy.x); |
||||||
|
detail::tvec2<T> g10(gx.y, gy.y); |
||||||
|
detail::tvec2<T> g01(gx.z, gy.z); |
||||||
|
detail::tvec2<T> g11(gx.w, gy.w); |
||||||
|
|
||||||
|
detail::tvec4<T> norm = taylorInvSqrt(detail::tvec4<T>(dot(g00, g00), dot(g01, g01), dot(g10, g10), dot(g11, g11))); |
||||||
|
g00 *= norm.x; |
||||||
|
g01 *= norm.y; |
||||||
|
g10 *= norm.z; |
||||||
|
g11 *= norm.w; |
||||||
|
|
||||||
|
T n00 = dot(g00, detail::tvec2<T>(fx.x, fy.x)); |
||||||
|
T n10 = dot(g10, detail::tvec2<T>(fx.y, fy.y)); |
||||||
|
T n01 = dot(g01, detail::tvec2<T>(fx.z, fy.z)); |
||||||
|
T n11 = dot(g11, detail::tvec2<T>(fx.w, fy.w)); |
||||||
|
|
||||||
|
detail::tvec2<T> fade_xy = fade(detail::tvec2<T>(Pf.x, Pf.y)); |
||||||
|
detail::tvec2<T> n_x = mix(detail::tvec2<T>(n00, n01), detail::tvec2<T>(n10, n11), fade_xy.x); |
||||||
|
T n_xy = mix(n_x.x, n_x.y, fade_xy.y); |
||||||
|
return T(2.3) * n_xy; |
||||||
|
} |
||||||
|
|
||||||
|
// Classic Perlin noise, periodic variant |
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T perlin(detail::tvec3<T> const & P, detail::tvec3<T> const & rep) |
||||||
|
{ |
||||||
|
detail::tvec3<T> Pi0 = mod(floor(P), rep); // Integer part, modulo period |
||||||
|
detail::tvec3<T> Pi1 = mod(Pi0 + detail::tvec3<T>(1.0), rep); // Integer part + 1, mod period |
||||||
|
Pi0 = mod(Pi0, T(289)); |
||||||
|
Pi1 = mod(Pi1, T(289)); |
||||||
|
detail::tvec3<T> Pf0 = fract(P); // Fractional part for interpolation |
||||||
|
detail::tvec3<T> Pf1 = Pf0 - detail::tvec3<T>(1.0); // Fractional part - 1.0 |
||||||
|
detail::tvec4<T> ix = detail::tvec4<T>(Pi0.x, Pi1.x, Pi0.x, Pi1.x); |
||||||
|
detail::tvec4<T> iy = detail::tvec4<T>(Pi0.y, Pi0.y, Pi1.y, Pi1.y); |
||||||
|
detail::tvec4<T> iz0(Pi0.z); |
||||||
|
detail::tvec4<T> iz1(Pi1.z); |
||||||
|
|
||||||
|
detail::tvec4<T> ixy = permute(permute(ix) + iy); |
||||||
|
detail::tvec4<T> ixy0 = permute(ixy + iz0); |
||||||
|
detail::tvec4<T> ixy1 = permute(ixy + iz1); |
||||||
|
|
||||||
|
detail::tvec4<T> gx0 = ixy0 / T(7); |
||||||
|
detail::tvec4<T> gy0 = fract(floor(gx0) / T(7)) - T(0.5); |
||||||
|
gx0 = fract(gx0); |
||||||
|
detail::tvec4<T> gz0 = detail::tvec4<T>(0.5) - abs(gx0) - abs(gy0); |
||||||
|
detail::tvec4<T> sz0 = step(gz0, detail::tvec4<T>(0)); |
||||||
|
gx0 -= sz0 * (step(0.0, gx0) - T(0.5)); |
||||||
|
gy0 -= sz0 * (step(0.0, gy0) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx1 = ixy1 / T(7); |
||||||
|
detail::tvec4<T> gy1 = fract(floor(gx1) / T(7)) - T(0.5); |
||||||
|
gx1 = fract(gx1); |
||||||
|
detail::tvec4<T> gz1 = detail::tvec4<T>(0.5) - abs(gx1) - abs(gy1); |
||||||
|
detail::tvec4<T> sz1 = step(gz1, detail::tvec4<T>(0.0)); |
||||||
|
gx1 -= sz1 * (step(0.0, gx1) - T(0.5)); |
||||||
|
gy1 -= sz1 * (step(0.0, gy1) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec3<T> g000 = detail::tvec3<T>(gx0.x, gy0.x, gz0.x); |
||||||
|
detail::tvec3<T> g100 = detail::tvec3<T>(gx0.y, gy0.y, gz0.y); |
||||||
|
detail::tvec3<T> g010 = detail::tvec3<T>(gx0.z, gy0.z, gz0.z); |
||||||
|
detail::tvec3<T> g110 = detail::tvec3<T>(gx0.w, gy0.w, gz0.w); |
||||||
|
detail::tvec3<T> g001 = detail::tvec3<T>(gx1.x, gy1.x, gz1.x); |
||||||
|
detail::tvec3<T> g101 = detail::tvec3<T>(gx1.y, gy1.y, gz1.y); |
||||||
|
detail::tvec3<T> g011 = detail::tvec3<T>(gx1.z, gy1.z, gz1.z); |
||||||
|
detail::tvec3<T> g111 = detail::tvec3<T>(gx1.w, gy1.w, gz1.w); |
||||||
|
|
||||||
|
detail::tvec4<T> norm0 = taylorInvSqrt(detail::tvec4<T>(dot(g000, g000), dot(g010, g010), dot(g100, g100), dot(g110, g110))); |
||||||
|
g000 *= norm0.x; |
||||||
|
g010 *= norm0.y; |
||||||
|
g100 *= norm0.z; |
||||||
|
g110 *= norm0.w; |
||||||
|
detail::tvec4<T> norm1 = taylorInvSqrt(detail::tvec4<T>(dot(g001, g001), dot(g011, g011), dot(g101, g101), dot(g111, g111))); |
||||||
|
g001 *= norm1.x; |
||||||
|
g011 *= norm1.y; |
||||||
|
g101 *= norm1.z; |
||||||
|
g111 *= norm1.w; |
||||||
|
|
||||||
|
T n000 = dot(g000, Pf0); |
||||||
|
T n100 = dot(g100, detail::tvec3<T>(Pf1.x, Pf0.y, Pf0.z)); |
||||||
|
T n010 = dot(g010, detail::tvec3<T>(Pf0.x, Pf1.y, Pf0.z)); |
||||||
|
T n110 = dot(g110, detail::tvec3<T>(Pf1.x, Pf1.y, Pf0.z)); |
||||||
|
T n001 = dot(g001, detail::tvec3<T>(Pf0.x, Pf0.y, Pf1.z)); |
||||||
|
T n101 = dot(g101, detail::tvec3<T>(Pf1.x, Pf0.y, Pf1.z)); |
||||||
|
T n011 = dot(g011, detail::tvec3<T>(Pf0.x, Pf1.y, Pf1.z)); |
||||||
|
T n111 = dot(g111, Pf1); |
||||||
|
|
||||||
|
detail::tvec3<T> fade_xyz = fade(Pf0); |
||||||
|
detail::tvec4<T> n_z = mix(detail::tvec4<T>(n000, n100, n010, n110), detail::tvec4<T>(n001, n101, n011, n111), fade_xyz.z); |
||||||
|
detail::tvec2<T> n_yz = mix(detail::tvec2<T>(n_z.x, n_z.y), detail::tvec2<T>(n_z.z, n_z.w), fade_xyz.y); |
||||||
|
T n_xyz = mix(n_yz.x, n_yz.y, fade_xyz.x); |
||||||
|
return T(2.2) * n_xyz; |
||||||
|
} |
||||||
|
|
||||||
|
// Classic Perlin noise, periodic version |
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T perlin(detail::tvec4<T> const & P, detail::tvec4<T> const & rep) |
||||||
|
{ |
||||||
|
detail::tvec4<T> Pi0 = mod(floor(P), rep); // Integer part modulo rep |
||||||
|
detail::tvec4<T> Pi1 = mod(Pi0 + T(1), rep); // Integer part + 1 mod rep |
||||||
|
detail::tvec4<T> Pf0 = fract(P); // Fractional part for interpolation |
||||||
|
detail::tvec4<T> Pf1 = Pf0 - T(1); // Fractional part - 1.0 |
||||||
|
detail::tvec4<T> ix = detail::tvec4<T>(Pi0.x, Pi1.x, Pi0.x, Pi1.x); |
||||||
|
detail::tvec4<T> iy = detail::tvec4<T>(Pi0.y, Pi0.y, Pi1.y, Pi1.y); |
||||||
|
detail::tvec4<T> iz0(Pi0.z); |
||||||
|
detail::tvec4<T> iz1(Pi1.z); |
||||||
|
detail::tvec4<T> iw0(Pi0.w); |
||||||
|
detail::tvec4<T> iw1(Pi1.w); |
||||||
|
|
||||||
|
detail::tvec4<T> ixy = permute(permute(ix) + iy); |
||||||
|
detail::tvec4<T> ixy0 = permute(ixy + iz0); |
||||||
|
detail::tvec4<T> ixy1 = permute(ixy + iz1); |
||||||
|
detail::tvec4<T> ixy00 = permute(ixy0 + iw0); |
||||||
|
detail::tvec4<T> ixy01 = permute(ixy0 + iw1); |
||||||
|
detail::tvec4<T> ixy10 = permute(ixy1 + iw0); |
||||||
|
detail::tvec4<T> ixy11 = permute(ixy1 + iw1); |
||||||
|
|
||||||
|
detail::tvec4<T> gx00 = ixy00 / T(7); |
||||||
|
detail::tvec4<T> gy00 = floor(gx00) / T(7); |
||||||
|
detail::tvec4<T> gz00 = floor(gy00) / T(6); |
||||||
|
gx00 = fract(gx00) - T(0.5); |
||||||
|
gy00 = fract(gy00) - T(0.5); |
||||||
|
gz00 = fract(gz00) - T(0.5); |
||||||
|
detail::tvec4<T> gw00 = detail::tvec4<T>(0.75) - abs(gx00) - abs(gy00) - abs(gz00); |
||||||
|
detail::tvec4<T> sw00 = step(gw00, detail::tvec4<T>(0)); |
||||||
|
gx00 -= sw00 * (step(0.0, gx00) - T(0.5)); |
||||||
|
gy00 -= sw00 * (step(0.0, gy00) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx01 = ixy01 / T(7); |
||||||
|
detail::tvec4<T> gy01 = floor(gx01) / T(7); |
||||||
|
detail::tvec4<T> gz01 = floor(gy01) / T(6); |
||||||
|
gx01 = fract(gx01) - T(0.5); |
||||||
|
gy01 = fract(gy01) - T(0.5); |
||||||
|
gz01 = fract(gz01) - T(0.5); |
||||||
|
detail::tvec4<T> gw01 = detail::tvec4<T>(0.75) - abs(gx01) - abs(gy01) - abs(gz01); |
||||||
|
detail::tvec4<T> sw01 = step(gw01, detail::tvec4<T>(0.0)); |
||||||
|
gx01 -= sw01 * (step(0.0, gx01) - T(0.5)); |
||||||
|
gy01 -= sw01 * (step(0.0, gy01) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx10 = ixy10 / T(7); |
||||||
|
detail::tvec4<T> gy10 = floor(gx10) / T(7); |
||||||
|
detail::tvec4<T> gz10 = floor(gy10) / T(6); |
||||||
|
gx10 = fract(gx10) - T(0.5); |
||||||
|
gy10 = fract(gy10) - T(0.5); |
||||||
|
gz10 = fract(gz10) - T(0.5); |
||||||
|
detail::tvec4<T> gw10 = detail::tvec4<T>(0.75) - abs(gx10) - abs(gy10) - abs(gz10); |
||||||
|
detail::tvec4<T> sw10 = step(gw10, detail::tvec4<T>(0.0)); |
||||||
|
gx10 -= sw10 * (step(0.0, gx10) - T(0.5)); |
||||||
|
gy10 -= sw10 * (step(0.0, gy10) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> gx11 = ixy11 / T(7); |
||||||
|
detail::tvec4<T> gy11 = floor(gx11) / T(7); |
||||||
|
detail::tvec4<T> gz11 = floor(gy11) / T(6); |
||||||
|
gx11 = fract(gx11) - T(0.5); |
||||||
|
gy11 = fract(gy11) - T(0.5); |
||||||
|
gz11 = fract(gz11) - T(0.5); |
||||||
|
detail::tvec4<T> gw11 = detail::tvec4<T>(0.75) - abs(gx11) - abs(gy11) - abs(gz11); |
||||||
|
detail::tvec4<T> sw11 = step(gw11, detail::tvec4<T>(0.0)); |
||||||
|
gx11 -= sw11 * (step(0.0, gx11) - T(0.5)); |
||||||
|
gy11 -= sw11 * (step(0.0, gy11) - T(0.5)); |
||||||
|
|
||||||
|
detail::tvec4<T> g0000(gx00.x, gy00.x, gz00.x, gw00.x); |
||||||
|
detail::tvec4<T> g1000(gx00.y, gy00.y, gz00.y, gw00.y); |
||||||
|
detail::tvec4<T> g0100(gx00.z, gy00.z, gz00.z, gw00.z); |
||||||
|
detail::tvec4<T> g1100(gx00.w, gy00.w, gz00.w, gw00.w); |
||||||
|
detail::tvec4<T> g0010(gx10.x, gy10.x, gz10.x, gw10.x); |
||||||
|
detail::tvec4<T> g1010(gx10.y, gy10.y, gz10.y, gw10.y); |
||||||
|
detail::tvec4<T> g0110(gx10.z, gy10.z, gz10.z, gw10.z); |
||||||
|
detail::tvec4<T> g1110(gx10.w, gy10.w, gz10.w, gw10.w); |
||||||
|
detail::tvec4<T> g0001(gx01.x, gy01.x, gz01.x, gw01.x); |
||||||
|
detail::tvec4<T> g1001(gx01.y, gy01.y, gz01.y, gw01.y); |
||||||
|
detail::tvec4<T> g0101(gx01.z, gy01.z, gz01.z, gw01.z); |
||||||
|
detail::tvec4<T> g1101(gx01.w, gy01.w, gz01.w, gw01.w); |
||||||
|
detail::tvec4<T> g0011(gx11.x, gy11.x, gz11.x, gw11.x); |
||||||
|
detail::tvec4<T> g1011(gx11.y, gy11.y, gz11.y, gw11.y); |
||||||
|
detail::tvec4<T> g0111(gx11.z, gy11.z, gz11.z, gw11.z); |
||||||
|
detail::tvec4<T> g1111(gx11.w, gy11.w, gz11.w, gw11.w); |
||||||
|
|
||||||
|
detail::tvec4<T> norm00 = taylorInvSqrt(detail::tvec4<T>(dot(g0000, g0000), dot(g0100, g0100), dot(g1000, g1000), dot(g1100, g1100))); |
||||||
|
g0000 *= norm00.x; |
||||||
|
g0100 *= norm00.y; |
||||||
|
g1000 *= norm00.z; |
||||||
|
g1100 *= norm00.w; |
||||||
|
|
||||||
|
detail::tvec4<T> norm01 = taylorInvSqrt(detail::tvec4<T>(dot(g0001, g0001), dot(g0101, g0101), dot(g1001, g1001), dot(g1101, g1101))); |
||||||
|
g0001 *= norm01.x; |
||||||
|
g0101 *= norm01.y; |
||||||
|
g1001 *= norm01.z; |
||||||
|
g1101 *= norm01.w; |
||||||
|
|
||||||
|
detail::tvec4<T> norm10 = taylorInvSqrt(detail::tvec4<T>(dot(g0010, g0010), dot(g0110, g0110), dot(g1010, g1010), dot(g1110, g1110))); |
||||||
|
g0010 *= norm10.x; |
||||||
|
g0110 *= norm10.y; |
||||||
|
g1010 *= norm10.z; |
||||||
|
g1110 *= norm10.w; |
||||||
|
|
||||||
|
detail::tvec4<T> norm11 = taylorInvSqrt(detail::tvec4<T>(dot(g0011, g0011), dot(g0111, g0111), dot(g1011, g1011), dot(g1111, g1111))); |
||||||
|
g0011 *= norm11.x; |
||||||
|
g0111 *= norm11.y; |
||||||
|
g1011 *= norm11.z; |
||||||
|
g1111 *= norm11.w; |
||||||
|
|
||||||
|
T n0000 = dot(g0000, Pf0); |
||||||
|
T n1000 = dot(g1000, detail::tvec4<T>(Pf1.x, Pf0.y, Pf0.z, Pf0.w)); |
||||||
|
T n0100 = dot(g0100, detail::tvec4<T>(Pf0.x, Pf1.y, Pf0.z, Pf0.w)); |
||||||
|
T n1100 = dot(g1100, detail::tvec4<T>(Pf1.x, Pf1.y, Pf0.z, Pf0.w)); |
||||||
|
T n0010 = dot(g0010, detail::tvec4<T>(Pf0.x, Pf0.y, Pf1.z, Pf0.w)); |
||||||
|
T n1010 = dot(g1010, detail::tvec4<T>(Pf1.x, Pf0.y, Pf1.z, Pf0.w)); |
||||||
|
T n0110 = dot(g0110, detail::tvec4<T>(Pf0.x, Pf1.y, Pf1.z, Pf0.w)); |
||||||
|
T n1110 = dot(g1110, detail::tvec4<T>(Pf1.x, Pf1.y, Pf1.z, Pf0.w)); |
||||||
|
T n0001 = dot(g0001, detail::tvec4<T>(Pf0.x, Pf0.y, Pf0.z, Pf1.w)); |
||||||
|
T n1001 = dot(g1001, detail::tvec4<T>(Pf1.x, Pf0.y, Pf0.z, Pf1.w)); |
||||||
|
T n0101 = dot(g0101, detail::tvec4<T>(Pf0.x, Pf1.y, Pf0.z, Pf1.w)); |
||||||
|
T n1101 = dot(g1101, detail::tvec4<T>(Pf1.x, Pf1.y, Pf0.z, Pf1.w)); |
||||||
|
T n0011 = dot(g0011, detail::tvec4<T>(Pf0.x, Pf0.y, Pf1.z, Pf1.w)); |
||||||
|
T n1011 = dot(g1011, detail::tvec4<T>(Pf1.x, Pf0.y, Pf1.z, Pf1.w)); |
||||||
|
T n0111 = dot(g0111, detail::tvec4<T>(Pf0.x, Pf1.y, Pf1.z, Pf1.w)); |
||||||
|
T n1111 = dot(g1111, Pf1); |
||||||
|
|
||||||
|
detail::tvec4<T> fade_xyzw = fade(Pf0); |
||||||
|
detail::tvec4<T> n_0w = mix(detail::tvec4<T>(n0000, n1000, n0100, n1100), detail::tvec4<T>(n0001, n1001, n0101, n1101), fade_xyzw.w); |
||||||
|
detail::tvec4<T> n_1w = mix(detail::tvec4<T>(n0010, n1010, n0110, n1110), detail::tvec4<T>(n0011, n1011, n0111, n1111), fade_xyzw.w); |
||||||
|
detail::tvec4<T> n_zw = mix(n_0w, n_1w, fade_xyzw.z); |
||||||
|
detail::tvec2<T> n_yzw = mix(detail::tvec2<T>(n_zw.x, n_zw.y), detail::tvec2<T>(n_zw.z, n_zw.w), fade_xyzw.y); |
||||||
|
T n_xyzw = mix(n_yzw.x, n_yzw.y, fade_xyzw.x); |
||||||
|
return T(2.2) * n_xyzw; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T simplex(glm::detail::tvec2<T> const & v) |
||||||
|
{ |
||||||
|
detail::tvec4<T> const C = detail::tvec4<T>( |
||||||
|
T( 0.211324865405187), // (3.0 - sqrt(3.0)) / 6.0 |
||||||
|
T( 0.366025403784439), // 0.5 * (sqrt(3.0) - 1.0) |
||||||
|
T(-0.577350269189626), // -1.0 + 2.0 * C.x |
||||||
|
T( 0.024390243902439)); // 1.0 / 41.0 |
||||||
|
|
||||||
|
// First corner |
||||||
|
detail::tvec2<T> i = floor(v + dot(v, detail::tvec2<T>(C[1]))); |
||||||
|
detail::tvec2<T> x0 = v - i + dot(i, detail::tvec2<T>(C[0])); |
||||||
|
|
||||||
|
// Other corners |
||||||
|
//i1.x = step( x0.y, x0.x ); // x0.x > x0.y ? 1.0 : 0.0 |
||||||
|
//i1.y = 1.0 - i1.x; |
||||||
|
detail::tvec2<T> i1 = (x0.x > x0.y) ? detail::tvec2<T>(1, 0) : detail::tvec2<T>(0, 1); |
||||||
|
// x0 = x0 - 0.0 + 0.0 * C.xx ; |
||||||
|
// x1 = x0 - i1 + 1.0 * C.xx ; |
||||||
|
// x2 = x0 - 1.0 + 2.0 * C.xx ; |
||||||
|
detail::tvec4<T> x12 = detail::tvec4<T>(x0.x, x0.y, x0.x, x0.y) + detail::tvec4<T>(C.x, C.x, C.z, C.z); |
||||||
|
x12 = detail::tvec4<T>(detail::tvec2<T>(x12) - i1, x12.z, x12.w); |
||||||
|
|
||||||
|
// Permutations |
||||||
|
i = mod(i, T(289)); // Avoid truncation effects in permutation |
||||||
|
detail::tvec3<T> p = permute( |
||||||
|
permute(i.y + detail::tvec3<T>(T(0), i1.y, T(1))) |
||||||
|
+ i.x + detail::tvec3<T>(T(0), i1.x, T(1))); |
||||||
|
|
||||||
|
detail::tvec3<T> m = max(T(0.5) - detail::tvec3<T>( |
||||||
|
dot(x0, x0), |
||||||
|
dot(detail::tvec2<T>(x12.x, x12.y), detail::tvec2<T>(x12.x, x12.y)), |
||||||
|
dot(detail::tvec2<T>(x12.z, x12.w), detail::tvec2<T>(x12.z, x12.w))), T(0)); |
||||||
|
m = m * m ; |
||||||
|
m = m * m ; |
||||||
|
|
||||||
|
// Gradients: 41 points uniformly over a line, mapped onto a diamond. |
||||||
|
// The ring size 17*17 = 289 is close to a multiple of 41 (41*7 = 287) |
||||||
|
|
||||||
|
detail::tvec3<T> x = T(2) * fract(p * C.w) - T(1); |
||||||
|
detail::tvec3<T> h = abs(x) - T(0.5); |
||||||
|
detail::tvec3<T> ox = floor(x + T(0.5)); |
||||||
|
detail::tvec3<T> a0 = x - ox; |
||||||
|
|
||||||
|
// Normalise gradients implicitly by scaling m |
||||||
|
// Inlined for speed: m *= taylorInvSqrt( a0*a0 + h*h ); |
||||||
|
m *= T(1.79284291400159) - T(0.85373472095314) * (a0 * a0 + h * h); |
||||||
|
|
||||||
|
// Compute final noise value at P |
||||||
|
detail::tvec3<T> g; |
||||||
|
g.x = a0.x * x0.x + h.x * x0.y; |
||||||
|
//g.yz = a0.yz * x12.xz + h.yz * x12.yw; |
||||||
|
g.y = a0.y * x12.x + h.y * x12.y; |
||||||
|
g.z = a0.z * x12.z + h.z * x12.w; |
||||||
|
return T(130) * dot(m, g); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T simplex(detail::tvec3<T> const & v) |
||||||
|
{ |
||||||
|
detail::tvec2<T> const C(1.0 / 6.0, 1.0 / 3.0); |
||||||
|
detail::tvec4<T> const D(0.0, 0.5, 1.0, 2.0); |
||||||
|
|
||||||
|
// First corner |
||||||
|
detail::tvec3<T> i(floor(v + dot(v, detail::tvec3<T>(C.y)))); |
||||||
|
detail::tvec3<T> x0(v - i + dot(i, detail::tvec3<T>(C.x))); |
||||||
|
|
||||||
|
// Other corners |
||||||
|
detail::tvec3<T> g(step(detail::tvec3<T>(x0.y, x0.z, x0.x), x0)); |
||||||
|
detail::tvec3<T> l(T(1) - g); |
||||||
|
detail::tvec3<T> i1(min(g, detail::tvec3<T>(l.z, l.x, l.y))); |
||||||
|
detail::tvec3<T> i2(max(g, detail::tvec3<T>(l.z, l.x, l.y))); |
||||||
|
|
||||||
|
// x0 = x0 - 0.0 + 0.0 * C.xxx; |
||||||
|
// x1 = x0 - i1 + 1.0 * C.xxx; |
||||||
|
// x2 = x0 - i2 + 2.0 * C.xxx; |
||||||
|
// x3 = x0 - 1.0 + 3.0 * C.xxx; |
||||||
|
detail::tvec3<T> x1(x0 - i1 + C.x); |
||||||
|
detail::tvec3<T> x2(x0 - i2 + C.y); // 2.0*C.x = 1/3 = C.y |
||||||
|
detail::tvec3<T> x3(x0 - D.y); // -1.0+3.0*C.x = -0.5 = -D.y |
||||||
|
|
||||||
|
// Permutations |
||||||
|
i = mod289(i); |
||||||
|
detail::tvec4<T> p(permute(permute(permute( |
||||||
|
i.z + detail::tvec4<T>(T(0), i1.z, i2.z, T(1))) + |
||||||
|
i.y + detail::tvec4<T>(T(0), i1.y, i2.y, T(1))) + |
||||||
|
i.x + detail::tvec4<T>(T(0), i1.x, i2.x, T(1)))); |
||||||
|
|
||||||
|
// Gradients: 7x7 points over a square, mapped onto an octahedron. |
||||||
|
// The ring size 17*17 = 289 is close to a multiple of 49 (49*6 = 294) |
||||||
|
T n_ = T(0.142857142857); // 1.0/7.0 |
||||||
|
detail::tvec3<T> ns(n_ * detail::tvec3<T>(D.w, D.y, D.z) - detail::tvec3<T>(D.x, D.z, D.x)); |
||||||
|
|
||||||
|
detail::tvec4<T> j(p - T(49) * floor(p * ns.z * ns.z)); // mod(p,7*7) |
||||||
|
|
||||||
|
detail::tvec4<T> x_(floor(j * ns.z)); |
||||||
|
detail::tvec4<T> y_(floor(j - T(7) * x_)); // mod(j,N) |
||||||
|
|
||||||
|
detail::tvec4<T> x(x_ * ns.x + ns.y); |
||||||
|
detail::tvec4<T> y(y_ * ns.x + ns.y); |
||||||
|
detail::tvec4<T> h(T(1) - abs(x) - abs(y)); |
||||||
|
|
||||||
|
detail::tvec4<T> b0(x.x, x.y, y.x, y.y); |
||||||
|
detail::tvec4<T> b1(x.z, x.w, y.z, y.w); |
||||||
|
|
||||||
|
// vec4 s0 = vec4(lessThan(b0,0.0))*2.0 - 1.0; |
||||||
|
// vec4 s1 = vec4(lessThan(b1,0.0))*2.0 - 1.0; |
||||||
|
detail::tvec4<T> s0(floor(b0) * T(2) + T(1)); |
||||||
|
detail::tvec4<T> s1(floor(b1) * T(2) + T(1)); |
||||||
|
detail::tvec4<T> sh(-step(h, detail::tvec4<T>(0.0))); |
||||||
|
|
||||||
|
detail::tvec4<T> a0 = detail::tvec4<T>(b0.x, b0.z, b0.y, b0.w) + detail::tvec4<T>(s0.x, s0.z, s0.y, s0.w) * detail::tvec4<T>(sh.x, sh.x, sh.y, sh.y); |
||||||
|
detail::tvec4<T> a1 = detail::tvec4<T>(b1.x, b1.z, b1.y, b1.w) + detail::tvec4<T>(s1.x, s1.z, s1.y, s1.w) * detail::tvec4<T>(sh.z, sh.z, sh.w, sh.w); |
||||||
|
|
||||||
|
detail::tvec3<T> p0(a0.x, a0.y, h.x); |
||||||
|
detail::tvec3<T> p1(a0.z, a0.w, h.y); |
||||||
|
detail::tvec3<T> p2(a1.x, a1.y, h.z); |
||||||
|
detail::tvec3<T> p3(a1.z, a1.w, h.w); |
||||||
|
|
||||||
|
// Normalise gradients |
||||||
|
detail::tvec4<T> norm = taylorInvSqrt(detail::tvec4<T>(dot(p0, p0), dot(p1, p1), dot(p2, p2), dot(p3, p3))); |
||||||
|
p0 *= norm.x; |
||||||
|
p1 *= norm.y; |
||||||
|
p2 *= norm.z; |
||||||
|
p3 *= norm.w; |
||||||
|
|
||||||
|
// Mix final noise value |
||||||
|
detail::tvec4<T> m = max(T(0.6) - detail::tvec4<T>(dot(x0, x0), dot(x1, x1), dot(x2, x2), dot(x3, x3)), T(0)); |
||||||
|
m = m * m; |
||||||
|
return T(42) * dot(m * m, detail::tvec4<T>(dot(p0, x0), dot(p1, x1), dot(p2, x2), dot(p3, x3))); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T simplex(detail::tvec4<T> const & v) |
||||||
|
{ |
||||||
|
detail::tvec4<T> const C( |
||||||
|
0.138196601125011, // (5 - sqrt(5))/20 G4 |
||||||
|
0.276393202250021, // 2 * G4 |
||||||
|
0.414589803375032, // 3 * G4 |
||||||
|
-0.447213595499958); // -1 + 4 * G4 |
||||||
|
|
||||||
|
// (sqrt(5) - 1)/4 = F4, used once below |
||||||
|
T const F4 = T(0.309016994374947451); |
||||||
|
|
||||||
|
// First corner |
||||||
|
detail::tvec4<T> i = floor(v + dot(v, vec4(F4))); |
||||||
|
detail::tvec4<T> x0 = v - i + dot(i, vec4(C.x)); |
||||||
|
|
||||||
|
// Other corners |
||||||
|
|
||||||
|
// Rank sorting originally contributed by Bill Licea-Kane, AMD (formerly ATI) |
||||||
|
detail::tvec4<T> i0; |
||||||
|
detail::tvec3<T> isX = step(detail::tvec3<T>(x0.y, x0.z, x0.w), detail::tvec3<T>(x0.x)); |
||||||
|
detail::tvec3<T> isYZ = step(detail::tvec3<T>(x0.z, x0.w, x0.w), detail::tvec3<T>(x0.y, x0.y, x0.z)); |
||||||
|
// i0.x = dot(isX, vec3(1.0)); |
||||||
|
//i0.x = isX.x + isX.y + isX.z; |
||||||
|
//i0.yzw = T(1) - isX; |
||||||
|
i0 = detail::tvec4<T>(isX.x + isX.y + isX.z, T(1) - isX); |
||||||
|
// i0.y += dot(isYZ.xy, vec2(1.0)); |
||||||
|
i0.y += isYZ.x + isYZ.y; |
||||||
|
//i0.zw += 1.0 - detail::tvec2<T>(isYZ.x, isYZ.y); |
||||||
|
i0.z += T(1) - isYZ.x; |
||||||
|
i0.w += T(1) - isYZ.y; |
||||||
|
i0.z += isYZ.z; |
||||||
|
i0.w += T(1) - isYZ.z; |
||||||
|
|
||||||
|
// i0 now contains the unique values 0,1,2,3 in each channel |
||||||
|
detail::tvec4<T> i3 = clamp(i0, 0.0, 1.0); |
||||||
|
detail::tvec4<T> i2 = clamp(i0 - 1.0, 0.0, 1.0); |
||||||
|
detail::tvec4<T> i1 = clamp(i0 - 2.0, 0.0, 1.0); |
||||||
|
|
||||||
|
// x0 = x0 - 0.0 + 0.0 * C.xxxx |
||||||
|
// x1 = x0 - i1 + 0.0 * C.xxxx |
||||||
|
// x2 = x0 - i2 + 0.0 * C.xxxx |
||||||
|
// x3 = x0 - i3 + 0.0 * C.xxxx |
||||||
|
// x4 = x0 - 1.0 + 4.0 * C.xxxx |
||||||
|
detail::tvec4<T> x1 = x0 - i1 + C.x; |
||||||
|
detail::tvec4<T> x2 = x0 - i2 + C.y; |
||||||
|
detail::tvec4<T> x3 = x0 - i3 + C.z; |
||||||
|
detail::tvec4<T> x4 = x0 + C.w; |
||||||
|
|
||||||
|
// Permutations |
||||||
|
i = mod(i, T(289)); |
||||||
|
T j0 = permute(permute(permute(permute(i.w) + i.z) + i.y) + i.x); |
||||||
|
detail::tvec4<T> j1 = permute(permute(permute(permute( |
||||||
|
i.w + detail::tvec4<T>(i1.w, i2.w, i3.w, T(1))) |
||||||
|
+ i.z + detail::tvec4<T>(i1.z, i2.z, i3.z, T(1))) |
||||||
|
+ i.y + detail::tvec4<T>(i1.y, i2.y, i3.y, T(1))) |
||||||
|
+ i.x + detail::tvec4<T>(i1.x, i2.x, i3.x, T(1))); |
||||||
|
|
||||||
|
// Gradients: 7x7x6 points over a cube, mapped onto a 4-cross polytope |
||||||
|
// 7*7*6 = 294, which is close to the ring size 17*17 = 289. |
||||||
|
detail::tvec4<T> ip = detail::tvec4<T>(T(1) / T(294), T(1) / T(49), T(1) / T(7), T(0)); |
||||||
|
|
||||||
|
detail::tvec4<T> p0 = grad4(j0, ip); |
||||||
|
detail::tvec4<T> p1 = grad4(j1.x, ip); |
||||||
|
detail::tvec4<T> p2 = grad4(j1.y, ip); |
||||||
|
detail::tvec4<T> p3 = grad4(j1.z, ip); |
||||||
|
detail::tvec4<T> p4 = grad4(j1.w, ip); |
||||||
|
|
||||||
|
// Normalise gradients |
||||||
|
detail::tvec4<T> norm = taylorInvSqrt(detail::tvec4<T>(dot(p0, p0), dot(p1, p1), dot(p2, p2), dot(p3, p3))); |
||||||
|
p0 *= norm.x; |
||||||
|
p1 *= norm.y; |
||||||
|
p2 *= norm.z; |
||||||
|
p3 *= norm.w; |
||||||
|
p4 *= taylorInvSqrt(dot(p4, p4)); |
||||||
|
|
||||||
|
// Mix contributions from the five corners |
||||||
|
detail::tvec3<T> m0 = max(T(0.6) - detail::tvec3<T>(dot(x0, x0), dot(x1, x1), dot(x2, x2)), T(0)); |
||||||
|
detail::tvec2<T> m1 = max(T(0.6) - detail::tvec2<T>(dot(x3, x3), dot(x4, x4) ), T(0)); |
||||||
|
m0 = m0 * m0; |
||||||
|
m1 = m1 * m1; |
||||||
|
return T(49) * |
||||||
|
(dot(m0 * m0, detail::tvec3<T>(dot(p0, x0), dot(p1, x1), dot(p2, x2))) + |
||||||
|
dot(m1 * m1, detail::tvec2<T>(dot(p3, x3), dot(p4, x4)))); |
||||||
|
} |
||||||
|
|
||||||
|
}//namespace glm |
||||||
@ -0,0 +1,177 @@ |
|||||||
|
///////////////////////////////////////////////////////////////////////////////////
|
||||||
|
/// OpenGL Mathematics (glm.g-truc.net)
|
||||||
|
///
|
||||||
|
/// Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net)
|
||||||
|
/// Permission is hereby granted, free of charge, to any person obtaining a copy
|
||||||
|
/// of this software and associated documentation files (the "Software"), to deal
|
||||||
|
/// in the Software without restriction, including without limitation the rights
|
||||||
|
/// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
||||||
|
/// copies of the Software, and to permit persons to whom the Software is
|
||||||
|
/// furnished to do so, subject to the following conditions:
|
||||||
|
///
|
||||||
|
/// The above copyright notice and this permission notice shall be included in
|
||||||
|
/// all copies or substantial portions of the Software.
|
||||||
|
///
|
||||||
|
/// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
||||||
|
/// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
||||||
|
/// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
||||||
|
/// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
||||||
|
/// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
||||||
|
/// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
|
||||||
|
/// THE SOFTWARE.
|
||||||
|
///
|
||||||
|
/// @ref constants
|
||||||
|
/// @file glm/gtx/constants.hpp
|
||||||
|
/// @date 2011-09-30 / 2011-09-30
|
||||||
|
/// @author Christophe Riccio
|
||||||
|
///
|
||||||
|
/// @see core (dependence)
|
||||||
|
/// @see gtc_half_float (dependence)
|
||||||
|
///
|
||||||
|
/// @defgroup gtx_constants GLM_GTX_constants: Provide build-in constants
|
||||||
|
/// @ingroup gtx
|
||||||
|
///
|
||||||
|
/// @brief Allow to perform bit operations on integer values
|
||||||
|
///
|
||||||
|
/// <glm/gtx/constants.hpp> need to be included to use these functionalities.
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////
|
||||||
|
|
||||||
|
#ifndef GLM_GTX_constants |
||||||
|
#define GLM_GTX_constants GLM_VERSION |
||||||
|
|
||||||
|
// Dependency:
|
||||||
|
#include "../glm.hpp" |
||||||
|
#include "../gtc/half_float.hpp" |
||||||
|
|
||||||
|
#if(defined(GLM_MESSAGES) && !defined(glm_ext)) |
||||||
|
# pragma message("GLM: GLM_GTX_constants extension included") |
||||||
|
#endif |
||||||
|
|
||||||
|
namespace glm |
||||||
|
{ |
||||||
|
/// @addtogroup gtx_constants
|
||||||
|
/// @{
|
||||||
|
|
||||||
|
template <typename T> |
||||||
|
T pi(); |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T pi() |
||||||
|
{ |
||||||
|
return T(3.14159265358979323846264338327950288); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T root_pi() |
||||||
|
{ |
||||||
|
return T(1.772453850905516027); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T half_pi() |
||||||
|
{ |
||||||
|
return T(1.57079632679489661923132169163975144); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T quarter_pi() |
||||||
|
{ |
||||||
|
return T(0.785398163397448309615660845819875721); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T one_over_pi() |
||||||
|
{ |
||||||
|
return T(0.318309886183790671537767526745028724); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T two_over_pi() |
||||||
|
{ |
||||||
|
return T(0.636619772367581343075535053490057448); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T two_over_root_pi() |
||||||
|
{ |
||||||
|
return T(1.12837916709551257389615890312154517); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T one_over_root_two() |
||||||
|
{ |
||||||
|
return T(0.707106781186547524400844362104849039); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T root_half_pi() |
||||||
|
{ |
||||||
|
return T(1.253314137315500251); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T root_two_pi() |
||||||
|
{ |
||||||
|
return T(2.506628274631000502); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T root_ln_four() |
||||||
|
{ |
||||||
|
return T(1.17741002251547469); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T e() |
||||||
|
{ |
||||||
|
return T(2.71828182845904523536); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T euler() |
||||||
|
{ |
||||||
|
return T(0.577215664901532860606); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T root_two() |
||||||
|
{ |
||||||
|
return T(1.41421356237309504880168872420969808); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T ln_two() |
||||||
|
{ |
||||||
|
return T(0.693147180559945309417232121458176568); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T ln_ten(2.30258509299404568401799145468436421) |
||||||
|
{ |
||||||
|
return T(); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T ln_ln_two() |
||||||
|
{ |
||||||
|
return T(-0.3665129205816643); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T third() |
||||||
|
{ |
||||||
|
return T(0.333333333333333333); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename T> |
||||||
|
GLM_FUNC_QUALIFIER T twothirds() |
||||||
|
{ |
||||||
|
return T(0.666666666666666666); |
||||||
|
} |
||||||
|
|
||||||
|
/// @}
|
||||||
|
} //namespace glm
|
||||||
|
|
||||||
|
#include "constants.inl" |
||||||
|
|
||||||
|
#endif//GLM_GTX_constants
|
||||||
@ -0,0 +1,792 @@ |
|||||||
|
/////////////////////////////////////////////////////////////////////////////////// |
||||||
|
/// OpenGL Mathematics (glm.g-truc.net) |
||||||
|
/// |
||||||
|
/// Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net) |
||||||
|
/// Permission is hereby granted, free of charge, to any person obtaining a copy |
||||||
|
/// of this software and associated documentation files (the "Software"), to deal |
||||||
|
/// in the Software without restriction, including without limitation the rights |
||||||
|
/// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell |
||||||
|
/// copies of the Software, and to permit persons to whom the Software is |
||||||
|
/// furnished to do so, subject to the following conditions: |
||||||
|
/// |
||||||
|
/// The above copyright notice and this permission notice shall be included in |
||||||
|
/// all copies or substantial portions of the Software. |
||||||
|
/// |
||||||
|
/// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR |
||||||
|
/// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, |
||||||
|
/// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE |
||||||
|
/// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER |
||||||
|
/// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, |
||||||
|
/// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN |
||||||
|
/// THE SOFTWARE. |
||||||
|
/// |
||||||
|
/// @ref gtc_half_float |
||||||
|
/// @file glm/gtc/half_float.inl |
||||||
|
/// @date 2009-04-29 / 2011-06-05 |
||||||
|
/// @author Christophe Riccio |
||||||
|
/////////////////////////////////////////////////////////////////////////////////// |
||||||
|
|
||||||
|
#include "../core/_detail.hpp" |
||||||
|
|
||||||
|
namespace glm{ |
||||||
|
|
||||||
|
template <typename genIType> |
||||||
|
GLM_FUNC_QUALIFIER genIType mask |
||||||
|
( |
||||||
|
genIType const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return ((genIType(1) << (count)) - genIType(1)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valIType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<valIType> mask |
||||||
|
( |
||||||
|
detail::tvec2<valIType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<valIType>( |
||||||
|
mask(count[0]), |
||||||
|
mask(count[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valIType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<valIType> mask |
||||||
|
( |
||||||
|
detail::tvec3<valIType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<valIType>( |
||||||
|
mask(count[0]), |
||||||
|
mask(count[1]), |
||||||
|
mask(count[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valIType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<valIType> mask |
||||||
|
( |
||||||
|
detail::tvec4<valIType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<valIType>( |
||||||
|
mask(count[0]), |
||||||
|
mask(count[1]), |
||||||
|
mask(count[2]), |
||||||
|
mask(count[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// extractField |
||||||
|
template <typename genIType> |
||||||
|
GLM_FUNC_QUALIFIER genIType extractField |
||||||
|
( |
||||||
|
half const & value, |
||||||
|
genIType const & first, |
||||||
|
genIType const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
assert(first + count < sizeof(half)); |
||||||
|
return (value._data() << first) >> ((sizeof(half) << 3) - count); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIType> |
||||||
|
GLM_FUNC_QUALIFIER genIType extractField |
||||||
|
( |
||||||
|
float const & value, |
||||||
|
genIType const & first, |
||||||
|
genIType const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
assert(first + count < sizeof(float)); |
||||||
|
return (detail::uif32(value).i << first) >> ((sizeof(float) << 3) - count); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIType> |
||||||
|
GLM_FUNC_QUALIFIER genIType extractField |
||||||
|
( |
||||||
|
double const & value, |
||||||
|
genIType const & first, |
||||||
|
genIType const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
assert(first + count < sizeof(double)); |
||||||
|
return (detail::uif64(value).i << first) >> ((sizeof(double) << genIType(3)) - count); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER genIUType extractField |
||||||
|
( |
||||||
|
genIUType const & Value, |
||||||
|
sizeType const & First, |
||||||
|
sizeType const & Count |
||||||
|
) |
||||||
|
{ |
||||||
|
sizeType GenSize = sizeof(genIUType) << 3; |
||||||
|
|
||||||
|
assert(First + Count <= GenSize); |
||||||
|
|
||||||
|
genIUType ShiftLeft = Count ? Value << (GenSize - (Count + First)) : 0; |
||||||
|
genIUType ShiftBack = ShiftLeft >> genIUType(GenSize - Count); |
||||||
|
|
||||||
|
return ShiftBack; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<genIUType> extractField |
||||||
|
( |
||||||
|
detail::tvec2<genIUType> const & value, |
||||||
|
sizeType const & first, |
||||||
|
sizeType const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<genIUType>( |
||||||
|
extractField(value[0], first, count), |
||||||
|
extractField(value[1], first, count)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<genIUType> extractField |
||||||
|
( |
||||||
|
detail::tvec3<genIUType> const & value, |
||||||
|
sizeType const & first, |
||||||
|
sizeType const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<genIUType>( |
||||||
|
extractField(value[0], first, count), |
||||||
|
extractField(value[1], first, count), |
||||||
|
extractField(value[2], first, count)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<genIUType> extractField |
||||||
|
( |
||||||
|
detail::tvec4<genIUType> const & value, |
||||||
|
sizeType const & first, |
||||||
|
sizeType const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<genIUType>( |
||||||
|
extractField(value[0], first, count), |
||||||
|
extractField(value[1], first, count), |
||||||
|
extractField(value[2], first, count), |
||||||
|
extractField(value[3], first, count)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<genIUType> extractField |
||||||
|
( |
||||||
|
detail::tvec2<genIUType> const & value, |
||||||
|
detail::tvec2<sizeType> const & first, |
||||||
|
detail::tvec2<sizeType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<genIUType>( |
||||||
|
extractField(value[0], first[0], count[0]), |
||||||
|
extractField(value[1], first[1], count[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<genIUType> extractField |
||||||
|
( |
||||||
|
detail::tvec3<genIUType> const & value, |
||||||
|
detail::tvec3<sizeType> const & first, |
||||||
|
detail::tvec3<sizeType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<genIUType>( |
||||||
|
extractField(value[0], first[0], count[0]), |
||||||
|
extractField(value[1], first[1], count[1]), |
||||||
|
extractField(value[2], first[2], count[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<genIUType> extractField |
||||||
|
( |
||||||
|
detail::tvec4<genIUType> const & value, |
||||||
|
detail::tvec4<sizeType> const & first, |
||||||
|
detail::tvec4<sizeType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<genIUType>( |
||||||
|
extractField(value[0], first[0], count[0]), |
||||||
|
extractField(value[1], first[1], count[1]), |
||||||
|
extractField(value[2], first[2], count[2]), |
||||||
|
extractField(value[3], first[3], count[3])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<genIUType> extractField |
||||||
|
( |
||||||
|
genIUType const & value, |
||||||
|
detail::tvec2<sizeType> const & first, |
||||||
|
detail::tvec2<sizeType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<genIUType>( |
||||||
|
extractField(value, first[0], count[0]), |
||||||
|
extractField(value, first[1], count[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<genIUType> extractField |
||||||
|
( |
||||||
|
genIUType const & value, |
||||||
|
detail::tvec3<sizeType> const & first, |
||||||
|
detail::tvec3<sizeType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<genIUType>( |
||||||
|
extractField(value, first[0], count[0]), |
||||||
|
extractField(value, first[1], count[1]), |
||||||
|
extractField(value, first[2], count[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType, typename sizeType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<genIUType> extractField |
||||||
|
( |
||||||
|
genIUType const & value, |
||||||
|
detail::tvec4<sizeType> const & first, |
||||||
|
detail::tvec4<sizeType> const & count |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<genIUType>( |
||||||
|
extractField(value, first[0], count[0]), |
||||||
|
extractField(value, first[1], count[1]), |
||||||
|
extractField(value, first[2], count[2]), |
||||||
|
extractField(value, first[3], count[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// lowestBit |
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER int lowestBit |
||||||
|
( |
||||||
|
genType const & Value |
||||||
|
) |
||||||
|
{ |
||||||
|
assert(Value != genType(0)); // not valid call |
||||||
|
|
||||||
|
genType Bit; |
||||||
|
for(Bit = genType(0); !(Value & (1 << Bit)); ++Bit){} |
||||||
|
return Bit; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<int> lowestBit |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<int>( |
||||||
|
lowestBit(value[0]), |
||||||
|
lowestBit(value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<int> lowestBit |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<int>( |
||||||
|
lowestBit(value[0]), |
||||||
|
lowestBit(value[1]), |
||||||
|
lowestBit(value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<int> lowestBit |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<int>( |
||||||
|
lowestBit(value[0]), |
||||||
|
lowestBit(value[1]), |
||||||
|
lowestBit(value[2]), |
||||||
|
lowestBit(value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// highestBit |
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER int highestBit |
||||||
|
( |
||||||
|
genType const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
assert(value != genType(0)); // not valid call |
||||||
|
|
||||||
|
genType bit = genType(-1); |
||||||
|
for(genType tmp = value; tmp; tmp >>= 1, ++bit){} |
||||||
|
return bit; |
||||||
|
} |
||||||
|
|
||||||
|
//template <> |
||||||
|
//GLM_FUNC_QUALIFIER int highestBit<int> |
||||||
|
//( |
||||||
|
// int value |
||||||
|
//) |
||||||
|
//{ |
||||||
|
// int bit = -1; |
||||||
|
// for(int tmp = value; tmp; tmp >>= 1, ++bit); |
||||||
|
// return bit; |
||||||
|
//} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<int> highestBit |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<int>( |
||||||
|
highestBit(value[0]), |
||||||
|
highestBit(value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<int> highestBit |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<int>( |
||||||
|
highestBit(value[0]), |
||||||
|
highestBit(value[1]), |
||||||
|
highestBit(value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<int> highestBit |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<int>( |
||||||
|
highestBit(value[0]), |
||||||
|
highestBit(value[1]), |
||||||
|
highestBit(value[2]), |
||||||
|
highestBit(value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// highestBitValue |
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER genType highestBitValue |
||||||
|
( |
||||||
|
genType const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
genType tmp = value; |
||||||
|
genType result = genType(0); |
||||||
|
while(tmp) |
||||||
|
{ |
||||||
|
result = (tmp & (~tmp + 1)); // grab lowest bit |
||||||
|
tmp &= ~result; // clear lowest bit |
||||||
|
} |
||||||
|
return result; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<int> highestBitValue |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<int>( |
||||||
|
highestBitValue(value[0]), |
||||||
|
highestBitValue(value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<int> highestBitValue |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<int>( |
||||||
|
highestBitValue(value[0]), |
||||||
|
highestBitValue(value[1]), |
||||||
|
highestBitValue(value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<int> highestBitValue |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<int>( |
||||||
|
highestBitValue(value[0]), |
||||||
|
highestBitValue(value[1]), |
||||||
|
highestBitValue(value[2]), |
||||||
|
highestBitValue(value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// isPowerOfTwo |
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER bool isPowerOfTwo(genType const & Value) |
||||||
|
{ |
||||||
|
//detail::If<std::numeric_limits<genType>::is_signed>::apply(abs, Value); |
||||||
|
//return !(Value & (Value - 1)); |
||||||
|
|
||||||
|
// For old complier? |
||||||
|
genType Result = Value; |
||||||
|
if(std::numeric_limits<genType>::is_signed) |
||||||
|
Result = abs(Result); |
||||||
|
return !(Result & (Result - 1)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<bool> isPowerOfTwo |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<bool>( |
||||||
|
isPowerOfTwo(value[0]), |
||||||
|
isPowerOfTwo(value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<bool> isPowerOfTwo |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<bool>( |
||||||
|
isPowerOfTwo(value[0]), |
||||||
|
isPowerOfTwo(value[1]), |
||||||
|
isPowerOfTwo(value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<bool> isPowerOfTwo |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<bool>( |
||||||
|
isPowerOfTwo(value[0]), |
||||||
|
isPowerOfTwo(value[1]), |
||||||
|
isPowerOfTwo(value[2]), |
||||||
|
isPowerOfTwo(value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// powerOfTwoAbove |
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER genType powerOfTwoAbove(genType const & value) |
||||||
|
{ |
||||||
|
return isPowerOfTwo(value) ? value : highestBitValue(value) << 1; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<valType> powerOfTwoAbove |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<valType>( |
||||||
|
powerOfTwoAbove(value[0]), |
||||||
|
powerOfTwoAbove(value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<valType> powerOfTwoAbove |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<valType>( |
||||||
|
powerOfTwoAbove(value[0]), |
||||||
|
powerOfTwoAbove(value[1]), |
||||||
|
powerOfTwoAbove(value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<valType> powerOfTwoAbove |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<valType>( |
||||||
|
powerOfTwoAbove(value[0]), |
||||||
|
powerOfTwoAbove(value[1]), |
||||||
|
powerOfTwoAbove(value[2]), |
||||||
|
powerOfTwoAbove(value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// powerOfTwoBelow |
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER genType powerOfTwoBelow |
||||||
|
( |
||||||
|
genType const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return isPowerOfTwo(value) ? value : highestBitValue(value); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<valType> powerOfTwoBelow |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<valType>( |
||||||
|
powerOfTwoBelow(value[0]), |
||||||
|
powerOfTwoBelow(value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<valType> powerOfTwoBelow |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<valType>( |
||||||
|
powerOfTwoBelow(value[0]), |
||||||
|
powerOfTwoBelow(value[1]), |
||||||
|
powerOfTwoBelow(value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<valType> powerOfTwoBelow |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<valType>( |
||||||
|
powerOfTwoBelow(value[0]), |
||||||
|
powerOfTwoBelow(value[1]), |
||||||
|
powerOfTwoBelow(value[2]), |
||||||
|
powerOfTwoBelow(value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
// powerOfTwoNearest |
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER genType powerOfTwoNearest |
||||||
|
( |
||||||
|
genType const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
if(isPowerOfTwo(value)) |
||||||
|
return value; |
||||||
|
|
||||||
|
genType prev = highestBitValue(value); |
||||||
|
genType next = prev << 1; |
||||||
|
return (next - value) < (value - prev) ? next : prev; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<valType> powerOfTwoNearest |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<valType>( |
||||||
|
powerOfTwoNearest(value[0]), |
||||||
|
powerOfTwoNearest(value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<valType> powerOfTwoNearest |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<valType>( |
||||||
|
powerOfTwoNearest(value[0]), |
||||||
|
powerOfTwoNearest(value[1]), |
||||||
|
powerOfTwoNearest(value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<valType> powerOfTwoNearest |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<valType>( |
||||||
|
powerOfTwoNearest(value[0]), |
||||||
|
powerOfTwoNearest(value[1]), |
||||||
|
powerOfTwoNearest(value[2]), |
||||||
|
powerOfTwoNearest(value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER genType bitRevert(genType const & In) |
||||||
|
{ |
||||||
|
GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_integer, "'bitRevert' only accept integer values"); |
||||||
|
|
||||||
|
genType Out = 0; |
||||||
|
std::size_t BitSize = sizeof(genType) * 8; |
||||||
|
for(std::size_t i = 0; i < BitSize; ++i) |
||||||
|
if(In & (genType(1) << i)) |
||||||
|
Out |= genType(1) << (BitSize - 1 - i); |
||||||
|
return Out; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<valType> bitRevert |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & Value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<valType>( |
||||||
|
bitRevert(Value[0]), |
||||||
|
bitRevert(Value[1])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<valType> bitRevert |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & Value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<valType>( |
||||||
|
bitRevert(Value[0]), |
||||||
|
bitRevert(Value[1]), |
||||||
|
bitRevert(Value[2])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<valType> bitRevert |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & Value |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<valType>( |
||||||
|
bitRevert(Value[0]), |
||||||
|
bitRevert(Value[1]), |
||||||
|
bitRevert(Value[2]), |
||||||
|
bitRevert(Value[3])); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER genType bitRotateRight(genType const & In, std::size_t Shift) |
||||||
|
{ |
||||||
|
GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_integer, "'bitRotateRight' only accept integer values"); |
||||||
|
|
||||||
|
std::size_t BitSize = sizeof(genType) * 8; |
||||||
|
return (In << Shift) | (In >> (BitSize - Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<valType> bitRotateRight |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & Value, |
||||||
|
std::size_t Shift |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<valType>( |
||||||
|
bitRotateRight(Value[0], Shift), |
||||||
|
bitRotateRight(Value[1], Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<valType> bitRotateRight |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & Value, |
||||||
|
std::size_t Shift |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<valType>( |
||||||
|
bitRotateRight(Value[0], Shift), |
||||||
|
bitRotateRight(Value[1], Shift), |
||||||
|
bitRotateRight(Value[2], Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<valType> bitRotateRight |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & Value, |
||||||
|
std::size_t Shift |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<valType>( |
||||||
|
bitRotateRight(Value[0], Shift), |
||||||
|
bitRotateRight(Value[1], Shift), |
||||||
|
bitRotateRight(Value[2], Shift), |
||||||
|
bitRotateRight(Value[3], Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genType> |
||||||
|
GLM_FUNC_QUALIFIER genType bitRotateLeft(genType const & In, std::size_t Shift) |
||||||
|
{ |
||||||
|
GLM_STATIC_ASSERT(std::numeric_limits<genType>::is_integer, "'bitRotateLeft' only accept integer values"); |
||||||
|
|
||||||
|
std::size_t BitSize = sizeof(genType) * 8; |
||||||
|
return (In >> Shift) | (In << (BitSize - Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec2<valType> bitRotateLeft |
||||||
|
( |
||||||
|
detail::tvec2<valType> const & Value, |
||||||
|
std::size_t Shift |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec2<valType>( |
||||||
|
bitRotateLeft(Value[0], Shift), |
||||||
|
bitRotateLeft(Value[1], Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec3<valType> bitRotateLeft |
||||||
|
( |
||||||
|
detail::tvec3<valType> const & Value, |
||||||
|
std::size_t Shift |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec3<valType>( |
||||||
|
bitRotateLeft(Value[0], Shift), |
||||||
|
bitRotateLeft(Value[1], Shift), |
||||||
|
bitRotateLeft(Value[2], Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename valType> |
||||||
|
GLM_FUNC_QUALIFIER detail::tvec4<valType> bitRotateLeft |
||||||
|
( |
||||||
|
detail::tvec4<valType> const & Value, |
||||||
|
std::size_t Shift |
||||||
|
) |
||||||
|
{ |
||||||
|
return detail::tvec4<valType>( |
||||||
|
bitRotateLeft(Value[0], Shift), |
||||||
|
bitRotateLeft(Value[1], Shift), |
||||||
|
bitRotateLeft(Value[2], Shift), |
||||||
|
bitRotateLeft(Value[3], Shift)); |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType> |
||||||
|
GLM_FUNC_QUALIFIER genIUType fillBitfieldWithOne |
||||||
|
( |
||||||
|
genIUType const & Value, |
||||||
|
int const & FromBit, |
||||||
|
int const & ToBit |
||||||
|
) |
||||||
|
{ |
||||||
|
assert(FromBit <= ToBit); |
||||||
|
assert(ToBit <= sizeof(genIUType) * std::size_t(8)); |
||||||
|
|
||||||
|
genIUType Result = Value; |
||||||
|
for(std::size_t i = 0; i <= ToBit; ++i) |
||||||
|
Result |= (1 << i); |
||||||
|
return Result; |
||||||
|
} |
||||||
|
|
||||||
|
template <typename genIUType> |
||||||
|
GLM_FUNC_QUALIFIER genIUType fillBitfieldWithZero |
||||||
|
( |
||||||
|
genIUType const & Value, |
||||||
|
int const & FromBit, |
||||||
|
int const & ToBit |
||||||
|
) |
||||||
|
{ |
||||||
|
assert(FromBit <= ToBit); |
||||||
|
assert(ToBit <= sizeof(genIUType) * std::size_t(8)); |
||||||
|
|
||||||
|
genIUType Result = Value; |
||||||
|
for(std::size_t i = 0; i <= ToBit; ++i) |
||||||
|
Result &= ~(1 << i); |
||||||
|
return Result; |
||||||
|
} |
||||||
|
|
||||||
|
}//namespace glm |
||||||
@ -0,0 +1,199 @@ |
|||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
// OpenGL Mathematics Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net)
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
// Created : 2011-04-21
|
||||||
|
// Updated : 2011-04-26
|
||||||
|
// Licence : This source is under MIT licence
|
||||||
|
// File : test/gtx/noise.cpp
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
|
||||||
|
#include <glm/glm.hpp> |
||||||
|
#include <glm/gtx/noise.hpp> |
||||||
|
#include <gli/gli.hpp> |
||||||
|
#include <gli/gtx/loader.hpp> |
||||||
|
#include <iostream> |
||||||
|
|
||||||
|
int test_simplex() |
||||||
|
{ |
||||||
|
std::size_t const Size = 256; |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::simplex(glm::vec2(x / 64.f, y / 64.f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_simplex2d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::simplex(glm::vec3(x / 64.f, y / 64.f, 0.5f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_simplex3d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::simplex(glm::vec4(x / 64.f, y / 64.f, 0.5f, 0.5f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_simplex4d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
return 0; |
||||||
|
} |
||||||
|
|
||||||
|
int test_perlin() |
||||||
|
{ |
||||||
|
std::size_t const Size = 256; |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::perlin(glm::vec2(x / 64.f, y / 64.f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_perlin2d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::perlin(glm::vec3(x / 64.f, y / 64.f, 0.5f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_perlin3d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::perlin(glm::vec4(x / 64.f, y / 64.f, 0.5f, 0.5f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_perlin4d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
return 0; |
||||||
|
} |
||||||
|
|
||||||
|
int test_perlin_pedioric() |
||||||
|
{ |
||||||
|
std::size_t const Size = 256; |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::perlin(glm::vec2(x / 64.f, y / 64.f), glm::vec2(2.0f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_perlin_pedioric_2d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::perlin(glm::vec3(x / 64.f, y / 64.f, 0.5f), glm::vec3(2.0f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_perlin_pedioric_3d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
{ |
||||||
|
std::vector<glm::byte> ImageData(Size * Size * 3); |
||||||
|
|
||||||
|
for(std::size_t y = 0; y < Size; ++y) |
||||||
|
for(std::size_t x = 0; x < Size; ++x) |
||||||
|
{ |
||||||
|
ImageData[(x + y * Size) * 3 + 0] = glm::byte(glm::perlin(glm::vec4(x / 64.f, y / 64.f, 0.5f, 0.5f), glm::vec4(2.0f)) * 128.f + 127.f); |
||||||
|
ImageData[(x + y * Size) * 3 + 1] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
ImageData[(x + y * Size) * 3 + 2] = ImageData[(x + y * Size) * 3 + 0]; |
||||||
|
} |
||||||
|
|
||||||
|
gli::texture2D Texture(1); |
||||||
|
Texture[0] = gli::image2D(glm::uvec2(Size), gli::RGB8U); |
||||||
|
memcpy(Texture[0].data(), &ImageData[0], ImageData.size()); |
||||||
|
gli::saveDDS9(Texture, "texture_perlin_pedioric_4d_256.dds"); |
||||||
|
} |
||||||
|
|
||||||
|
return 0; |
||||||
|
} |
||||||
|
|
||||||
|
int main() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
Error += test_simplex(); |
||||||
|
Error += test_perlin(); |
||||||
|
Error += test_perlin_pedioric(); |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
@ -0,0 +1,41 @@ |
|||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
// OpenGL Mathematics Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net)
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
// Created : 2011-10-13
|
||||||
|
// Updated : 2011-10-13
|
||||||
|
// Licence : This source is under MIT licence
|
||||||
|
// File : test/gtx/gradient_paint.cpp
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
|
||||||
|
#include <glm/glm.hpp> |
||||||
|
#include <glm/gtx/gradient_paint.hpp> |
||||||
|
|
||||||
|
int test_radialGradient() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
float Gradient = glm::radialGradient(glm::vec2(0), 1.0f, glm::vec2(1), glm::vec2(0.5)); |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
|
|
||||||
|
int test_linearGradient() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
float Gradient = glm::linearGradient(glm::vec2(0), glm::vec2(1), glm::vec2(0.5)); |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
|
|
||||||
|
int main() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
Error += test_radialGradient(); |
||||||
|
Error += test_linearGradient(); |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
|
|
||||||
|
|
||||||
@ -0,0 +1,66 @@ |
|||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
// OpenGL Mathematics Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net)
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
// Created : 2011-10-11
|
||||||
|
// Updated : 2011-10-11
|
||||||
|
// Licence : This source is under MIT licence
|
||||||
|
// File : test/gtx/gtx_integer.cpp
|
||||||
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
||||||
|
|
||||||
|
#include <glm/glm.hpp> |
||||||
|
#include <glm/gtx/integer.hpp> |
||||||
|
#include <glm/gtx/epsilon.hpp> |
||||||
|
#include <iostream> |
||||||
|
|
||||||
|
int test_floor_log2() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
for(std::size_t i = 1; i < 1000000; ++i) |
||||||
|
{ |
||||||
|
glm::uint A = glm::floor_log2(glm::uint(i)); |
||||||
|
glm::uint B = glm::uint(glm::log2(double(i))); // Will fail with float, lack of accuracy
|
||||||
|
|
||||||
|
Error += A == B ? 0 : 1; |
||||||
|
assert(!Error); |
||||||
|
} |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
|
|
||||||
|
int test_log2() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
for(std::size_t i = 1; i < 1000000; ++i) |
||||||
|
{ |
||||||
|
glm::uint A = glm::log2(glm::uint(i)); |
||||||
|
double B = glm::log2(double(i)); |
||||||
|
|
||||||
|
Error += glm::equalEpsilon(double(A), B, 1.0) ? 0 : 1; |
||||||
|
//assert(!Error);
|
||||||
|
} |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
|
|
||||||
|
int test_nlz() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
for(std::size_t i = 1; i < 33; ++i) |
||||||
|
printf("%d, %d\n", glm::nlz(i), 31u - glm::findMSB(i)); |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
|
|
||||||
|
int main() |
||||||
|
{ |
||||||
|
int Error = 0; |
||||||
|
|
||||||
|
Error += test_nlz(); |
||||||
|
Error += test_floor_log2(); |
||||||
|
Error += test_log2(); |
||||||
|
|
||||||
|
return Error; |
||||||
|
} |
||||||
Loading…
Reference in New Issue