summaryrefslogtreecommitdiff
path: root/test/randomtest.cc
blob: a1869f770efd1a15bcc4e153ec0df0a2b5dbcaae (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
/* -*- Mode: C++; tab-width: 2; indent-tabs-mode: nil; c-basic-offset: 2 -*- */
/***************************************************************************
 *            randomtest.cc
 *
 *  Sat Mar 26 08:48:53 CET 2016
 *  Copyright 2016 André Nusser
 *  andre.nusser@googlemail.com
 ****************************************************************************/

/*
 *  This file is part of DrumGizmo.
 *
 *  DrumGizmo is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU Lesser General Public License as published by
 *  the Free Software Foundation; either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  DrumGizmo is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with DrumGizmo; if not, write to the Free Software
 *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA.
 */
#include <uunit.h>

#include <random.h>

#include <vector>
#include <cmath>

class RandomTest
	: public uUnit
{
public:
	RandomTest()
	{
		uUNIT_TEST(RandomTest::rangeTest);
		uUNIT_TEST(RandomTest::normalTest);
		uUNIT_TEST(RandomTest::chooseTest);
	}

	void rangeTest()
	{
		// check if random numbers are in the wanted range
		Random rand;

		float float_lb = -42.23;
		float float_ub = 666.666;
		int int_lb = -42;
		int int_ub = 23;
		for (int i = 0; i<10000; i++)
		{
			float rand_float = rand.floatInRange(float_lb, float_ub);
			float rand_int = rand.intInRange(int_lb, int_ub);
			uUNIT_ASSERT(rand_float >= float_lb && rand_float <= float_ub);
			uUNIT_ASSERT(rand_int >= int_lb && rand_int <= int_ub);
		}

		// check if the series of random numbers is the one we expect
		// for a certain seed.
		rand = Random(666);
		uUNIT_ASSERT_EQUAL(0, rand.intInRange(0,100));
		uUNIT_ASSERT_EQUAL(61, rand.intInRange(0,100));
		uUNIT_ASSERT_EQUAL(23, rand.intInRange(0,100));
	}

	void normalTest()
	{
		// This test might theoretically fail but it is extremely unlikely to happen.
		Random rand;

		float real_mean = 42.0;
		float real_stddev = 2.0;
		int nr_of_samples = 50000;

		float sum = 0.;
		float sum_of_squares = 0.;
		for (int i=0; i<nr_of_samples; i++)
		{
			float rand_float = rand.normalDistribution(real_mean, real_stddev);
			sum += rand_float;
			sum_of_squares += rand_float*rand_float;
		}

		float estimated_mean = sum/nr_of_samples;
		float estimated_stddev = sqrt(sum_of_squares/nr_of_samples - estimated_mean*estimated_mean);

		float epsilon = 0.1;
		uUNIT_ASSERT(estimated_mean >= real_mean-epsilon && estimated_mean <= real_mean+epsilon);
		uUNIT_ASSERT(estimated_stddev >= real_stddev-epsilon && estimated_stddev <= real_stddev+epsilon);
	}

	void chooseTest()
	{
		Random rand;

		std::vector<int> vec = { 42, 42, 42 };
		int nr_of_samples = 1000;

		for (int i=0; i<nr_of_samples; i++)
		{
			uUNIT_ASSERT_EQUAL(42, rand.choose(vec));
		}
	}
};

// Registers the fixture into the 'registry'
static RandomTest test;