?  search:
summary | log | graphiclog1 | graphiclog2 | commit | commitdiff | tree | refs | edit | fork
blame | history | raw | HEAD
barvinok_e.cc: directly include required headers
[barvinok.git] / iscc.c
blob00c1c93f714f736ce374d61b63bc3763abd452cb
1 #include <assert.h>
2 #include <ctype.h>
3 #include <stdio.h>
4 #include <string.h>
5 #include <unistd.h>
6 #include <isl/ctx.h>
7 #include <isl/val.h>
8 #include <isl/space.h>
9 #include <isl/aff.h>
10 #include <isl/obj.h>
11 #include <isl/stream.h>
12 #include <isl/set.h>
13 #include <isl/map.h>
14 #include <isl/union_set.h>
15 #include <isl/union_map.h>
16 #include <isl/vertices.h>
17 #include <isl/flow.h>
18 #include <isl/schedule.h>
19 #include <isl/ast_build.h>
20 #include <isl/printer.h>
21 #include <isl_obj_list.h>
22 #include <isl_obj_str.h>
23 #include <barvinok/isl.h>
24 #include <barvinok/options.h>
25 #include "lattice_width.h"
26
27 #include "config.h"
28
29 #ifdef HAVE_SIGACTION
30 #include <signal.h>
31
32 static isl_ctx *main_ctx;
33
34 static void handler(int signum)
35 {
36 if (isl_ctx_aborted(main_ctx))
37 exit(EXIT_FAILURE);
38 isl_ctx_abort(main_ctx);
39 }
40
41 static struct sigaction sa_old;
42
43 static void install_signal_handler(isl_ctx *ctx)
44 {
45 struct sigaction sa;
46
47 main_ctx = ctx;
48
49 memset(&sa, 0, sizeof(struct sigaction));
50 sa.sa_handler = &handler;
51 sa.sa_flags = SA_RESTART;
52 sigaction(SIGINT, &sa, &sa_old);
53 }
54
55 static void remove_signal_handler(isl_ctx *ctx)
56 {
57 sigaction(SIGINT, &sa_old, NULL);
58 }
59
60 #else
61
62 static void install_signal_handler(isl_ctx *ctx)
63 {
64 }
65
66 static void remove_signal_handler(isl_ctx *ctx)
67 {
68 }
69
70 #endif
71
72 #ifdef HAVE_PET
73 #include <pet.h>
74 #else
75 struct pet_options;
76 int pet_options_set_autodetect(isl_ctx *ctx, int val)
77 {
78 return -1;
79 }
80 int pet_options_set_encapsulate_dynamic_control(isl_ctx *ctx, int val)
81 {
82 return -1;
83 }
84 #endif
85
86 static int iscc_bool_false = 0;
87 static int iscc_bool_true = 1;
88 static int iscc_bool_error = -1;
89
90 enum iscc_op { ISCC_READ, ISCC_WRITE, ISCC_SOURCE, ISCC_VERTICES,
91 ISCC_LAST, ISCC_ANY, ISCC_BEFORE, ISCC_UNDER,
92 ISCC_SCHEDULE,
93 ISCC_MINIMIZING, ISCC_RESPECTING,
94 ISCC_CODEGEN, ISCC_USING,
95 ISCC_TYPEOF, ISCC_PRINT, ISCC_ASSERT,
96 ISCC_N_OP };
97 static const char *op_name[ISCC_N_OP] = {
98 [ISCC_ASSERT] = "assert",
99 [ISCC_READ] = "read",
100 [ISCC_WRITE] = "write",
101 [ISCC_PRINT] = "print",
102 [ISCC_SOURCE] = "source",
103 [ISCC_VERTICES] = "vertices",
104 [ISCC_LAST] = "last",
105 [ISCC_ANY] = "any",
106 [ISCC_BEFORE] = "before",
107 [ISCC_UNDER] = "under",
108 [ISCC_SCHEDULE] = "schedule",
109 [ISCC_MINIMIZING] = "minimizing",
110 [ISCC_RESPECTING] = "respecting",
111 [ISCC_CODEGEN] = "codegen",
112 [ISCC_USING] = "using",
113 [ISCC_TYPEOF] = "typeof"
114 };
115 static enum isl_token_type iscc_op[ISCC_N_OP];
116
117 struct isl_arg_choice iscc_format[] = {
118 {"isl", ISL_FORMAT_ISL},
119 {"omega", ISL_FORMAT_OMEGA},
120 {"polylib", ISL_FORMAT_POLYLIB},
121 {"ext-polylib", ISL_FORMAT_EXT_POLYLIB},
122 {"latex", ISL_FORMAT_LATEX},
123 {"C", ISL_FORMAT_C},
124 {0}
125 };
126
127 struct iscc_options {
128 struct barvinok_options *barvinok;
129 struct pet_options *pet;
130 unsigned format;
131 int io;
132 };
133
134 ISL_ARGS_START(struct iscc_options, iscc_options_args)
135 ISL_ARG_CHILD(struct iscc_options, barvinok, "barvinok", &barvinok_options_args,
136 "barvinok options")
137 #ifdef HAVE_PET
138 ISL_ARG_CHILD(struct iscc_options, pet, "pet", &pet_options_args, "pet options")
139 #endif
140 ISL_ARG_CHOICE(struct iscc_options, format, 0, "format", \
141 iscc_format, ISL_FORMAT_ISL, "output format")
142 ISL_ARG_BOOL(struct iscc_options, io, 0, "io", 1,
143 "allow read and write operations")
144 ISL_ARGS_END
145
146 ISL_ARG_DEF(iscc_options, struct iscc_options, iscc_options_args)
147 ISL_ARG_CTX_DEF(iscc_options, struct iscc_options, iscc_options_args)
148
149 static void *isl_obj_bool_copy(void *v)
150 {
151 return v;
152 }
153
154 static void isl_obj_bool_free(void *v)
155 {
156 }
157
158 static __isl_give isl_printer *isl_obj_bool_print(__isl_take isl_printer *p,
159 void *v)
160 {
161 if (v == &iscc_bool_true)
162 return isl_printer_print_str(p, "True");
163 else if (v == &iscc_bool_false)
164 return isl_printer_print_str(p, "False");
165 else
166 return isl_printer_print_str(p, "Error");
167 }
168
169 static void *isl_obj_bool_add(void *v1, void *v2)
170 {
171 return v1;
172 }
173
174 struct isl_obj_vtable isl_obj_bool_vtable = {
175 isl_obj_bool_copy,
176 isl_obj_bool_add,
177 isl_obj_bool_print,
178 isl_obj_bool_free
179 };
180 #define isl_obj_bool (&isl_obj_bool_vtable)
181
182 int *iscc_bool_from_int(int res)
183 {
184 return res < 0 ? &iscc_bool_error :
185 res ? &iscc_bool_true : &iscc_bool_false;
186 }
187
188 /* Conjunction of "b1" and "b2".
189 * The result is returned as an integer because it is post-processed by
190 * iscc_bool_from_int.
191 */
192 static int isl_bool_and(isl_bool *b1, __isl_take isl_bool *b2)
193 {
194 if (b1 == &iscc_bool_error || b2 == &iscc_bool_error)
195 return -1;
196 return b1 == &iscc_bool_true && b2 == &iscc_bool_true;
197 }
198
199 /* Disjunction of "b1" and "b2".
200 * The result is returned as an integer because it is post-processed by
201 * iscc_bool_from_int.
202 */
203 static int isl_bool_or(isl_bool *b1, __isl_take isl_bool *b2)
204 {
205 if (b1 == &iscc_bool_error || b2 == &iscc_bool_error)
206 return -1;
207 return b1 == &iscc_bool_true || b2 == &iscc_bool_true;
208 }
209
210 static int isl_union_map_is_superset(__isl_take isl_union_map *map1,
211 __isl_take isl_union_map *map2)
212 {
213 return isl_union_map_is_subset(map2, map1);
214 }
215 static int isl_union_set_is_superset(__isl_take isl_union_set *set1,
216 __isl_take isl_union_set *set2)
217 {
218 return isl_union_set_is_subset(set2, set1);
219 }
220
221 static int isl_union_map_is_strict_superset(__isl_take isl_union_map *map1,
222 __isl_take isl_union_map *map2)
223 {
224 return isl_union_map_is_strict_subset(map2, map1);
225 }
226 static int isl_union_set_is_strict_superset(__isl_take isl_union_set *set1,
227 __isl_take isl_union_set *set2)
228 {
229 return isl_union_set_is_strict_subset(set2, set1);
230 }
231
232 extern struct isl_obj_vtable isl_obj_list_vtable;
233 #define isl_obj_list (&isl_obj_list_vtable)
234
235 typedef void *(*isc_bin_op_fn)(void *lhs, void *rhs);
236 typedef int (*isc_bin_test_fn)(void *lhs, void *rhs);
237 struct isc_bin_op {
238 enum isl_token_type op;
239 isl_obj_type lhs;
240 isl_obj_type rhs;
241 isl_obj_type res;
242 union {
243 isc_bin_op_fn fn;
244 isc_bin_test_fn test;
245 } o;
246 };
247 struct isc_named_bin_op {
248 char *name;
249 struct isc_bin_op op;
250 };
251 /* Compound binary operator.
252 * "full" is only used to generate a unique token number for op.op
253 * in register_named_ops.
254 * "op1" is the first part of the compound operator.
255 * "op2" is the second part of the compound operator.
256 */
257 struct iscc_compound_bin_op {
258 char *full;
259 enum isl_token_type op1;
260 enum isl_token_type op2;
261 struct isc_bin_op op;
262 };
263
264 struct iscc_at {
265 isl_union_pw_qpolynomial *upwqp;
266 isl_union_pw_qpolynomial *res;
267 };
268
269 static isl_stat eval_at(__isl_take isl_point *pnt, void *user)
270 {
271 struct iscc_at *at = (struct iscc_at *) user;
272 isl_val *v;
273 isl_qpolynomial *qp;
274 isl_set *set;
275
276 set = isl_set_from_point(isl_point_copy(pnt));
277 v = isl_union_pw_qpolynomial_eval(
278 isl_union_pw_qpolynomial_copy(at->upwqp), pnt);
279 qp = isl_qpolynomial_val_on_domain(isl_set_get_space(set), v);
280
281 at->res = isl_union_pw_qpolynomial_add(at->res,
282 isl_union_pw_qpolynomial_from_pw_qpolynomial(
283 isl_pw_qpolynomial_alloc(set, qp)));
284
285 return isl_stat_ok;
286 }
287
288 __isl_give isl_union_pw_qpolynomial *isl_union_pw_qpolynomial_at(
289 __isl_take isl_union_pw_qpolynomial *upwqp,
290 __isl_take isl_union_set *uset)
291 {
292 struct iscc_at at;
293
294 at.upwqp = upwqp;
295 at.res = isl_union_pw_qpolynomial_zero(isl_union_set_get_space(uset));
296
297 isl_union_set_foreach_point(uset, eval_at, &at);
298
299 isl_union_pw_qpolynomial_free(upwqp);
300 isl_union_set_free(uset);
301
302 return at.res;
303 }
304
305 struct iscc_fold_at {
306 isl_union_pw_qpolynomial_fold *upwf;
307 isl_union_pw_qpolynomial *res;
308 };
309
310 static isl_stat eval_fold_at(__isl_take isl_point *pnt, void *user)
311 {
312 struct iscc_fold_at *at = (struct iscc_fold_at *) user;
313 isl_val *v;
314 isl_qpolynomial *qp;
315 isl_set *set;
316
317 set = isl_set_from_point(isl_point_copy(pnt));
318 v = isl_union_pw_qpolynomial_fold_eval(
319 isl_union_pw_qpolynomial_fold_copy(at->upwf), pnt);
320 qp = isl_qpolynomial_val_on_domain(isl_set_get_space(set), v);
321
322 at->res = isl_union_pw_qpolynomial_add(at->res,
323 isl_union_pw_qpolynomial_from_pw_qpolynomial(
324 isl_pw_qpolynomial_alloc(set, qp)));
325
326 return isl_stat_ok;
327 }
328
329 __isl_give isl_union_pw_qpolynomial *isl_union_pw_qpolynomial_fold_at(
330 __isl_take isl_union_pw_qpolynomial_fold *upwf,
331 __isl_take isl_union_set *uset)
332 {
333 struct iscc_fold_at at;
334
335 at.upwf = upwf;
336 at.res = isl_union_pw_qpolynomial_zero(isl_union_set_get_space(uset));
337
338 isl_union_set_foreach_point(uset, eval_fold_at, &at);
339
340 isl_union_pw_qpolynomial_fold_free(upwf);
341 isl_union_set_free(uset);
342
343 return at.res;
344 }
345
346 static __isl_give isl_union_pw_qpolynomial_fold *union_pw_qpolynomial_add_union_pw_qpolynomial_fold(
347 __isl_take isl_union_pw_qpolynomial *upwqp,
348 __isl_take isl_union_pw_qpolynomial_fold *upwf)
349 {
350 return isl_union_pw_qpolynomial_fold_add_union_pw_qpolynomial(upwf,
351 upwqp);
352 }
353
354 static __isl_give struct isl_list *union_map_apply_union_pw_qpolynomial_fold(
355 __isl_take isl_union_map *umap,
356 __isl_take isl_union_pw_qpolynomial_fold *upwf)
357 {
358 isl_ctx *ctx;
359 struct isl_list *list;
360 int tight;
361
362 ctx = isl_union_map_get_ctx(umap);
363 list = isl_list_alloc(ctx, 2);
364 if (!list)
365 goto error2;
366
367 list->obj[0].type = isl_obj_union_pw_qpolynomial_fold;
368 list->obj[0].v = isl_union_map_apply_union_pw_qpolynomial_fold(umap,
369 upwf, &tight);
370 list->obj[1].type = isl_obj_bool;
371 list->obj[1].v = tight ? &iscc_bool_true : &iscc_bool_false;
372 if (tight < 0 || !list->obj[0].v)
373 goto error;
374
375 return list;
376 error2:
377 isl_union_map_free(umap);
378 isl_union_pw_qpolynomial_fold_free(upwf);
379 error:
380 isl_list_free(list);
381 return NULL;
382 }
383
384 static __isl_give struct isl_list *union_set_apply_union_pw_qpolynomial_fold(
385 __isl_take isl_union_set *uset,
386 __isl_take isl_union_pw_qpolynomial_fold *upwf)
387 {
388 isl_ctx *ctx;
389 struct isl_list *list;
390 int tight;
391
392 ctx = isl_union_set_get_ctx(uset);
393 list = isl_list_alloc(ctx, 2);
394 if (!list)
395 goto error2;
396
397 list->obj[0].type = isl_obj_union_pw_qpolynomial_fold;
398 list->obj[0].v = isl_union_set_apply_union_pw_qpolynomial_fold(uset,
399 upwf, &tight);
400 list->obj[1].type = isl_obj_bool;
401 list->obj[1].v = tight ? &iscc_bool_true : &iscc_bool_false;
402 if (tight < 0 || !list->obj[0].v)
403 goto error;
404
405 return list;
406 error2:
407 isl_union_set_free(uset);
408 isl_union_pw_qpolynomial_fold_free(upwf);
409 error:
410 isl_list_free(list);
411 return NULL;
412 }
413
414 static __isl_give isl_union_pw_qpolynomial *isl_val_mul_union_pw_qpolynomial(
415 __isl_take isl_val *v, __isl_take isl_union_pw_qpolynomial *upwqp)
416 {
417 return isl_union_pw_qpolynomial_scale_val(upwqp, v);
418 }
419
420 static __isl_give isl_union_pw_qpolynomial_fold *
421 int_val_mul_union_pw_qpolynomial_fold(__isl_take isl_val *v,
422 __isl_take isl_union_pw_qpolynomial_fold *upwf)
423 {
424 return isl_union_pw_qpolynomial_fold_scale_val(upwf, v);
425 }
426
427 /* Are the two strings "str1" and "str2" equal to each other?
428 */
429 static int str_eq(__isl_keep isl_str *str1, __isl_keep isl_str *str2)
430 {
431 if (!str1 || !str2)
432 return -1;
433
434 return !strcmp(str1->s, str2->s);
435 }
436
437 struct isc_bin_op bin_ops[] = {
438 { '+', isl_obj_bool, isl_obj_bool, isl_obj_bool,
439 (isc_bin_op_fn) &isl_bool_or },
440 { '*', isl_obj_bool, isl_obj_bool, isl_obj_bool,
441 (isc_bin_op_fn) &isl_bool_and },
442 { '+', isl_obj_val, isl_obj_val, isl_obj_val,
443 (isc_bin_op_fn) &isl_val_add },
444 { '-', isl_obj_val, isl_obj_val, isl_obj_val,
445 (isc_bin_op_fn) &isl_val_sub },
446 { '*', isl_obj_val, isl_obj_val, isl_obj_val,
447 (isc_bin_op_fn) &isl_val_mul },
448 { '+', isl_obj_pw_multi_aff, isl_obj_pw_multi_aff,
449 isl_obj_pw_multi_aff,
450 (isc_bin_op_fn) &isl_pw_multi_aff_add },
451 { '+', isl_obj_union_set, isl_obj_union_set,
452 isl_obj_union_set,
453 (isc_bin_op_fn) &isl_union_set_union },
454 { '+', isl_obj_union_map, isl_obj_union_map,
455 isl_obj_union_map,
456 (isc_bin_op_fn) &isl_union_map_union },
457 { '-', isl_obj_union_set, isl_obj_union_set,
458 isl_obj_union_set,
459 (isc_bin_op_fn) &isl_union_set_subtract },
460 { '-', isl_obj_union_map, isl_obj_union_map,
461 isl_obj_union_map,
462 (isc_bin_op_fn) &isl_union_map_subtract },
463 { '-', isl_obj_union_map, isl_obj_union_set,
464 isl_obj_union_map,
465 (isc_bin_op_fn) &isl_union_map_subtract_domain },
466 { '*', isl_obj_union_set, isl_obj_union_set,
467 isl_obj_union_set,
468 (isc_bin_op_fn) &isl_union_set_intersect },
469 { '*', isl_obj_union_map, isl_obj_union_map,
470 isl_obj_union_map,
471 (isc_bin_op_fn) &isl_union_map_intersect },
472 { '*', isl_obj_union_map, isl_obj_union_set,
473 isl_obj_union_map,
474 (isc_bin_op_fn) &isl_union_map_intersect_domain },
475 { '.', isl_obj_union_map, isl_obj_union_map,
476 isl_obj_union_map,
477 (isc_bin_op_fn) &isl_union_map_apply_range },
478 { '.', isl_obj_union_map, isl_obj_union_pw_qpolynomial,
479 isl_obj_union_pw_qpolynomial,
480 (isc_bin_op_fn) &isl_union_map_apply_union_pw_qpolynomial },
481 { '.', isl_obj_union_map, isl_obj_union_pw_qpolynomial_fold,
482 isl_obj_list,
483 (isc_bin_op_fn) &union_map_apply_union_pw_qpolynomial_fold },
484 { ISL_TOKEN_TO, isl_obj_union_set, isl_obj_union_set,
485 isl_obj_union_map,
486 (isc_bin_op_fn) &isl_union_map_from_domain_and_range },
487 { '=', isl_obj_union_set, isl_obj_union_set, isl_obj_bool,
488 { .test = (isc_bin_test_fn) &isl_union_set_is_equal } },
489 { '=', isl_obj_union_map, isl_obj_union_map, isl_obj_bool,
490 { .test = (isc_bin_test_fn) &isl_union_map_is_equal } },
491 { ISL_TOKEN_LE, isl_obj_union_set, isl_obj_union_set,
492 isl_obj_bool,
493 { .test = (isc_bin_test_fn) &isl_union_set_is_subset } },
494 { ISL_TOKEN_LE, isl_obj_union_map, isl_obj_union_map,
495 isl_obj_bool,
496 { .test = (isc_bin_test_fn) &isl_union_map_is_subset } },
497 { ISL_TOKEN_LT, isl_obj_union_set, isl_obj_union_set,
498 isl_obj_bool,
499 { .test = (isc_bin_test_fn) &isl_union_set_is_strict_subset } },
500 { ISL_TOKEN_LT, isl_obj_union_map, isl_obj_union_map,
501 isl_obj_bool,
502 { .test = (isc_bin_test_fn) &isl_union_map_is_strict_subset } },
503 { ISL_TOKEN_GE, isl_obj_union_set, isl_obj_union_set,
504 isl_obj_bool,
505 { .test = (isc_bin_test_fn) &isl_union_set_is_superset } },
506 { ISL_TOKEN_GE, isl_obj_union_map, isl_obj_union_map,
507 isl_obj_bool,
508 { .test = (isc_bin_test_fn) &isl_union_map_is_superset } },
509 { ISL_TOKEN_GT, isl_obj_union_set, isl_obj_union_set,
510 isl_obj_bool,
511 { .test =
512 (isc_bin_test_fn) &isl_union_set_is_strict_superset } },
513 { ISL_TOKEN_GT, isl_obj_union_map, isl_obj_union_map,
514 isl_obj_bool,
515 { .test =
516 (isc_bin_test_fn) &isl_union_map_is_strict_superset } },
517 { ISL_TOKEN_LEX_LE, isl_obj_union_set, isl_obj_union_set,
518 isl_obj_union_map,
519 (isc_bin_op_fn) &isl_union_set_lex_le_union_set },
520 { ISL_TOKEN_LEX_LT, isl_obj_union_set, isl_obj_union_set,
521 isl_obj_union_map,
522 (isc_bin_op_fn) &isl_union_set_lex_lt_union_set },
523 { ISL_TOKEN_LEX_GE, isl_obj_union_set, isl_obj_union_set,
524 isl_obj_union_map,
525 (isc_bin_op_fn) &isl_union_set_lex_ge_union_set },
526 { ISL_TOKEN_LEX_GT, isl_obj_union_set, isl_obj_union_set,
527 isl_obj_union_map,
528 (isc_bin_op_fn) &isl_union_set_lex_gt_union_set },
529 { ISL_TOKEN_LEX_LE, isl_obj_union_map, isl_obj_union_map,
530 isl_obj_union_map,
531 (isc_bin_op_fn) &isl_union_map_lex_le_union_map },
532 { ISL_TOKEN_LEX_LT, isl_obj_union_map, isl_obj_union_map,
533 isl_obj_union_map,
534 (isc_bin_op_fn) &isl_union_map_lex_lt_union_map },
535 { ISL_TOKEN_LEX_GE, isl_obj_union_map, isl_obj_union_map,
536 isl_obj_union_map,
537 (isc_bin_op_fn) &isl_union_map_lex_ge_union_map },
538 { ISL_TOKEN_LEX_GT, isl_obj_union_map, isl_obj_union_map,
539 isl_obj_union_map,
540 (isc_bin_op_fn) &isl_union_map_lex_gt_union_map },
541 { '.', isl_obj_union_pw_qpolynomial_fold,
542 isl_obj_union_pw_qpolynomial_fold,
543 isl_obj_union_pw_qpolynomial_fold,
544 (isc_bin_op_fn) &isl_union_pw_qpolynomial_fold_fold },
545 { '+', isl_obj_union_pw_qpolynomial, isl_obj_union_pw_qpolynomial,
546 isl_obj_union_pw_qpolynomial,
547 (isc_bin_op_fn) &isl_union_pw_qpolynomial_add },
548 { '+', isl_obj_union_pw_qpolynomial,
549 isl_obj_union_pw_qpolynomial_fold,
550 isl_obj_union_pw_qpolynomial_fold,
551 (isc_bin_op_fn) &union_pw_qpolynomial_add_union_pw_qpolynomial_fold },
552 { '+', isl_obj_union_pw_qpolynomial_fold,
553 isl_obj_union_pw_qpolynomial,
554 isl_obj_union_pw_qpolynomial_fold,
555 (isc_bin_op_fn) &isl_union_pw_qpolynomial_fold_add_union_pw_qpolynomial },
556 { '-', isl_obj_union_pw_qpolynomial, isl_obj_union_pw_qpolynomial,
557 isl_obj_union_pw_qpolynomial,
558 (isc_bin_op_fn) &isl_union_pw_qpolynomial_sub },
559 { '*', isl_obj_val, isl_obj_union_pw_qpolynomial,
560 isl_obj_union_pw_qpolynomial,
561 (isc_bin_op_fn) &isl_val_mul_union_pw_qpolynomial },
562 { '*', isl_obj_union_pw_qpolynomial, isl_obj_val,
563 isl_obj_union_pw_qpolynomial,
564 (isc_bin_op_fn) &isl_union_pw_qpolynomial_scale_val },
565 { '*', isl_obj_val, isl_obj_union_pw_qpolynomial_fold,
566 isl_obj_union_pw_qpolynomial_fold,
567 (isc_bin_op_fn) &int_val_mul_union_pw_qpolynomial_fold },
568 { '*', isl_obj_union_pw_qpolynomial_fold, isl_obj_val,
569 isl_obj_union_pw_qpolynomial_fold,
570 (isc_bin_op_fn) &isl_union_pw_qpolynomial_fold_scale_val },
571 { '*', isl_obj_union_pw_qpolynomial, isl_obj_union_pw_qpolynomial,
572 isl_obj_union_pw_qpolynomial,
573 (isc_bin_op_fn) &isl_union_pw_qpolynomial_mul },
574 { '*', isl_obj_union_pw_qpolynomial, isl_obj_union_set,
575 isl_obj_union_pw_qpolynomial,
576 (isc_bin_op_fn) &isl_union_pw_qpolynomial_intersect_domain },
577 { '*', isl_obj_union_pw_qpolynomial_fold, isl_obj_union_set,
578 isl_obj_union_pw_qpolynomial_fold,
579 (isc_bin_op_fn) &isl_union_pw_qpolynomial_fold_intersect_domain },
580 { '@', isl_obj_union_pw_qpolynomial, isl_obj_union_set,
581 isl_obj_union_pw_qpolynomial,
582 (isc_bin_op_fn) &isl_union_pw_qpolynomial_at },
583 { '@', isl_obj_union_pw_qpolynomial_fold, isl_obj_union_set,
584 isl_obj_union_pw_qpolynomial,
585 (isc_bin_op_fn) &isl_union_pw_qpolynomial_fold_at },
586 { '%', isl_obj_union_set, isl_obj_union_set,
587 isl_obj_union_set,
588 (isc_bin_op_fn) &isl_union_set_gist },
589 { '%', isl_obj_union_map, isl_obj_union_map,
590 isl_obj_union_map,
591 (isc_bin_op_fn) &isl_union_map_gist },
592 { '%', isl_obj_union_map, isl_obj_union_set,
593 isl_obj_union_map,
594 (isc_bin_op_fn) &isl_union_map_gist_domain },
595 { '%', isl_obj_union_pw_qpolynomial, isl_obj_union_set,
596 isl_obj_union_pw_qpolynomial,
597 (isc_bin_op_fn) &isl_union_pw_qpolynomial_gist },
598 { '%', isl_obj_union_pw_qpolynomial_fold, isl_obj_union_set,
599 isl_obj_union_pw_qpolynomial_fold,
600 (isc_bin_op_fn) &isl_union_pw_qpolynomial_fold_gist },
601 { ISL_TOKEN_EQ_EQ, isl_obj_union_pw_qpolynomial,
602 isl_obj_union_pw_qpolynomial, isl_obj_bool,
603 { .test = (isc_bin_test_fn)
604 &isl_union_pw_qpolynomial_plain_is_equal } },
605 { ISL_TOKEN_EQ_EQ, isl_obj_union_pw_qpolynomial_fold,
606 isl_obj_union_pw_qpolynomial_fold, isl_obj_bool,
607 { .test = (isc_bin_test_fn)
608 &isl_union_pw_qpolynomial_fold_plain_is_equal } },
609 { '+', isl_obj_str, isl_obj_str, isl_obj_str,
610 (isc_bin_op_fn) &isl_str_concat },
611 { '=', isl_obj_str, isl_obj_str, isl_obj_bool,
612 { .test = (isc_bin_test_fn) &str_eq } },
613 0
614 };
615
616 struct iscc_compound_bin_op compound_bin_ops[] = {
617 { "->*", ISL_TOKEN_TO, '*',
618 { -1, isl_obj_union_map, isl_obj_union_set,
619 isl_obj_union_map,
620 (isc_bin_op_fn) &isl_union_map_intersect_range } },
621 { "->-", ISL_TOKEN_TO, '-',
622 { -1, isl_obj_union_map, isl_obj_union_set,
623 isl_obj_union_map,
624 (isc_bin_op_fn) &isl_union_map_subtract_range } },
625 0
626 };
627
628 static __isl_give isl_union_map *map_after_map(__isl_take isl_union_map *umap1,
629 __isl_take isl_union_map *umap2)
630 {
631 return isl_union_map_apply_range(umap2, umap1);
632 }
633
634 static __isl_give isl_union_pw_qpolynomial *qpolynomial_after_map(
635 __isl_take isl_union_pw_qpolynomial *upwqp,
636 __isl_take isl_union_map *umap)
637 {
638 return isl_union_map_apply_union_pw_qpolynomial(umap, upwqp);
639 }
640
641 static __isl_give struct isl_list *qpolynomial_fold_after_map(
642 __isl_take isl_union_pw_qpolynomial_fold *upwf,
643 __isl_take isl_union_map *umap)
644 {
645 return union_map_apply_union_pw_qpolynomial_fold(umap, upwf);
646 }
647
648 struct isc_named_bin_op named_bin_ops[] = {
649 { "after", { -1, isl_obj_union_map, isl_obj_union_map,
650 isl_obj_union_map,
651 (isc_bin_op_fn) &map_after_map } },
652 { "after", { -1, isl_obj_union_pw_qpolynomial,
653 isl_obj_union_map, isl_obj_union_pw_qpolynomial,
654 (isc_bin_op_fn) &qpolynomial_after_map } },
655 { "after", { -1, isl_obj_union_pw_qpolynomial_fold,
656 isl_obj_union_map, isl_obj_list,
657 (isc_bin_op_fn) &qpolynomial_fold_after_map } },
658 { "before", { -1, isl_obj_union_map, isl_obj_union_map,
659 isl_obj_union_map,
660 (isc_bin_op_fn) &isl_union_map_apply_range } },
661 { "before", { -1, isl_obj_union_map,
662 isl_obj_union_pw_qpolynomial, isl_obj_union_pw_qpolynomial,
663 (isc_bin_op_fn) &isl_union_map_apply_union_pw_qpolynomial } },
664 { "before", { -1, isl_obj_union_map,
665 isl_obj_union_pw_qpolynomial_fold, isl_obj_list,
666 (isc_bin_op_fn) &union_map_apply_union_pw_qpolynomial_fold } },
667 { "cross", { -1, isl_obj_union_set, isl_obj_union_set,
668 isl_obj_union_set,
669 (isc_bin_op_fn) &isl_union_set_product } },
670 { "cross", { -1, isl_obj_union_map, isl_obj_union_map,
671 isl_obj_union_map,
672 (isc_bin_op_fn) &isl_union_map_product } },
673 { "cross_range", { -1, isl_obj_union_map,
674 isl_obj_union_map, isl_obj_union_map,
675 (isc_bin_op_fn) &isl_union_map_range_product } },
676 NULL
677 };
678
679 __isl_give isl_set *union_set_sample(__isl_take isl_union_set *uset)
680 {
681 return isl_set_from_basic_set(isl_union_set_sample(uset));
682 }
683
684 __isl_give isl_map *union_map_sample(__isl_take isl_union_map *umap)
685 {
686 return isl_map_from_basic_map(isl_union_map_sample(umap));
687 }
688
689 static __isl_give struct isl_list *union_map_power(
690 __isl_take isl_union_map *umap)
691 {
692 isl_ctx *ctx;
693 struct isl_list *list;
694 int exact;
695
696 ctx = isl_union_map_get_ctx(umap);
697 list = isl_list_alloc(ctx, 2);
698 if (!list)
699 goto error2;
700
701 list->obj[0].type = isl_obj_union_map;
702 list->obj[0].v = isl_union_map_power(umap, &exact);
703 list->obj[1].type = isl_obj_bool;
704 list->obj[1].v = exact ? &iscc_bool_true : &iscc_bool_false;
705 if (exact < 0 || !list->obj[0].v)
706 goto error;
707
708 return list;
709 error2:
710 isl_union_map_free(umap);
711 error:
712 isl_list_free(list);
713 return NULL;
714 }
715
716 /* Compute a lower or upper bound on "upwqp" depending on "type" and
717 * return a list containing two elements, the bound and a boolean
718 * indicating whether the result is tight.
719 */
720 static __isl_give struct isl_list *union_pw_qpolynomial_bound(
721 __isl_take isl_union_pw_qpolynomial *upwqp, enum isl_fold type)
722 {
723 isl_ctx *ctx;
724 struct isl_list *list;
725 int tight;
726
727 ctx = isl_union_pw_qpolynomial_get_ctx(upwqp);
728 list = isl_list_alloc(ctx, 2);
729 if (!list)
730 goto error2;
731
732 list->obj[0].type = isl_obj_union_pw_qpolynomial_fold;
733 list->obj[0].v = isl_union_pw_qpolynomial_bound(upwqp, type, &tight);
734 list->obj[1].type = isl_obj_bool;
735 list->obj[1].v = tight ? &iscc_bool_true : &iscc_bool_false;
736 if (tight < 0 || !list->obj[0].v)
737 goto error;
738
739 return list;
740 error2:
741 isl_union_pw_qpolynomial_free(upwqp);
742 error:
743 isl_list_free(list);
744 return NULL;
745 }
746
747 /* Compute a lower bound on "upwqp" and return a list containing
748 * two elements, the bound and a booleanindicating whether
749 * the result is tight.
750 */
751 static __isl_give struct isl_list *union_pw_qpolynomial_lower_bound(
752 __isl_take isl_union_pw_qpolynomial *upwqp)
753 {
754 return union_pw_qpolynomial_bound(upwqp, isl_fold_min);
755 }
756
757 /* Compute a upper bound on "upwqp" and return a list containing
758 * two elements, the bound and a booleanindicating whether
759 * the result is tight.
760 */
761 static __isl_give struct isl_list *union_pw_qpolynomial_upper_bound(
762 __isl_take isl_union_pw_qpolynomial *upwqp)
763 {
764 return union_pw_qpolynomial_bound(upwqp, isl_fold_max);
765 }
766
767 #ifdef HAVE_PET
768 /* Collect all statement instances, except those of kill statements.
769 */
770 static __isl_give isl_union_set *collect_non_kill_instances(
771 struct pet_scop *scop)
772 {
773 int i;
774 isl_set *domain_i;
775 isl_union_set *domain;
776
777 if (!scop)
778 return NULL;
779
780 domain = isl_union_set_empty(isl_set_get_space(scop->context));
781
782 for (i = 0; i < scop->n_stmt; ++i) {
783 struct pet_stmt *stmt = scop->stmts[i];
784
785 if (pet_stmt_is_kill(stmt))
786 continue;
787 domain_i = isl_set_copy(stmt->domain);
788 if (scop->stmts[i]->n_arg > 0)
789 domain_i = isl_map_domain(isl_set_unwrap(domain_i));
790 domain = isl_union_set_add_set(domain, domain_i);
791 }
792
793 return domain;
794 }
795
796 static __isl_give isl_list *parse(__isl_take isl_str *str)
797 {
798 isl_ctx *ctx;
799 struct isl_list *list;
800 struct pet_scop *scop;
801 isl_schedule *sched;
802 isl_union_map *may_reads, *must_writes, *may_writes;
803 isl_union_set *domain;
804 struct iscc_options *options;
805
806 if (!str)
807 return NULL;
808 ctx = str->ctx;
809
810 options = isl_ctx_peek_iscc_options(ctx);
811 if (!options || !options->io) {
812 isl_str_free(str);
813 isl_die(ctx, isl_error_invalid,
814 "parse_file operation not allowed", return NULL);
815 }
816
817 list = isl_list_alloc(ctx, 5);
818 if (!list)
819 goto error;
820
821 scop = pet_scop_extract_from_C_source(ctx, str->s, NULL);
822 domain = collect_non_kill_instances(scop);
823 sched = pet_scop_get_schedule(scop);
824 sched = isl_schedule_intersect_domain(sched,
825 isl_union_set_copy(domain));
826 may_reads = pet_scop_get_may_reads(scop);
827 may_writes = pet_scop_get_may_writes(scop);
828 must_writes = pet_scop_get_must_writes(scop);
829 pet_scop_free(scop);
830
831 list->obj[0].type = isl_obj_union_set;
832 list->obj[0].v = domain;
833 list->obj[1].type = isl_obj_union_map;
834 list->obj[1].v = must_writes;
835 list->obj[2].type = isl_obj_union_map;
836 list->obj[2].v = may_writes;
837 list->obj[3].type = isl_obj_union_map;
838 list->obj[3].v = may_reads;
839 list->obj[4].type = isl_obj_schedule;
840 list->obj[4].v = sched;
841
842 if (!list->obj[0].v || !list->obj[1].v ||
843 !list->obj[2].v || !list->obj[3].v || !list->obj[4].v)
844 goto error;
845
846 isl_str_free(str);
847 return list;
848 error:
849 isl_list_free(list);
850 isl_str_free(str);
851 return NULL;
852 }
853 #endif
854
855 static isl_stat add_point(__isl_take isl_point *pnt, void *user)
856 {
857 isl_union_set **scan = (isl_union_set **) user;
858
859 *scan = isl_union_set_add_set(*scan, isl_set_from_point(pnt));
860
861 return isl_stat_ok;
862 }
863
864 static __isl_give isl_union_set *union_set_scan(__isl_take isl_union_set *uset)
865 {
866 isl_union_set *scan;
867
868 scan = isl_union_set_empty(isl_union_set_get_space(uset));
869
870 if (isl_union_set_foreach_point(uset, add_point, &scan) < 0) {
871 isl_union_set_free(scan);
872 return uset;
873 }
874
875 isl_union_set_free(uset);
876 return scan;
877 }
878
879 static __isl_give isl_union_map *union_map_scan(__isl_take isl_union_map *umap)
880 {
881 return isl_union_set_unwrap(union_set_scan(isl_union_map_wrap(umap)));
882 }
883
884 static __isl_give isl_union_pw_qpolynomial *union_pw_qpolynomial_poly(
885 __isl_take isl_union_pw_qpolynomial *upwqp)
886 {
887 return isl_union_pw_qpolynomial_to_polynomial(upwqp, 0);
888 }
889
890 static __isl_give isl_union_pw_qpolynomial *union_pw_qpolynomial_lpoly(
891 __isl_take isl_union_pw_qpolynomial *upwqp)
892 {
893 return isl_union_pw_qpolynomial_to_polynomial(upwqp, -1);
894 }
895
896 static __isl_give isl_union_pw_qpolynomial *union_pw_qpolynomial_upoly(
897 __isl_take isl_union_pw_qpolynomial *upwqp)
898 {
899 return isl_union_pw_qpolynomial_to_polynomial(upwqp, 1);
900 }
901
902 /* Return the domain of "schedule".
903 */
904 static __isl_give isl_union_set *schedule_domain(
905 __isl_take isl_schedule *schedule)
906 {
907 isl_union_set *domain;
908
909 domain = isl_schedule_get_domain(schedule);
910 isl_schedule_free(schedule);
911
912 return domain;
913 }
914
915 /* Convert "schedule" to a union map representation.
916 */
917 static __isl_give isl_union_map *schedule_map(__isl_take isl_schedule *schedule)
918 {
919 isl_union_map *map;
920
921 map = isl_schedule_get_map(schedule);
922 isl_schedule_free(schedule);
923
924 return map;
925 }
926
927 typedef void *(*isc_un_op_fn)(void *arg);
928 struct isc_un_op {
929 enum isl_token_type op;
930 isl_obj_type arg;
931 isl_obj_type res;
932 isc_un_op_fn fn;
933 };
934 struct isc_named_un_op {
935 char *name;
936 struct isc_un_op op;
937 };
938 struct isc_named_un_op named_un_ops[] = {
939 {"aff", { -1, isl_obj_union_map, isl_obj_union_map,
940 (isc_un_op_fn) &isl_union_map_affine_hull } },
941 {"aff", { -1, isl_obj_union_set, isl_obj_union_set,
942 (isc_un_op_fn) &isl_union_set_affine_hull } },
943 {"card", { -1, isl_obj_union_set,
944 isl_obj_union_pw_qpolynomial,
945 (isc_un_op_fn) &isl_union_set_card } },
946 {"card", { -1, isl_obj_union_map,
947 isl_obj_union_pw_qpolynomial,
948 (isc_un_op_fn) &isl_union_map_card } },
949 {"coalesce", { -1, isl_obj_union_set, isl_obj_union_set,
950 (isc_un_op_fn) &isl_union_set_coalesce } },
951 {"coalesce", { -1, isl_obj_union_map, isl_obj_union_map,
952 (isc_un_op_fn) &isl_union_map_coalesce } },
953 {"coalesce", { -1, isl_obj_union_pw_qpolynomial,
954 isl_obj_union_pw_qpolynomial,
955 (isc_un_op_fn) &isl_union_pw_qpolynomial_coalesce } },
956 {"coalesce", { -1, isl_obj_union_pw_qpolynomial_fold,
957 isl_obj_union_pw_qpolynomial_fold,
958 (isc_un_op_fn) &isl_union_pw_qpolynomial_fold_coalesce } },
959 {"coefficients", { -1, isl_obj_union_set,
960 isl_obj_union_set,
961 (isc_un_op_fn) &isl_union_set_coefficients } },
962 {"solutions", { -1, isl_obj_union_set, isl_obj_union_set,
963 (isc_un_op_fn) &isl_union_set_solutions } },
964 {"deltas", { -1, isl_obj_union_map, isl_obj_union_set,
965 (isc_un_op_fn) &isl_union_map_deltas } },
966 {"deltas_map", { -1, isl_obj_union_map, isl_obj_union_map,
967 (isc_un_op_fn) &isl_union_map_deltas_map } },
968 {"dom", { -1, isl_obj_schedule, isl_obj_union_set,
969 (isc_un_op_fn) &schedule_domain } },
970 {"dom", { -1, isl_obj_union_map, isl_obj_union_set,
971 (isc_un_op_fn) &isl_union_map_domain } },
972 {"dom", { -1, isl_obj_union_pw_qpolynomial, isl_obj_union_set,
973 (isc_un_op_fn) &isl_union_pw_qpolynomial_domain } },
974 {"dom", { -1, isl_obj_union_pw_qpolynomial_fold,
975 isl_obj_union_set,
976 (isc_un_op_fn) &isl_union_pw_qpolynomial_fold_domain } },
977 {"domain", { -1, isl_obj_schedule, isl_obj_union_set,
978 (isc_un_op_fn) &schedule_domain } },
979 {"domain", { -1, isl_obj_union_map, isl_obj_union_set,
980 (isc_un_op_fn) &isl_union_map_domain } },
981 {"domain", { -1, isl_obj_union_pw_qpolynomial,
982 isl_obj_union_set,
983 (isc_un_op_fn) &isl_union_pw_qpolynomial_domain } },
984 {"domain", { -1, isl_obj_union_pw_qpolynomial_fold,
985 isl_obj_union_set,
986 (isc_un_op_fn) &isl_union_pw_qpolynomial_fold_domain } },
987 {"domain_map", { -1, isl_obj_union_map, isl_obj_union_map,
988 (isc_un_op_fn) &isl_union_map_domain_map } },
989 {"ran", { -1, isl_obj_union_map, isl_obj_union_set,
990 (isc_un_op_fn) &isl_union_map_range } },
991 {"range", { -1, isl_obj_union_map, isl_obj_union_set,
992 (isc_un_op_fn) &isl_union_map_range } },
993 {"range_map", { -1, isl_obj_union_map, isl_obj_union_map,
994 (isc_un_op_fn) &isl_union_map_range_map } },
995 {"identity", { -1, isl_obj_union_set, isl_obj_union_map,
996 (isc_un_op_fn) &isl_union_set_identity } },
997 {"lattice_width", { -1, isl_obj_union_set,
998 isl_obj_union_pw_qpolynomial,
999 (isc_un_op_fn) &isl_union_set_lattice_width } },
1000 {"lb", { -1, isl_obj_union_pw_qpolynomial, isl_obj_list,
1001 (isc_un_op_fn) &union_pw_qpolynomial_lower_bound } },
1002 {"lexmin", { -1, isl_obj_union_map, isl_obj_union_map,
1003 (isc_un_op_fn) &isl_union_map_lexmin } },
1004 {"lexmax", { -1, isl_obj_union_map, isl_obj_union_map,
1005 (isc_un_op_fn) &isl_union_map_lexmax } },
1006 {"lexmin", { -1, isl_obj_union_set, isl_obj_union_set,
1007 (isc_un_op_fn) &isl_union_set_lexmin } },
1008 {"lexmax", { -1, isl_obj_union_set, isl_obj_union_set,
1009 (isc_un_op_fn) &isl_union_set_lexmax } },
1010 {"lift", { -1, isl_obj_union_set, isl_obj_union_set,
1011 (isc_un_op_fn) &isl_union_set_lift } },
1012 {"map", { -1, isl_obj_schedule, isl_obj_union_map,
1013 (isc_un_op_fn) &schedule_map } },
1014 {"params", { -1, isl_obj_union_map, isl_obj_set,
1015 (isc_un_op_fn) &isl_union_map_params } },
1016 {"params", { -1, isl_obj_union_set, isl_obj_set,
1017 (isc_un_op_fn) &isl_union_set_params } },
1018 {"poly", { -1, isl_obj_union_map, isl_obj_union_map,
1019 (isc_un_op_fn) &isl_union_map_polyhedral_hull } },
1020 {"poly", { -1, isl_obj_union_set, isl_obj_union_set,
1021 (isc_un_op_fn) &isl_union_set_polyhedral_hull } },
1022 {"poly", { -1, isl_obj_union_pw_qpolynomial,
1023 isl_obj_union_pw_qpolynomial,
1024 (isc_un_op_fn) &union_pw_qpolynomial_poly } },
1025 {"lpoly", { -1, isl_obj_union_pw_qpolynomial,
1026 isl_obj_union_pw_qpolynomial,
1027 (isc_un_op_fn) &union_pw_qpolynomial_lpoly } },
1028 {"upoly", { -1, isl_obj_union_pw_qpolynomial,
1029 isl_obj_union_pw_qpolynomial,
1030 (isc_un_op_fn) &union_pw_qpolynomial_upoly } },
1031 #ifdef HAVE_PET
1032 {"parse_file", { -1, isl_obj_str, isl_obj_list,
1033 (isc_un_op_fn) &parse } },
1034 #endif
1035 {"pow", { -1, isl_obj_union_map, isl_obj_list,
1036 (isc_un_op_fn) &union_map_power } },
1037 {"sample", { -1, isl_obj_union_set, isl_obj_set,
1038 (isc_un_op_fn) &union_set_sample } },
1039 {"sample", { -1, isl_obj_union_map, isl_obj_map,
1040 (isc_un_op_fn) &union_map_sample } },
1041 {"scan", { -1, isl_obj_union_set, isl_obj_union_set,
1042 (isc_un_op_fn) &union_set_scan } },
1043 {"scan", { -1, isl_obj_union_map, isl_obj_union_map,
1044 (isc_un_op_fn) &union_map_scan } },
1045 {"sum", { -1, isl_obj_union_pw_qpolynomial,
1046 isl_obj_union_pw_qpolynomial,
1047 (isc_un_op_fn) &isl_union_pw_qpolynomial_sum } },
1048 {"ub", { -1, isl_obj_union_pw_qpolynomial, isl_obj_list,
1049 (isc_un_op_fn) &union_pw_qpolynomial_upper_bound } },
1050 {"unwrap", { -1, isl_obj_union_set, isl_obj_union_map,
1051 (isc_un_op_fn) &isl_union_set_unwrap } },
1052 {"wrap", { -1, isl_obj_union_map, isl_obj_union_set,
1053 (isc_un_op_fn) &isl_union_map_wrap } },
1054 {"zip", { -1, isl_obj_union_map, isl_obj_union_map,
1055 (isc_un_op_fn) &isl_union_map_zip } },
1056 NULL
1057 };
1058
1059 struct isl_named_obj {
1060 char *name;
1061 struct isl_obj obj;
1062 };
1063
1064 static void free_obj(struct isl_obj obj)
1065 {
1066 obj.type->free(obj.v);
1067 }
1068
1069 static int same_name(const void *entry, const void *val)
1070 {
1071 const struct isl_named_obj *named = (const struct isl_named_obj *)entry;
1072
1073 return !strcmp(named->name, val);
1074 }
1075
1076 static int do_assign(struct isl_ctx *ctx, struct isl_hash_table *table,
1077 char *name, struct isl_obj obj)
1078 {
1079 struct isl_hash_table_entry *entry;
1080 uint32_t name_hash;
1081 struct isl_named_obj *named;
1082
1083 name_hash = isl_hash_string(isl_hash_init(), name);
1084 entry = isl_hash_table_find(ctx, table, name_hash, same_name, name, 1);
1085 if (!entry)
1086 goto error;
1087 if (entry->data) {
1088 named = entry->data;
1089 free_obj(named->obj);
1090 free(name);
1091 } else {
1092 named = isl_alloc_type(ctx, struct isl_named_obj);
1093 if (!named)
1094 goto error;
1095 named->name = name;
1096 entry->data = named;
1097 }
1098 named->obj = obj;
1099
1100 return 0;
1101 error:
1102 free_obj(obj);
1103 free(name);
1104 return -1;
1105 }
1106
1107 static struct isl_obj stored_obj(struct isl_ctx *ctx,
1108 struct isl_hash_table *table, char *name)
1109 {
1110 struct isl_obj obj = { isl_obj_none, NULL };
1111 struct isl_hash_table_entry *entry;
1112 uint32_t name_hash;
1113
1114 name_hash = isl_hash_string(isl_hash_init(), name);
1115 entry = isl_hash_table_find(ctx, table, name_hash, same_name, name, 0);
1116 if (entry) {
1117 struct isl_named_obj *named;
1118 named = entry->data;
1119 obj = named->obj;
1120 } else if (isdigit(name[0]))
1121 fprintf(stderr, "unknown identifier '$%s'\n", name);
1122 else
1123 fprintf(stderr, "unknown identifier '%s'\n", name);
1124
1125 free(name);
1126 obj.v = obj.type->copy(obj.v);
1127 return obj;
1128 }
1129
1130 static int is_subtype(struct isl_obj obj, isl_obj_type super)
1131 {
1132 if (obj.type == super)
1133 return 1;
1134 if (obj.type == isl_obj_map && super == isl_obj_union_map)
1135 return 1;
1136 if (obj.type == isl_obj_set && super == isl_obj_union_set)
1137 return 1;
1138 if (obj.type == isl_obj_schedule && super == isl_obj_union_map)
1139 return 1;
1140 if (obj.type == isl_obj_pw_multi_aff && super == isl_obj_union_set) {
1141 isl_space *space = isl_pw_multi_aff_get_space(obj.v);
1142 int is_set = isl_space_is_set(space);
1143 isl_space_free(space);
1144 return is_set;
1145 }
1146 if (obj.type == isl_obj_pw_qpolynomial &&
1147 super == isl_obj_union_pw_qpolynomial)
1148 return 1;
1149 if (obj.type == isl_obj_pw_qpolynomial_fold &&
1150 super == isl_obj_union_pw_qpolynomial_fold)
1151 return 1;
1152 if (obj.type == isl_obj_union_set && isl_union_set_is_empty(obj.v))
1153 return 1;
1154 if (obj.type == isl_obj_list) {
1155 struct isl_list *list = obj.v;
1156 if (list->n == 2 && list->obj[1].type == isl_obj_bool)
1157 return is_subtype(list->obj[0], super);
1158 }
1159 if (super == isl_obj_str)
1160 return 1;
1161 return 0;
1162 }
1163
1164 static struct isl_obj obj_at(struct isl_obj obj, int i)
1165 {
1166 struct isl_list *list = obj.v;
1167
1168 obj = list->obj[i];
1169 obj.v = obj.type->copy(obj.v);
1170
1171 isl_list_free(list);
1172
1173 return obj;
1174 }
1175
1176 static struct isl_obj convert(isl_ctx *ctx, struct isl_obj obj,
1177 isl_obj_type type)
1178 {
1179 if (obj.type == type)
1180 return obj;
1181 if (obj.type == isl_obj_pw_multi_aff && type == isl_obj_union_set) {
1182 isl_set *set = isl_set_from_pw_multi_aff(obj.v);
1183 obj.type = isl_obj_union_set;
1184 obj.v = isl_union_set_from_set(set);
1185 return obj;
1186 }
1187 if (obj.type == isl_obj_map && type == isl_obj_union_map) {
1188 obj.type = isl_obj_union_map;
1189 obj.v = isl_union_map_from_map(obj.v);
1190 return obj;
1191 }
1192 if (obj.type == isl_obj_schedule && type == isl_obj_union_map) {
1193 obj.type = isl_obj_union_map;
1194 obj.v = schedule_map(obj.v);
1195 return obj;
1196 }
1197 if (obj.type == isl_obj_set && type == isl_obj_union_set) {
1198 obj.type = isl_obj_union_set;
1199 obj.v = isl_union_set_from_set(obj.v);
1200 return obj;
1201 }
1202 if (obj.type == isl_obj_pw_qpolynomial &&
1203 type == isl_obj_union_pw_qpolynomial) {
1204 obj.type = isl_obj_union_pw_qpolynomial;
1205 obj.v = isl_union_pw_qpolynomial_from_pw_qpolynomial(obj.v);
1206 return obj;
1207 }
1208 if (obj.type == isl_obj_pw_qpolynomial_fold &&
1209 type == isl_obj_union_pw_qpolynomial_fold) {
1210 obj.type = isl_obj_union_pw_qpolynomial_fold;
1211 obj.v = isl_union_pw_qpolynomial_fold_from_pw_qpolynomial_fold(obj.v);
1212 return obj;
1213 }
1214 if (obj.type == isl_obj_union_set && isl_union_set_is_empty(obj.v)) {
1215 if (type == isl_obj_union_map) {
1216 obj.type = isl_obj_union_map;
1217 return obj;
1218 }
1219 if (type == isl_obj_union_pw_qpolynomial) {
1220 isl_space *dim = isl_union_set_get_space(obj.v);
1221 isl_union_set_free(obj.v);
1222 obj.v = isl_union_pw_qpolynomial_zero(dim);
1223 obj.type = isl_obj_union_pw_qpolynomial;
1224 return obj;
1225 }
1226 if (type == isl_obj_union_pw_qpolynomial_fold) {
1227 isl_space *dim = isl_union_set_get_space(obj.v);
1228 isl_union_set_free(obj.v);
1229 obj.v = isl_union_pw_qpolynomial_fold_zero(dim,
1230 isl_fold_list);
1231 obj.type = isl_obj_union_pw_qpolynomial_fold;
1232 return obj;
1233 }
1234 }
1235 if (obj.type == isl_obj_list) {
1236 struct isl_list *list = obj.v;
1237 if (list->n == 2 && list->obj[1].type == isl_obj_bool)
1238 return convert(ctx, obj_at(obj, 0), type);
1239 }
1240 if (type == isl_obj_str) {
1241 isl_str *str;
1242 isl_printer *p;
1243 char *s;
1244
1245 p = isl_printer_to_str(ctx);
1246 if (!p)
1247 goto error;
1248 p = obj.type->print(p, obj.v);
1249 s = isl_printer_get_str(p);
1250 isl_printer_free(p);
1251
1252 str = isl_str_from_string(ctx, s);
1253 if (!str)
1254 goto error;
1255 free_obj(obj);
1256 obj.v = str;
1257 obj.type = isl_obj_str;
1258 return obj;
1259
1260 }
1261 error:
1262 free_obj(obj);
1263 obj.type = isl_obj_none;
1264 obj.v = NULL;
1265 return obj;
1266 }
1267
1268 static struct isc_bin_op *read_bin_op_if_available(struct isl_stream *s,
1269 struct isl_obj lhs)
1270 {
1271 int i;
1272 int read_tok2 = 0;
1273 struct isl_token *tok, *tok2;
1274
1275 tok = isl_stream_next_token(s);
1276 if (!tok)
1277 return NULL;
1278
1279 for (i = 0; ; ++i) {
1280 if (!bin_ops[i].op)
1281 break;
1282 if (bin_ops[i].op != isl_token_get_type(tok))
1283 continue;
1284 if (!is_subtype(lhs, bin_ops[i].lhs))
1285 continue;
1286
1287 isl_token_free(tok);
1288 return &bin_ops[i];
1289 }
1290
1291 for (i = 0; ; ++i) {
1292 if (!named_bin_ops[i].name)
1293 break;
1294 if (named_bin_ops[i].op.op != isl_token_get_type(tok))
1295 continue;
1296 if (!is_subtype(lhs, named_bin_ops[i].op.lhs))
1297 continue;
1298
1299 isl_token_free(tok);
1300 return &named_bin_ops[i].op;
1301 }
1302
1303 for (i = 0; ; ++i) {
1304 if (!compound_bin_ops[i].full)
1305 break;
1306 if (compound_bin_ops[i].op1 != isl_token_get_type(tok))
1307 continue;
1308 if (!read_tok2)
1309 tok2 = isl_stream_next_token(s);
1310 read_tok2 = 1;
1311 if (compound_bin_ops[i].op2 != isl_token_get_type(tok2))
1312 continue;
1313 if (!is_subtype(lhs, compound_bin_ops[i].op.lhs))
1314 continue;
1315
1316 isl_token_free(tok2);
1317 isl_token_free(tok);
1318 return &compound_bin_ops[i].op;
1319 }
1320
1321 if (read_tok2)
1322 isl_stream_push_token(s, tok2);
1323 isl_stream_push_token(s, tok);
1324
1325 return NULL;
1326 }
1327
1328 static struct isc_un_op *read_prefix_un_op_if_available(struct isl_stream *s)
1329 {
1330 int i;
1331 struct isl_token *tok;
1332
1333 tok = isl_stream_next_token(s);
1334 if (!tok)
1335 return NULL;
1336
1337 for (i = 0; ; ++i) {
1338 if (!named_un_ops[i].name)
1339 break;
1340 if (named_un_ops[i].op.op != isl_token_get_type(tok))
1341 continue;
1342
1343 isl_token_free(tok);
1344 return &named_un_ops[i].op;
1345 }
1346
1347 isl_stream_push_token(s, tok);
1348
1349 return NULL;
1350 }
1351
1352 static struct isc_un_op *find_matching_un_op(struct isc_un_op *like,
1353 struct isl_obj arg)
1354 {
1355 int i;
1356
1357 for (i = 0; ; ++i) {
1358 if (!named_un_ops[i].name)
1359 break;
1360 if (named_un_ops[i].op.op != like->op)
1361 continue;
1362 if (!is_subtype(arg, named_un_ops[i].op.arg))
1363 continue;
1364
1365 return &named_un_ops[i].op;
1366 }
1367
1368 return NULL;
1369 }
1370
1371 static int is_assign(struct isl_stream *s)
1372 {
1373 struct isl_token *tok;
1374 struct isl_token *tok2;
1375 int assign;
1376
1377 tok = isl_stream_next_token(s);
1378 if (!tok)
1379 return 0;
1380 if (isl_token_get_type(tok) != ISL_TOKEN_IDENT) {
1381 isl_stream_push_token(s, tok);
1382 return 0;
1383 }
1384
1385 tok2 = isl_stream_next_token(s);
1386 if (!tok2) {
1387 isl_stream_push_token(s, tok);
1388 return 0;
1389 }
1390 assign = isl_token_get_type(tok2) == ISL_TOKEN_DEF;
1391 isl_stream_push_token(s, tok2);
1392 isl_stream_push_token(s, tok);
1393
1394 return assign;
1395 }
1396
1397 static struct isl_obj read_obj(struct isl_stream *s,
1398 struct isl_hash_table *table);
1399 static struct isl_obj read_expr(struct isl_stream *s,
1400 struct isl_hash_table *table);
1401
1402 static struct isl_obj read_un_op_expr(struct isl_stream *s,
1403 struct isl_hash_table *table, struct isc_un_op *op)
1404 {
1405 isl_ctx *ctx;
1406 struct isl_obj obj = { isl_obj_none, NULL };
1407
1408 obj = read_obj(s, table);
1409 if (!obj.v)
1410 goto error;
1411
1412 op = find_matching_un_op(op, obj);
1413
1414 ctx = isl_stream_get_ctx(s);
1415 if (!op)
1416 isl_die(ctx, isl_error_invalid,
1417 "no such unary operator defined on given operand",
1418 goto error);
1419
1420 obj = convert(ctx, obj, op->arg);
1421 obj.v = op->fn(obj.v);
1422 obj.type = op->res;
1423
1424 return obj;
1425 error:
1426 free_obj(obj);
1427 obj.type = isl_obj_none;
1428 obj.v = NULL;
1429 return obj;
1430 }
1431
1432 static struct isl_obj transitive_closure(struct isl_ctx *ctx, struct isl_obj obj)
1433 {
1434 struct isl_list *list;
1435 int exact;
1436
1437 if (obj.type != isl_obj_union_map)
1438 obj = convert(ctx, obj, isl_obj_union_map);
1439 isl_assert(ctx, obj.type == isl_obj_union_map, goto error);
1440 list = isl_list_alloc(ctx, 2);
1441 if (!list)
1442 goto error;
1443
1444 list->obj[0].type = isl_obj_union_map;
1445 list->obj[0].v = isl_union_map_transitive_closure(obj.v, &exact);
1446 list->obj[1].type = isl_obj_bool;
1447 list->obj[1].v = exact ? &iscc_bool_true : &iscc_bool_false;
1448 obj.v = list;
1449 obj.type = isl_obj_list;
1450 if (exact < 0 || !list->obj[0].v)
1451 goto error;
1452
1453 return obj;
1454 error:
1455 free_obj(obj);
1456 obj.type = isl_obj_none;
1457 obj.v = NULL;
1458 return obj;
1459 }
1460
1461 static struct isl_obj obj_at_index(struct isl_stream *s, struct isl_obj obj)
1462 {
1463 struct isl_list *list = obj.v;
1464 struct isl_token *tok;
1465 isl_ctx *ctx;
1466 isl_val *v;
1467 int i;
1468
1469 tok = isl_stream_next_token(s);
1470 if (!tok || isl_token_get_type(tok) != ISL_TOKEN_VALUE) {
1471 isl_stream_error(s, tok, "expecting index");
1472 if (tok)
1473 isl_stream_push_token(s, tok);
1474 goto error;
1475 }
1476 ctx = isl_stream_get_ctx(s);
1477 v = isl_token_get_val(ctx, tok);
1478 i = isl_val_get_num_si(v);
1479 isl_val_free(v);
1480 isl_token_free(tok);
1481 isl_assert(ctx, i < list->n, goto error);
1482 if (isl_stream_eat(s, ']'))
1483 goto error;
1484
1485 return obj_at(obj, i);
1486 error:
1487 free_obj(obj);
1488 obj.type = isl_obj_none;
1489 obj.v = NULL;
1490 return obj;
1491 }
1492
1493 static struct isl_obj apply(struct isl_stream *s, __isl_take isl_union_map *umap,
1494 struct isl_hash_table *table)
1495 {
1496 isl_ctx *ctx;
1497 struct isl_obj obj;
1498
1499 obj = read_expr(s, table);
1500 ctx = isl_stream_get_ctx(s);
1501 isl_assert(ctx, is_subtype(obj, isl_obj_union_set) ||
1502 is_subtype(obj, isl_obj_union_map), goto error);
1503
1504 if (obj.type == isl_obj_list) {
1505 struct isl_list *list = obj.v;
1506 if (list->n == 2 && list->obj[1].type == isl_obj_bool)
1507 obj = obj_at(obj, 0);
1508 }
1509 if (obj.type == isl_obj_set)
1510 obj = convert(ctx, obj, isl_obj_union_set);
1511 else if (obj.type == isl_obj_map)
1512 obj = convert(ctx, obj, isl_obj_union_map);
1513 if (obj.type == isl_obj_union_set) {
1514 obj.v = isl_union_set_apply(obj.v, umap);
1515 } else
1516 obj.v = isl_union_map_apply_range(obj.v, umap);
1517 if (!obj.v)
1518 goto error2;
1519
1520 if (isl_stream_eat(s, ')'))
1521 goto error2;
1522
1523 return obj;
1524 error:
1525 isl_union_map_free(umap);
1526 error2:
1527 free_obj(obj);
1528 obj.type = isl_obj_none;
1529 obj.v = NULL;
1530 return obj;
1531 }
1532
1533 static struct isl_obj apply_fun_set(struct isl_obj obj,
1534 __isl_take isl_union_set *uset)
1535 {
1536 if (obj.type == isl_obj_union_pw_qpolynomial) {
1537 obj.v = isl_union_set_apply_union_pw_qpolynomial(uset, obj.v);
1538 } else {
1539 obj.type = isl_obj_list;
1540 obj.v = union_set_apply_union_pw_qpolynomial_fold(uset, obj.v);
1541 }
1542 return obj;
1543 }
1544
1545 static struct isl_obj apply_fun_map(struct isl_obj obj,
1546 __isl_take isl_union_map *umap)
1547 {
1548 if (obj.type == isl_obj_union_pw_qpolynomial) {
1549 obj.v = isl_union_map_apply_union_pw_qpolynomial(umap, obj.v);
1550 } else {
1551 obj.type = isl_obj_list;
1552 obj.v = union_map_apply_union_pw_qpolynomial_fold(umap, obj.v);
1553 }
1554 return obj;
1555 }
1556
1557 static struct isl_obj apply_fun(struct isl_stream *s,
1558 struct isl_obj obj, struct isl_hash_table *table)
1559 {
1560 struct isl_obj arg;
1561 isl_ctx *ctx;
1562
1563 arg = read_expr(s, table);
1564 ctx = isl_stream_get_ctx(s);
1565 if (!is_subtype(arg, isl_obj_union_map) &&
1566 !is_subtype(arg, isl_obj_union_set))
1567 isl_die(ctx, isl_error_invalid,
1568 "expecting set of map argument", goto error);
1569
1570 if (arg.type == isl_obj_list) {
1571 struct isl_list *list = arg.v;
1572 if (list->n == 2 && list->obj[1].type == isl_obj_bool)
1573 arg = obj_at(arg, 0);
1574 }
1575 if (arg.type == isl_obj_set)
1576 arg = convert(ctx, arg, isl_obj_union_set);
1577 else if (arg.type == isl_obj_map)
1578 arg = convert(ctx, arg, isl_obj_union_map);
1579 if (arg.type == isl_obj_union_set)
1580 obj = apply_fun_set(obj, arg.v);
1581 else
1582 obj = apply_fun_map(obj, arg.v);
1583 if (!obj.v)
1584 goto error2;
1585
1586 if (isl_stream_eat(s, ')'))
1587 goto error2;
1588
1589 return obj;
1590 error:
1591 free_obj(arg);
1592 error2:
1593 free_obj(obj);
1594 obj.type = isl_obj_none;
1595 obj.v = NULL;
1596 return obj;
1597 }
1598
1599 struct add_vertex_data {
1600 struct isl_list *list;
1601 int i;
1602 };
1603
1604 static isl_stat add_vertex(__isl_take isl_vertex *vertex, void *user)
1605 {
1606 struct add_vertex_data *data = (struct add_vertex_data *)user;
1607 isl_multi_aff *ma;
1608 isl_set *dom;
1609
1610 ma = isl_vertex_get_expr(vertex);
1611 dom = isl_set_from_basic_set(isl_vertex_get_domain(vertex));
1612
1613 data->list->obj[data->i].type = isl_obj_pw_multi_aff;
1614 data->list->obj[data->i].v = isl_pw_multi_aff_alloc(dom, ma);
1615 data->i++;
1616
1617 isl_vertex_free(vertex);
1618
1619 return isl_stat_ok;
1620 }
1621
1622 static isl_stat set_vertices(__isl_take isl_set *set, void *user)
1623 {
1624 isl_ctx *ctx;
1625 isl_basic_set *hull;
1626 isl_vertices *vertices = NULL;
1627 struct isl_list *list = NULL;
1628 isl_stat r;
1629 struct add_vertex_data *data = (struct add_vertex_data *)user;
1630
1631 set = isl_set_remove_divs(set);
1632 hull = isl_set_convex_hull(set);
1633 vertices = isl_basic_set_compute_vertices(hull);
1634 isl_basic_set_free(hull);
1635
1636 list = data->list;
1637
1638 ctx = isl_vertices_get_ctx(vertices);
1639 data->list = isl_list_alloc(ctx, isl_vertices_get_n_vertices(vertices));
1640 if (!data->list)
1641 goto error;
1642
1643 data->i = 0;
1644 r = isl_vertices_foreach_vertex(vertices, &add_vertex, user);
1645
1646 data->list = isl_list_concat(list, data->list);
1647
1648 isl_vertices_free(vertices);
1649
1650 return r;
1651 error:
1652 data->list = list;
1653 isl_vertices_free(vertices);
1654 return isl_stat_error;
1655 }
1656
1657 static struct isl_obj vertices(struct isl_stream *s,
1658 struct isl_hash_table *table)
1659 {
1660 isl_ctx *ctx;
1661 struct isl_obj obj;
1662 struct isl_list *list = NULL;
1663 isl_union_set *uset = NULL;
1664 struct add_vertex_data data = { NULL };
1665
1666 obj = read_expr(s, table);
1667 ctx = isl_stream_get_ctx(s);
1668 obj = convert(ctx, obj, isl_obj_union_set);
1669 isl_assert(ctx, obj.type == isl_obj_union_set, goto error);
1670 uset = obj.v;
1671 obj.v = NULL;
1672
1673 list = isl_list_alloc(ctx, 0);
1674 if (!list)
1675 goto error;
1676
1677 data.list = list;
1678
1679 if (isl_union_set_foreach_set(uset, &set_vertices, &data) < 0)
1680 goto error;
1681
1682 isl_union_set_free(uset);
1683
1684 obj.type = isl_obj_list;
1685 obj.v = data.list;
1686
1687 return obj;
1688 error:
1689 isl_union_set_free(uset);
1690 isl_list_free(data.list);
1691 free_obj(obj);
1692 obj.type = isl_obj_none;
1693 obj.v = NULL;
1694 return obj;
1695 }
1696
1697 static struct isl_obj type_of(struct isl_stream *s,
1698 struct isl_hash_table *table)
1699 {
1700 struct isl_obj obj;
1701 const char *type = "unknown";
1702
1703 obj = read_expr(s, table);
1704
1705 if (obj.type == isl_obj_map ||
1706 obj.type == isl_obj_union_map)
1707 type = "map";
1708 if (obj.type == isl_obj_set ||
1709 obj.type == isl_obj_union_set)
1710 type = "set";
1711 if (obj.type == isl_obj_pw_multi_aff)
1712 type = "piecewise multi-quasiaffine expression";
1713 if (obj.type == isl_obj_pw_qpolynomial ||
1714 obj.type == isl_obj_union_pw_qpolynomial)
1715 type = "piecewise quasipolynomial";
1716 if (obj.type == isl_obj_pw_qpolynomial_fold ||
1717 obj.type == isl_obj_union_pw_qpolynomial_fold)
1718 type = "piecewise quasipolynomial fold";
1719 if (obj.type == isl_obj_list)
1720 type = "list";
1721 if (obj.type == isl_obj_bool)
1722 type = "boolean";
1723 if (obj.type == isl_obj_str)
1724 type = "string";
1725 if (obj.type == isl_obj_val)
1726 type = "value";
1727 if (obj.type == isl_obj_schedule)
1728 type = "schedule";
1729
1730 free_obj(obj);
1731 obj.type = isl_obj_str;
1732 obj.v = isl_str_from_string(isl_stream_get_ctx(s), strdup(type));
1733
1734 return obj;
1735 }
1736
1737 static __isl_give isl_union_set *read_set(struct isl_stream *s,
1738 struct isl_hash_table *table)
1739 {
1740 struct isl_obj obj;
1741 isl_ctx *ctx;
1742
1743 obj = read_obj(s, table);
1744 ctx = isl_stream_get_ctx(s);
1745 obj = convert(ctx, obj, isl_obj_union_set);
1746 isl_assert(ctx, obj.type == isl_obj_union_set, goto error);
1747 return obj.v;
1748 error:
1749 free_obj(obj);
1750 return NULL;
1751 }
1752
1753 static __isl_give isl_union_map *read_map(struct isl_stream *s,
1754 struct isl_hash_table *table)
1755 {
1756 struct isl_obj obj;
1757 isl_ctx *ctx;
1758
1759 obj = read_obj(s, table);
1760 ctx = isl_stream_get_ctx(s);
1761 obj = convert(ctx, obj, isl_obj_union_map);
1762 isl_assert(ctx, obj.type == isl_obj_union_map, goto error);
1763 return obj.v;
1764 error:
1765 free_obj(obj);
1766 return NULL;
1767 }
1768
1769 /* Read a schedule in the form of either a schedule (tree) or a union map
1770 * from "s" and store the schedule in "access".
1771 */
1772 static __isl_give isl_union_access_info *access_info_set_schedule(
1773 __isl_take isl_union_access_info *access, struct isl_stream *s,
1774 struct isl_hash_table *table)
1775 {
1776 struct isl_obj obj;
1777 isl_ctx *ctx;
1778
1779 obj = read_obj(s, table);
1780 if (obj.type == isl_obj_schedule)
1781 return isl_union_access_info_set_schedule(access, obj.v);
1782 ctx = isl_stream_get_ctx(s);
1783 obj = convert(ctx, obj, isl_obj_union_map);
1784
1785 return isl_union_access_info_set_schedule_map(access, obj.v);
1786 }
1787
1788 static struct isl_obj last_any(struct isl_stream *s,
1789 struct isl_hash_table *table, __isl_take isl_union_map *must_source,
1790 __isl_take isl_union_map *may_source)
1791 {
1792 struct isl_obj obj = { isl_obj_none, NULL };
1793 isl_union_access_info *access;
1794 isl_union_flow *flow;
1795 isl_union_map *sink = NULL;
1796 isl_union_map *may_dep;
1797
1798 if (isl_stream_eat(s, iscc_op[ISCC_BEFORE]))
1799 goto error;
1800
1801 sink = read_map(s, table);
1802 if (!sink)
1803 goto error;
1804
1805 if (isl_stream_eat(s, iscc_op[ISCC_UNDER]))
1806 goto error;
1807
1808 access = isl_union_access_info_from_sink(sink);
1809 access = isl_union_access_info_set_must_source(access, must_source);
1810 access = isl_union_access_info_set_may_source(access, may_source);
1811 access = access_info_set_schedule(access, s, table);
1812 flow = isl_union_access_info_compute_flow(access);
1813 may_dep = isl_union_flow_get_may_dependence(flow);
1814 isl_union_flow_free(flow);
1815
1816 if (!may_dep)
1817 return obj;
1818
1819 obj.type = isl_obj_union_map;
1820 obj.v = may_dep;
1821
1822 return obj;
1823 error:
1824 isl_union_map_free(may_source);
1825 isl_union_map_free(must_source);
1826 isl_union_map_free(sink);
1827 free_obj(obj);
1828 obj.type = isl_obj_none;
1829 obj.v = NULL;
1830 return obj;
1831 }
1832
1833 static struct isl_obj any(struct isl_stream *s, struct isl_hash_table *table)
1834 {
1835 struct isl_obj obj = { isl_obj_none, NULL };
1836 isl_union_access_info *access;
1837 isl_union_flow *flow;
1838 isl_union_map *may_source = NULL;
1839 isl_union_map *sink = NULL;
1840 isl_union_map *may_dep;
1841
1842 may_source = read_map(s, table);
1843 if (!may_source)
1844 goto error;
1845
1846 if (isl_stream_eat_if_available(s, iscc_op[ISCC_LAST])) {
1847 isl_union_map *must_source;
1848 must_source = read_map(s, table);
1849 if (!must_source)
1850 goto error;
1851 return last_any(s, table, must_source, may_source);
1852 }
1853
1854 if (isl_stream_eat(s, iscc_op[ISCC_BEFORE]))
1855 goto error;
1856
1857 sink = read_map(s, table);
1858 if (!sink)
1859 goto error;
1860
1861 if (isl_stream_eat(s, iscc_op[ISCC_UNDER]))
1862 goto error;
1863
1864 access = isl_union_access_info_from_sink(sink);
1865 access = isl_union_access_info_set_may_source(access, may_source);
1866 access = access_info_set_schedule(access, s, table);
1867 flow = isl_union_access_info_compute_flow(access);
1868 may_dep = isl_union_flow_get_may_dependence(flow);
1869 isl_union_flow_free(flow);
1870
1871 if (!may_dep)
1872 return obj;
1873
1874 obj.type = isl_obj_union_map;
1875 obj.v = may_dep;
1876
1877 return obj;
1878 error:
1879 isl_union_map_free(may_source);
1880 isl_union_map_free(sink);
1881 free_obj(obj);
1882 obj.type = isl_obj_none;
1883 obj.v = NULL;
1884 return obj;
1885 }
1886
1887 static struct isl_obj last(struct isl_stream *s, struct isl_hash_table *table)
1888 {
1889 struct isl_obj obj = { isl_obj_none, NULL };
1890 struct isl_list *list = NULL;
1891 isl_union_access_info *access;
1892 isl_union_flow *flow;
1893 isl_union_map *must_source = NULL;
1894 isl_union_map *sink = NULL;
1895 isl_union_map *must_dep;
1896 isl_union_map *must_no_source;
1897
1898 must_source = read_map(s, table);
1899 if (!must_source)
1900 goto error;
1901
1902 if (isl_stream_eat_if_available(s, iscc_op[ISCC_ANY])) {
1903 isl_union_map *may_source;
1904 may_source = read_map(s, table);
1905 if (!may_source)
1906 goto error;
1907 return last_any(s, table, must_source, may_source);
1908 }
1909
1910 list = isl_list_alloc(isl_stream_get_ctx(s), 2);
1911 if (!list)
1912 goto error;
1913
1914 if (isl_stream_eat(s, iscc_op[ISCC_BEFORE]))
1915 goto error;
1916
1917 sink = read_map(s, table);
1918 if (!sink)
1919 goto error;
1920
1921 if (isl_stream_eat(s, iscc_op[ISCC_UNDER]))
1922 goto error;
1923
1924 access = isl_union_access_info_from_sink(sink);
1925 access = isl_union_access_info_set_must_source(access, must_source);
1926 access = access_info_set_schedule(access, s, table);
1927 flow = isl_union_access_info_compute_flow(access);
1928 must_dep = isl_union_flow_get_must_dependence(flow);
1929 must_no_source = isl_union_flow_get_must_no_source(flow);
1930 isl_union_flow_free(flow);
1931
1932 list->obj[0].type = isl_obj_union_map;
1933 list->obj[0].v = must_dep;
1934 list->obj[1].type = isl_obj_union_map;
1935 list->obj[1].v = must_no_source;
1936
1937 if (!must_dep || !must_no_source) {
1938 isl_list_free(list);
1939 return obj;
1940 }
1941
1942 obj.v = list;
1943 obj.type = isl_obj_list;
1944
1945 return obj;
1946 error:
1947 isl_list_free(list);
1948 isl_union_map_free(must_source);
1949 isl_union_map_free(sink);
1950 free_obj(obj);
1951 obj.type = isl_obj_none;
1952 obj.v = NULL;
1953 return obj;
1954 }
1955
1956 static __isl_give isl_schedule *get_schedule(struct isl_stream *s,
1957 struct isl_hash_table *table)
1958 {
1959 isl_union_set *domain;
1960 isl_union_map *validity;
1961 isl_union_map *proximity;
1962
1963 domain = read_set(s, table);
1964 if (!domain)
1965 return NULL;
1966
1967 validity = isl_union_map_empty(isl_union_set_get_space(domain));
1968 proximity = isl_union_map_empty(isl_union_set_get_space(domain));
1969
1970 for (;;) {
1971 isl_union_map *umap;
1972 if (isl_stream_eat_if_available(s, iscc_op[ISCC_RESPECTING])) {
1973 umap = read_map(s, table);
1974 validity = isl_union_map_union(validity, umap);
1975 } else if (isl_stream_eat_if_available(s, iscc_op[ISCC_MINIMIZING])) {
1976 umap = read_map(s, table);
1977 proximity = isl_union_map_union(proximity, umap);
1978 } else
1979 break;
1980 }
1981
1982 return isl_union_set_compute_schedule(domain, validity, proximity);
1983 }
1984
1985 static struct isl_obj schedule(struct isl_stream *s,
1986 struct isl_hash_table *table)
1987 {
1988 struct isl_obj obj = { isl_obj_none, NULL };
1989 isl_schedule *schedule;
1990
1991 schedule = get_schedule(s, table);
1992
1993 obj.v = schedule;
1994 obj.type = isl_obj_schedule;
1995
1996 return obj;
1997 }
1998
1999 /* Read a schedule for code generation in the form of either
2000 * a schedule tree or a union map.
2001 */
2002 static struct isl_obj get_codegen_schedule(struct isl_stream *s,
2003 struct isl_hash_table *table)
2004 {
2005 struct isl_obj obj;
2006 isl_ctx *ctx;
2007
2008 obj = read_obj(s, table);
2009 ctx = isl_stream_get_ctx(s);
2010
2011 if (obj.type == isl_obj_schedule)
2012 return obj;
2013 if (is_subtype(obj, isl_obj_union_map))
2014 return convert(ctx, obj, isl_obj_union_map);
2015
2016 free_obj(obj);
2017 obj.v = NULL;
2018 obj.type = isl_obj_none;
2019 isl_die(ctx, isl_error_invalid, "expecting schedule or map",
2020 return obj);
2021 }
2022
2023 /* Generate an AST for the given schedule and options and return the AST.
2024 */
2025 static __isl_give isl_ast_node *get_ast_from_union_map(
2026 __isl_take isl_union_map *schedule, __isl_take isl_union_map *options)
2027 {
2028 isl_space *space;
2029 isl_set *context;
2030 isl_ast_build *build;
2031 isl_ast_node *tree;
2032
2033 space = isl_union_map_get_space(schedule);
2034 context = isl_set_universe(isl_space_params(space));
2035
2036 build = isl_ast_build_from_context(context);
2037 build = isl_ast_build_set_options(build, options);
2038 tree = isl_ast_build_ast_from_schedule(build, schedule);
2039 isl_ast_build_free(build);
2040
2041 return tree;
2042 }
2043
2044 /* Generate an AST for the given schedule and return the AST.
2045 */
2046 static __isl_give isl_ast_node *get_ast_from_schedule(
2047 __isl_take isl_schedule *schedule)
2048 {
2049 isl_ast_build *build;
2050 isl_ast_node *tree;
2051
2052 build = isl_ast_build_alloc(isl_schedule_get_ctx(schedule));
2053 tree = isl_ast_build_node_from_schedule(build, schedule);
2054 isl_ast_build_free(build);
2055
2056 return tree;
2057 }
2058
2059 /* Print the AST "tree" on the printer "p".
2060 */
2061 static __isl_give isl_printer *print_ast(__isl_take isl_printer *p,
2062 __isl_take isl_ast_node *tree)
2063 {
2064 int format;
2065
2066 format = isl_printer_get_output_format(p);
2067 p = isl_printer_set_output_format(p, ISL_FORMAT_C);
2068 p = isl_printer_print_ast_node(p, tree);
2069 p = isl_printer_set_output_format(p, format);
2070
2071 isl_ast_node_free(tree);
2072
2073 return p;
2074 }
2075
2076 /* Perform the codegen operation.
2077 * In particular, read a schedule, check if the user has specified any options
2078 * and then generate an AST from the schedule (and options) and print it.
2079 * In case the schedule is specified as a schedule tree, the AST generation
2080 * options are embedded in the schedule, so they are not read in separately.
2081 */
2082 static __isl_give isl_printer *codegen(struct isl_stream *s,
2083 struct isl_hash_table *table, __isl_take isl_printer *p)
2084 {
2085 struct isl_obj obj;
2086 isl_ast_node *tree;
2087
2088 obj = get_codegen_schedule(s, table);
2089 if (!obj.v)
2090 return p;
2091
2092 if (obj.type == isl_obj_schedule) {
2093 isl_schedule *schedule = obj.v;
2094
2095 tree = get_ast_from_schedule(schedule);
2096 } else {
2097 isl_union_map *schedule = obj.v;
2098 isl_union_map *options;
2099
2100 if (isl_stream_eat_if_available(s, iscc_op[ISCC_USING]))
2101 options = read_map(s, table);
2102 else
2103 options = isl_union_map_empty(
2104 isl_union_map_get_space(schedule));
2105
2106 tree = get_ast_from_union_map(schedule, options);
2107 }
2108
2109 if (tree)
2110 p = print_ast(p, tree);
2111
2112 isl_stream_eat(s, ';');
2113
2114 return p;
2115 }
2116
2117 static struct isl_obj power(struct isl_stream *s, struct isl_obj obj)
2118 {
2119 struct isl_token *tok;
2120 isl_ctx *ctx;
2121 isl_val *v;
2122
2123 ctx = isl_stream_get_ctx(s);
2124 if (isl_stream_eat_if_available(s, '+'))
2125 return transitive_closure(ctx, obj);
2126
2127 isl_assert(ctx, is_subtype(obj, isl_obj_union_map), goto error);
2128 if (obj.type != isl_obj_union_map)
2129 obj = convert(ctx, obj, isl_obj_union_map);
2130
2131 tok = isl_stream_next_token(s);
2132 if (!tok || isl_token_get_type(tok) != ISL_TOKEN_VALUE) {
2133 isl_stream_error(s, tok, "expecting integer exponent");
2134 if (tok)
2135 isl_stream_push_token(s, tok);
2136 goto error;
2137 }
2138
2139 v = isl_token_get_val(ctx, tok);
2140 if (isl_val_is_zero(v)) {
2141 isl_stream_error(s, tok, "expecting non-zero exponent");
2142 isl_val_free(v);
2143 if (tok)
2144 isl_stream_push_token(s, tok);
2145 goto error;
2146 }
2147
2148 obj.v = isl_union_map_fixed_power_val(obj.v, v);
2149 isl_token_free(tok);
2150 if (!obj.v)
2151 goto error;
2152
2153 return obj;
2154 error:
2155 free_obj(obj);
2156 obj.type = isl_obj_none;
2157 obj.v = NULL;
2158 return obj;
2159 }
2160
2161 static struct isl_obj check_assert(struct isl_stream *s,
2162 struct isl_hash_table *table)
2163 {
2164 struct isl_obj obj;
2165 isl_ctx *ctx;
2166
2167 obj = read_expr(s, table);
2168 ctx = isl_stream_get_ctx(s);
2169 if (obj.type != isl_obj_bool)
2170 isl_die(ctx, isl_error_invalid,
2171 "expecting boolean expression", goto error);
2172 if (obj.v != &iscc_bool_true)
2173 isl_die(ctx, isl_error_unknown,
2174 "assertion failed", abort());
2175 error:
2176 free_obj(obj);
2177 obj.type = isl_obj_none;
2178 obj.v = NULL;
2179 return obj;
2180 }
2181
2182 static struct isl_obj read_from_file(struct isl_stream *s)
2183 {
2184 isl_ctx *ctx;
2185 struct isl_obj obj;
2186 struct isl_token *tok;
2187 struct isl_stream *s_file;
2188 struct iscc_options *options;
2189 char *name;
2190 FILE *file;
2191
2192 tok = isl_stream_next_token(s);
2193 if (!tok || isl_token_get_type(tok) != ISL_TOKEN_STRING) {
2194 isl_stream_error(s, tok, "expecting filename");
2195 isl_token_free(tok);
2196 goto error;
2197 }
2198
2199 ctx = isl_stream_get_ctx(s);
2200 options = isl_ctx_peek_iscc_options(ctx);
2201 if (!options || !options->io) {
2202 isl_token_free(tok);
2203 isl_die(ctx, isl_error_invalid,
2204 "read operation not allowed", goto error);
2205 }
2206
2207 name = isl_token_get_str(ctx, tok);
2208 isl_token_free(tok);
2209 file = fopen(name, "r");
2210 free(name);
2211 isl_assert(ctx, file, goto error);
2212
2213 s_file = isl_stream_new_file(ctx, file);
2214 if (!s_file) {
2215 fclose(file);
2216 goto error;
2217 }
2218
2219 obj = isl_stream_read_obj(s_file);
2220
2221 isl_stream_free(s_file);
2222 fclose(file);
2223
2224 return obj;
2225 error:
2226 obj.type = isl_obj_none;
2227 obj.v = NULL;
2228 return obj;
2229 }
2230
2231 static struct isl_obj write_to_file(struct isl_stream *s,
2232 struct isl_hash_table *table)
2233 {
2234 struct isl_obj obj = { isl_obj_none, NULL };
2235 struct isl_token *tok;
2236 struct iscc_options *options;
2237 char *name;
2238 FILE *file;
2239 isl_ctx *ctx;
2240 isl_printer *p;
2241
2242 tok = isl_stream_next_token(s);
2243 if (!tok || isl_token_get_type(tok) != ISL_TOKEN_STRING) {
2244 isl_stream_error(s, tok, "expecting filename");
2245 isl_token_free(tok);
2246 goto error;
2247 }
2248
2249 obj = read_expr(s, table);
2250
2251 ctx = isl_stream_get_ctx(s);
2252 options = isl_ctx_peek_iscc_options(ctx);
2253 if (!options || !options->io) {
2254 isl_token_free(tok);
2255 isl_die(ctx, isl_error_invalid,
2256 "write operation not allowed", goto error);
2257 }
2258
2259 name = isl_token_get_str(ctx, tok);
2260 isl_token_free(tok);
2261 file = fopen(name, "w");
2262 free(name);
2263 if (!file)
2264 isl_die(ctx, isl_error_unknown,
2265 "could not open file for writing", goto error);
2266
2267 p = isl_printer_to_file(ctx, file);
2268 p = isl_printer_set_output_format(p, options->format);
2269 p = obj.type->print(p, obj.v);
2270 p = isl_printer_end_line(p);
2271 isl_printer_free(p);
2272
2273 fclose(file);
2274 error:
2275 free_obj(obj);
2276 obj.type = isl_obj_none;
2277 obj.v = NULL;
2278 return obj;
2279 }
2280
2281 static struct isl_obj read_string_if_available(struct isl_stream *s)
2282 {
2283 struct isl_token *tok;
2284 struct isl_obj obj = { isl_obj_none, NULL };
2285
2286 tok = isl_stream_next_token(s);
2287 if (!tok)
2288 return obj;
2289 if (isl_token_get_type(tok) == ISL_TOKEN_STRING) {
2290 isl_str *str;
2291 str = isl_str_alloc(isl_stream_get_ctx(s));
2292 if (!str)
2293 goto error;
2294 str->s = isl_token_get_str(isl_stream_get_ctx(s), tok);
2295 isl_token_free(tok);
2296 obj.v = str;
2297 obj.type = isl_obj_str;
2298 } else
2299 isl_stream_push_token(s, tok);
2300 return obj;
2301 error:
2302 isl_token_free(tok);
2303 return obj;
2304 }
2305
2306 static struct isl_obj read_bool_if_available(struct isl_stream *s)
2307 {
2308 struct isl_token *tok;
2309 struct isl_obj obj = { isl_obj_none, NULL };
2310 int type;
2311
2312 tok = isl_stream_next_token(s);
2313 if (!tok)
2314 return obj;
2315 type = isl_token_get_type(tok);
2316 if (type == ISL_TOKEN_FALSE || type == ISL_TOKEN_TRUE) {
2317 int is_true = type == ISL_TOKEN_TRUE;
2318 isl_token_free(tok);
2319 obj.v = is_true ? &iscc_bool_true : &iscc_bool_false;
2320 obj.type = isl_obj_bool;
2321 } else
2322 isl_stream_push_token(s, tok);
2323 return obj;
2324 }
2325
2326 static __isl_give char *read_ident(struct isl_stream *s)
2327 {
2328 char *name;
2329 isl_val *v;
2330 struct isl_token *tok, *tok2;
2331
2332 name = isl_stream_read_ident_if_available(s);
2333 if (name)
2334 return name;
2335
2336 tok = isl_stream_next_token(s);
2337 if (!tok)
2338 return NULL;
2339 if (isl_token_get_type(tok) != '$') {
2340 isl_stream_push_token(s, tok);
2341 return NULL;
2342 }
2343 tok2 = isl_stream_next_token(s);
2344 if (!tok2 || isl_token_get_type(tok2) != ISL_TOKEN_VALUE) {
2345 if (tok2)
2346 isl_stream_push_token(s, tok2);
2347 isl_stream_push_token(s, tok);
2348 return NULL;
2349 }
2350
2351 v = isl_token_get_val(isl_stream_get_ctx(s), tok2);
2352 name = isl_val_to_str(v);
2353 isl_val_free(v);
2354 isl_token_free(tok);
2355 isl_token_free(tok2);
2356
2357 return name;
2358 }
2359
2360 static struct isl_obj read_list(struct isl_stream *s,
2361 struct isl_hash_table *table, struct isl_obj obj)
2362 {
2363 struct isl_list *list;
2364
2365 list = isl_list_alloc(isl_stream_get_ctx(s), 2);
2366 if (!list)
2367 goto error;
2368 list->obj[0] = obj;
2369 list->obj[1] = read_obj(s, table);
2370 obj.v = list;
2371 obj.type = isl_obj_list;
2372
2373 if (!list->obj[1].v)
2374 goto error;
2375
2376 while (isl_stream_eat_if_available(s, ',')) {
2377 obj.v = list = isl_list_add_obj(list, read_obj(s, table));
2378 if (!obj.v)
2379 goto error;
2380 }
2381
2382 return obj;
2383 error:
2384 free_obj(obj);
2385 obj.type = isl_obj_none;
2386 obj.v = NULL;
2387 return obj;
2388 }
2389
2390 static struct isl_obj read_obj(struct isl_stream *s,
2391 struct isl_hash_table *table)
2392 {
2393 isl_ctx *ctx;
2394 struct isl_obj obj = { isl_obj_none, NULL };
2395 char *name = NULL;
2396 struct isc_un_op *op = NULL;
2397
2398 obj = read_string_if_available(s);
2399 if (obj.v)
2400 return obj;
2401 obj = read_bool_if_available(s);
2402 if (obj.v)
2403 return obj;
2404 ctx = isl_stream_get_ctx(s);
2405 if (isl_stream_eat_if_available(s, '(')) {
2406 if (isl_stream_next_token_is(s, ')')) {
2407 obj.type = isl_obj_list;
2408 obj.v = isl_list_alloc(ctx, 0);
2409 } else {
2410 obj = read_expr(s, table);
2411 if (obj.v && isl_stream_eat_if_available(s, ','))
2412 obj = read_list(s, table, obj);
2413 }
2414 if (!obj.v || isl_stream_eat(s, ')'))
2415 goto error;
2416 } else {
2417 op = read_prefix_un_op_if_available(s);
2418 if (op)
2419 return read_un_op_expr(s, table, op);
2420
2421 if (isl_stream_eat_if_available(s, iscc_op[ISCC_ASSERT]))
2422 return check_assert(s, table);
2423 if (isl_stream_eat_if_available(s, iscc_op[ISCC_READ]))
2424 return read_from_file(s);
2425 if (isl_stream_eat_if_available(s, iscc_op[ISCC_WRITE]))
2426 return write_to_file(s, table);
2427 if (isl_stream_eat_if_available(s, iscc_op[ISCC_VERTICES]))
2428 return vertices(s, table);
2429 if (isl_stream_eat_if_available(s, iscc_op[ISCC_ANY]))
2430 return any(s, table);
2431 if (isl_stream_eat_if_available(s, iscc_op[ISCC_LAST]))
2432 return last(s, table);
2433 if (isl_stream_eat_if_available(s, iscc_op[ISCC_SCHEDULE]))
2434 return schedule(s, table);
2435 if (isl_stream_eat_if_available(s, iscc_op[ISCC_TYPEOF]))
2436 return type_of(s, table);
2437
2438 name = read_ident(s);
2439 if (name)
2440 obj = stored_obj(ctx, table, name);
2441 else
2442 obj = isl_stream_read_obj(s);
2443 if (!obj.v)
2444 goto error;
2445 }
2446
2447 if (isl_stream_eat_if_available(s, '^'))
2448 obj = power(s, obj);
2449 else if (obj.type == isl_obj_list && isl_stream_eat_if_available(s, '['))
2450 obj = obj_at_index(s, obj);
2451 else if (is_subtype(obj, isl_obj_union_map) &&
2452 isl_stream_eat_if_available(s, '(')) {
2453 obj = convert(ctx, obj, isl_obj_union_map);
2454 obj = apply(s, obj.v, table);
2455 } else if (is_subtype(obj, isl_obj_union_pw_qpolynomial) &&
2456 isl_stream_eat_if_available(s, '(')) {
2457 obj = convert(ctx, obj, isl_obj_union_pw_qpolynomial);
2458 obj = apply_fun(s, obj, table);
2459 } else if (is_subtype(obj, isl_obj_union_pw_qpolynomial_fold) &&
2460 isl_stream_eat_if_available(s, '(')) {
2461 obj = convert(ctx, obj, isl_obj_union_pw_qpolynomial_fold);
2462 obj = apply_fun(s, obj, table);
2463 }
2464
2465 return obj;
2466 error:
2467 free_obj(obj);
2468 obj.type = isl_obj_none;
2469 obj.v = NULL;
2470 return obj;
2471 }
2472
2473 static struct isc_bin_op *find_matching_bin_op(struct isc_bin_op *like,
2474 struct isl_obj lhs, struct isl_obj rhs)
2475 {
2476 int i;
2477
2478 for (i = 0; ; ++i) {
2479 if (!bin_ops[i].op)
2480 break;
2481 if (bin_ops[i].op != like->op)
2482 continue;
2483 if (!is_subtype(lhs, bin_ops[i].lhs))
2484 continue;
2485 if (!is_subtype(rhs, bin_ops[i].rhs))
2486 continue;
2487
2488 return &bin_ops[i];
2489 }
2490
2491 for (i = 0; ; ++i) {
2492 if (!named_bin_ops[i].name)
2493 break;
2494 if (named_bin_ops[i].op.op != like->op)
2495 continue;
2496 if (!is_subtype(lhs, named_bin_ops[i].op.lhs))
2497 continue;
2498 if (!is_subtype(rhs, named_bin_ops[i].op.rhs))
2499 continue;
2500
2501 return &named_bin_ops[i].op;
2502 }
2503
2504 for (i = 0; ; ++i) {
2505 if (!compound_bin_ops[i].full)
2506 break;
2507 if (compound_bin_ops[i].op.op != like->op)
2508 continue;
2509 if (!is_subtype(lhs, compound_bin_ops[i].op.lhs))
2510 continue;
2511 if (!is_subtype(rhs, compound_bin_ops[i].op.rhs))
2512 continue;
2513
2514 return &compound_bin_ops[i].op;
2515 }
2516
2517 return NULL;
2518 }
2519
2520 static int next_is_neg_int(struct isl_stream *s)
2521 {
2522 struct isl_token *tok;
2523 int ret;
2524
2525 tok = isl_stream_next_token(s);
2526 if (tok && isl_token_get_type(tok) == ISL_TOKEN_VALUE) {
2527 isl_val *v;
2528 v = isl_token_get_val(isl_stream_get_ctx(s), tok);
2529 ret = isl_val_is_neg(v);
2530 isl_val_free(v);
2531 } else
2532 ret = 0;
2533 isl_stream_push_token(s, tok);
2534
2535 return ret;
2536 }
2537
2538 static struct isl_obj call_bin_op(isl_ctx *ctx, struct isc_bin_op *op,
2539 struct isl_obj lhs, struct isl_obj rhs)
2540 {
2541 struct isl_obj obj;
2542
2543 lhs = convert(ctx, lhs, op->lhs);
2544 rhs = convert(ctx, rhs, op->rhs);
2545 if (op->res != isl_obj_bool)
2546 obj.v = op->o.fn(lhs.v, rhs.v);
2547 else {
2548 int res = op->o.test(lhs.v, rhs.v);
2549 free_obj(lhs);
2550 free_obj(rhs);
2551 obj.v = iscc_bool_from_int(res);
2552 }
2553 obj.type = op->res;
2554
2555 return obj;
2556 }
2557
2558 static struct isl_obj read_expr(struct isl_stream *s,
2559 struct isl_hash_table *table)
2560 {
2561 isl_ctx *ctx;
2562 struct isl_obj obj = { isl_obj_none, NULL };
2563 struct isl_obj right_obj = { isl_obj_none, NULL };
2564
2565 obj = read_obj(s, table);
2566 ctx = isl_stream_get_ctx(s);
2567 for (; obj.v;) {
2568 struct isc_bin_op *op = NULL;
2569
2570 op = read_bin_op_if_available(s, obj);
2571 if (!op)
2572 break;
2573
2574 right_obj = read_obj(s, table);
2575
2576 op = find_matching_bin_op(op, obj, right_obj);
2577
2578 if (!op)
2579 isl_die(ctx, isl_error_invalid,
2580 "no such binary operator defined on given operands",
2581 goto error);
2582
2583 obj = call_bin_op(ctx, op, obj, right_obj);
2584 }
2585
2586 if (obj.type == isl_obj_val && next_is_neg_int(s)) {
2587 right_obj = read_obj(s, table);
2588 obj.v = isl_val_add(obj.v, right_obj.v);
2589 }
2590
2591 return obj;
2592 error:
2593 free_obj(right_obj);
2594 free_obj(obj);
2595 obj.type = isl_obj_none;
2596 obj.v = NULL;
2597 return obj;
2598 }
2599
2600 static __isl_give isl_printer *source_file(struct isl_stream *s,
2601 struct isl_hash_table *table, __isl_take isl_printer *p);
2602
2603 /* Print "obj" to the printer "p".
2604 * If the object is a schedule, then print it in block format.
2605 */
2606 static __isl_give isl_printer *print_obj(__isl_take isl_printer *p,
2607 struct isl_obj obj)
2608 {
2609 if (obj.type != isl_obj_schedule)
2610 return obj.type->print(p, obj.v);
2611
2612 p = isl_printer_set_yaml_style(p, ISL_YAML_STYLE_BLOCK);
2613 p = obj.type->print(p, obj.v);
2614 p = isl_printer_set_yaml_style(p, ISL_YAML_STYLE_FLOW);
2615
2616 return p;
2617 }
2618
2619 static __isl_give isl_printer *read_line(struct isl_stream *s,
2620 struct isl_hash_table *table, __isl_take isl_printer *p, int tty)
2621 {
2622 isl_ctx *ctx;
2623 struct isl_obj obj = { isl_obj_none, NULL };
2624 char *lhs = NULL;
2625 int assign = 0;
2626 int only_print = 0;
2627 char buf[30];
2628
2629 if (!p)
2630 return NULL;
2631 if (isl_stream_is_empty(s))
2632 return p;
2633
2634 if (isl_stream_eat_if_available(s, iscc_op[ISCC_SOURCE]))
2635 return source_file(s, table, p);
2636 if (isl_stream_eat_if_available(s, iscc_op[ISCC_CODEGEN]))
2637 return codegen(s, table, p);
2638
2639 assign = is_assign(s);
2640 if (assign) {
2641 lhs = isl_stream_read_ident_if_available(s);
2642 if (isl_stream_eat(s, ISL_TOKEN_DEF))
2643 goto error;
2644 } else if (isl_stream_eat_if_available(s, iscc_op[ISCC_PRINT]))
2645 only_print = 1;
2646 else if (!tty)
2647 only_print = 1;
2648
2649 obj = read_expr(s, table);
2650 ctx = isl_stream_get_ctx(s);
2651 if (isl_stream_eat(s, ';'))
2652 goto error;
2653
2654 if (only_print) {
2655 if (obj.type != isl_obj_none && obj.v != NULL) {
2656 p = print_obj(p, obj);
2657 p = isl_printer_end_line(p);
2658 }
2659 free_obj(obj);
2660 return p;
2661 }
2662 if (!assign && obj.type != isl_obj_none && obj.v != NULL) {
2663 static int count = 0;
2664 snprintf(buf, sizeof(buf), "$%d", count++);
2665 lhs = strdup(buf + 1);
2666
2667 p = isl_printer_print_str(p, buf);
2668 p = isl_printer_print_str(p, " := ");
2669 p = obj.type->print(p, obj.v);
2670 p = isl_printer_end_line(p);
2671 }
2672 if (lhs && do_assign(ctx, table, lhs, obj))
2673 return p;
2674
2675 return p;
2676 error:
2677 isl_stream_flush_tokens(s);
2678 isl_stream_skip_line(s);
2679 free(lhs);
2680 free_obj(obj);
2681 return p;
2682 }
2683
2684 static isl_stat free_cb(void **entry, void *user)
2685 {
2686 struct isl_named_obj *named = *entry;
2687
2688 free_obj(named->obj);
2689 free(named->name);
2690 free(named);
2691
2692 return isl_stat_ok;
2693 }
2694
2695 static void register_named_ops(struct isl_stream *s)
2696 {
2697 int i;
2698
2699 for (i = 0; i < ISCC_N_OP; ++i) {
2700 iscc_op[i] = isl_stream_register_keyword(s, op_name[i]);
2701 assert(iscc_op[i] != ISL_TOKEN_ERROR);
2702 }
2703
2704 for (i = 0; ; ++i) {
2705 if (!named_un_ops[i].name)
2706 break;
2707 named_un_ops[i].op.op = isl_stream_register_keyword(s,
2708 named_un_ops[i].name);
2709 assert(named_un_ops[i].op.op != ISL_TOKEN_ERROR);
2710 }
2711
2712 for (i = 0; ; ++i) {
2713 if (!named_bin_ops[i].name)
2714 break;
2715 named_bin_ops[i].op.op = isl_stream_register_keyword(s,
2716 named_bin_ops[i].name);
2717 assert(named_bin_ops[i].op.op != ISL_TOKEN_ERROR);
2718 }
2719
2720 for (i = 0; ; ++i) {
2721 if (!compound_bin_ops[i].full)
2722 break;
2723 compound_bin_ops[i].op.op = isl_stream_register_keyword(s,
2724 compound_bin_ops[i].full);
2725 assert(compound_bin_ops[i].op.op != ISL_TOKEN_ERROR);
2726 }
2727 }
2728
2729 static __isl_give isl_printer *source_file(struct isl_stream *s,
2730 struct isl_hash_table *table, __isl_take isl_printer *p)
2731 {
2732 isl_ctx *ctx;
2733 struct isl_token *tok;
2734 struct isl_stream *s_file;
2735 struct iscc_options *options;
2736 char *name;
2737 FILE *file;
2738
2739 tok = isl_stream_next_token(s);
2740 if (!tok || isl_token_get_type(tok) != ISL_TOKEN_STRING) {
2741 isl_stream_error(s, tok, "expecting filename");
2742 isl_token_free(tok);
2743 return p;
2744 }
2745
2746 isl_stream_eat(s, ';');
2747
2748 ctx = isl_stream_get_ctx(s);
2749 options = isl_ctx_peek_iscc_options(ctx);
2750 if (!options || !options->io) {
2751 isl_token_free(tok);
2752 isl_die(ctx, isl_error_invalid,
2753 "source operation not allowed", return p);
2754 }
2755
2756 name = isl_token_get_str(ctx, tok);
2757 isl_token_free(tok);
2758 file = fopen(name, "r");
2759 free(name);
2760 isl_assert(ctx, file, return p);
2761
2762 s_file = isl_stream_new_file(ctx, file);
2763 if (!s_file) {
2764 fclose(file);
2765 return p;
2766 }
2767
2768 register_named_ops(s_file);
2769
2770 while (!isl_stream_is_empty(s_file))
2771 p = read_line(s_file, table, p, 0);
2772
2773 isl_stream_free(s_file);
2774 fclose(file);
2775
2776 return p;
2777 }
2778
2779 int main(int argc, char **argv)
2780 {
2781 struct isl_ctx *ctx;
2782 struct isl_stream *s;
2783 struct isl_hash_table *table;
2784 struct iscc_options *options;
2785 isl_printer *p;
2786 int tty = isatty(0);
2787
2788 options = iscc_options_new_with_defaults();
2789 assert(options);
2790
2791 ctx = isl_ctx_alloc_with_options(&iscc_options_args, options);
2792 pet_options_set_autodetect(ctx, 1);
2793 pet_options_set_encapsulate_dynamic_control(ctx, 1);
2794 argc = isl_ctx_parse_options(ctx, argc, argv, ISL_ARG_ALL);
2795 s = isl_stream_new_file(ctx, stdin);
2796 assert(s);
2797 table = isl_hash_table_alloc(ctx, 10);
2798 assert(table);
2799 p = isl_printer_to_file(ctx, stdout);
2800 p = isl_printer_set_output_format(p, options->format);
2801 assert(p);
2802
2803 register_named_ops(s);
2804
2805 install_signal_handler(ctx);
2806
2807 while (p) {
2808 int empty;
2809
2810 empty = isl_stream_is_empty(s);
2811 if (empty && isl_ctx_last_error(ctx) != isl_error_abort)
2812 break;
2813 if (!empty)
2814 p = read_line(s, table, p, tty);
2815 if (isl_ctx_last_error(ctx) == isl_error_abort) {
2816 fprintf(stderr, "Interrupted\n");
2817 isl_ctx_resume(ctx);
2818 isl_ctx_reset_error(ctx);
2819 }
2820 }
2821
2822 remove_signal_handler(ctx);
2823
2824 isl_printer_free(p);
2825 isl_hash_table_foreach(ctx, table, free_cb, NULL);
2826 isl_hash_table_free(ctx, table);
2827 isl_stream_free(s);
2828 isl_ctx_free(ctx);
2829
2830 return 0;
2831 }
lattice point enumeration library