Rev 7692 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 7692 | Rev 8775 | ||
---|---|---|---|
Line 16... | Line 16... | ||
16 | :proc |
16 | :proc |
17 | !distribute item $wims_read_parm into slib_m,slib_mu,slib_sigma |
17 | !distribute item $wims_read_parm into slib_m,slib_mu,slib_sigma |
18 | !default slib_mu=0 |
18 | !default slib_mu=0 |
19 | !default slib_sigma=1 |
19 | !default slib_sigma=1 |
20 | 20 | ||
21 | slib_m2= $[floor($slib_m/2)] |
21 | slib_m2= $[floor(($slib_m)/2)] |
22 | slib_U1=!random 10^-20,1 repeat $[$slib_m2+1] |
22 | slib_U1=!random 10^-20,1 repeat $[$slib_m2+1] |
23 | slib_U2=!random 0,1 repeat $[$slib_m2+1] |
23 | slib_U2=!random 0,1 repeat $[$slib_m2+1] |
24 | 24 | ||
25 | slib_w=!values sqrt(-2*log(x)) for x in $slib_U1 |
25 | slib_w=!values sqrt(-2*log(x)) for x in $slib_U1 |
26 | slib_w1=!values cos(2*Pi*x) for x in $slib_U2 |
26 | slib_w1=!values cos(2*Pi*x) for x in $slib_U2 |