- Timestamp:
- 11/19/09 11:29:41 (14 years ago)
- Branches:
- master
- Children:
- 38d77eb
- Parents:
- 6650a56
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
twister.f
r6650a56 r32289cd 47 47 * May 03, 2000 48 48 * 49 * compile with: f90 -i32 -dp -o mt19937 mt19937.f 49 * compile with: f90 -i32 -dp -o mt19937 mt19937.f 50 50 * Changes: 51 51 * Do not use UMASK as parameter due to restrictions of CF90 … … 55 55 * 56 56 * Now the T3E produces the same random number as a RS6000 running AIX 57 * 57 * 58 58 * from other programs sgrnd() must be called with integer*4 59 59 * … … 64 64 subroutine sgrnd(seed) 65 65 * 66 implicit integer(a-z) 66 integer seed, tmaskb, tmaskc, tshftu, tshfts, tshftt 67 integer tshftl, y 68 69 integer N, N1, mti, mt, m, mata, lmask, mag01 70 integer kk 67 71 * 68 72 * Period parameters … … 88 92 89 93 block data twbloks 94 integer n,n1,mt,mti 90 95 parameter(N = 624) 91 96 parameter(N1 = N+1) … … 99 104 double precision function grnd() 100 105 * 101 implicit integer(a-z) 106 integer tmaskb, tmaskc, tshftu, tshfts, tshftt, tshftl, y 107 108 integer n, n1, mti, mt, m, mata, lmask, mag01, kk 109 102 110 * 103 111 * Period parameters
Note:
See TracChangeset
for help on using the changeset viewer.