1 /* Graphite polyhedral representation.
2 Copyright (C) 2009, 2010 Free Software Foundation, Inc.
3 Contributed by Sebastian Pop <sebastian.pop@amd.com> and
4 Tobias Grosser <grosser@fim.uni-passau.de>.
6 This file is part of GCC.
8 GCC is free software; you can redistribute it and/or modify
9 it under the terms of the GNU General Public License as published by
10 the Free Software Foundation; either version 3, or (at your option)
13 GCC is distributed in the hope that it will be useful,
14 but WITHOUT ANY WARRANTY; without even the implied warranty of
15 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 GNU General Public License for more details.
18 You should have received a copy of the GNU General Public License
19 along with GCC; see the file COPYING3. If not see
20 <http://www.gnu.org/licenses/>. */
23 #include "coretypes.h"
29 #include "basic-block.h"
30 #include "diagnostic.h"
31 #include "tree-pretty-print.h"
32 #include "gimple-pretty-print.h"
33 #include "tree-flow.h"
35 #include "tree-dump.h"
38 #include "tree-chrec.h"
39 #include "tree-data-ref.h"
40 #include "tree-scalar-evolution.h"
41 #include "tree-pass.h"
43 #include "value-prof.h"
44 #include "pointer-set.h"
51 #include "graphite-ppl.h"
53 #include "graphite-poly.h"
54 #include "graphite-dependences.h"
55 #include "graphite-cloog-util.h"
57 #define OPENSCOP_MAX_STRING 256
59 /* Return the maximal loop depth in SCOP. */
62 scop_max_loop_depth (scop_p scop
)
68 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
70 int nb_loops
= pbb_dim_iter_domain (pbb
);
71 if (max_nb_loops
< nb_loops
)
72 max_nb_loops
= nb_loops
;
78 /* Extend the scattering matrix of PBB to MAX_SCATTERING scattering
82 extend_scattering (poly_bb_p pbb
, int max_scattering
)
84 ppl_dimension_type nb_old_dims
, nb_new_dims
;
86 ppl_Coefficient_t coef
;
89 nb_added_dims
= max_scattering
- pbb_nb_scattering_transform (pbb
);
92 ppl_new_Coefficient (&coef
);
93 ppl_assign_Coefficient_from_mpz_t (coef
, one
);
95 gcc_assert (nb_added_dims
>= 0);
97 nb_old_dims
= pbb_nb_scattering_transform (pbb
) + pbb_dim_iter_domain (pbb
)
98 + scop_nb_params (PBB_SCOP (pbb
));
99 nb_new_dims
= nb_old_dims
+ nb_added_dims
;
101 ppl_insert_dimensions (PBB_TRANSFORMED_SCATTERING (pbb
),
102 pbb_nb_scattering_transform (pbb
), nb_added_dims
);
103 PBB_NB_SCATTERING_TRANSFORM (pbb
) += nb_added_dims
;
105 /* Add identity matrix for the added dimensions. */
106 for (i
= max_scattering
- nb_added_dims
; i
< max_scattering
; i
++)
108 ppl_Constraint_t cstr
;
109 ppl_Linear_Expression_t expr
;
111 ppl_new_Linear_Expression_with_dimension (&expr
, nb_new_dims
);
112 ppl_Linear_Expression_add_to_coefficient (expr
, i
, coef
);
113 ppl_new_Constraint (&cstr
, expr
, PPL_CONSTRAINT_TYPE_EQUAL
);
114 ppl_Polyhedron_add_constraint (PBB_TRANSFORMED_SCATTERING (pbb
), cstr
);
115 ppl_delete_Constraint (cstr
);
116 ppl_delete_Linear_Expression (expr
);
119 ppl_delete_Coefficient (coef
);
123 /* All scattering matrices in SCOP will have the same number of scattering
127 unify_scattering_dimensions (scop_p scop
)
131 graphite_dim_t max_scattering
= 0;
133 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
134 max_scattering
= MAX (pbb_nb_scattering_transform (pbb
), max_scattering
);
136 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
137 extend_scattering (pbb
, max_scattering
);
139 return max_scattering
;
142 /* Print to FILE the pdr PH in OpenScop format. NB_SUBSCRIPTS is the number
143 of subscripts in PH, ALIAS_SET_DIM is the dimension of the alias set and
144 NB_PARAMS is the number of parameters in PH. */
147 openscop_print_pdr_polyhedron (FILE *file
, ppl_const_Polyhedron_t ph
,
148 int nb_subscripts
, int alias_set_dimension
,
151 int input
, locals
, output
;
152 ppl_dimension_type alias_set_dim
= (ppl_dimension_type
) alias_set_dimension
;
153 ppl_dimension_type sub_dim_last
= alias_set_dim
+ nb_subscripts
;
154 ppl_dimension_type
*map
, i
, ph_space_dim
= sub_dim_last
+ 1;
155 ppl_Polyhedron_t pph
;
157 ppl_new_C_Polyhedron_from_C_Polyhedron (&pph
, ph
);
159 map
= (ppl_dimension_type
*) XNEWVEC (ppl_dimension_type
, ph_space_dim
);
161 for (i
= 0; i
< alias_set_dim
- 1; i
++)
162 map
[i
] = nb_subscripts
+ 1 + i
;
164 for (i
= alias_set_dim
- 1; i
< sub_dim_last
; i
++)
165 map
[i
] = i
- alias_set_dim
+ 1;
167 ppl_Polyhedron_map_space_dimensions (pph
, map
, ph_space_dim
- 1);
170 input
= alias_set_dim
- nb_params
- 1;
172 /* According to OpenScop specification, the alias set column is a part of
173 the output columns. */
174 output
= nb_subscripts
+ 1;
176 openscop_print_polyhedron_matrix (file
, pph
, output
, input
, locals
, nb_params
);
179 /* Print to FILE the powerset PDR. NB_SUBSCRIPTS is the number of subscripts
180 in PDR, ALIAS_SET_DIM is the dimension of the alias set in PDR and
181 NB_PARAMS is the number of parameters in PDR. */
184 openscop_print_pdr_powerset (FILE *file
,
185 ppl_Pointset_Powerset_C_Polyhedron_t ps
,
191 ppl_Pointset_Powerset_C_Polyhedron_iterator_t it
, end
;
193 ppl_new_Pointset_Powerset_C_Polyhedron_iterator (&it
);
194 ppl_new_Pointset_Powerset_C_Polyhedron_iterator (&end
);
196 ppl_Pointset_Powerset_C_Polyhedron_size (ps
, &nb_disjuncts
);
197 fprintf (file
, "%d\n", (int) nb_disjuncts
);
199 for (ppl_Pointset_Powerset_C_Polyhedron_iterator_begin (ps
, it
),
200 ppl_Pointset_Powerset_C_Polyhedron_iterator_end (ps
, end
);
201 !ppl_Pointset_Powerset_C_Polyhedron_iterator_equal_test (it
, end
);
202 ppl_Pointset_Powerset_C_Polyhedron_iterator_increment (it
))
204 ppl_const_Polyhedron_t ph
;
206 ppl_Pointset_Powerset_C_Polyhedron_iterator_dereference (it
, &ph
);
207 openscop_print_pdr_polyhedron (file
, ph
, nb_subscripts
, alias_set_dim
,
211 ppl_delete_Pointset_Powerset_C_Polyhedron_iterator (it
);
212 ppl_delete_Pointset_Powerset_C_Polyhedron_iterator (end
);
215 /* Print to FILE the powerset PS in its OpenScop matrix form. */
218 openscop_print_powerset_matrix (FILE *file
,
219 ppl_Pointset_Powerset_C_Polyhedron_t ps
,
220 int output
, int input
, int locals
,
224 ppl_Pointset_Powerset_C_Polyhedron_iterator_t it
, end
;
226 ppl_new_Pointset_Powerset_C_Polyhedron_iterator (&it
);
227 ppl_new_Pointset_Powerset_C_Polyhedron_iterator (&end
);
229 ppl_Pointset_Powerset_C_Polyhedron_size (ps
, &nb_disjuncts
);
230 fprintf (file
, "%d\n", (int) nb_disjuncts
);
232 for (ppl_Pointset_Powerset_C_Polyhedron_iterator_begin (ps
, it
),
233 ppl_Pointset_Powerset_C_Polyhedron_iterator_end (ps
, end
);
234 !ppl_Pointset_Powerset_C_Polyhedron_iterator_equal_test (it
, end
);
235 ppl_Pointset_Powerset_C_Polyhedron_iterator_increment (it
))
237 ppl_const_Polyhedron_t ph
;
239 ppl_Pointset_Powerset_C_Polyhedron_iterator_dereference (it
, &ph
);
240 openscop_print_polyhedron_matrix (file
, ph
, output
, input
, locals
,
244 ppl_delete_Pointset_Powerset_C_Polyhedron_iterator (it
);
245 ppl_delete_Pointset_Powerset_C_Polyhedron_iterator (end
);
248 /* Prints to FILE the scattering function of PBB in OpenScop format, at some
252 openscop_print_scattering_function_1 (FILE *file
, poly_bb_p pbb
, int verbosity
)
255 ppl_const_Polyhedron_t ph
;
259 fprintf (file
, "# scattering bb_%d (\n", pbb_index (pbb
));
260 fprintf (file
, "#eq");
262 for (i
= 0; i
< pbb_nb_scattering_transform (pbb
); i
++)
263 fprintf (file
, " s%d", (int) i
);
265 for (i
= 0; i
< pbb_nb_local_vars (pbb
); i
++)
266 fprintf (file
, " lv%d", (int) i
);
268 for (i
= 0; i
< pbb_dim_iter_domain (pbb
); i
++)
269 fprintf (file
, " i%d", (int) i
);
271 for (i
= 0; i
< pbb_nb_params (pbb
); i
++)
272 fprintf (file
, " p%d", (int) i
);
274 fprintf (file
, " cst\n");
277 /* Number of disjunct components. Remove this when
278 PBB_TRANSFORMED_SCATTERING will be a pointset_powerset. */
279 fprintf (file
, "1\n");
281 ph
= PBB_TRANSFORMED_SCATTERING (pbb
)
282 ? PBB_TRANSFORMED_SCATTERING (pbb
)
283 : PBB_ORIGINAL_SCATTERING (pbb
);
285 openscop_print_polyhedron_matrix (file
, ph
,
286 pbb_nb_scattering_transform (pbb
),
287 pbb_dim_iter_domain (pbb
),
288 pbb_nb_local_vars (pbb
),
289 pbb_nb_params (pbb
));
292 fprintf (file
, "#)\n");
295 /* Prints to FILE the scattering function of PBB, at some VERBOSITY
299 print_scattering_function_1 (FILE *file
, poly_bb_p pbb
, int verbosity
)
305 fprintf (file
, "# scattering bb_%d (\n", pbb_index (pbb
));
306 fprintf (file
, "#eq");
308 for (i
= 0; i
< pbb_nb_scattering_transform (pbb
); i
++)
309 fprintf (file
, " s%d", (int) i
);
311 for (i
= 0; i
< pbb_nb_local_vars (pbb
); i
++)
312 fprintf (file
, " lv%d", (int) i
);
314 for (i
= 0; i
< pbb_dim_iter_domain (pbb
); i
++)
315 fprintf (file
, " i%d", (int) i
);
317 for (i
= 0; i
< pbb_nb_params (pbb
); i
++)
318 fprintf (file
, " p%d", (int) i
);
320 fprintf (file
, " cst\n");
323 /* Number of disjunct components. Remove this when
324 PBB_TRANSFORMED_SCATTERING will be a pointset_powerset. */
325 fprintf (file
, "1\n");
326 ppl_print_polyhedron_matrix (file
, PBB_TRANSFORMED_SCATTERING (pbb
)
327 ? PBB_TRANSFORMED_SCATTERING (pbb
)
328 : PBB_ORIGINAL_SCATTERING (pbb
));
331 fprintf (file
, "#)\n");
334 /* Prints to FILE the scattering function of PBB, at some VERBOSITY
338 print_scattering_function (FILE *file
, poly_bb_p pbb
, int verbosity
)
340 if (!PBB_TRANSFORMED (pbb
))
343 if (PBB_TRANSFORMED_SCATTERING (pbb
)
344 || PBB_ORIGINAL_SCATTERING (pbb
))
347 fprintf (file
, "# Scattering function is provided\n");
349 fprintf (file
, "1\n");
354 fprintf (file
, "# Scattering function is not provided\n");
356 fprintf (file
, "0\n");
360 openscop_print_scattering_function_1 (file
, pbb
, verbosity
);
363 fprintf (file
, "# Scattering names are not provided\n");
365 fprintf (file
, "0\n");
369 /* Prints to FILE the iteration domain of PBB, at some VERBOSITY
373 print_iteration_domain (FILE *file
, poly_bb_p pbb
, int verbosity
)
375 print_pbb_domain (file
, pbb
, verbosity
);
378 /* Prints to FILE the scattering functions of every PBB of SCOP. */
381 print_scattering_functions (FILE *file
, scop_p scop
, int verbosity
)
386 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
387 print_scattering_function (file
, pbb
, verbosity
);
390 /* Prints to FILE the iteration domains of every PBB of SCOP, at some
394 print_iteration_domains (FILE *file
, scop_p scop
, int verbosity
)
399 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
400 print_iteration_domain (file
, pbb
, verbosity
);
403 /* Prints to STDERR the scattering function of PBB, at some VERBOSITY
407 debug_scattering_function (poly_bb_p pbb
, int verbosity
)
409 print_scattering_function (stderr
, pbb
, verbosity
);
412 /* Prints to STDERR the iteration domain of PBB, at some VERBOSITY
416 debug_iteration_domain (poly_bb_p pbb
, int verbosity
)
418 print_iteration_domain (stderr
, pbb
, verbosity
);
421 /* Prints to STDERR the scattering functions of every PBB of SCOP, at
422 some VERBOSITY level. */
425 debug_scattering_functions (scop_p scop
, int verbosity
)
427 print_scattering_functions (stderr
, scop
, verbosity
);
430 /* Prints to STDERR the iteration domains of every PBB of SCOP, at
431 some VERBOSITY level. */
434 debug_iteration_domains (scop_p scop
, int verbosity
)
436 print_iteration_domains (stderr
, scop
, verbosity
);
439 /* Read N integer from FILE. */
442 openscop_read_N_int (FILE *file
, int N
)
444 char s
[OPENSCOP_MAX_STRING
];
446 int i
, *res
= (int *) xmalloc (OPENSCOP_MAX_STRING
* sizeof (int));
448 /* Skip blank and commented lines. */
449 while (fgets (s
, sizeof s
, file
) == (char *) 0
456 for (i
= 0; i
< N
; i
++)
458 sscanf (str
, "%d", &res
[i
]);
460 /* Jump the integer that was read. */
461 while ((*str
) && !ISSPACE (*str
) && (*str
!= '#'))
465 while ((*str
) && ISSPACE (*str
) && (*str
!= '#'))
472 /* Read one integer from FILE. */
475 openscop_read_one_int (FILE *file
)
477 int *x
= openscop_read_N_int (file
, 1);
484 /* Read N string from FILE. */
487 openscop_read_N_string (FILE *file
, int N
)
490 char str
[OPENSCOP_MAX_STRING
];
491 char *tmp
= (char *) xmalloc (sizeof (char) * OPENSCOP_MAX_STRING
);
494 /* Skip blank and commented lines. */
495 while (fgets (str
, sizeof str
, file
) == (char *) 0
503 for (i
= 0; i
< N
; i
++)
505 /* Read the first word. */
506 for (; (*s
) && (!ISSPACE (*s
)) && (*s
!= '#'); ++count
)
513 while ((*s
) && ISSPACE (*s
) && (*s
!= '#'))
522 /* Read one string from FILE. */
525 openscop_read_one_string (FILE *file
)
527 return openscop_read_N_string (file
, 1);
530 /* Read from FILE the powerset PS in its OpenScop matrix form. OUTPUT is the
531 number of output dimensions, INPUT is the number of input dimensions,
532 LOCALS is the number of existentially quantified variables and PARAMS is
533 the number of parameters. */
536 openscop_read_powerset_matrix (FILE *file
,
537 ppl_Pointset_Powerset_C_Polyhedron_t
*ps
,
538 int *output
, int *input
, int *locals
,
543 nb_disjuncts
= openscop_read_one_int (file
);
545 for (i
= 0; i
< nb_disjuncts
; i
++)
549 openscop_read_polyhedron_matrix (file
, &ph
, output
, input
, locals
,
554 ppl_new_Pointset_Powerset_C_Polyhedron_from_C_Polyhedron (ps
, ph
);
556 ppl_Pointset_Powerset_C_Polyhedron_add_disjunct (*ps
, ph
);
560 /* Read a scattering function from FILE and save it to PBB. Return whether
561 the scattering function was provided or not. */
564 graphite_read_scatt (FILE *file
, poly_bb_p pbb
)
566 bool scattering_provided
= false;
567 int output
, input
, locals
, params
;
568 ppl_Polyhedron_t newp
;
570 if (openscop_read_one_int (file
) > 0)
572 /* Read number of disjunct components. */
573 openscop_read_one_int (file
);
575 /* Read scattering function. */
576 openscop_read_polyhedron_matrix (file
, &newp
, &output
, &input
,
578 store_scattering (PBB_SCOP (pbb
));
579 PBB_TRANSFORMED (pbb
) = poly_scattering_new ();
580 PBB_TRANSFORMED_SCATTERING (pbb
) = newp
;
581 PBB_NB_LOCAL_VARIABLES (pbb
) = locals
;
583 /* New scattering dimension. */
584 PBB_NB_SCATTERING_TRANSFORM (pbb
) = output
;
586 scattering_provided
= true;
589 return scattering_provided
;
592 /* Read a scop file. Return true if the scop is transformed. */
595 graphite_read_scop_file (FILE *file
, scop_p scop
)
597 char *tmp
, *language
;
598 size_t i
, j
, nb_statements
, nbr
, nbw
;
599 int input
, output
, locals
, params
;
600 ppl_Pointset_Powerset_C_Polyhedron_t ps
;
604 /* Ensure that the file is in OpenScop format. */
605 tmp
= openscop_read_N_string (file
, 2);
607 if (strcmp (tmp
, "SCoP 1"))
609 error ("The file is not in OpenScop format.\n");
615 /* Read the language. */
616 language
= openscop_read_one_string (file
);
618 if (strcmp (language
, "Gimple"))
620 error ("The language is not recognized\n");
626 /* Read the context but do not use it. */
627 openscop_read_powerset_matrix (file
, &ps
, &input
, &output
, &locals
, ¶ms
);
629 if ((size_t) params
!= scop
->nb_params
)
631 error ("Parameters number in the scop file is different from the"
632 " internal scop parameter number.");
636 /* Read parameter names if provided. */
637 if (openscop_read_one_int (file
))
638 openscop_read_N_string (file
, scop
->nb_params
);
640 nb_statements
= openscop_read_one_int (file
);
642 if (nb_statements
!= VEC_length (poly_bb_p
, SCOP_BBS (scop
)))
644 error ("Number of statements in the OpenScop file does not match"
645 " the graphite internal statements number.");
649 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
651 /* Read iteration domain. */
652 openscop_read_powerset_matrix (file
, &ps
, &input
, &output
, &locals
,
655 /* Read scattering. */
656 transform_done
= graphite_read_scatt (file
, pbb
);
658 /* Scattering names. */
659 openscop_read_one_int (file
);
661 /* Read access functions. */
662 if (openscop_read_one_int (file
) > 0)
664 nbr
= openscop_read_one_int (file
);
666 /* Read access functions. */
667 for (j
= 0; j
< nbr
; j
++)
668 openscop_read_powerset_matrix (file
, &ps
, &input
, &output
, &locals
,
671 nbw
= openscop_read_one_int (file
);
673 /* Write access functions. */
674 for (j
= 0; j
< nbw
; j
++)
675 openscop_read_powerset_matrix (file
, &ps
, &input
, &output
, &locals
,
679 /* Statement body. */
680 openscop_read_one_int (file
);
683 return transform_done
;
686 /* Initialize and return a file that will be used to write a scop. SCOP_NUMBER
687 is a sequential number (identifier) used to differentiate scop files.
688 Examples of the generated file names: dump_base_name.0.graphite,
689 dump_base_name.1.graphite, dump_base_name.2.graphite, etc. */
692 init_graphite_out_file (int scop_number
)
694 FILE *graphite_out_file
;
695 int len
= strlen (dump_base_name
);
696 char *dumpname
= XNEWVEC (char, len
+ 25);
697 char *s_scop_number
= XNEWVEC (char, 15);
699 memcpy (dumpname
, dump_base_name
, len
+ 1);
700 strip_off_ending (dumpname
, len
);
701 sprintf (s_scop_number
, ".%d", scop_number
);
702 strcat (dumpname
, s_scop_number
);
703 strcat (dumpname
, ".graphite");
704 graphite_out_file
= fopen (dumpname
, "w+b");
706 if (graphite_out_file
== 0)
707 fatal_error ("can%'t open %s for writing: %m", dumpname
);
711 return graphite_out_file
;
714 /* Open and return a file used for scop reading. SCOP_NUMBER is a sequential
715 number (identifier) used to differentiate scop files. Examples of the
716 generated file names: dump_base_name.0.graphite, dump_base_name.1.graphite,
717 dump_base_name.2.graphite, etc. */
720 init_graphite_in_file (int scop_number
)
722 FILE *graphite_in_file
;
723 int len
= strlen (dump_base_name
);
724 char *dumpname
= XNEWVEC (char, len
+ 25);
725 char *s_scop_number
= XNEWVEC (char, 15);
727 memcpy (dumpname
, dump_base_name
, len
+ 1);
728 strip_off_ending (dumpname
, len
);
729 sprintf (s_scop_number
, ".%d", scop_number
);
730 strcat (dumpname
, s_scop_number
);
731 strcat (dumpname
, ".graphite");
732 graphite_in_file
= fopen (dumpname
, "r+b");
734 if (graphite_in_file
== 0)
735 fatal_error ("can%'t open %s for reading: %m", dumpname
);
739 return graphite_in_file
;
742 /* Apply graphite transformations to all the basic blocks of SCOP. */
745 apply_poly_transforms (scop_p scop
)
747 bool transform_done
= false;
749 static size_t file_scop_number
= 0;
751 /* This feature is only enabled in the Graphite branch. */
754 graphite_file
= init_graphite_in_file (file_scop_number
);
755 transform_done
|= graphite_read_scop_file (graphite_file
, scop
);
756 gcc_assert (graphite_legal_transform (scop
));
760 /* Generate code even if we did not apply any real transformation.
761 This also allows to check the performance for the identity
762 transformation: GIMPLE -> GRAPHITE -> GIMPLE
763 Keep in mind that CLooG optimizes in control, so the loop structure
764 may change, even if we only use -fgraphite-identity. */
765 if (flag_graphite_identity
)
766 transform_done
= true;
768 if (flag_loop_parallelize_all
)
769 transform_done
= true;
772 transform_done
|= scop_do_block (scop
);
775 if (flag_loop_strip_mine
)
776 transform_done
|= scop_do_strip_mine (scop
);
778 if (flag_loop_interchange
)
779 transform_done
|= scop_do_interchange (scop
);
782 /* This feature is only enabled in the Graphite branch. */
785 graphite_file
= init_graphite_out_file (file_scop_number
);
786 print_scop (graphite_file
, scop
, 1);
790 return transform_done
;
793 /* Returns true when it PDR1 is a duplicate of PDR2: same PBB, and
794 their ACCESSES, TYPE, and NB_SUBSCRIPTS are the same. */
797 can_collapse_pdrs (poly_dr_p pdr1
, poly_dr_p pdr2
)
800 ppl_Pointset_Powerset_C_Polyhedron_t af1
, af2
, diff
;
802 if (PDR_PBB (pdr1
) != PDR_PBB (pdr2
)
803 || PDR_NB_SUBSCRIPTS (pdr1
) != PDR_NB_SUBSCRIPTS (pdr2
)
804 || PDR_TYPE (pdr1
) != PDR_TYPE (pdr2
))
807 af1
= PDR_ACCESSES (pdr1
);
808 af2
= PDR_ACCESSES (pdr2
);
809 ppl_new_Pointset_Powerset_C_Polyhedron_from_Pointset_Powerset_C_Polyhedron
811 ppl_Pointset_Powerset_C_Polyhedron_difference_assign (diff
, af2
);
813 res
= ppl_Pointset_Powerset_C_Polyhedron_is_empty (diff
);
814 ppl_delete_Pointset_Powerset_C_Polyhedron (diff
);
818 /* Removes duplicated data references in PBB. */
821 pbb_remove_duplicate_pdrs (poly_bb_p pbb
)
824 poly_dr_p pdr1
, pdr2
;
825 unsigned n
= VEC_length (poly_dr_p
, PBB_DRS (pbb
));
826 VEC (poly_dr_p
, heap
) *collapsed
= VEC_alloc (poly_dr_p
, heap
, n
);
828 FOR_EACH_VEC_ELT (poly_dr_p
, PBB_DRS (pbb
), i
, pdr1
)
829 FOR_EACH_VEC_ELT (poly_dr_p
, collapsed
, j
, pdr2
)
830 if (!can_collapse_pdrs (pdr1
, pdr2
))
831 VEC_quick_push (poly_dr_p
, collapsed
, pdr1
);
833 VEC_free (poly_dr_p
, heap
, collapsed
);
834 PBB_PDR_DUPLICATES_REMOVED (pbb
) = true;
837 /* Create a new polyhedral data reference and add it to PBB. It is
838 defined by its ACCESSES, its TYPE, and the number of subscripts
842 new_poly_dr (poly_bb_p pbb
, int dr_base_object_set
,
843 ppl_Pointset_Powerset_C_Polyhedron_t accesses
,
844 enum poly_dr_type type
, void *cdr
, graphite_dim_t nb_subscripts
)
847 poly_dr_p pdr
= XNEW (struct poly_dr
);
850 PDR_BASE_OBJECT_SET (pdr
) = dr_base_object_set
;
851 PDR_NB_REFS (pdr
) = 1;
853 PDR_ACCESSES (pdr
) = accesses
;
854 PDR_TYPE (pdr
) = type
;
856 PDR_NB_SUBSCRIPTS (pdr
) = nb_subscripts
;
857 VEC_safe_push (poly_dr_p
, heap
, PBB_DRS (pbb
), pdr
);
860 /* Free polyhedral data reference PDR. */
863 free_poly_dr (poly_dr_p pdr
)
865 ppl_delete_Pointset_Powerset_C_Polyhedron (PDR_ACCESSES (pdr
));
869 /* Create a new polyhedral black box. */
872 new_poly_bb (scop_p scop
, void *black_box
, bool reduction
)
874 poly_bb_p pbb
= XNEW (struct poly_bb
);
876 PBB_DOMAIN (pbb
) = NULL
;
877 PBB_SCOP (pbb
) = scop
;
878 pbb_set_black_box (pbb
, black_box
);
879 PBB_TRANSFORMED (pbb
) = NULL
;
880 PBB_SAVED (pbb
) = NULL
;
881 PBB_ORIGINAL (pbb
) = NULL
;
882 PBB_DRS (pbb
) = VEC_alloc (poly_dr_p
, heap
, 3);
883 PBB_IS_REDUCTION (pbb
) = reduction
;
884 PBB_PDR_DUPLICATES_REMOVED (pbb
) = false;
885 VEC_safe_push (poly_bb_p
, heap
, SCOP_BBS (scop
), pbb
);
888 /* Free polyhedral black box. */
891 free_poly_bb (poly_bb_p pbb
)
896 ppl_delete_Pointset_Powerset_C_Polyhedron (PBB_DOMAIN (pbb
));
898 if (PBB_TRANSFORMED (pbb
))
899 poly_scattering_free (PBB_TRANSFORMED (pbb
));
902 poly_scattering_free (PBB_SAVED (pbb
));
904 if (PBB_ORIGINAL (pbb
))
905 poly_scattering_free (PBB_ORIGINAL (pbb
));
908 FOR_EACH_VEC_ELT (poly_dr_p
, PBB_DRS (pbb
), i
, pdr
)
911 VEC_free (poly_dr_p
, heap
, PBB_DRS (pbb
));
916 print_pdr_access_layout (FILE *file
, poly_bb_p pbb
, poly_dr_p pdr
)
920 fprintf (file
, "# eq");
922 fprintf (file
, " alias");
924 for (i
= 0; i
< PDR_NB_SUBSCRIPTS (pdr
); i
++)
925 fprintf (file
, " sub%d", (int) i
);
927 for (i
= 0; i
< pbb_dim_iter_domain (pbb
); i
++)
928 fprintf (file
, " i%d", (int) i
);
930 for (i
= 0; i
< pbb_nb_params (pbb
); i
++)
931 fprintf (file
, " p%d", (int) i
);
933 fprintf (file
, " cst\n");
936 /* Prints to FILE the polyhedral data reference PDR, at some VERBOSITY
940 print_pdr (FILE *file
, poly_dr_p pdr
, int verbosity
)
946 fprintf (file
, "# pdr_%d (", PDR_ID (pdr
));
948 switch (PDR_TYPE (pdr
))
951 fprintf (file
, "read \n");
955 fprintf (file
, "write \n");
959 fprintf (file
, "may_write \n");
966 dump_data_reference (file
, (data_reference_p
) PDR_CDR (pdr
));
971 fprintf (file
, "# data accesses (\n");
972 print_pdr_access_layout (file
, PDR_PBB (pdr
), pdr
);
975 alias_set_dim
= pdr_alias_set_dim (pdr
) + 1;
977 openscop_print_pdr_powerset (file
,
979 PDR_NB_SUBSCRIPTS (pdr
),
981 pbb_nb_params (PDR_PBB (pdr
)));
984 fprintf (file
, "#)\n");
987 fprintf (file
, "#)\n");
990 /* Prints to STDERR the polyhedral data reference PDR, at some
994 debug_pdr (poly_dr_p pdr
, int verbosity
)
996 print_pdr (stderr
, pdr
, verbosity
);
999 /* Creates a new SCOP containing REGION. */
1002 new_scop (void *region
)
1004 scop_p scop
= XNEW (struct scop
);
1006 SCOP_CONTEXT (scop
) = NULL
;
1007 scop_set_region (scop
, region
);
1008 SCOP_BBS (scop
) = VEC_alloc (poly_bb_p
, heap
, 3);
1009 SCOP_ORIGINAL_PDDRS (scop
) = htab_create (10, hash_poly_ddr_p
,
1010 eq_poly_ddr_p
, free_poly_ddr
);
1011 SCOP_ORIGINAL_SCHEDULE (scop
) = NULL
;
1012 SCOP_TRANSFORMED_SCHEDULE (scop
) = NULL
;
1013 SCOP_SAVED_SCHEDULE (scop
) = NULL
;
1014 POLY_SCOP_P (scop
) = false;
1022 free_scop (scop_p scop
)
1027 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
1030 VEC_free (poly_bb_p
, heap
, SCOP_BBS (scop
));
1032 if (SCOP_CONTEXT (scop
))
1033 ppl_delete_Pointset_Powerset_C_Polyhedron (SCOP_CONTEXT (scop
));
1035 htab_delete (SCOP_ORIGINAL_PDDRS (scop
));
1036 free_lst (SCOP_ORIGINAL_SCHEDULE (scop
));
1037 free_lst (SCOP_TRANSFORMED_SCHEDULE (scop
));
1038 free_lst (SCOP_SAVED_SCHEDULE (scop
));
1042 /* Print to FILE the domain of PBB in OpenScop format, at some VERBOSITY
1046 openscop_print_pbb_domain (FILE *file
, poly_bb_p pbb
, int verbosity
)
1049 gimple_bb_p gbb
= PBB_BLACK_BOX (pbb
);
1051 if (!PBB_DOMAIN (pbb
))
1056 fprintf (file
, "\n# Iteration domain of bb_%d (\n", GBB_BB (gbb
)->index
);
1057 fprintf (file
, "#eq");
1059 for (i
= 0; i
< pbb_dim_iter_domain (pbb
); i
++)
1060 fprintf (file
, " i%d", (int) i
);
1062 for (i
= 0; i
< pbb_nb_params (pbb
); i
++)
1063 fprintf (file
, " p%d", (int) i
);
1065 fprintf (file
, " cst\n");
1068 if (PBB_DOMAIN (pbb
))
1069 openscop_print_powerset_matrix (file
, PBB_DOMAIN (pbb
),
1070 pbb_dim_iter_domain (pbb
),
1073 pbb_nb_params (pbb
));
1075 fprintf (file
, "0\n");
1078 fprintf (file
, "#)\n");
1081 /* Print to FILE the domain of PBB, at some VERBOSITY level. */
1084 print_pbb_domain (FILE *file
, poly_bb_p pbb
, int verbosity
)
1087 gimple_bb_p gbb
= PBB_BLACK_BOX (pbb
);
1089 if (!PBB_DOMAIN (pbb
))
1094 fprintf (file
, "# Iteration domain of bb_%d (\n", GBB_BB (gbb
)->index
);
1095 fprintf (file
, "# eq");
1097 for (i
= 0; i
< pbb_dim_iter_domain (pbb
); i
++)
1098 fprintf (file
, " i%d", (int) i
);
1100 for (i
= 0; i
< pbb_nb_params (pbb
); i
++)
1101 fprintf (file
, " p%d", (int) i
);
1103 fprintf (file
, " cst\n");
1106 if (PBB_DOMAIN (pbb
))
1107 ppl_print_powerset_matrix (file
, PBB_DOMAIN (pbb
));
1109 fprintf (file
, "0\n");
1112 fprintf (file
, "#)\n");
1115 /* Dump the cases of a graphite basic block GBB on FILE. */
1118 dump_gbb_cases (FILE *file
, gimple_bb_p gbb
)
1122 VEC (gimple
, heap
) *cases
;
1127 cases
= GBB_CONDITION_CASES (gbb
);
1128 if (VEC_empty (gimple
, cases
))
1131 fprintf (file
, "# cases bb_%d (\n", GBB_BB (gbb
)->index
);
1133 FOR_EACH_VEC_ELT (gimple
, cases
, i
, stmt
)
1135 fprintf (file
, "# ");
1136 print_gimple_stmt (file
, stmt
, 0, 0);
1139 fprintf (file
, "#)\n");
1142 /* Dump conditions of a graphite basic block GBB on FILE. */
1145 dump_gbb_conditions (FILE *file
, gimple_bb_p gbb
)
1149 VEC (gimple
, heap
) *conditions
;
1154 conditions
= GBB_CONDITIONS (gbb
);
1155 if (VEC_empty (gimple
, conditions
))
1158 fprintf (file
, "# conditions bb_%d (\n", GBB_BB (gbb
)->index
);
1160 FOR_EACH_VEC_ELT (gimple
, conditions
, i
, stmt
)
1162 fprintf (file
, "# ");
1163 print_gimple_stmt (file
, stmt
, 0, 0);
1166 fprintf (file
, "#)\n");
1169 /* Print to FILE all the data references of PBB, at some VERBOSITY
1173 print_pdrs (FILE *file
, poly_bb_p pbb
, int verbosity
)
1180 if (VEC_length (poly_dr_p
, PBB_DRS (pbb
)) == 0)
1183 fprintf (file
, "# Access informations are not provided\n");\
1184 fprintf (file
, "0\n");
1189 fprintf (file
, "# Data references (\n");
1192 fprintf (file
, "# Access informations are provided\n");
1193 fprintf (file
, "1\n");
1195 FOR_EACH_VEC_ELT (poly_dr_p
, PBB_DRS (pbb
), i
, pdr
)
1196 if (PDR_TYPE (pdr
) == PDR_READ
)
1202 fprintf (file
, "# Read data references (\n");
1205 fprintf (file
, "# Read access informations\n");
1206 fprintf (file
, "%d\n", nb_reads
);
1208 FOR_EACH_VEC_ELT (poly_dr_p
, PBB_DRS (pbb
), i
, pdr
)
1209 if (PDR_TYPE (pdr
) == PDR_READ
)
1210 print_pdr (file
, pdr
, verbosity
);
1213 fprintf (file
, "#)\n");
1216 fprintf (file
, "# Write data references (\n");
1219 fprintf (file
, "# Write access informations\n");
1220 fprintf (file
, "%d\n", nb_writes
);
1222 FOR_EACH_VEC_ELT (poly_dr_p
, PBB_DRS (pbb
), i
, pdr
)
1223 if (PDR_TYPE (pdr
) != PDR_READ
)
1224 print_pdr (file
, pdr
, verbosity
);
1227 fprintf (file
, "#)\n");
1230 fprintf (file
, "#)\n");
1233 /* Print to STDERR all the data references of PBB. */
1236 debug_pdrs (poly_bb_p pbb
, int verbosity
)
1238 print_pdrs (stderr
, pbb
, verbosity
);
1241 /* Print to FILE the body of PBB, at some VERBOSITY level.
1242 If statement_body_provided is false statement body is not printed. */
1245 print_pbb_body (FILE *file
, poly_bb_p pbb
, int verbosity
,
1246 bool statement_body_provided
)
1249 fprintf (file
, "# Body (\n");
1251 if (!statement_body_provided
)
1254 fprintf (file
, "# Statement body is not provided\n");
1256 fprintf (file
, "0\n");
1261 fprintf (file
, "# Statement body is provided\n");
1262 fprintf (file
, "1\n");
1265 fprintf (file
, "# Original iterator names\n# Iterator names are not provided yet.\n");
1268 fprintf (file
, "# Statement body\n");
1270 fprintf (file
, "{\n");
1271 dump_bb (pbb_bb (pbb
), file
, 0);
1272 fprintf (file
, "}\n");
1275 fprintf (file
, "#)\n");
1278 /* Print to FILE the domain and scattering function of PBB, at some
1282 print_pbb (FILE *file
, poly_bb_p pbb
, int verbosity
)
1286 fprintf (file
, "# pbb_%d (\n", pbb_index (pbb
));
1287 dump_gbb_conditions (file
, PBB_BLACK_BOX (pbb
));
1288 dump_gbb_cases (file
, PBB_BLACK_BOX (pbb
));
1291 openscop_print_pbb_domain (file
, pbb
, verbosity
);
1292 print_scattering_function (file
, pbb
, verbosity
);
1293 print_pdrs (file
, pbb
, verbosity
);
1294 print_pbb_body (file
, pbb
, verbosity
, false);
1297 fprintf (file
, "#)\n");
1300 /* Print to FILE the parameters of SCOP, at some VERBOSITY level. */
1303 print_scop_params (FILE *file
, scop_p scop
, int verbosity
)
1309 fprintf (file
, "# parameters (\n");
1311 if (VEC_length (tree
, SESE_PARAMS (SCOP_REGION (scop
))))
1314 fprintf (file
, "# Parameter names are provided\n");
1316 fprintf (file
, "1\n");
1319 fprintf (file
, "# Parameter names\n");
1324 fprintf (file
, "# Parameter names are not provided\n");
1325 fprintf (file
, "0\n");
1328 FOR_EACH_VEC_ELT (tree
, SESE_PARAMS (SCOP_REGION (scop
)), i
, t
)
1330 print_generic_expr (file
, t
, 0);
1331 fprintf (file
, " ");
1334 fprintf (file
, "\n");
1337 fprintf (file
, "#)\n");
1340 /* Print to FILE the context of SCoP in OpenScop format, at some VERBOSITY
1344 openscop_print_scop_context (FILE *file
, scop_p scop
, int verbosity
)
1350 fprintf (file
, "# Context (\n");
1351 fprintf (file
, "#eq");
1353 for (i
= 0; i
< scop_nb_params (scop
); i
++)
1354 fprintf (file
, " p%d", (int) i
);
1356 fprintf (file
, " cst\n");
1359 if (SCOP_CONTEXT (scop
))
1360 openscop_print_powerset_matrix (file
, SCOP_CONTEXT (scop
), 0, 0, 0,
1361 scop_nb_params (scop
));
1363 fprintf (file
, "0 %d 0 0 0 %d\n", (int) scop_nb_params (scop
) + 2,
1364 (int) scop_nb_params (scop
));
1367 fprintf (file
, "# )\n");
1370 /* Print to FILE the context of SCoP, at some VERBOSITY level. */
1373 print_scop_context (FILE *file
, scop_p scop
, int verbosity
)
1379 fprintf (file
, "# Context (\n");
1380 fprintf (file
, "#eq");
1382 for (i
= 0; i
< scop_nb_params (scop
); i
++)
1383 fprintf (file
, " p%d", (int) i
);
1385 fprintf (file
, " cst\n");
1388 if (SCOP_CONTEXT (scop
))
1389 ppl_print_powerset_matrix (file
, SCOP_CONTEXT (scop
));
1391 fprintf (file
, "0 %d\n", (int) scop_nb_params (scop
) + 2);
1394 fprintf (file
, "# )\n");
1397 /* Print to FILE the SCOP header: context, parameters, and statements
1401 print_scop_header (FILE *file
, scop_p scop
, int verbosity
)
1403 fprintf (file
, "SCoP 1\n#(\n");
1404 fprintf (file
, "# Language\nGimple\n");
1405 openscop_print_scop_context (file
, scop
, verbosity
);
1406 print_scop_params (file
, scop
, verbosity
);
1409 fprintf (file
, "# Number of statements\n");
1411 fprintf (file
, "%d\n",VEC_length (poly_bb_p
, SCOP_BBS (scop
)));
1414 /* Print to FILE the SCOP, at some VERBOSITY level. */
1417 print_scop (FILE *file
, scop_p scop
, int verbosity
)
1422 print_scop_header (file
, scop
, verbosity
);
1424 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
1425 print_pbb (file
, pbb
, verbosity
);
1429 fprintf (file
, "# original_lst (\n");
1430 print_lst (file
, SCOP_ORIGINAL_SCHEDULE (scop
), 0);
1431 fprintf (file
, "\n#)\n");
1433 fprintf (file
, "# transformed_lst (\n");
1434 print_lst (file
, SCOP_TRANSFORMED_SCHEDULE (scop
), 0);
1435 fprintf (file
, "\n#)\n");
1438 fprintf (file
, "#)\n");
1441 /* Print to FILE the input file that CLooG would expect as input, at
1442 some VERBOSITY level. */
1445 print_cloog (FILE *file
, scop_p scop
, int verbosity
)
1450 fprintf (file
, "# SCoP (generated by GCC/Graphite\n");
1452 fprintf (file
, "# CLooG output language\n");
1453 fprintf (file
, "c\n");
1455 print_scop_context (file
, scop
, verbosity
);
1456 print_scop_params (file
, scop
, verbosity
);
1459 fprintf (file
, "# Number of statements\n");
1461 fprintf (file
, "%d\n", VEC_length (poly_bb_p
, SCOP_BBS (scop
)));
1463 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
1466 fprintf (file
, "# pbb_%d (\n", pbb_index (pbb
));
1468 print_pbb_domain (file
, pbb
, verbosity
);
1469 fprintf (file
, "0 0 0");
1472 fprintf (file
, "# For future CLooG options.\n");
1474 fprintf (file
, "\n");
1477 fprintf (file
, "#)\n");
1480 fprintf (file
, "0");
1482 fprintf (file
, "# Don't set the iterator names.\n");
1484 fprintf (file
, "\n");
1487 fprintf (file
, "# Number of scattering functions\n");
1489 fprintf (file
, "%d\n", VEC_length (poly_bb_p
, SCOP_BBS (scop
)));
1490 unify_scattering_dimensions (scop
);
1492 FOR_EACH_VEC_ELT (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
)
1494 if (!PBB_TRANSFORMED (pbb
)
1495 || !(PBB_TRANSFORMED_SCATTERING (pbb
)
1496 || PBB_ORIGINAL_SCATTERING (pbb
)))
1500 fprintf (file
, "# pbb_%d (\n", pbb_index (pbb
));
1502 print_scattering_function_1 (file
, pbb
, verbosity
);
1505 fprintf (file
, "#)\n");
1508 fprintf (file
, "0");
1510 fprintf (file
, "# Don't set the scattering dimension names.\n");
1512 fprintf (file
, "\n");
1514 fprintf (file
, "#)\n");
1517 /* Print to STDERR the domain of PBB, at some VERBOSITY level. */
1520 debug_pbb_domain (poly_bb_p pbb
, int verbosity
)
1522 print_pbb_domain (stderr
, pbb
, verbosity
);
1525 /* Print to FILE the domain and scattering function of PBB, at some
1529 debug_pbb (poly_bb_p pbb
, int verbosity
)
1531 print_pbb (stderr
, pbb
, verbosity
);
1534 /* Print to STDERR the context of SCOP, at some VERBOSITY level. */
1537 debug_scop_context (scop_p scop
, int verbosity
)
1539 print_scop_context (stderr
, scop
, verbosity
);
1542 /* Print to STDERR the SCOP, at some VERBOSITY level. */
1545 debug_scop (scop_p scop
, int verbosity
)
1547 print_scop (stderr
, scop
, verbosity
);
1550 /* Print to STDERR the SCOP under CLooG format, at some VERBOSITY
1554 debug_cloog (scop_p scop
, int verbosity
)
1556 print_cloog (stderr
, scop
, verbosity
);
1559 /* Print to STDERR the parameters of SCOP, at some VERBOSITY
1563 debug_scop_params (scop_p scop
, int verbosity
)
1565 print_scop_params (stderr
, scop
, verbosity
);
1569 /* The dimension in the transformed scattering polyhedron of PBB
1570 containing the scattering iterator for the loop at depth LOOP_DEPTH. */
1573 psct_scattering_dim_for_loop_depth (poly_bb_p pbb
, graphite_dim_t loop_depth
)
1575 ppl_const_Constraint_System_t pcs
;
1576 ppl_Constraint_System_const_iterator_t cit
, cend
;
1577 ppl_const_Constraint_t cstr
;
1578 ppl_Polyhedron_t ph
= PBB_TRANSFORMED_SCATTERING (pbb
);
1579 ppl_dimension_type iter
= psct_iterator_dim (pbb
, loop_depth
);
1580 ppl_Linear_Expression_t expr
;
1581 ppl_Coefficient_t coef
;
1586 ppl_new_Coefficient (&coef
);
1587 ppl_Polyhedron_get_constraints (ph
, &pcs
);
1588 ppl_new_Constraint_System_const_iterator (&cit
);
1589 ppl_new_Constraint_System_const_iterator (&cend
);
1591 for (ppl_Constraint_System_begin (pcs
, cit
),
1592 ppl_Constraint_System_end (pcs
, cend
);
1593 !ppl_Constraint_System_const_iterator_equal_test (cit
, cend
);
1594 ppl_Constraint_System_const_iterator_increment (cit
))
1596 ppl_Constraint_System_const_iterator_dereference (cit
, &cstr
);
1597 ppl_new_Linear_Expression_from_Constraint (&expr
, cstr
);
1598 ppl_Linear_Expression_coefficient (expr
, iter
, coef
);
1599 ppl_Coefficient_to_mpz_t (coef
, val
);
1601 if (mpz_sgn (val
) == 0)
1603 ppl_delete_Linear_Expression (expr
);
1607 for (i
= 0; i
< pbb_nb_scattering_transform (pbb
); i
++)
1609 ppl_dimension_type scatter
= psct_scattering_dim (pbb
, i
);
1611 ppl_Linear_Expression_coefficient (expr
, scatter
, coef
);
1612 ppl_Coefficient_to_mpz_t (coef
, val
);
1614 if (mpz_sgn (val
) != 0)
1617 ppl_delete_Linear_Expression (expr
);
1618 ppl_delete_Coefficient (coef
);
1619 ppl_delete_Constraint_System_const_iterator (cit
);
1620 ppl_delete_Constraint_System_const_iterator (cend
);
1630 /* Returns the number of iterations NITER of the loop around PBB at
1631 depth LOOP_DEPTH. */
1634 pbb_number_of_iterations (poly_bb_p pbb
,
1635 graphite_dim_t loop_depth
,
1638 ppl_Linear_Expression_t le
;
1639 ppl_dimension_type dim
;
1641 ppl_Pointset_Powerset_C_Polyhedron_space_dimension (PBB_DOMAIN (pbb
), &dim
);
1642 ppl_new_Linear_Expression_with_dimension (&le
, dim
);
1643 ppl_set_coef (le
, pbb_iterator_dim (pbb
, loop_depth
), 1);
1644 mpz_set_si (niter
, -1);
1645 ppl_max_for_le_pointset (PBB_DOMAIN (pbb
), le
, niter
);
1646 ppl_delete_Linear_Expression (le
);
1649 /* Returns the number of iterations NITER of the loop around PBB at
1650 time(scattering) dimension TIME_DEPTH. */
1653 pbb_number_of_iterations_at_time (poly_bb_p pbb
,
1654 graphite_dim_t time_depth
,
1657 ppl_Pointset_Powerset_C_Polyhedron_t ext_domain
, sctr
;
1658 ppl_Linear_Expression_t le
;
1659 ppl_dimension_type dim
;
1661 /* Takes together domain and scattering polyhedrons, and composes
1662 them into the bigger polyhedron that has the following format:
1664 t0..t_{n-1} | l0..l_{nlcl-1} | i0..i_{niter-1} | g0..g_{nparm-1}
1667 | t0..t_{n-1} are time dimensions (scattering dimensions)
1668 | l0..l_{nclc-1} are local variables in scattering function
1669 | i0..i_{niter-1} are original iteration variables
1670 | g0..g_{nparam-1} are global parameters. */
1672 ppl_new_Pointset_Powerset_C_Polyhedron_from_C_Polyhedron (&sctr
,
1673 PBB_TRANSFORMED_SCATTERING (pbb
));
1675 /* Extend the iteration domain with the scattering dimensions:
1676 0..0 | 0..0 | i0..i_{niter-1} | g0..g_{nparm-1}. */
1677 ppl_new_Pointset_Powerset_C_Polyhedron_from_Pointset_Powerset_C_Polyhedron
1678 (&ext_domain
, PBB_DOMAIN (pbb
));
1679 ppl_insert_dimensions_pointset (ext_domain
, 0,
1680 pbb_nb_scattering_transform (pbb
)
1681 + pbb_nb_local_vars (pbb
));
1683 /* Add to sctr the extended domain. */
1684 ppl_Pointset_Powerset_C_Polyhedron_intersection_assign (sctr
, ext_domain
);
1686 /* Extract the number of iterations. */
1687 ppl_Pointset_Powerset_C_Polyhedron_space_dimension (sctr
, &dim
);
1688 ppl_new_Linear_Expression_with_dimension (&le
, dim
);
1689 ppl_set_coef (le
, time_depth
, 1);
1690 mpz_set_si (niter
, -1);
1691 ppl_max_for_le_pointset (sctr
, le
, niter
);
1693 ppl_delete_Linear_Expression (le
);
1694 ppl_delete_Pointset_Powerset_C_Polyhedron (sctr
);
1695 ppl_delete_Pointset_Powerset_C_Polyhedron (ext_domain
);
1698 /* Translates LOOP to LST. */
1701 loop_to_lst (loop_p loop
, VEC (poly_bb_p
, heap
) *bbs
, int *i
)
1704 VEC (lst_p
, heap
) *seq
= VEC_alloc (lst_p
, heap
, 5);
1706 for (; VEC_iterate (poly_bb_p
, bbs
, *i
, pbb
); (*i
)++)
1709 basic_block bb
= GBB_BB (PBB_BLACK_BOX (pbb
));
1711 if (bb
->loop_father
== loop
)
1712 stmt
= new_lst_stmt (pbb
);
1713 else if (flow_bb_inside_loop_p (loop
, bb
))
1715 loop_p next
= loop
->inner
;
1717 while (next
&& !flow_bb_inside_loop_p (next
, bb
))
1720 stmt
= loop_to_lst (next
, bbs
, i
);
1725 return new_lst_loop (seq
);
1728 VEC_safe_push (lst_p
, heap
, seq
, stmt
);
1731 return new_lst_loop (seq
);
1734 /* Reads the original scattering of the SCOP and returns an LST
1738 scop_to_lst (scop_p scop
)
1741 int i
, n
= VEC_length (poly_bb_p
, SCOP_BBS (scop
));
1742 VEC (lst_p
, heap
) *seq
= VEC_alloc (lst_p
, heap
, 5);
1743 sese region
= SCOP_REGION (scop
);
1745 for (i
= 0; i
< n
; i
++)
1747 poly_bb_p pbb
= VEC_index (poly_bb_p
, SCOP_BBS (scop
), i
);
1748 loop_p loop
= outermost_loop_in_sese (region
, GBB_BB (PBB_BLACK_BOX (pbb
)));
1750 if (loop_in_sese_p (loop
, region
))
1751 res
= loop_to_lst (loop
, SCOP_BBS (scop
), &i
);
1753 res
= new_lst_stmt (pbb
);
1755 VEC_safe_push (lst_p
, heap
, seq
, res
);
1758 res
= new_lst_loop (seq
);
1759 SCOP_ORIGINAL_SCHEDULE (scop
) = res
;
1760 SCOP_TRANSFORMED_SCHEDULE (scop
) = copy_lst (res
);
1763 /* Print to FILE on a new line COLUMN white spaces. */
1766 lst_indent_to (FILE *file
, int column
)
1771 fprintf (file
, "\n#");
1773 for (i
= 0; i
< column
; i
++)
1774 fprintf (file
, " ");
1777 /* Print LST to FILE with INDENT spaces of indentation. */
1780 print_lst (FILE *file
, lst_p lst
, int indent
)
1785 lst_indent_to (file
, indent
);
1787 if (LST_LOOP_P (lst
))
1792 if (LST_LOOP_FATHER (lst
))
1793 fprintf (file
, "%d (loop", lst_dewey_number (lst
));
1795 fprintf (file
, "#(root");
1797 FOR_EACH_VEC_ELT (lst_p
, LST_SEQ (lst
), i
, l
)
1798 print_lst (file
, l
, indent
+ 2);
1800 fprintf (file
, ")");
1803 fprintf (file
, "%d stmt_%d", lst_dewey_number (lst
), pbb_index (LST_PBB (lst
)));
1806 /* Print LST to STDERR. */
1809 debug_lst (lst_p lst
)
1811 print_lst (stderr
, lst
, 0);
1814 /* Pretty print to FILE the loop statement tree LST in DOT format. */
1817 dot_lst_1 (FILE *file
, lst_p lst
)
1822 if (LST_LOOP_P (lst
))
1827 if (!LST_LOOP_FATHER (lst
))
1828 fprintf (file
, "L -> L_%d_%d\n",
1830 lst_dewey_number (lst
));
1832 fprintf (file
, "L_%d_%d -> L_%d_%d\n",
1833 lst_depth (LST_LOOP_FATHER (lst
)),
1834 lst_dewey_number (LST_LOOP_FATHER (lst
)),
1836 lst_dewey_number (lst
));
1838 FOR_EACH_VEC_ELT (lst_p
, LST_SEQ (lst
), i
, l
)
1839 dot_lst_1 (file
, l
);
1843 fprintf (file
, "L_%d_%d -> S_%d\n",
1844 lst_depth (LST_LOOP_FATHER (lst
)),
1845 lst_dewey_number (LST_LOOP_FATHER (lst
)),
1846 pbb_index (LST_PBB (lst
)));
1850 /* Display the LST using dotty. */
1855 /* When debugging, enable the following code. This cannot be used
1856 in production compilers because it calls "system". */
1859 FILE *stream
= fopen ("/tmp/lst.dot", "w");
1860 gcc_assert (stream
);
1862 fputs ("digraph all {\n", stream
);
1863 dot_lst_1 (stream
, lst
);
1864 fputs ("}\n\n", stream
);
1867 x
= system ("dotty /tmp/lst.dot &");
1869 fputs ("digraph all {\n", stderr
);
1870 dot_lst_1 (stderr
, lst
);
1871 fputs ("}\n\n", stderr
);