aboutsummaryrefslogtreecommitdiff
path: root/shared/common/rand.c
blob: 4a54b3a247b2ebeba573f78eab72717bfa59d1e1 (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
/* Copyright (c) 1997-2004 Miller Puckette, krzYszcz, and others.
 * For information on usage and redistribution, and for a DISCLAIMER OF ALL
 * WARRANTIES, see the file, "LICENSE.txt," in this distribution.  */

#include <time.h>
#include "m_pd.h"
EXTERN double sys_getrealtime(void);  /* used to be in m_imp.h */
#include "rand.h"

/* borrowed from x_misc.c, LATER rethink */
void rand_seed(unsigned int *statep, unsigned int seed)
{
    if (seed) *statep = (seed & 0x7fffffff);
    else
    {
	/* LATER consider using time elapsed from system startup,
	   (or login -- in linux we might call getutent) */
	static unsigned int failsafe = 1489853723;
	static int shift = 0;
	static unsigned int lastticks = 0;
	/* LATER rethink -- it might fail on faster machine than mine
	   (but does it matter?) */
	unsigned int newticks = (unsigned int)(sys_getrealtime() * 1000000.);
	if (newticks == lastticks)
	{
	    failsafe = failsafe * 435898247 + 938284287;
	    *statep = (failsafe & 0x7fffffff);
#ifdef RAND_DEBUG
	    post("rand_seed failed (newticks %d)", newticks);
#endif
	}
	else
	{
	    if (!shift)
		shift = (int)time(0);  /* LATER deal with error return (-1) */
	    *statep = ((newticks + shift) & 0x7fffffff);
#if 0
	    post("rand_seed: newticks %d, shift %d", newticks, shift);
#endif
	}
	lastticks = newticks;
    }
}

/* borrowed from x_misc.c, LATER rethink */
int rand_int(unsigned int *statep, int range)
{
    int result;
    *statep = *statep * 472940017 + 832416023;
    result = ((double)range) * ((double)*statep) * (1./4294967296.);
    return (result < range ? result : range - 1);
}

float rand_unipolar(unsigned int *statep)
{
    float result;
    *statep = *statep * 472940017 + 832416023;
    result = (float)((double)*statep * (1./4294967296.));
    return (result);
}

/* borrowed from d_osc.c, LATER rethink */
float rand_bipolar(unsigned int *statep)
{
    float result = ((float)(((int)*statep & 0x7fffffff) - 0x40000000))
	* (float)(1.0 / 0x40000000);
    *statep = (unsigned)((int)*statep * 435898247 + 382842987);
    return (result);
}