3 void RandomSelection_Init()
5 RandomSelection_totalweight = 0;
6 RandomSelection_chosen_ent = NULL;
7 RandomSelection_chosen_float = 0;
8 RandomSelection_chosen_string = string_null;
9 RandomSelection_best_priority = -1;
12 void RandomSelection_Add(entity e, float f, string s, float weight, float priority)
14 if (priority > RandomSelection_best_priority)
16 RandomSelection_best_priority = priority;
17 RandomSelection_chosen_ent = e;
18 RandomSelection_chosen_float = f;
19 RandomSelection_chosen_string = s;
20 RandomSelection_totalweight = weight;
22 else if (priority == RandomSelection_best_priority)
24 RandomSelection_totalweight += weight;
25 if (random() * RandomSelection_totalweight <= weight)
27 RandomSelection_chosen_ent = e;
28 RandomSelection_chosen_float = f;
29 RandomSelection_chosen_string = s;
34 float DistributeEvenly_amount;
35 float DistributeEvenly_totalweight;
37 void DistributeEvenly_Init(float amount, float totalweight)
39 if (DistributeEvenly_amount)
41 LOG_TRACE("DistributeEvenly_Init: UNFINISHED DISTRIBUTION (", ftos(DistributeEvenly_amount), " for ", ftos(DistributeEvenly_totalweight), " left!)");
43 if (totalweight == 0) DistributeEvenly_amount = 0;
44 else DistributeEvenly_amount = amount;
45 DistributeEvenly_totalweight = totalweight;
48 float DistributeEvenly_Get(float weight)
51 if (weight <= 0) return 0;
52 f = floor(0.5 + DistributeEvenly_amount * weight / DistributeEvenly_totalweight);
53 DistributeEvenly_totalweight -= weight;
54 DistributeEvenly_amount -= f;
58 float DistributeEvenly_GetRandomized(float weight)
61 if (weight <= 0) return 0;
62 f = floor(random() + DistributeEvenly_amount * weight / DistributeEvenly_totalweight);
63 DistributeEvenly_totalweight -= weight;
64 DistributeEvenly_amount -= f;
68 // from the GNU Scientific Library
69 float gsl_ran_gaussian_lastvalue;
70 float gsl_ran_gaussian_lastvalue_set;
71 float gsl_ran_gaussian(float sigma)
73 if (gsl_ran_gaussian_lastvalue_set)
75 gsl_ran_gaussian_lastvalue_set = 0;
76 return sigma * gsl_ran_gaussian_lastvalue;
80 float a = random() * 2 * M_PI;
81 float b = sqrt(-2 * log(random()));
82 gsl_ran_gaussian_lastvalue = cos(a) * b;
83 gsl_ran_gaussian_lastvalue_set = 1;
84 return sigma * sin(a) * b;
88 // prandom - PREDICTABLE random number generator (not seeded yet)
95 c = crc16(false, strcat(ftos(prandom_seed), ftos(prandom_seed + M_PI)));
98 #ifdef USE_PRANDOM_DEBUG
99 LOG_TRACE("RANDOM -> ", ftos(c));
102 return c / 65536; // in [0..1[
120 void psrandom(float seed)
123 #ifdef USE_PRANDOM_DEBUG
124 LOG_TRACE("SRANDOM ", ftos(seed));
128 #ifdef USE_PRANDOM_DEBUG
131 LOG_TRACE("Current random seed = ", ftos(prandom_seed));