home *** CD-ROM | disk | FTP | other *** search
- ; Program: ULTRA.ASM (80386 IBM Fortran/2 version)
- ;
- ; Description: The greatest random number generator that ever was
- ; or ever will be. Way beyond Super-Duper.
- ; (Just kidding, but we think its a good one.)
- ; Authors: Arif Zaman (arif@stat.fsu.edu) and
- ; George Marsaglia (geo@stat.fsu.edu).
- ; Date: 26 March 1992
- ; Version: 1.01
- ; Copyright: To obtain permission to incorporate this program into
- ; any commercial product, please contact the authors at
- ; the e-mail address given above or at
- ;
- ; Department of Statistics and
- ; Supercomputer Computations Research Institute
- ; Florida State University
- ; Tallahassee, FL 32306.
- ;
- ; See Also: The file ULTRA.DOC contains detailed comments
- ;
- DOSSEG
- .MODEL huge,fortran
-
- ;===== B: GLOBAL FUNCTIONS, SUBROUTINES and VARIABLES =================
- ;
- public i31bit, i15bit, i7bit, i1bit, uni, duni,
- public i32bit, i16bit, i8bit, rinit, vni, dvni
- public swbfill,swbx, swbn
-
- ;===== D. MACRO DEFINITIONS ===========================================
- ;
- ; RinitProcStart should take two 32-bit arguments conx and shrx
- ; and place them in eax and ebx.
- ; es and ds should both point to the data segment.
- ; conx must be odd (so we or with 1 just to make sure).
- ; FillProc
- ; DS is already the data segment. ES should also be made the same.
-
- EnterProcedure macro
- mov dx,seg x
- push ds
- mov ds,dx
- endm
-
- ExitProcedure macro
- pop ds
- ret
- endm
-
- RinitProcStart macro ; fortran has call by reference.
- rinit proc conx:dword, shrx:dword
- EnterProcedure
- mov ax,ds
- mov es,ax
- lds si,ShrX
- lodsd ; load shrx in ebx
- mov ebx,eax
- lds si,ConX ; and conx in eax
- lodsd
- or al,1 ; ( congx must be odd )
- endm
-
- RinitProcEnd macro
- ExitProcedure
- rinit endp
- endm
-
- FillProcStart macro
- fillswb proc near
- mov ax,ds
- mov es,ax
- endm
-
- FillProcEnd macro
- ret
- fillswb endp
- endm
-
- DwordFn macro
- endm
-
- WordFn macro
- endm
-
- ByteFn macro
- cbw ; fortran doesn't have bytes so return integer*2
- endm
-
- RealFn macro
- endm
-
- DoubleFn macro
- endm
-
- DATA SEGMENT WORD PUBLIC 'F@DATA'
- INCLUDE ULT32DAT.INC
- DATA ENDS
-
- F@ICODE SEGMENT WORD PUBLIC 'CODE'
- ASSUME CS:F@ICODE, DS:DATA
- INCLUDE ULT32COD.INC
- F@ICODE ENDS
- END