home *** CD-ROM | disk | FTP | other *** search
/ CP/M / CPM_CDROM.iso / simtel / sigm / vols000 / vol021 / random.src < prev    next >
Encoding:
Text File  |  1984-04-29  |  3.0 KB  |  181 lines
  1. ;==================================================
  2. ;  PROGRAM TITLE:    Random Number Generator
  3. ; 
  4. ;  WRITTEN BY:        Raymond E. Penley
  5. ;  DATE WRITTEN:    27 June 1980
  6. ;
  7. ;  WRITTEN FOR:        Pascal/Z Users Group
  8. ;
  9. ;  SUMMARY:
  10. ;    Implements a Fibonacci series Random number generator.
  11. ;     RANDOM will return numbers from 0 to x
  12. ;
  13. ;         Call as              Returns
  14. ;         -------            ------------
  15. ;         real := RANDOM(10);        0.0 to 10.0
  16. ;         real := RANDOM(112);    0.0 to 112.0
  17. ;          I := TRUNC(RANDOM(10));      0 to 10
  18. ;
  19. ; **
  20. ;   Add these lines to your PASCAL source program:
  21. ;
  22. ;     Procedure SEEDRAND; EXTERNAL;
  23. ;     Functin RANDOM(X: Integer): REAL; external;
  24. ;
  25. ;   Also within the body of the main program
  26. ;   but BEFORE calling RANDOM(X);
  27. ; 
  28. ;       SEEDRAND;
  29. ;
  30. ;*** What would happen if you did not call SEEDRAND ? ***
  31. ; 
  32. ;===============================================
  33. ;
  34. ; PROCEDURE SEEDRAND;
  35. ; (* INITIAL VALUES FOR SEED1 AND SEED2 ARE HERE  *)
  36. ;
  37.      NAME RANDOM
  38.      ENTRY SEEDRAND,RANDOM
  39. ;
  40. SEEDRAND:
  41.     ENTR    D,2,0
  42. ;    SEED1 := 10946;
  43.     MVI    0(IY),42
  44.     MVI    -1(IY),194
  45. ;    SEED2 := 17711
  46. ; END;
  47.     MVI    -2(IY),69
  48.     MVI    -3(IY),47
  49.     EXIT    D,0
  50. ;
  51. ;
  52. ;
  53. ; Function Random(x: integer): real;
  54. ; (*
  55. ; GLOBAL
  56. ;   SEED1, SEED2 : INTEGER        *)
  57. ; CONST
  58. ;   factor = Maxint;
  59. ;   HALFINT = 16383; (* 1/2 OF MAXINT *)
  60. ; VAR
  61. ;   x1 : real;
  62. ;   temp1, temp2, HALF_ADDER : INTEGER;
  63.  
  64. RANDOM:
  65.     ENTR    D,2,10
  66. ; (* Take 1/2 of the seeds for the comparison test *)
  67. ;   temp1 := SEED1 DIV 2;
  68.     MOV    L,-1(IY)
  69.     MOV    H,0(IY)
  70.     LXI    D,2
  71.     DIVD    D,0
  72.     MOV    -6(IX),H
  73.     MOV    -7(IX),L
  74. ;   temp2 := SEED2 DIV 2;
  75.     MOV    L,-3(IY)
  76.     MOV    H,-2(IY)
  77.     LXI    D,2
  78.     DIVD    D,0
  79.     MOV    -8(IX),H
  80.     MOV    -9(IX),L
  81. ;   IF (temp1+temp2) >= HALFINT then{the number is too big -}
  82.     MOV    L,-7(IX)
  83.     MOV    H,-6(IX)
  84.     MOV    E,-9(IX)
  85.     MOV    D,-8(IX)
  86.     DADD    D,0
  87.     LXI    D,16383
  88.     GE    D,0
  89. ;     { scale it down }
  90. ;     HALF_ADDER := temp1 + temp2 - HALFINT
  91.     JNC    L177
  92.     MOV    L,-7(IX)
  93.     MOV    H,-6(IX)
  94.     MOV    E,-9(IX)
  95.     MOV    D,-8(IX)
  96.     DADD    D,0
  97. ;   ELSE
  98.     LXI    D,-16383
  99.     DADD    D,0
  100.     MOV    -4(IX),H
  101.     MOV    -5(IX),L
  102. ;     HALF_ADDER := temp1 + temp2;
  103.     JMP    L197
  104. L177
  105.     MOV    L,-7(IX)
  106.     MOV    H,-6(IX)
  107.     MOV    E,-9(IX)
  108.     MOV    D,-8(IX)
  109.     DADD    D,0
  110.     MOV    -4(IX),H
  111.     MOV    -5(IX),L
  112. L197
  113. ;   SEED1 := SEED2;
  114.     MOV    L,-3(IY)
  115.     MOV    H,-2(IY)
  116.     MOV    0(IY),H
  117.     MOV    -1(IY),L
  118. ;   (* Restore from previous DIVision *)
  119. ;   SEED2 := HALF_ADDER * 2;
  120.     MOV    L,-5(IX)
  121.     MOV    H,-4(IX)
  122.     DADD    C
  123.     MOV    -2(IY),H
  124.     MOV    -3(IY),L
  125. ;   (*---Convert X to real and divide by factor---*)
  126. ;   x1 := ((X*1.0)/factor);
  127.     MOV    L,8(IX)
  128.     MOV    H,9(IX)
  129.     PUSH    H
  130.     LXI    H,320
  131.     MOV    D,A
  132.     MOV    E,A
  133.     PUSH    H
  134.     PUSH    D
  135.     CVTF    C
  136.     MULT    D,-4
  137.     CVTF    A,32767
  138.     FDVD    D,-4
  139.     LXI    H,3
  140.     DADD    S
  141.     XCHG
  142.     PUSH    IX
  143.     POP    H
  144.     XCHG
  145.     LXI    B,4
  146.     LDDR
  147.     POP    H
  148.     POP    H
  149. ;   (*---Return random number scaled by factor---*)
  150. ;   RANDOM := ( SEED2 * x1 );
  151.     MOV    L,-3(IY)
  152.     MOV    H,-2(IY)
  153.     PUSH    H
  154.     LXI    H,-4
  155.     DADD    S
  156.     SPHL
  157.     XCHG
  158.     PUSH    IX
  159.     POP    H
  160.     DCX    H
  161.     DCX    H
  162.     DCX    H
  163.     LXI    B,4
  164.     LDIR
  165.     CVTF    C
  166.     MULT    D,-4
  167.     LXI    H,3
  168.     DADD    S
  169.     XCHG
  170.     PUSH    IX
  171.     POP    H
  172.     LXI    B,13
  173.     DADD    B
  174.     XCHG
  175.     LXI    B,4
  176.     LDDR
  177.     POP    H
  178.     POP    H
  179. ; End{ of RANDOM(X) };
  180.     EXIT    D,2
  181.