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
112
113
114
115
116
|
/**
* $Id$
*
* ***** BEGIN GPL/BL DUAL LICENSE BLOCK *****
*
* This program is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License
* as published by the Free Software Foundation; either version 2
* of the License, or (at your option) any later version. The Blender
* Foundation also sells licenses for use in proprietary software under
* the Blender License. See http://www.blender.org/BL/ for information
* about this.
*
* 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 GNU General Public License
* along with this program; if not, write to the Free Software Foundation,
* Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
*
* The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
* All rights reserved.
*
* The Original Code is: all of this file.
*
* Contributor(s): none yet.
*
* ***** END GPL/BL DUAL LICENSE BLOCK *****
*/
#include <stdlib.h>
#include <string.h>
#include "PIL_time.h"
#include "BLI_rand.h"
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#if defined(WIN32) && !defined(FREE_WINDOWS)
typedef unsigned __int64 r_uint64;
#else
typedef unsigned long long r_uint64;
#endif
#define MULTIPLIER 0x5DEECE66D
#define ADDEND 0xB
#define LOWSEED 0x330E
static r_uint64 X= 0;
void BLI_srand(unsigned int seed) {
X= (((r_uint64) seed)<<16) | LOWSEED;
}
int BLI_rand(void) {
X= (MULTIPLIER*X + ADDEND)&0x0000FFFFFFFFFFFF;
return (int) (X>>17);
}
double BLI_drand(void) {
return (double) BLI_rand()/0x80000000;
}
float BLI_frand(void) {
return (float) BLI_rand()/0x80000000;
}
void BLI_storerand(unsigned int loc_r[2]) {
loc_r[0]= (unsigned int) (X>>32);
loc_r[1]= (unsigned int) (X&0xFFFFFFFF);
}
void BLI_restorerand(unsigned int loc[2]) {
X= ((r_uint64) loc[0])<<32;
X|= loc[1];
}
void BLI_fillrand(void *addr, int len) {
unsigned char *p= addr;
unsigned int save[2];
BLI_storerand(save);
BLI_srand((unsigned int) (PIL_check_seconds_timer()*0x7FFFFFFF));
while (len--) *p++= BLI_rand()&0xFF;
BLI_restorerand(save);
}
void BLI_array_randomize(void *data, int elemSize, int numElems, unsigned int seed)
{
unsigned int oldrand[2];
int i = numElems;
void *temp = malloc(elemSize);
BLI_storerand(oldrand);
BLI_srand(seed);
while (--i) {
int j = BLI_rand()%i;
void *iElem = (unsigned char*)data + i*elemSize;
void *jElem = (unsigned char*)data + j*elemSize;
memcpy(temp, iElem, elemSize);
memcpy(iElem, jElem, elemSize);
memcpy(jElem, temp, elemSize);
}
BLI_restorerand(oldrand);
free(temp);
}
|