1
2 /*
3 * M_APM - mapm_rnd.c
4 *
5 * Copyright (C) 1999 - 2007 Michael C. Ring
6 *
7 * Permission to use, copy, and distribute this software and its
8 * documentation for any purpose with or without fee is hereby granted,
9 * provided that the above copyright notice appear in all copies and
10 * that both that copyright notice and this permission notice appear
11 * in supporting documentation.
12 *
13 * Permission to modify the software is granted. Permission to distribute
14 * the modified code is granted. Modifications are to be distributed by
15 * using the file 'license.txt' as a template to modify the file header.
16 * 'license.txt' is available in the official MAPM distribution.
17 *
18 * This software is provided "as is" without express or implied warranty.
19 */
20
21 /*
22 * $Id: mapm_rnd.c,v 1.12 2007/12/03 01:47:17 mike Exp $
23 *
24 * This file contains the Random Number Generator function.
25 *
26 * $Log: mapm_rnd.c,v $
27 * Revision 1.12 2007/12/03 01:47:17 mike
28 * Update license
29 *
30 * Revision 1.11 2003/10/25 22:55:43 mike
31 * add support for National Instruments LabWindows CVI
32 *
33 * Revision 1.10 2002/11/03 22:41:03 mike
34 * Updated function parameters to use the modern style
35 *
36 * Revision 1.9 2002/02/14 21:50:45 mike
37 * add _set_random_seed
38 *
39 * Revision 1.8 2001/07/16 19:30:32 mike
40 * add function M_free_all_rnd
41 *
42 * Revision 1.7 2001/03/20 17:19:45 mike
43 * use a new multiplier
44 *
45 * Revision 1.6 2000/08/20 23:46:07 mike
46 * add more possible multupliers (no code changes)
47 *
48 * Revision 1.5 1999/09/19 23:32:14 mike
49 * added comments
50 *
51 * Revision 1.4 1999/09/18 03:49:25 mike
52 * *** empty log message ***
53 *
54 * Revision 1.3 1999/09/18 03:35:36 mike
55 * only prototype get_microsec for non-DOS
56 *
57 * Revision 1.2 1999/09/18 02:35:36 mike
58 * delete debug printf's
59 *
60 * Revision 1.1 1999/09/18 02:26:52 mike
61 * Initial revision
62 */
63
64 #include "m_apm_lc.h"
65
66 #ifndef _HAVE_NI_LABWIN_CVI_
67 #ifdef MSDOS
68 #include <time.h>
69 #include <sys/timeb.h>
70 #else
71 #include <sys/time.h>
72 extern void M_get_microsec(unsigned long *, long *);
73 #endif
74 #endif
75
76 #ifdef _HAVE_NI_LABWIN_CVI_
77 #include <time.h>
78 #include <utility.h>
79 #include <ansi/math.h>
80 #endif
81
82 extern void M_reverse_string(char *);
83 extern void M_get_rnd_seed(M_APM);
84
85 static M_APM M_rnd_aa;
86 static M_APM M_rnd_mm;
87 static M_APM M_rnd_XX;
88 static M_APM M_rtmp0;
89 static M_APM M_rtmp1;
90
91 static int M_firsttime2 = TRUE;
92
93 /*
94 Used Knuth's The Art of Computer Programming, Volume 2 as
95 the basis. Assuming the random number is X, compute
96 (where all the math is performed on integers) :
97
98 X = (a * X + c) MOD m
99
100 From Knuth:
101
102 'm' should be large, at least 2^30 : we use 1.0E+15
103
104 'a' should be between .01m and .99m and not have a simple
105 pattern. 'a' should not have any large factors in common
106 with 'm' and (since 'm' is a power of 10) if 'a' MOD 200
107 = 21 then all 'm' different possible values will be
108 generated before 'X' starts to repeat.
109
110 We use 'a' = 716805947629621.
111
112 This is a prime number and also meets 'a' MOD 200 = 21.
113 Commented out below are many potential multipliers that
114 are all prime and meet 'a' MOD 200 = 21.
115
116 There are few restrictions on 'c' except 'c' can have no
117 factor in common with 'm', hence we set 'c' = 'a'.
118
119 On the first call, the system time is used to initialize X.
120 */
121
122 /*
123 * the following constants are all potential multipliers. they are
124 * all prime numbers that also meet the criteria of NUM mod 200 = 21.
125 */
126
127 /*
128 439682071525421 439682071528421 439682071529221 439682071529821
129 439682071530421 439682071532021 439682071538821 439682071539421
130 439682071540021 439682071547021 439682071551221 439682071553821
131 439682071555421 439682071557221 439682071558021 439682071558621
132 439682071559821 439652381461621 439652381465221 439652381465621
133 439652381466421 439652381467421 439652381468621 439652381470021
134 439652381471221 439652381477021 439652381484221 439652381488421
135 439652381491021 439652381492021 439652381494021 439652381496821
136 617294387035621 617294387038621 617294387039221 617294387044421
137 617294387045221 617294387048621 617294387051621 617294387051821
138 617294387053621 617294387058421 617294387064221 617294387065621
139 617294387068621 617294387069221 617294387069821 617294387070421
140 617294387072021 617294387072621 617294387073821 617294387076821
141 649378126517621 649378126517821 649378126518221 649378126520821
142 649378126523821 649378126525621 649378126526621 649378126528421
143 649378126529621 649378126530821 649378126532221 649378126533221
144 649378126535221 649378126539421 649378126543621 649378126546021
145 649378126546421 649378126549421 649378126550821 649378126555021
146 649378126557421 649378126560221 649378126561621 649378126562021
147 649378126564621 649378126565821 672091582360421 672091582364221
148 672091582364621 672091582367021 672091582368421 672091582369021
149 672091582370821 672091582371421 672091582376821 672091582380821
150 716805243983221 716805243984821 716805947623621 716805947624621
151 716805947629021 716805947629621 716805947630621 716805947633621
152 716805947634221 716805947635021 716805947635621 716805947642221
153 */
154
155 /****************************************************************************/
M_free_all_rnd()156 void M_free_all_rnd()
157 {
158 if (M_firsttime2 == FALSE)
159 {
160 m_apm_free(M_rnd_aa);
161 m_apm_free(M_rnd_mm);
162 m_apm_free(M_rnd_XX);
163 m_apm_free(M_rtmp0);
164 m_apm_free(M_rtmp1);
165
166 M_firsttime2 = TRUE;
167 }
168 }
169 /****************************************************************************/
m_apm_set_random_seed(char * ss)170 void m_apm_set_random_seed(char *ss)
171 {
172 M_APM btmp;
173
174 if (M_firsttime2)
175 {
176 btmp = M_get_stack_var();
177 m_apm_get_random(btmp);
178 M_restore_stack(1);
179 }
180
181 m_apm_set_string(M_rnd_XX, ss);
182 }
183 /****************************************************************************/
184 /*
185 * compute X = (a * X + c) MOD m where c = a
186 */
m_apm_get_random(M_APM mrnd)187 void m_apm_get_random(M_APM mrnd)
188 {
189
190 if (M_firsttime2) /* use the system time as the initial seed value */
191 {
192 M_firsttime2 = FALSE;
193
194 M_rnd_aa = m_apm_init();
195 M_rnd_XX = m_apm_init();
196 M_rnd_mm = m_apm_init();
197 M_rtmp0 = m_apm_init();
198 M_rtmp1 = m_apm_init();
199
200 /* set the multiplier M_rnd_aa and M_rnd_mm */
201
202 m_apm_set_string(M_rnd_aa, "716805947629621");
203 m_apm_set_string(M_rnd_mm, "1.0E15");
204
205 M_get_rnd_seed(M_rnd_XX);
206 }
207
208 m_apm_multiply(M_rtmp0, M_rnd_XX, M_rnd_aa);
209 m_apm_add(M_rtmp1, M_rtmp0, M_rnd_aa);
210 m_apm_integer_div_rem(M_rtmp0, M_rnd_XX, M_rtmp1, M_rnd_mm);
211 m_apm_copy(mrnd, M_rnd_XX);
212 mrnd->m_apm_exponent -= 15;
213 }
214 /****************************************************************************/
M_reverse_string(char * s)215 void M_reverse_string(char *s)
216 {
217 int ct;
218 char ch, *p1, *p2;
219
220 if ((ct = strlen(s)) <= 1)
221 return;
222
223 p1 = s;
224 p2 = s + ct - 1;
225 ct /= 2;
226
227 while (TRUE)
228 {
229 ch = *p1;
230 *p1++ = *p2;
231 *p2-- = ch;
232
233 if (--ct == 0)
234 break;
235 }
236 }
237 /****************************************************************************/
238
239 #ifndef _HAVE_NI_LABWIN_CVI_
240
241 #ifdef MSDOS
242
243 /****************************************************************************/
244 /*
245 * for DOS / Win 9x/NT systems : use 'ftime'
246 */
M_get_rnd_seed(M_APM mm)247 void M_get_rnd_seed(M_APM mm)
248 {
249 int millisec;
250 time_t timestamp;
251 unsigned long ul;
252 char ss[32], buf1[48], buf2[32];
253 struct timeb timebuffer;
254 M_APM atmp;
255
256 atmp = M_get_stack_var();
257
258 ftime(&timebuffer);
259
260 millisec = (int)timebuffer.millitm;
261 timestamp = timebuffer.time;
262 ul = (unsigned long)(timestamp / 7);
263 ul += timestamp + 537;
264 strcpy(ss,ctime(×tamp)); /* convert to string and copy to ss */
265
266 sprintf(buf1,"%d",(millisec / 10));
267 sprintf(buf2,"%lu",ul);
268
269 ss[0] = ss[18];
270 ss[1] = ss[17];
271 ss[2] = ss[15];
272 ss[3] = ss[14];
273 ss[4] = ss[12];
274 ss[5] = ss[11];
275 ss[6] = ss[9];
276 ss[7] = ss[23];
277 ss[8] = ss[20];
278 ss[9] = '\0';
279
280 M_reverse_string(buf2);
281 strcat(buf1,buf2);
282 strcat(buf1,ss);
283
284 m_apm_set_string(atmp, buf1);
285 atmp->m_apm_exponent = 15;
286 m_apm_integer_divide(mm, atmp, MM_One);
287
288 M_restore_stack(1);
289 }
290 /****************************************************************************/
291
292 #else
293
294 /****************************************************************************/
295 /*
296 * for unix systems : use 'gettimeofday'
297 */
M_get_rnd_seed(M_APM mm)298 void M_get_rnd_seed(M_APM mm)
299 {
300 unsigned long sec3;
301 long usec3;
302 char buf1[32], buf2[32];
303 M_APM atmp;
304
305 atmp = M_get_stack_var();
306 M_get_microsec(&sec3,&usec3);
307
308 sprintf(buf1,"%ld",usec3);
309 sprintf(buf2,"%lu",sec3);
310 M_reverse_string(buf2);
311 strcat(buf1,buf2);
312
313 m_apm_set_string(atmp, buf1);
314 atmp->m_apm_exponent = 15;
315 m_apm_integer_divide(mm, atmp, MM_One);
316
317 M_restore_stack(1);
318 }
319 /****************************************************************************/
M_get_microsec(unsigned long * sec,long * usec)320 void M_get_microsec(unsigned long *sec, long *usec)
321 {
322 struct timeval time_now; /* current time for elapsed time check */
323 struct timezone time_zone; /* time zone for gettimeofday call */
324
325 gettimeofday(&time_now, &time_zone); /* get current time */
326
327 *sec = time_now.tv_sec;
328 *usec = time_now.tv_usec;
329 }
330 /****************************************************************************/
331
332 #endif
333 #endif
334
335 #ifdef _HAVE_NI_LABWIN_CVI_
336
337 /****************************************************************************/
338 /*
339 * for National Instruments LabWindows CVI
340 */
341
M_get_rnd_seed(M_APM mm)342 void M_get_rnd_seed(M_APM mm)
343 {
344 double timer0;
345 int millisec;
346 char *cvi_time, *cvi_date, buf1[64], buf2[32];
347 M_APM atmp;
348
349 atmp = M_get_stack_var();
350
351 cvi_date = DateStr();
352 cvi_time = TimeStr();
353 timer0 = Timer();
354
355 /*
356 * note that Timer() is not syncronized to TimeStr(),
357 * but we don't care here since we are just looking
358 * for a random source of digits.
359 */
360
361 millisec = (int)(0.01 + 1000.0 * (timer0 - floor(timer0)));
362
363 sprintf(buf1, "%d", millisec);
364
365 buf2[0] = cvi_time[6]; /* time format: "HH:MM:SS" */
366 buf2[1] = cvi_time[7];
367 buf2[2] = cvi_time[3];
368 buf2[3] = cvi_time[4];
369 buf2[4] = cvi_time[0];
370 buf2[5] = cvi_time[1];
371
372 buf2[6] = cvi_date[3]; /* date format: "MM-DD-YYYY" */
373 buf2[7] = cvi_date[4];
374 buf2[8] = cvi_date[0];
375 buf2[9] = cvi_date[1];
376 buf2[10] = cvi_date[8];
377 buf2[11] = cvi_date[9];
378 buf2[12] = cvi_date[7];
379
380 buf2[13] = '4';
381 buf2[14] = '7';
382 buf2[15] = '\0';
383
384 strcat(buf1, buf2);
385
386 m_apm_set_string(atmp, buf1);
387 atmp->m_apm_exponent = 15;
388 m_apm_integer_divide(mm, atmp, MM_One);
389
390 M_restore_stack(1);
391 }
392
393 #endif
394
395 /****************************************************************************/
396
397