Line | |
---|
1 | subroutine loduni(nbase,n,y) |
---|
2 | C----------------------------------------------------------------------- |
---|
3 | C Quiet start |
---|
4 | C Load an uniform distribution using the Hammersley's sequence. |
---|
5 | C (NBASE=0 ==> Random sampling !) |
---|
6 | C----------------------------------------------------------------------- |
---|
7 | dimension y(*) |
---|
8 | C |
---|
9 | C |
---|
10 | C----Random and Quasi-Random Loading |
---|
11 | C |
---|
12 | if(nbase.eq.0) then |
---|
13 | c |
---|
14 | c Random |
---|
15 | do 100 i=1,n |
---|
16 | 100 y(i)=ranf() |
---|
17 | c |
---|
18 | c First element of Hammersley's sequence |
---|
19 | else if(nbase.eq.1) then |
---|
20 | do 110 i=1,n |
---|
21 | 110 y(i) = (i-0.5)/n |
---|
22 | c |
---|
23 | c Radical-inverse Function in base NBASE |
---|
24 | else if(nbase.gt.1) then |
---|
25 | do 120 i=1,n |
---|
26 | xs = 0. |
---|
27 | xsi = 1.0 |
---|
28 | j2 = i |
---|
29 | 1 xsi = xsi/nbase |
---|
30 | j1 = j2/nbase |
---|
31 | xs = xs + (j2-nbase*j1)*xsi |
---|
32 | j2 = j1 |
---|
33 | if( j2.gt.0 ) goto 1 |
---|
34 | 120 y(i) = xs |
---|
35 | end if |
---|
36 | c |
---|
37 | return |
---|
38 | end |
---|
39 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|
Note: See
TracBrowser
for help on using the repository browser.