Rev 4330 | Rev 7193 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed
Rev 4330 | Rev 7192 | ||
---|---|---|---|
Line 3... | Line 3... | ||
3 | !endif |
3 | !endif |
4 | 4 | ||
5 | slib_parms=2\ |
5 | slib_parms=2\ |
6 | ,[x_1,...,x_n] \ |
6 | ,[x_1,...,x_n] \ |
7 | ,[a_1,...,a_m] |
7 | ,[a_1,...,a_m] |
8 | slib_author=Adeline, Grelot; Bernadette, |
8 | slib_author=Adeline, Grelot; Bernadette, Perrin-Riou |
9 | 9 | ||
10 | slib_example=[1,1,1,3,4,5,6,2],[ |
10 | slib_example=[1,1,1,3,4,5,6,2],[1,4]\ |
11 | [ |
11 | [1,1,1,3,4,5,6,2],[1,4,6]\ |
12 | [ |
12 | [1,1,1,3,4,5,6,2,7,4,7],[2,3,6]\ |
13 | [1,1,1,3,4,5,6,2] |
13 | [1,1,1,3,4,5,6,2]\ |
- | 14 | [1,1,2] |
|
14 | 15 | ||
15 | !exit |
16 | !exit |
16 | :proc |
17 | :proc |
17 | slib_out= |
18 | slib_out= |
18 | 19 | ||
19 | !readproc slib/stat/dataproc $wims_read_parm |
20 | !readproc slib/stat/dataproc $wims_read_parm |
20 | 21 | ||
21 | slib_v=!declosing $slib_data |
22 | slib_v=!declosing $slib_data |
22 | slib_w=!declosing $slib_weight |
23 | slib_w=!declosing $slib_weight |
23 | !if $slib_v$slib_w=$empty |
24 | !if $slib_v$slib_w=$empty |
24 |
|
25 | slib_out=0 |
25 |
|
26 | !exit |
26 | !else |
27 | !else |
27 |
|
28 | slib_v=!sort numeric item $slib_v |
28 |
|
29 | slib_max=!item -1 of $slib_v |
29 |
|
30 | slib_min=!item 1 of $slib_v |
30 |
|
31 | slib_cnt= !itemcnt $slib_v |
31 |
|
32 | slib_cntw=!itemcnt $slib_w |
32 |
|
33 | !if $slib_cntw=2 |
33 |
|
34 | !distribute item $slib_w into slib_b,slib_p |
34 |
|
35 | slib_c=$[ceil(($slib_max-$slib_b)/$slib_p)] |
35 |
|
36 | slib_d=$[min(0,floor(($slib_min-($slib_b))/$slib_p))] |
36 |
|
37 | slib_w=!values $slib_b+x*$slib_p for x=$slib_d to $slib_c |
37 |
|
38 | !endif |
38 |
|
39 | slib_w=!sort numeric item $slib_w |
39 |
|
40 | !if $slib_w=$empty |
40 |
|
41 | slib_N=$[ceil(4*($slib_cnt)^(1/4))] |
41 |
|
42 | slib_w=!values $slib_min+x*$[(($slib_max-($slib_min))/$slib_N)] for x=0 to $slib_N |
42 |
|
43 | !else |
43 |
|
44 | slib_m=!item -1 of $slib_w |
44 |
|
45 | slib_max=$[max($slib_max,$slib_m)] |
45 |
|
46 | slib_m=!item 1 of $slib_w |
46 |
|
47 | slib_min=$[min($slib_min,$slib_m)] |
47 |
|
48 | slib_w=$slib_min,$slib_w,$slib_max |
48 |
|
49 | slib_w= !listuniq $slib_w |
49 |
|
50 | !endif |
50 |
|
51 | slib_cpt=0 |
51 |
|
52 | slib_cntw=!itemcnt $slib_w |
52 |
|
53 | slib_wpointer=2 |
53 | # end est le point de separation |
54 | # end est le point de separation |
54 |
|
55 | slib_end=!item $slib_wpointer of $slib_w |
55 |
|
56 | slib_result= |
56 |
|
57 | !for slib_i=1 to $slib_cnt |
57 |
|
58 | slib_dat=!item $slib_i of $slib_v |
58 |
|
59 | !if $slib_dat<$slib_end or ($slib_dat=$slib_max and $slib_wpointer=$slib_cntw) |
59 | # On est toujours dans le meme intervalle |
60 | # On est toujours dans le meme intervalle |
60 |
|
61 | !advance slib_cpt |
61 |
|
62 | !else |
62 | # Le point de separation franchi |
63 | # Le point de separation franchi |
63 | # creer un nouveau point de separation |
64 | # creer un nouveau point de separation |
64 | :point |
65 | :point |
65 |
|
66 | !advance slib_wpointer |
66 |
|
67 | slib_end=!item $slib_wpointer of $slib_w |
67 |
|
68 | !if $slib_wpointer < $slib_cntw |
68 |
|
69 | slib_end1=!item $[$slib_wpointer+1] of $slib_w |
69 |
|
70 | !else |
70 |
|
71 | slib_end1=$[$slib_end+2] |
71 |
|
72 | !endif |
72 |
|
73 | slib_result=!append item $[$slib_cpt] to $slib_result |
73 |
|
74 | !ifval $slib_dat>=$slib_end |
74 |
|
75 | slib_cpt=0 |
75 |
|
76 | !goto point |
76 |
|
77 | !else |
77 |
|
78 | slib_cpt=1 |
78 |
|
79 | !endif |
79 |
|
80 | !endif |
80 |
|
81 | !next slib_i |
81 | #Le dernier intervalle |
82 | #Le dernier intervalle |
82 |
|
83 | slib_result=!append item $[$slib_cpt] to $slib_result |
83 |
|
84 | !if $[$slib_wpointer+1] <=$slib_cntw |
84 |
|
85 | slib_zero=!values 0 for x=$[$slib_wpointer+1] to $slib_cntw |
85 |
|
86 | slib_result=$slib_result,$slib_zero |
86 |
|
87 | !endif |
87 | !endif |
88 | !endif |
88 | slib_out=!trim $slib_result |
89 | slib_out=!trim $slib_result |
89 | slib_out=!line -1 of $slib_out |
90 | slib_out=!line -1 of $slib_out |
90 | slib_out=[$slib_out],[$slib_w] |
91 | slib_out=[$slib_out],[$slib_w] |