Subversion Repositories wimsdev

Rev

Rev 8775 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 8775 Rev 13586
Line 1... Line 1...
1 
!if $wims_read_parm!=slib_header
 1 
!if $wims_read_parm!=slib_header
2 
 !goto proc
 2 
  !goto proc
3 
!endif
 3 
!endif
4 
 
 4 
 
5 
 
 5 
 
6 
slib_parms=3\
 6 
slib_parms=3\
7 
,m\
 7 
,m\
Line 26... Line 26...
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
27 
slib_w2=!values sin(2*Pi*x) for x in $slib_U2
 27 
slib_w2=!values sin(2*Pi*x) for x in $slib_U2
28 
 
 28 
 
29 
slib_V=$empty
 29 
slib_V=$empty
30 
!for slib_i=1 to $[$slib_m2 +1]
 30 
!for slib_i=1 to $[$slib_m2 +1]
31 
slib_v=!item  $slib_i of $slib_w
 31 
  slib_v=!item  $slib_i of $slib_w
32 
slib_v1=!item  $slib_i of $slib_w1
 32 
  slib_v1=!item  $slib_i of $slib_w1
33 
slib_v2=!item  $slib_i of $slib_w2
 33 
  slib_v2=!item  $slib_i of $slib_w2
34 
slib_v1=$[$slib_sigma*$slib_v1*$slib_v+$slib_mu]
 34 
  slib_v1=$[$slib_sigma*$slib_v1*$slib_v+$slib_mu]
35 
slib_v2=$[$slib_sigma*$slib_v2*$slib_v+$slib_mu]
 35 
  slib_v2=$[$slib_sigma*$slib_v2*$slib_v+$slib_mu]
36 
 
 36 
 
37 
!if  $slib_i<=$slib_m2
 37 
  !if  $slib_i<=$slib_m2
38 
slib_V=!append item $slib_v1 to $slib_V
 38 
    slib_V=!append item $slib_v1 to $slib_V
39 
slib_V=!append item $slib_v2 to $slib_V
 39 
    slib_V=!append item $slib_v2 to $slib_V
40 
!endif
 40 
  !endif
41 
!if $slib_i>$slib_m2 and   $[2*$slib_m2]!=$slib_m
 41 
  !if $slib_i>$slib_m2 and   $[2*$slib_m2]!=$slib_m
42 
slib_V=!append item $slib_v1 to $slib_V
 42 
    slib_V=!append item $slib_v1 to $slib_V
43 
!endif
 43 
  !endif
44 
!next slib_i
 44 
!next slib_i
45 
 
 45 
 
46 
!if $slib_m >1
 46 
!if $slib_m >1
47 
slib_V=!shuffle $slib_V
 47 
  slib_V=!shuffle $slib_V
48 
!endif
 48 
!endif
49 
 
 49 
 
50 
slib_out=!trim $slib_V
 50 
slib_out=!trim $slib_V