Rev 4158 | Rev 7692 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
20 | reyssat | 1 | !if $wims_read_parm!=slib_header |
2 | !goto proc |
||
3 | !endif |
||
4158 | bpr | 4 | |
4349 | bpr | 5 | slib_author=Bernadette, Perrin-Riou |
4158 | bpr | 6 | |
7 | slib_example=x^2+sin(x*y)+exp(y^2),[x,y],0,3,-3,2,50 |
||
20 | reyssat | 8 | |
9 | !exit |
||
10 | |||
11 | :proc |
||
12 | !distribute items $wims_read_parm into slib_fn,slib_var,slib_x1,slib_x2,slib_y1,slib_y2,slib_steps |
||
13 | slib_fn=!rawmath $slib_fn |
||
14 | slib_var=!declosing $slib_var |
||
15 | slib_var1=!item 1 of $slib_var |
||
16 | slib_var2=!item 2 of $slib_var |
||
17 | !default slib_var1=x |
||
18 | !default slib_var2=y |
||
19 | slib_x1=$[$slib_x1] |
||
20 | !if $slib_x1=$empty or $slib_x1 isitemof NaN,Inf |
||
21 | slib_x1=0 |
||
22 | !endif |
||
23 | slib_x2=$[$slib_x2] |
||
24 | !if $slib_x2=$empty or $slib_x2 isitemof NaN,Inf |
||
25 | slib_x2=1 |
||
26 | !endif |
||
27 | slib_y1=$[$slib_y1] |
||
28 | !if $slib_y1=$empty or $slib_y1 isitemof NaN,Inf |
||
29 | slib_y1=0 |
||
30 | !endif |
||
31 | slib_y2=$[$slib_y2] |
||
32 | !if $slib_y2=$empty or $slib_y2 isitemof NaN,Inf |
||
33 | slib_x2=1 |
||
34 | !endif |
||
35 | !bound slib_steps between integer 1 and 1000 default 50 |
||
36 | !bound slib_steps2 between integer 1 and 1000 default 50 |
||
37 | slib_step=$[($slib_x2-($slib_x1))/$slib_steps] |
||
38 | slib_step2=$[($slib_y2-($slib_y1))/$slib_steps] |
||
39 | slib_min=$empty |
||
40 | slib_max=$empty |
||
41 | |||
42 | !for slib_i=0 to $slib_steps |
||
43 | slib_fxn=!mathsubst $slib_var1= $[$slib_x1+$slib_i*$slib_step] in $slib_fn |
||
44 | slib_fxn=!replace item NaN by $ in $slib_fxn |
||
45 | slib_vals=!values $slib_fxn for $slib_var2=$slib_y1 to $slib_y2 step $slib_step |
||
46 | slib_vals=!replace item NaN by $ in $slib_vals |
||
47 | slib_vals=!nonempty items $slib_vals |
||
48 | !if $slib_i>1 |
||
49 | slib_vals=!sort numeric items $slib_min,$slib_vals,$slib_max |
||
50 | !else |
||
51 | slib_vals=!sort numeric items $slib_vals |
||
52 | !endif |
||
53 | slib_vals=!trim $slib_vals |
||
54 | slib_min=!item 1 of $slib_vals |
||
55 | slib_max=!item -1 of $slib_vals |
||
56 | !next slib_i |
||
57 | |||
58 | slib_out=$slib_min,$slib_max |