First commit : 0.14.0 version (with roadmap in doc instead of
[cloog/uuh.git] / test / published / Web / web5.cloog
blob381339b298f83d20408911cb4b233bb9ed69eabe
1 # CLooG example file #5.
2 # Please read the first and second examples which are fully documented to
3 # understand the different parts of the input file.
5 ################################################################################
6 #  do i=1,n                            The problem here is to regenerate a     #
7 #  | do j =1,i-1                       real-life Cholesau kernel according to  #
8 #  | | if (j.EQ.1) then                the original scheduling (see the user's #
9 #S1| | | s1(i,j)=a(i,j)s4(j,i)**2      manual for more details). The original  #
10 #  | | else                            program is given on the left. For each  #
11 #S2| | | s1(i,j)=s1(i,j-1)-s4(j,i)**2  statement the original schedule is:     #
12 #  | if (i .EQ. 1) then                        T_S1(i,j)  =(i,0,j,0,0,0)       #
13 #S3| | s2(i)=sqrt(a(i,i))                      T_S2(i,j)  =(i,0,j,1,0,0)       #
14 #  | else                                      T_S3(i)    =(i,1,0,0,0,0)       #
15 #S4| | s2(i)=sqrt (s1(i,i-1))                  T_S4(i)    =(i,2,0,0,0,0)       #
16 #  | do k=i+1,n                                T_S5(i,j,k)=(i,3,j,0,k,0)       #
17 #  | | do l=1,i-1                              T_S6(i,j,k)=(i,3,j,0,k,1)       #
18 #  | | | if (l .EQ. 1) then                    T_S7(i,j)  =(i,3,j,1,0,0)       #
19 #S5| | | | s3(i,k,l)=a(k,i)-(s4(l,k)*s4(l,i))  T_S8(i,j)  =(i,3,j,2,0,0)       #
20 #  | | | else                                                                  #
21 #S6| | | | s3(i,k,l)=s3(i,k,l-1)-(s4(l,k)*s4(l,i))                             #
22 #  | | if (i .EQ.1) then                                                       #
23 #S7| | | s4(i,k)=a(k,i)/s2(i)          Note that in the generated code there   #
24 #  | | else                            are no more conditions.                 #
25 #S8| | | s4(i,k)=s3(i,k,i-1)/s2(i)                                             #
26 ################################################################################
28 #------------------------------------CONTEXT------------------------------------
30 # 1. language: FORTRAN
33 # 2. Parameters {n | n>=10}
34 1 3
35 #  n  1
36 1  1 -10 # n>=10
38 # 3. We set manually the parameter name: n
42 #-----------------------------------POLYHEDRA-----------------------------------
44 # 4. Number of polyhedra:
47 # Polyhedron #1
49 # {i, j | 1<=i<=n; 1<=j<=i-1; j=1}
50 5 5
51 #  i  j  n  1
52 1  1  0  0 -1 # 1<=i
53 1 -1  0  1  0 # i<=n
54 1  0  1  0 -1 # 1<=j
55 1  1 -1  0 -1 # j<=i-1
56 0  0  1  0 -1 # j=1
57 0  0  0       # 3 zeroes !
59 # Polyhedron #2
61 # {i, j | 1<=i<=n; 1<=j<=i-1; j!=1}
62 5 5
63 #  i  j  n  1
64 1  1  0  0 -1 # 1<=i
65 1 -1  0  1  0 # i<=n
66 1  0  1  0 -1 # 1<=j
67 1  1 -1  0 -1 # j<=i-1
68 1  0  1  0 -2 # j>=2
69 5 5
70 #  i  j  n  1
71 1  1  0  0 -1 # 1<=i
72 1 -1  0  1  0 # i<=n
73 1  0  1  0 -1 # 1<=j
74 1  1 -1  0 -1 # j<=i-1
75 1  0 -1  0  0 # j<=0
76 0  0  0       # 3 zeroes !
78 # Polyhedron #3
80 # {i | 1<=i<=n; i=1}
81 3 4
82 #  i  n  1
83 1  1  0 -1 # 1<=i
84 1 -1  1  0 # i<=n
85 0  1  0 -1 # i=1
86 0  0  0    # 3 zeroes !
88 # Polyhedron #4
90 # {i | 1<=i<=n; i!=1}
91 3 4
92 #  i  n  1
93 1  1  0 -1 # 1<=i
94 1 -1  1  0 # i<=n
95 1  1  0 -2 # i>=2
96 3 4
97 #  i  n  1
98 1  1  0 -1 # 1<=i
99 1 -1  1  0 # i<=n
100 1 -1  0  0 # i<=0
101 0  0  0    # 3 zeroes !
103 # Polyhedron #5
105 # {i, j | 1<=i<=n; i+1<=j<=n; 1<=k<=i-1; k=1}
106 7 6
107 #  i  j  k  n  1
108 1  1  0  0  0 -1 # 1<=i
109 1 -1  0  0  1  0 # i<=n
110 1 -1  1  0  0 -1 # i+1<=j
111 1  0 -1  0  1  0 # j<=n
112 1  0  0  1  0 -1 # 1<=k
113 1  1  0 -1  0 -1 # k<=i-1
114 0  0  0  1  0 -1 # k=1
115 0  0  0          # 3 zeroes !
117 # Polyhedron #6
119 # {i, j | 1<=i<=n; i+1<=j<=n; 1<=k<=i-1; k!=1}
120 7 6
121 #  i  j  k  n  1
122 1  1  0  0  0 -1 # 1<=i
123 1 -1  0  0  1  0 # i<=n
124 1 -1  1  0  0 -1 # i+1<=j
125 1  0 -1  0  1  0 # j<=n
126 1  0  0  1  0 -1 # 1<=k
127 1  1  0 -1  0 -1 # k<=i-1
128 1  0  0  1  0 -2 # k>=2
129 7 6
130 #  i  j  k  n  1
131 1  1  0  0  0 -1 # 1<=i
132 1 -1  0  0  1  0 # i<=n
133 1 -1  1  0  0 -1 # i+1<=j
134 1  0 -1  0  1  0 # j<=n
135 1  0  0  1  0 -1 # 1<=k
136 1  1  0 -1  0 -1 # k<=i-1
137 1  0  0 -1  0  0 # k<=0
138 0  0  0          # 3 zeroes !
140 # Polyhedron #7
142 # {i, j | 1<=i<=n; i+1<=j<=n; i=1}
143 5 5
144 #  i  j  n  1
145 1  1  0  0 -1 # 1<=i
146 1 -1  0  1  0 # i<=n
147 1 -1  1  0 -1 # i+1<=j
148 1  0 -1  1  0 # j<=n
149 0  1  0  0 -1 # i=1
150 0  0  0       # 3 zeroes !
152 # Polyhedron #8
154 # {i, j | 1<=i<=n; i+1<=j<=n; i!=1}
155 5 5
156 #  i  j  n  1
157 1  1  0  0 -1 # 1<=i
158 1 -1  0  1  0 # i<=n
159 1 -1  1  0 -1 # i+1<=j
160 1  0 -1  1  0 # j<=n
161 1  1  0  0 -2 # i>=2
162 5 5
163 #  i  j  n  1
164 1  1  0  0 -1 # 1<=i
165 1 -1  0  1  0 # i<=n
166 1 -1  1  0 -1 # i+1<=j
167 1  0 -1  1  0 # j<=n
168 1 -1  0  0  0 # i<=0
169 0  0  0       # 3 zeroes !
171 # 6. We let CLooG choose the iterator names
174 #----------------------------------SCATTERING-----------------------------------
176 # 7. Scattering functions ORIGINAL SCHEDULING
179 # Scattering function for polyhedron #1: T_S1(i,j)  =(i,0,j,0,0,0)
180 6 11
181 # c1 c2 c3 c4 c5 c6  i  j  n  1
182 0  1  0  0  0  0  0 -1  0  0  0 # i
183 0  0  1  0  0  0  0  0  0  0  0 # 0
184 0  0  0  1  0  0  0  0 -1  0  0 # j
185 0  0  0  0  1  0  0  0  0  0  0 # 0
186 0  0  0  0  0  1  0  0  0  0  0 # 0
187 0  0  0  0  0  0  1  0  0  0  0 # 0
189 # Scattering function for polyhedron #2: T_S2(i,j)  =(i,0,j,1,0,0)
190 6 11
191 # c1 c2 c3 c4 c5 c6  i  j  n  1
192 0  1  0  0  0  0  0 -1  0  0  0 # i
193 0  0  1  0  0  0  0  0  0  0  0 # 0
194 0  0  0  1  0  0  0  0 -1  0  0 # j
195 0  0  0  0  1  0  0  0  0  0 -1 # 1
196 0  0  0  0  0  1  0  0  0  0  0 # 0
197 0  0  0  0  0  0  1  0  0  0  0 # 0
199 # Scattering function for polyhedron #3: T_S3(i)    =(i,1,0,0,0,0)
200 6 10
201 # c1 c2 c3 c4 c5 c6  i  n  1
202 0  1  0  0  0  0  0 -1  0  0 # i
203 0  0  1  0  0  0  0  0  0 -1 # 1
204 0  0  0  1  0  0  0  0  0  0 # 0
205 0  0  0  0  1  0  0  0  0  0 # 0
206 0  0  0  0  0  1  0  0  0  0 # 0
207 0  0  0  0  0  0  1  0  0  0 # 0
209 # Scattering function for polyhedron #4: T_S4(i)    =(i,2,0,0,0,0)
210 6 10
211 # c1 c2 c3 c4 c5 c6  i  n  1
212 0  1  0  0  0  0  0 -1  0  0 # i
213 0  0  1  0  0  0  0  0  0 -2 # 2
214 0  0  0  1  0  0  0  0  0  0 # 0
215 0  0  0  0  1  0  0  0  0  0 # 0
216 0  0  0  0  0  1  0  0  0  0 # 0
217 0  0  0  0  0  0  1  0  0  0 # 0
219 # Scattering function for polyhedron #5: T_S5(i,j,k)=(i,3,j,0,k,0)
220 6 12
221 # c1 c2 c3 c4 c5 c6  i  j  k  n  1
222 0  1  0  0  0  0  0 -1  0  0  0  0 # i
223 0  0  1  0  0  0  0  0  0  0  0 -3 # 3
224 0  0  0  1  0  0  0  0 -1  0  0  0 # j
225 0  0  0  0  1  0  0  0  0  0  0  0 # 0
226 0  0  0  0  0  1  0  0  0 -1  0  0 # k
227 0  0  0  0  0  0  1  0  0  0  0  0 # 0
229 # Scattering function for polyhedron #6: T_S6(i,j,k)=(i,3,j,0,k,1)
230 6 12
231 # c1 c2 c3 c4 c5 c6  i  j  k  n  1
232 0  1  0  0  0  0  0 -1  0  0  0  0 # i
233 0  0  1  0  0  0  0  0  0  0  0 -3 # 3
234 0  0  0  1  0  0  0  0 -1  0  0  0 # j
235 0  0  0  0  1  0  0  0  0  0  0  0 # 0
236 0  0  0  0  0  1  0  0  0 -1  0  0 # k
237 0  0  0  0  0  0  1  0  0  0  0 -1 # 1
239 # Scattering function for polyhedron #7: T_S7(i,j)  =(i,3,j,1,0,0)
240 6 11
241 # c1 c2 c3 c4 c5 c6  i  j  n  1
242 0  1  0  0  0  0  0 -1  0  0  0 # i
243 0  0  1  0  0  0  0  0  0  0 -3 # 3
244 0  0  0  1  0  0  0  0 -1  0  0 # j
245 0  0  0  0  1  0  0  0  0  0 -1 # 1
246 0  0  0  0  0  1  0  0  0  0  0 # 0
247 0  0  0  0  0  0  1  0  0  0  0 # 0
249 # Scattering function for polyhedron #8: T_S8(i,j)  =(i,3,j,2,0,0)
250 6 11
251 # c1 c2 c3 c4 c5 c6  i  j  n  1
252 0  1  0  0  0  0  0 -1  0  0  0 # i
253 0  0  1  0  0  0  0  0  0  0 -3 # 3
254 0  0  0  1  0  0  0  0 -1  0  0 # j
255 0  0  0  0  1  0  0  0  0  0 -2 # 2
256 0  0  0  0  0  1  0  0  0  0  0 # 0
257 0  0  0  0  0  0  1  0  0  0  0 # 0
259 # We want to set manually the scattering dimension names.
261 c1 c2 c3 c4 c5 c6