server/common/rng.C

/*
 * This file is part of Deliantra, the Roguelike Realtime MMORPG.
 *
 * Copyright (©) 2017,2018 Marc Alexander Lehmann / the Deliantra team
 * Copyright (©) 2010,2011,2012,2013,2014,2015,2016 Marc Alexander Lehmann / Robin Redeker / the Deliantra team
 *
 * Deliantra is free software: you can redistribute it and/or modify it under
 * the terms of the Affero GNU General Public License as published by the
 * Free Software Foundation, either version 3 of the License, or (at your
 * option) any later version.
 *
 * This program 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 General Public License for more details.
 *
 * You should have received a copy of the Affero GNU General Public License
 * and the GNU General Public License along with this program. If not, see
 * <http://www.gnu.org/licenses/>.
 *
 * The authors can be reached via e-mail to <support@deliantra.net>
 */

#include <cmath>

#include <global.h>

rand_gen rndm, rmg_rndm;

ecb_noinline void
tausworthe_rng::seed (uint32_t seed)
{
  state [0] = seed      * 69069U; if (state [0] <   2U) state [0] +=   2U;
  state [1] = state [0] * 69069U; if (state [1] <   8U) state [1] +=   8U;
  state [2] = state [1] * 69069U; if (state [2] <  16U) state [2] +=  16U;
  state [3] = state [2] * 69069U; if (state [3] < 128U) state [3] += 128U;

  for (int i = 11; --i; )
    next ();
}

uint32_t
tausworthe_rng::next ()
{
  state [0] = ((state [0] & 0xFFFFFFFEU) << 18U) ^ (((state [0] <<  6U) ^ state [0]) >> 13U);
  state [1] = ((state [1] & 0xFFFFFFF8U) <<  2U) ^ (((state [1] <<  2U) ^ state [1]) >> 27U);
  state [2] = ((state [2] & 0xFFFFFFF0U) <<  7U) ^ (((state [2] << 13U) ^ state [2]) >> 21U);
  state [3] = ((state [3] & 0xFFFFFF80U) << 13U) ^ (((state [3] <<  3U) ^ state [3]) >> 12U);

  return state [0] ^ state [1] ^ state [2] ^ state [3];
}

ecb_noinline void
r250521_rng::seed (uint32_t seed)
{
  xorshift_rng rng;

  rng.seed (seed);

  r250.seed_rng (rng);
  r521.seed_rng (rng);
}

uint32_t
r250521_rng::next ()
{
  return r250.next () ^ r521.next ();
}

template<class generator>
uint32_t
random_number_generator<generator>::get_range (uint32_t num)
{
  return (this->next () * (uint64_t)num) >> 32U;
}

template<class generator>
ecb_noinline uint32_t
random_number_generator<generator>::get_u32 ()
{
  return generator::next ();
}

template<class generator>
ecb_noinline uint64_t
random_number_generator<generator>::get_u64 ()
{
  return (uint64_t (get_u32 ()) << 32U) | get_u32 ();
}

template<class generator>
ecb_noinline float
random_number_generator<generator>::get_float ()
{
  return this->next () / (float)0x100000000ULL;
}

template<class generator>
ecb_noinline double
random_number_generator<generator>::get_double ()
{
  return this->next () / (double)0x100000000ULL;
}

// return a number within (min .. max)
template<class generator>
int
random_number_generator<generator>::get_range (int r_min, int r_max)
{
  return r_min + get_range (max (r_max - r_min + 1, 0));
}

template<class generator>
uint32_t
random_number_generator<generator>::poisson (double mean)
{
  double g = std::exp (-mean);
  uint32_t em = 0;

  double t = get_double ();

  while (t > g)
    {
      ++em;
      t *= get_double ();
    }

  return em;
}

template struct random_number_generator<bcpl_rng>;
template struct random_number_generator<borosh_niederreiter_rng>;
template struct random_number_generator<tausworthe_rng>;
template struct random_number_generator<xorshift_rng>;
template struct random_number_generator<freeciv_rng>;
template struct random_number_generator<r250_rng>;
template struct random_number_generator<r521_rng>;
template struct random_number_generator<r250521_rng>;

Revision:	1.9
Committed:	Sat Nov 17 23:40:00 2018 UTC (5 years, 6 months ago) by root
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Changes since 1.8:	+1 -0 lines
Log Message:	copyright update 2018
#	User	Rev	Content
1	root	1.1	/*
2			* This file is part of Deliantra, the Roguelike Realtime MMORPG.
3	root	1.5	*
4	root	1.9	* Copyright (©) 2017,2018 Marc Alexander Lehmann / the Deliantra team
5	root	1.6	* Copyright (©) 2010,2011,2012,2013,2014,2015,2016 Marc Alexander Lehmann / Robin Redeker / the Deliantra team
6	root	1.5	*
7	root	1.1	* Deliantra is free software: you can redistribute it and/or modify it under
8			* the terms of the Affero GNU General Public License as published by the
9			* Free Software Foundation, either version 3 of the License, or (at your
10			* option) any later version.
11	root	1.5	*
12	root	1.1	* This program is distributed in the hope that it will be useful,
13			* but WITHOUT ANY WARRANTY; without even the implied warranty of
14			* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15			* GNU General Public License for more details.
16	root	1.5	*
17	root	1.1	* You should have received a copy of the Affero GNU General Public License
18			* and the GNU General Public License along with this program. If not, see
19			* <http://www.gnu.org/licenses/>.
20	root	1.5	*
21	root	1.1	* The authors can be reached via e-mail to <support@deliantra.net>
22			*/
23
24	root	1.2	#include <cmath>
25
26	root	1.1	#include <global.h>
27
28			rand_gen rndm, rmg_rndm;
29
30	root	1.8	ecb_noinline void
31	root	1.2	tausworthe_rng::seed (uint32_t seed)
32	root	1.1	{
33			state [0] = seed * 69069U; if (state [0] < 2U) state [0] += 2U;
34			state [1] = state [0] * 69069U; if (state [1] < 8U) state [1] += 8U;
35			state [2] = state [1] * 69069U; if (state [2] < 16U) state [2] += 16U;
36			state [3] = state [2] * 69069U; if (state [3] < 128U) state [3] += 128U;
37
38			for (int i = 11; --i; )
39			next ();
40			}
41
42			uint32_t
43	root	1.2	tausworthe_rng::next ()
44	root	1.1	{
45			state [0] = ((state [0] & 0xFFFFFFFEU) << 18U) ^ (((state [0] << 6U) ^ state [0]) >> 13U);
46			state [1] = ((state [1] & 0xFFFFFFF8U) << 2U) ^ (((state [1] << 2U) ^ state [1]) >> 27U);
47			state [2] = ((state [2] & 0xFFFFFFF0U) << 7U) ^ (((state [2] << 13U) ^ state [2]) >> 21U);
48			state [3] = ((state [3] & 0xFFFFFF80U) << 13U) ^ (((state [3] << 3U) ^ state [3]) >> 12U);
49
50			return state [0] ^ state [1] ^ state [2] ^ state [3];
51			}
52
53	root	1.8	ecb_noinline void
54	root	1.2	r250521_rng::seed (uint32_t seed)
55	root	1.1	{
56	root	1.2	xorshift_rng rng;
57	root	1.1
58			rng.seed (seed);
59
60			r250.seed_rng (rng);
61			r521.seed_rng (rng);
62			}
63
64			uint32_t
65	root	1.2	r250521_rng::next ()
66	root	1.1	{
67			return r250.next () ^ r521.next ();
68			}
69
70			template<class generator>
71			uint32_t
72			random_number_generator<generator>::get_range (uint32_t num)
73			{
74			return (this->next () * (uint64_t)num) >> 32U;
75			}
76
77			template<class generator>
78	root	1.8	ecb_noinline uint32_t
79	root	1.1	random_number_generator<generator>::get_u32 ()
80			{
81			return generator::next ();
82			}
83
84			template<class generator>
85	root	1.8	ecb_noinline uint64_t
86	root	1.1	random_number_generator<generator>::get_u64 ()
87			{
88			return (uint64_t (get_u32 ()) << 32U) \| get_u32 ();
89			}
90
91	root	1.2	template<class generator>
92	root	1.8	ecb_noinline float
93	root	1.2	random_number_generator<generator>::get_float ()
94			{
95			return this->next () / (float)0x100000000ULL;
96			}
97
98			template<class generator>
99	root	1.8	ecb_noinline double
100	root	1.2	random_number_generator<generator>::get_double ()
101			{
102			return this->next () / (double)0x100000000ULL;
103			}
104
105	root	1.1	// return a number within (min .. max)
106			template<class generator>
107	root	1.7	int
108	root	1.1	random_number_generator<generator>::get_range (int r_min, int r_max)
109			{
110			return r_min + get_range (max (r_max - r_min + 1, 0));
111			}
112
113	root	1.2	template<class generator>
114			uint32_t
115			random_number_generator<generator>::poisson (double mean)
116			{
117			double g = std::exp (-mean);
118			uint32_t em = 0;
119
120			double t = get_double ();
121
122			while (t > g)
123			{
124			++em;
125			t *= get_double ();
126			}
127
128			return em;
129			}
130
131			template struct random_number_generator<bcpl_rng>;
132			template struct random_number_generator<borosh_niederreiter_rng>;
133			template struct random_number_generator<tausworthe_rng>;
134			template struct random_number_generator<xorshift_rng>;
135			template struct random_number_generator<freeciv_rng>;
136			template struct random_number_generator<r250_rng>;
137			template struct random_number_generator<r521_rng>;
138			template struct random_number_generator<r250521_rng>;
139	root	1.1