Rev 7692 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 7692 | 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 | slib_parms=2\ |
5 | slib_parms=2\ |
6 | ,[data]\ |
6 | ,[data]\ |
7 | 1,p |
7 | 1,p |
Line 16... | Line 16... | ||
16 | 16 | ||
17 | !distribute item $wims_read_parm into slib_data,slib_p |
17 | !distribute item $wims_read_parm into slib_data,slib_p |
18 | 18 | ||
19 | !default slib_N=1000 |
19 | !default slib_N=1000 |
20 | !if $slib_p<=0 |
20 | !if $slib_p<=0 |
21 | slib_out=0 |
21 | slib_out=0 |
22 | !exit |
22 | !exit |
23 | !endif |
23 | !endif |
24 | 24 | ||
25 | slib_data=!declosing $slib_data |
25 | slib_data=!declosing $slib_data |
26 | slib_out= |
26 | slib_out= |
27 | !for slib_x in $slib_data |
27 | !for slib_x in $slib_data |
28 | !if $slib_x<0 |
28 | !if $slib_x<0 |
29 | slib_out=!append item 0 to $slib_out |
29 | slib_out=!append item 0 to $slib_out |
30 | !else |
30 | !else |
31 | slib_out=!append item $[$slib_p *exp(-$slib_p*$slib_x)] to $slib_out |
31 | slib_out=!append item $[$slib_p *exp(-$slib_p*$slib_x)] to $slib_out |
32 | !endif |
32 | !endif |
33 | !next slib_x |
33 | !next slib_x |
34 | slib_out=!trim $slib_out |
34 | slib_out=!trim $slib_out |