5070880428
This allows building against the same version of GLM across all targets and not deal with API changes of newer GLM versions which introduce regressions. Long term we want to migrate away from GLM and stick to a minimal self-written subset.
408 lines
11 KiB
C++
408 lines
11 KiB
C++
///////////////////////////////////////////////////////////////////////////////////
|
|
/// OpenGL Mathematics (glm.g-truc.net)
|
|
///
|
|
/// Copyright (c) 2005 - 2015 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.
|
|
///
|
|
/// Restrictions:
|
|
/// By making use of the Software for military purposes, you choose to make
|
|
/// a Bunny unhappy.
|
|
///
|
|
/// 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.
|
|
///
|
|
/// @file test/gtc/gtc_random.cpp
|
|
/// @date 2011-09-19 / 2014-11-25
|
|
/// @author Christophe Riccio
|
|
///////////////////////////////////////////////////////////////////////////////////
|
|
|
|
#include <glm/gtc/random.hpp>
|
|
#include <glm/gtc/epsilon.hpp>
|
|
#if(GLM_LANG & GLM_LANG_CXX0X_FLAG)
|
|
# include <array>
|
|
#endif
|
|
|
|
int test_linearRand()
|
|
{
|
|
int Error = 0;
|
|
|
|
glm::int32 const Min = 16;
|
|
glm::int32 const Max = 32;
|
|
|
|
{
|
|
glm::u8vec2 AMin(std::numeric_limits<glm::u8>::max());
|
|
glm::u8vec2 AMax(std::numeric_limits<glm::u8>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::u8vec2 A = glm::linearRand(glm::u8vec2(Min), glm::u8vec2(Max));
|
|
AMin = glm::min(AMin, A);
|
|
AMax = glm::max(AMax, A);
|
|
|
|
if(!glm::all(glm::lessThanEqual(A, glm::u8vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(A, glm::u8vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(AMin, glm::u8vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(AMax, glm::u8vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
glm::u16vec2 BMin(std::numeric_limits<glm::u16>::max());
|
|
glm::u16vec2 BMax(std::numeric_limits<glm::u16>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::u16vec2 B = glm::linearRand(glm::u16vec2(Min), glm::u16vec2(Max));
|
|
BMin = glm::min(BMin, B);
|
|
BMax = glm::max(BMax, B);
|
|
|
|
if(!glm::all(glm::lessThanEqual(B, glm::u16vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(B, glm::u16vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(BMin, glm::u16vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(BMax, glm::u16vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
glm::u32vec2 CMin(std::numeric_limits<glm::u32>::max());
|
|
glm::u32vec2 CMax(std::numeric_limits<glm::u32>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::u32vec2 C = glm::linearRand(glm::u32vec2(Min), glm::u32vec2(Max));
|
|
CMin = glm::min(CMin, C);
|
|
CMax = glm::max(CMax, C);
|
|
|
|
if(!glm::all(glm::lessThanEqual(C, glm::u32vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(C, glm::u32vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(CMin, glm::u32vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(CMax, glm::u32vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
glm::u64vec2 DMin(std::numeric_limits<glm::u64>::max());
|
|
glm::u64vec2 DMax(std::numeric_limits<glm::u64>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::u64vec2 D = glm::linearRand(glm::u64vec2(Min), glm::u64vec2(Max));
|
|
DMin = glm::min(DMin, D);
|
|
DMax = glm::max(DMax, D);
|
|
|
|
if(!glm::all(glm::lessThanEqual(D, glm::u64vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(D, glm::u64vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(DMin, glm::u64vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(DMax, glm::u64vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
}
|
|
|
|
{
|
|
glm::i8vec2 AMin(std::numeric_limits<glm::i8>::max());
|
|
glm::i8vec2 AMax(std::numeric_limits<glm::i8>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::i8vec2 A = glm::linearRand(glm::i8vec2(Min), glm::i8vec2(Max));
|
|
AMin = glm::min(AMin, A);
|
|
AMax = glm::max(AMax, A);
|
|
|
|
if(!glm::all(glm::lessThanEqual(A, glm::i8vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(A, glm::i8vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(AMin, glm::i8vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(AMax, glm::i8vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
glm::i16vec2 BMin(std::numeric_limits<glm::i16>::max());
|
|
glm::i16vec2 BMax(std::numeric_limits<glm::i16>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::i16vec2 B = glm::linearRand(glm::i16vec2(Min), glm::i16vec2(Max));
|
|
BMin = glm::min(BMin, B);
|
|
BMax = glm::max(BMax, B);
|
|
|
|
if(!glm::all(glm::lessThanEqual(B, glm::i16vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(B, glm::i16vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(BMin, glm::i16vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(BMax, glm::i16vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
glm::i32vec2 CMin(std::numeric_limits<glm::i32>::max());
|
|
glm::i32vec2 CMax(std::numeric_limits<glm::i32>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::i32vec2 C = glm::linearRand(glm::i32vec2(Min), glm::i32vec2(Max));
|
|
CMin = glm::min(CMin, C);
|
|
CMax = glm::max(CMax, C);
|
|
|
|
if(!glm::all(glm::lessThanEqual(C, glm::i32vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(C, glm::i32vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(CMin, glm::i32vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(CMax, glm::i32vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
glm::i64vec2 DMin(std::numeric_limits<glm::i64>::max());
|
|
glm::i64vec2 DMax(std::numeric_limits<glm::i64>::min());
|
|
{
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::i64vec2 D = glm::linearRand(glm::i64vec2(Min), glm::i64vec2(Max));
|
|
DMin = glm::min(DMin, D);
|
|
DMax = glm::max(DMax, D);
|
|
|
|
if(!glm::all(glm::lessThanEqual(D, glm::i64vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(D, glm::i64vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
Error += glm::all(glm::equal(DMin, glm::i64vec2(Min))) ? 0 : 1;
|
|
Error += glm::all(glm::equal(DMax, glm::i64vec2(Max))) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
}
|
|
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
glm::f32vec2 const A(glm::linearRand(glm::f32vec2(static_cast<float>(Min)), glm::f32vec2(static_cast<float>(Max))));
|
|
if(!glm::all(glm::lessThanEqual(A, glm::f32vec2(static_cast<float>(Max)))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(A, glm::f32vec2(static_cast<float>(Min)))))
|
|
++Error;
|
|
|
|
glm::f64vec2 const B(glm::linearRand(glm::f64vec2(Min), glm::f64vec2(Max)));
|
|
if(!glm::all(glm::lessThanEqual(B, glm::f64vec2(Max))))
|
|
++Error;
|
|
if(!glm::all(glm::greaterThanEqual(B, glm::f64vec2(Min))))
|
|
++Error;
|
|
assert(!Error);
|
|
}
|
|
|
|
{
|
|
float ResultFloat = 0.0f;
|
|
double ResultDouble = 0.0f;
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
ResultFloat += glm::linearRand(-1.0f, 1.0f);
|
|
ResultDouble += glm::linearRand(-1.0, 1.0);
|
|
}
|
|
|
|
Error += glm::epsilonEqual(ResultFloat, 0.0f, 0.0001f);
|
|
Error += glm::epsilonEqual(ResultDouble, 0.0, 0.0001);
|
|
assert(!Error);
|
|
}
|
|
|
|
return Error;
|
|
}
|
|
|
|
int test_circularRand()
|
|
{
|
|
int Error = 0;
|
|
|
|
{
|
|
std::size_t Max = 100000;
|
|
float ResultFloat = 0.0f;
|
|
double ResultDouble = 0.0f;
|
|
double Radius = 2.0f;
|
|
|
|
for(std::size_t i = 0; i < Max; ++i)
|
|
{
|
|
ResultFloat += glm::length(glm::circularRand(1.0f));
|
|
ResultDouble += glm::length(glm::circularRand(Radius));
|
|
}
|
|
|
|
Error += glm::epsilonEqual(ResultFloat, float(Max), 0.01f) ? 0 : 1;
|
|
Error += glm::epsilonEqual(ResultDouble, double(Max) * double(Radius), 0.01) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
return Error;
|
|
}
|
|
|
|
int test_sphericalRand()
|
|
{
|
|
int Error = 0;
|
|
|
|
{
|
|
std::size_t Max = 100000;
|
|
float ResultFloatA = 0.0f;
|
|
float ResultFloatB = 0.0f;
|
|
float ResultFloatC = 0.0f;
|
|
double ResultDoubleA = 0.0f;
|
|
double ResultDoubleB = 0.0f;
|
|
double ResultDoubleC = 0.0f;
|
|
|
|
for(std::size_t i = 0; i < Max; ++i)
|
|
{
|
|
ResultFloatA += glm::length(glm::sphericalRand(1.0f));
|
|
ResultDoubleA += glm::length(glm::sphericalRand(1.0));
|
|
ResultFloatB += glm::length(glm::sphericalRand(2.0f));
|
|
ResultDoubleB += glm::length(glm::sphericalRand(2.0));
|
|
ResultFloatC += glm::length(glm::sphericalRand(3.0f));
|
|
ResultDoubleC += glm::length(glm::sphericalRand(3.0));
|
|
}
|
|
|
|
Error += glm::epsilonEqual(ResultFloatA, float(Max), 0.01f) ? 0 : 1;
|
|
Error += glm::epsilonEqual(ResultDoubleA, double(Max), 0.0001) ? 0 : 1;
|
|
Error += glm::epsilonEqual(ResultFloatB, float(Max * 2), 0.01f) ? 0 : 1;
|
|
Error += glm::epsilonEqual(ResultDoubleB, double(Max * 2), 0.0001) ? 0 : 1;
|
|
Error += glm::epsilonEqual(ResultFloatC, float(Max * 3), 0.01f) ? 0 : 1;
|
|
Error += glm::epsilonEqual(ResultDoubleC, double(Max * 3), 0.01) ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
return Error;
|
|
}
|
|
|
|
int test_diskRand()
|
|
{
|
|
int Error = 0;
|
|
|
|
{
|
|
float ResultFloat = 0.0f;
|
|
double ResultDouble = 0.0f;
|
|
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
ResultFloat += glm::length(glm::diskRand(2.0f));
|
|
ResultDouble += glm::length(glm::diskRand(2.0));
|
|
}
|
|
|
|
Error += ResultFloat < 200000.f ? 0 : 1;
|
|
Error += ResultDouble < 200000.0 ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
return Error;
|
|
}
|
|
|
|
int test_ballRand()
|
|
{
|
|
int Error = 0;
|
|
|
|
{
|
|
float ResultFloat = 0.0f;
|
|
double ResultDouble = 0.0f;
|
|
|
|
for(std::size_t i = 0; i < 100000; ++i)
|
|
{
|
|
ResultFloat += glm::length(glm::ballRand(2.0f));
|
|
ResultDouble += glm::length(glm::ballRand(2.0));
|
|
}
|
|
|
|
Error += ResultFloat < 200000.f ? 0 : 1;
|
|
Error += ResultDouble < 200000.0 ? 0 : 1;
|
|
assert(!Error);
|
|
}
|
|
|
|
return Error;
|
|
}
|
|
/*
|
|
#if(GLM_LANG & GLM_LANG_CXX0X_FLAG)
|
|
int test_grid()
|
|
{
|
|
int Error = 0;
|
|
|
|
typedef std::array<int, 8> colors;
|
|
typedef std::array<int, 8 * 8> grid;
|
|
|
|
grid Grid;
|
|
colors Colors;
|
|
|
|
grid GridBest;
|
|
colors ColorsBest;
|
|
|
|
while(true)
|
|
{
|
|
for(std::size_t i = 0; i < Grid.size(); ++i)
|
|
Grid[i] = int(glm::linearRand(0.0, 8.0 * 8.0 * 8.0 - 1.0) / 64.0);
|
|
|
|
for(std::size_t i = 0; i < Grid.size(); ++i)
|
|
++Colors[Grid[i]];
|
|
|
|
bool Exit = true;
|
|
for(std::size_t i = 0; i < Colors.size(); ++i)
|
|
{
|
|
if(Colors[i] == 8)
|
|
continue;
|
|
|
|
Exit = false;
|
|
break;
|
|
}
|
|
|
|
if(Exit == true)
|
|
break;
|
|
}
|
|
|
|
return Error;
|
|
}
|
|
#endif
|
|
*/
|
|
int main()
|
|
{
|
|
int Error = 0;
|
|
|
|
Error += test_linearRand();
|
|
Error += test_circularRand();
|
|
Error += test_sphericalRand();
|
|
Error += test_diskRand();
|
|
Error += test_ballRand();
|
|
/*
|
|
#if(GLM_LANG & GLM_LANG_CXX0X_FLAG)
|
|
Error += test_grid();
|
|
#endif
|
|
*/
|
|
return Error;
|
|
}
|