Subversion Repositories wimsdev

Rev

Rev 4349 | 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
7692 bpr 8
 
20 reyssat 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
7692 bpr 48
	!if $slib_i>1
20 reyssat 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