head     56.3;
access   paws bayes jws quist brad dew jwh;
symbols  ;
locks    ; strict;
comment  @# @;


56.3
date     93.01.27.13.30.00;  author jwh;  state Exp;
branches ;
next     56.2;

56.2
date     93.01.27.12.07.23;  author jwh;  state Exp;
branches ;
next     56.1;

56.1
date     91.11.05.09.51.51;  author jwh;  state Exp;
branches ;
next     55.1;

55.1
date     91.08.25.10.27.19;  author jwh;  state Exp;
branches ;
next     54.1;

54.1
date     91.03.18.15.30.17;  author jwh;  state Exp;
branches ;
next     53.1;

53.1
date     91.03.11.19.30.14;  author jwh;  state Exp;
branches ;
next     52.1;

52.1
date     91.02.19.09.14.13;  author jwh;  state Exp;
branches ;
next     51.1;

51.1
date     91.01.30.16.14.09;  author jwh;  state Exp;
branches ;
next     50.1;

50.1
date     90.10.29.16.28.58;  author jwh;  state Exp;
branches ;
next     49.1;

49.1
date     90.08.14.14.12.51;  author jwh;  state Exp;
branches ;
next     48.1;

48.1
date     90.07.26.11.19.11;  author jwh;  state Exp;
branches ;
next     47.1;

47.1
date     90.05.14.11.02.54;  author dew;  state Exp;
branches ;
next     46.1;

46.1
date     90.05.07.08.50.02;  author jwh;  state Exp;
branches ;
next     45.1;

45.1
date     90.04.19.15.57.58;  author jwh;  state Exp;
branches ;
next     44.1;

44.1
date     90.04.01.22.15.04;  author jwh;  state Exp;
branches ;
next     43.1;

43.1
date     90.03.20.14.07.06;  author jwh;  state Exp;
branches ;
next     42.1;

42.1
date     90.01.23.17.51.49;  author jwh;  state Exp;
branches ;
next     41.1;

41.1
date     89.12.22.11.34.07;  author jwh;  state Exp;
branches ;
next     40.1;

40.1
date     89.09.29.11.55.25;  author jwh;  state Exp;
branches ;
next     39.1;

39.1
date     89.09.26.16.40.03;  author dew;  state Exp;
branches ;
next     38.1;

38.1
date     89.08.29.11.32.01;  author jwh;  state Exp;
branches ;
next     37.1;

37.1
date     89.05.12.11.45.25;  author dew;  state Exp;
branches ;
next     36.1;

36.1
date     89.02.06.10.23.11;  author dew;  state Exp;
branches ;
next     35.1;

35.1
date     89.02.02.13.38.26;  author dew;  state Exp;
branches ;
next     34.1;

34.1
date     89.01.23.16.13.41;  author jwh;  state Exp;
branches ;
next     33.1;

33.1
date     89.01.16.11.45.23;  author dew;  state Exp;
branches ;
next     32.1;

32.1
date     89.01.10.11.54.19;  author bayes;  state Exp;
branches ;
next     31.1;

31.1
date     88.12.14.18.15.23;  author bayes;  state Exp;
branches ;
next     30.1;

30.1
date     88.12.09.13.52.18;  author dew;  state Exp;
branches ;
next     29.1;

29.1
date     88.10.31.15.37.12;  author bayes;  state Exp;
branches ;
next     28.1;

28.1
date     88.10.06.11.03.36;  author dew;  state Exp;
branches ;
next     27.1;

27.1
date     88.09.29.11.44.24;  author bayes;  state Exp;
branches ;
next     26.1;

26.1
date     88.09.28.13.22.09;  author bayes;  state Exp;
branches ;
next     25.1;

25.1
date     88.03.02.09.38.14;  author bayes;  state Exp;
branches ;
next     24.1;

24.1
date     87.08.31.10.04.37;  author jws;  state Exp;
branches ;
next     23.1;

23.1
date     87.08.26.10.42.41;  author bayes;  state Exp;
branches ;
next     22.1;

22.1
date     87.08.17.11.28.50;  author bayes;  state Exp;
branches ;
next     21.1;

21.1
date     87.08.12.14.11.31;  author bayes;  state Exp;
branches ;
next     20.1;

20.1
date     87.07.30.11.24.44;  author bayes;  state Exp;
branches ;
next     19.1;

19.1
date     87.06.01.08.37.43;  author jws;  state Exp;
branches ;
next     18.1;

18.1
date     87.05.20.15.41.14;  author bayes;  state Exp;
branches ;
next     17.1;

17.1
date     87.04.30.10.49.17;  author jws;  state Exp;
branches ;
next     16.1;

16.1
date     87.04.26.16.00.14;  author jws;  state Exp;
branches ;
next     15.1;

15.1
date     87.04.13.09.38.27;  author jws;  state Exp;
branches ;
next     14.1;

14.1
date     87.04.01.15.45.30;  author jws;  state Exp;
branches ;
next     13.1;

13.1
date     87.02.28.18.44.20;  author jws;  state Exp;
branches ;
next     12.1;

12.1
date     87.02.02.13.36.21;  author jws;  state Exp;
branches ;
next     11.1;

11.1
date     87.01.19.10.02.41;  author jws;  state Exp;
branches ;
next     10.1;

10.1
date     86.12.24.11.15.15;  author jws;  state Exp;
branches ;
next     9.1;

9.1
date     86.12.12.14.58.27;  author bayes;  state Exp;
branches ;
next     8.1;

8.1
date     86.11.27.12.11.21;  author jws;  state Exp;
branches ;
next     7.1;

7.1
date     86.11.20.14.17.15;  author hal;  state Exp;
branches ;
next     6.1;

6.1
date     86.11.04.18.15.22;  author paws;  state Exp;
branches ;
next     5.1;

5.1
date     86.10.28.17.03.27;  author hal;  state Exp;
branches ;
next     4.1;

4.1
date     86.09.30.19.59.42;  author hal;  state Exp;
branches ;
next     3.1;

3.1
date     86.09.01.12.09.57;  author hal;  state Exp;
branches ;
next     2.1;

2.1
date     86.07.30.14.58.41;  author hal;  state Exp;
branches ;
next     1.1;

1.1
date     86.06.30.15.59.34;  author danm;  state tmp;
branches ;
next     ;


desc
@Base file for PWS 3.2 release.

@


56.3
log
@
pws2rcs automatic delta on Wed Jan 27 13:14:25 MST 1993
@
text
@	mname   rnd

	src     module rnd;
	src     import sysglobals;
	src     export
	src       procedure random(var seed : integer);
	src       function rand( var seed : integer;
	src                      range : shortint) : shortint;
	src     end;

	def     rnd_rnd
	def     rnd_random
	def     rnd_rand
	def     random
	def     rand

*******************************************************************************
*
*       Procedure RANDOM(VAR SEED: INTEGER)
*
*       Description:
*               Generate a pseudo-random number with the formula
*               Xn <- (16807 * Xn-1) MOD (2^31 - 1), where Xn-1 is the
*               previous random number. A shortcut computation is:
*                   C <- 16807 * Xn-1.
*                   Xn <- C MOD 2^31 + C DIV 2^31.
*                   If Xn > 2^31 - 1, then Xn <- Xn - (2^31 - 1)
*
*       Parameters:
*               rndseed - the previous random number
*
*       Error conditions:
*               There are none.
*************************************************
rnd_rnd rts
rnd_random equ *
random  movea.l (sp)+,a0        return address
	movea.l (sp)+,a1        address of seed
	move.l  (a1),d0         get previous random seed Xn
	move.l  d0,d1           leave bottom 16 bits in d0
	swap    d1              get top 16 bits into d1
	mulu    #16807,d0       get one partial product in d0
	mulu    #16807,d1       high order partial product in d1
	swap    d1              align middle 16 bits of product in high d1
	add.w   d1,d1           most of (product div 2^31) is in low d1
	add.l   d1,d0           compute (product mod 2^31) + (product div 2^31)
	bcc.s   rnd1            any carries out of 32nd bit are part of the div
	addq.l  #2,d0           (so propagate into appropriate position)
rnd1    bpl.s   rnd2            bit 31 is also part of the div
	sub.l   #$7FFFFFFF,d0        so remove it and add it back to bit 0
rnd2    move.l  d0,(a1)
	jmp     (a0)

*************************************************
*
*       Function RAND(VAR SEED: INTEGER;
*               RANGE: SHORTINT): SHORTINT
*
*       Returns a 16 bit integer which is scaled
*       to the range 0..RANGE-1
*       (RANGE is treated as unsigned!)
*************************************************

rnd_rand equ *
rand    movea.l (sp)+,a2        return address
	move.w  (sp)+,d2        range parameter
	bsr.s   rnd_random      compute into d0
	asl.l   #1,d0           normalize
	swap    d0                   to 16 bits
	mulu    d2,d0           scale to range
	swap    d0
	move.w  d0,(sp)         return result
	jmp     (a2)

	end
@


56.2
log
@
pws2rcs automatic delta on Wed Jan 27 11:57:27 MST 1993
@
text
@d1 75
@


56.1
log
@Automatic bump of revision number for PWS version 3.25
@
text
@a0 75
	mname   rnd

	src     module rnd;
	src     import sysglobals;
	src     export
	src       procedure random(var seed : integer);
	src       function rand( var seed : integer;
	src                      range : shortint) : shortint;
	src     end;

	def     rnd_rnd
	def     rnd_random
	def     rnd_rand
	def     random
	def     rand

*******************************************************************************
*
*       Procedure RANDOM(VAR SEED: INTEGER)
*
*       Description:
*               Generate a pseudo-random number with the formula
*               Xn <- (16807 * Xn-1) MOD (2^31 - 1), where Xn-1 is the
*               previous random number. A shortcut computation is:
*                   C <- 16807 * Xn-1.
*                   Xn <- C MOD 2^31 + C DIV 2^31.
*                   If Xn > 2^31 - 1, then Xn <- Xn - (2^31 - 1)
*
*       Parameters:
*               rndseed - the previous random number
*
*       Error conditions:
*               There are none.
*************************************************
rnd_rnd rts
rnd_random equ *
random  movea.l (sp)+,a0        return address
	movea.l (sp)+,a1        address of seed
	move.l  (a1),d0         get previous random seed Xn
	move.l  d0,d1           leave bottom 16 bits in d0
	swap    d1              get top 16 bits into d1
	mulu    #16807,d0       get one partial product in d0
	mulu    #16807,d1       high order partial product in d1
	swap    d1              align middle 16 bits of product in high d1
	add.w   d1,d1           most of (product div 2^31) is in low d1
	add.l   d1,d0           compute (product mod 2^31) + (product div 2^31)
	bcc.s   rnd1            any carries out of 32nd bit are part of the div
	addq.l  #2,d0           (so propagate into appropriate position)
rnd1    bpl.s   rnd2            bit 31 is also part of the div
	sub.l   #$7FFFFFFF,d0        so remove it and add it back to bit 0
rnd2    move.l  d0,(a1)
	jmp     (a0)

*************************************************
*
*       Function RAND(VAR SEED: INTEGER;
*               RANGE: SHORTINT): SHORTINT
*
*       Returns a 16 bit integer which is scaled
*       to the range 0..RANGE-1
*       (RANGE is treated as unsigned!)
*************************************************

rnd_rand equ *
rand    movea.l (sp)+,a2        return address
	move.w  (sp)+,d2        range parameter
	bsr.s   rnd_random      compute into d0
	asl.l   #1,d0           normalize
	swap    d0                   to 16 bits
	mulu    d2,d0           scale to range
	swap    d0
	move.w  d0,(sp)         return result
	jmp     (a2)

	end
@


55.1
log
@Automatic bump of revision number for PWS version 3.25A
@
text
@@


54.1
log
@Automatic bump of revision number for PWS version 3.24
@
text
@@


53.1
log
@Automatic bump of revision number for PWS version 3.24B
@
text
@@


52.1
log
@Automatic bump of revision number for PWS version 3.24A
@
text
@@


51.1
log
@Automatic bump of revision number for PWS version 3.24d
@
text
@@


50.1
log
@Automatic bump of revision number for PWS version 3.23c
@
text
@@


49.1
log
@Automatic bump of revision number for PWS version 3.24b
@
text
@@


48.1
log
@Automatic bump of revision number for PWS version 3.24a
@
text
@@


47.1
log
@Automatic bump of revision number for PWS version 3.23
@
text
@@


46.1
log
@Automatic bump of revision number for PWS version 3.23
@
text
@@


45.1
log
@Automatic bump of revision number for PWS version 3.23C
@
text
@@


44.1
log
@Automatic bump of revision number for PWS version 3.23B
@
text
@@


43.1
log
@Automatic bump of revision number for PWS version 3.23aA
@
text
@@


42.1
log
@Automatic bump of revision number for PWS version 3.23e
@
text
@@


41.1
log
@Automatic bump of revision number for PWS version 3.23d
@
text
@@


40.1
log
@Automatic bump of revision number for PWS version 3.23c
@
text
@@


39.1
log
@Automatic bump of revision number for PWS version 3.23b
@
text
@@


38.1
log
@Automatic bump of revision number for PWS version 3.23a
@
text
@@


37.1
log
@Automatic bump of revision number for PWS version 3.3a
@
text
@@


36.1
log
@Automatic bump of revision number for PWS version 3.22
@
text
@@


35.1
log
@Automatic bump of revision number for PWS version 3.22
@
text
@@


34.1
log
@Automatic bump of revision number for PWS version 3.22
@
text
@@


33.1
log
@Automatic bump of revision number for PWS version 3.22D
@
text
@@


32.1
log
@Automatic bump of revision number for PWS version 3.22C
@
text
@@


31.1
log
@Automatic bump of revision number for PWS version 3.22B
@
text
@@


30.1
log
@Automatic bump of revision number for PWS version 3.22A
@
text
@@


29.1
log
@Automatic bump of revision number for PWS version 3.22b
@
text
@@


28.1
log
@Automatic bump of revision number for PWS version 3.3b
@
text
@@


27.1
log
@Automatic bump of revision number for PWS version 3.3a
@
text
@@


26.1
log
@Automatic bump of revision number for PWS version 3.3 Synch
@
text
@@


25.1
log
@Automatic bump of revision number for PWS version 3.2Y
@
text
@@


24.1
log
@Automatic bump of revision number for PWS version 3.2
@
text
@@


23.1
log
@Automatic bump of revision number for PWS version 3.2P
@
text
@@


22.1
log
@Automatic bump of revision number for PWS version 3.2N
@
text
@@


21.1
log
@Automatic bump of revision number for PWS version 3.2M
@
text
@@


20.1
log
@Automatic bump of revision number for PWS version 3.2L
@
text
@@


19.1
log
@Automatic bump of revision number for PWS version 3.2K
@
text
@@


18.1
log
@Automatic bump of revision number for PWS version 3.2J
@
text
@@


17.1
log
@Automatic bump of revision number for PWS version 3.2I+
@
text
@@


16.1
log
@Automatic bump of revision number for PWS version 3.2I
@
text
@@


15.1
log
@Automatic bump of revision number for PWS version 3.2H
@
text
@@


14.1
log
@Automatic bump of revision number for PWS version 3.2G
@
text
@@


13.1
log
@Automatic bump of revision number for PWS version 3.2F
@
text
@@


12.1
log
@Automatic bump of revision number for PWS version 3.2E
@
text
@@


11.1
log
@Automatic bump of revision number for PWS version 3.2D
@
text
@@


10.1
log
@Automatic bump of revision number for PWS version 3.2C
@
text
@@


9.1
log
@Automatic bump of revision number for PWS version 3.2B
@
text
@@


8.1
log
@Automatic bump of revision number for PWS version 3.2A
@
text
@@


7.1
log
@Automatic bump of revision number for PWS version 3.2l
@
text
@@


6.1
log
@Automatic bump of revision number for PWS version 3.2k
@
text
@@


5.1
log
@Automatic bump of revision number for PWS version 3.2j
@
text
@@


4.1
log
@Automatic bump of revision number for PWS version 3.2i
@
text
@@


3.1
log
@Automatic bump of revision number for PWS version 3.2h
@
text
@@


2.1
log
@Auto bump rev number to 2.1 for sys 3.2e.
@
text
@@


1.1
log
@Initial revision
@
text
@@
