1 /* Graphite polyhedral representation.
2 Copyright (C) 2009 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-flow.h"
33 #include "tree-dump.h"
36 #include "tree-chrec.h"
37 #include "tree-data-ref.h"
38 #include "tree-scalar-evolution.h"
39 #include "tree-pass.h"
41 #include "value-prof.h"
42 #include "pointer-set.h"
47 #include "cloog/cloog.h"
50 #include "graphite-ppl.h"
52 #include "graphite-poly.h"
53 #include "graphite-dependences.h"
55 /* Return the maximal loop depth in SCOP. */
58 scop_max_loop_depth (scop_p scop
)
64 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
66 int nb_loops
= pbb_dim_iter_domain (pbb
);
67 if (max_nb_loops
< nb_loops
)
68 max_nb_loops
= nb_loops
;
74 /* Extend the scattering matrix of PBB to MAX_SCATTERING scattering
78 extend_scattering (poly_bb_p pbb
, int max_scattering
)
80 ppl_dimension_type nb_old_dims
, nb_new_dims
;
82 ppl_Coefficient_t coef
;
85 nb_added_dims
= max_scattering
- pbb_nb_scattering_transform (pbb
);
87 value_set_si (one
, 1);
88 ppl_new_Coefficient (&coef
);
89 ppl_assign_Coefficient_from_mpz_t (coef
, one
);
91 gcc_assert (nb_added_dims
>= 0);
93 nb_old_dims
= pbb_nb_scattering_transform (pbb
) + pbb_dim_iter_domain (pbb
)
94 + scop_nb_params (PBB_SCOP (pbb
));
95 nb_new_dims
= nb_old_dims
+ nb_added_dims
;
97 ppl_insert_dimensions (PBB_TRANSFORMED_SCATTERING (pbb
),
98 pbb_nb_scattering_transform (pbb
), nb_added_dims
);
99 PBB_NB_SCATTERING_TRANSFORM (pbb
) += nb_added_dims
;
101 /* Add identity matrix for the added dimensions. */
102 for (i
= max_scattering
- nb_added_dims
; i
< max_scattering
; i
++)
104 ppl_Constraint_t cstr
;
105 ppl_Linear_Expression_t expr
;
107 ppl_new_Linear_Expression_with_dimension (&expr
, nb_new_dims
);
108 ppl_Linear_Expression_add_to_coefficient (expr
, i
, coef
);
109 ppl_new_Constraint (&cstr
, expr
, PPL_CONSTRAINT_TYPE_EQUAL
);
110 ppl_Polyhedron_add_constraint (PBB_TRANSFORMED_SCATTERING (pbb
), cstr
);
111 ppl_delete_Constraint (cstr
);
112 ppl_delete_Linear_Expression (expr
);
115 ppl_delete_Coefficient (coef
);
119 /* All scattering matrices in SCOP will have the same number of scattering
123 unify_scattering_dimensions (scop_p scop
)
127 graphite_dim_t max_scattering
= 0;
129 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
130 max_scattering
= MAX (pbb_nb_scattering_transform (pbb
), max_scattering
);
132 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
133 extend_scattering (pbb
, max_scattering
);
135 return max_scattering
;
138 /* Prints to FILE the scattering function of PBB. */
141 print_scattering_function (FILE *file
, poly_bb_p pbb
)
145 if (!PBB_TRANSFORMED_SCATTERING (pbb
))
148 fprintf (file
, "scattering bb_%d (\n", GBB_BB (PBB_BLACK_BOX (pbb
))->index
);
149 fprintf (file
, "# eq");
151 for (i
= 0; i
< pbb_nb_scattering_transform (pbb
); i
++)
152 fprintf (file
, " s%d", (int) i
);
154 for (i
= 0; i
< pbb_nb_local_vars (pbb
); i
++)
155 fprintf (file
, " lv%d", (int) i
);
157 for (i
= 0; i
< pbb_dim_iter_domain (pbb
); i
++)
158 fprintf (file
, " i%d", (int) i
);
160 for (i
= 0; i
< pbb_nb_params (pbb
); i
++)
161 fprintf (file
, " p%d", (int) i
);
163 fprintf (file
, " cst\n");
165 ppl_print_polyhedron_matrix (file
, PBB_TRANSFORMED_SCATTERING (pbb
));
167 fprintf (file
, ")\n");
170 /* Prints to FILE the iteration domain of PBB. */
173 print_iteration_domain (FILE *file
, poly_bb_p pbb
)
175 print_pbb_domain (file
, pbb
);
178 /* Prints to FILE the scattering functions of every PBB of SCOP. */
181 print_scattering_functions (FILE *file
, scop_p scop
)
186 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
187 print_scattering_function (file
, pbb
);
190 /* Prints to FILE the iteration domains of every PBB of SCOP. */
193 print_iteration_domains (FILE *file
, scop_p scop
)
198 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
199 print_iteration_domain (file
, pbb
);
202 /* Prints to STDERR the scattering function of PBB. */
205 debug_scattering_function (poly_bb_p pbb
)
207 print_scattering_function (stderr
, pbb
);
210 /* Prints to STDERR the iteration domain of PBB. */
213 debug_iteration_domain (poly_bb_p pbb
)
215 print_iteration_domain (stderr
, pbb
);
218 /* Prints to STDERR the scattering functions of every PBB of SCOP. */
221 debug_scattering_functions (scop_p scop
)
223 print_scattering_functions (stderr
, scop
);
226 /* Prints to STDERR the iteration domains of every PBB of SCOP. */
229 debug_iteration_domains (scop_p scop
)
231 print_iteration_domains (stderr
, scop
);
234 /* Apply graphite transformations to all the basic blocks of SCOP. */
237 apply_poly_transforms (scop_p scop
)
239 bool transform_done
= false;
241 gcc_assert (graphite_legal_transform (scop
));
243 /* Generate code even if we did not apply any real transformation.
244 This also allows to check the performance for the identity
245 transformation: GIMPLE -> GRAPHITE -> GIMPLE
246 Keep in mind that CLooG optimizes in control, so the loop structure
247 may change, even if we only use -fgraphite-identity. */
248 if (flag_graphite_identity
)
249 transform_done
= true;
251 if (flag_loop_parallelize_all
)
252 transform_done
= true;
255 gcc_unreachable (); /* Not yet supported. */
257 if (flag_loop_strip_mine
)
259 transform_done
|= scop_do_strip_mine (scop
);
260 gcc_assert (graphite_legal_transform (scop
));
263 if (flag_loop_interchange
)
265 transform_done
|= scop_do_interchange (scop
);
266 gcc_assert (graphite_legal_transform (scop
));
269 return transform_done
;
272 /* Create a new polyhedral data reference and add it to PBB. It is defined by
273 its ACCESSES, its TYPE*/
276 new_poly_dr (poly_bb_p pbb
,
277 ppl_Pointset_Powerset_C_Polyhedron_t accesses
,
278 ppl_Pointset_Powerset_C_Polyhedron_t data_container
,
279 enum POLY_DR_TYPE type
, void *cdr
)
281 poly_dr_p pdr
= XNEW (struct poly_dr
);
284 PDR_ACCESSES (pdr
) = accesses
;
285 PDR_DATA_CONTAINER (pdr
) = data_container
;
286 PDR_TYPE (pdr
) = type
;
288 VEC_safe_push (poly_dr_p
, heap
, PBB_DRS (pbb
), pdr
);
291 /* Free polyhedral data reference PDR. */
294 free_poly_dr (poly_dr_p pdr
)
296 ppl_delete_Pointset_Powerset_C_Polyhedron (PDR_ACCESSES (pdr
));
297 ppl_delete_Pointset_Powerset_C_Polyhedron (PDR_DATA_CONTAINER (pdr
));
302 /* Create a new polyhedral black box. */
305 new_poly_bb (scop_p scop
, void *black_box
)
307 poly_bb_p pbb
= XNEW (struct poly_bb
);
309 PBB_DOMAIN (pbb
) = NULL
;
310 PBB_SCOP (pbb
) = scop
;
311 pbb_set_black_box (pbb
, black_box
);
312 PBB_TRANSFORMED_SCATTERING (pbb
) = NULL
;
313 PBB_ORIGINAL_SCATTERING (pbb
) = NULL
;
314 PBB_DRS (pbb
) = VEC_alloc (poly_dr_p
, heap
, 3);
315 PBB_NB_SCATTERING_TRANSFORM (pbb
) = 0;
316 PBB_NB_LOCAL_VARIABLES (pbb
) = 0;
317 VEC_safe_push (poly_bb_p
, heap
, SCOP_BBS (scop
), pbb
);
320 /* Free polyhedral black box. */
323 free_poly_bb (poly_bb_p pbb
)
328 ppl_delete_Pointset_Powerset_C_Polyhedron (PBB_DOMAIN (pbb
));
330 if (PBB_TRANSFORMED_SCATTERING (pbb
))
331 ppl_delete_Polyhedron (PBB_TRANSFORMED_SCATTERING (pbb
));
333 if (PBB_ORIGINAL_SCATTERING (pbb
))
334 ppl_delete_Polyhedron (PBB_ORIGINAL_SCATTERING (pbb
));
337 for (i
= 0; VEC_iterate (poly_dr_p
, PBB_DRS (pbb
), i
, pdr
); i
++)
340 VEC_free (poly_dr_p
, heap
, PBB_DRS (pbb
));
345 print_pdr_access_layout (FILE *file
, poly_dr_p pdr
)
349 fprintf (file
, "# eq");
351 for (i
= 0; i
< pdr_dim_iter_domain (pdr
); i
++)
352 fprintf (file
, " i%d", (int) i
);
354 for (i
= 0; i
< pdr_nb_params (pdr
); i
++)
355 fprintf (file
, " p%d", (int) i
);
357 fprintf (file
, " alias");
359 for (i
= 0; i
< pdr_nb_subscripts (pdr
); i
++)
360 fprintf (file
, " sub%d", (int) i
);
362 fprintf (file
, " cst\n");
365 /* Prints to FILE the polyhedral data reference PDR. */
368 print_pdr (FILE *file
, poly_dr_p pdr
)
370 fprintf (file
, "pdr (");
372 switch (PDR_TYPE (pdr
))
375 fprintf (file
, "read \n");
379 fprintf (file
, "write \n");
383 fprintf (file
, "may_write \n");
390 dump_data_reference (file
, (data_reference_p
) PDR_CDR (pdr
));
392 fprintf (file
, "data accesses (\n");
393 print_pdr_access_layout (file
, pdr
);
394 ppl_print_powerset_matrix (file
, PDR_ACCESSES (pdr
));
395 fprintf (file
, ")\n");
397 fprintf (file
, "data container (\n");
398 print_pdr_access_layout (file
, pdr
);
399 ppl_print_powerset_matrix (file
, PDR_DATA_CONTAINER (pdr
));
400 fprintf (file
, ")\n");
402 fprintf (file
, ")\n");
405 /* Prints to STDERR the polyhedral data reference PDR. */
408 debug_pdr (poly_dr_p pdr
)
410 print_pdr (stderr
, pdr
);
413 /* Creates a new SCOP containing REGION. */
416 new_scop (void *region
)
418 scop_p scop
= XNEW (struct scop
);
420 SCOP_DEP_GRAPH (scop
) = NULL
;
421 SCOP_CONTEXT (scop
) = NULL
;
422 scop_set_region (scop
, region
);
423 SCOP_BBS (scop
) = VEC_alloc (poly_bb_p
, heap
, 3);
424 SCOP_ORIGINAL_PDR_PAIRS (scop
) = htab_create (10, hash_poly_dr_pair_p
,
425 eq_poly_dr_pair_p
, free
);
432 free_scop (scop_p scop
)
437 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
440 VEC_free (poly_bb_p
, heap
, SCOP_BBS (scop
));
442 if (SCOP_CONTEXT (scop
))
443 ppl_delete_Pointset_Powerset_C_Polyhedron (SCOP_CONTEXT (scop
));
445 htab_delete (SCOP_ORIGINAL_PDR_PAIRS (scop
));
449 /* Print to FILE the domain of PBB. */
452 print_pbb_domain (FILE *file
, poly_bb_p pbb
)
455 gimple_bb_p gbb
= PBB_BLACK_BOX (pbb
);
457 if (!PBB_DOMAIN (pbb
))
460 fprintf (file
, "domains bb_%d (\n", GBB_BB (gbb
)->index
);
461 fprintf (file
, "# eq");
463 for (i
= 0; i
< pbb_dim_iter_domain (pbb
); i
++)
464 fprintf (file
, " i%d", (int) i
);
466 for (i
= 0; i
< pbb_nb_params (pbb
); i
++)
467 fprintf (file
, " p%d", (int) i
);
469 fprintf (file
, " cst\n");
471 if (PBB_DOMAIN (pbb
))
472 ppl_print_powerset_matrix (file
, PBB_DOMAIN (pbb
));
474 fprintf (file
, ")\n");
477 /* Dump the cases of a graphite basic block GBB on FILE. */
480 dump_gbb_cases (FILE *file
, gimple_bb_p gbb
)
484 VEC (gimple
, heap
) *cases
;
489 cases
= GBB_CONDITION_CASES (gbb
);
490 if (VEC_empty (gimple
, cases
))
493 fprintf (file
, "cases bb_%d (", GBB_BB (gbb
)->index
);
495 for (i
= 0; VEC_iterate (gimple
, cases
, i
, stmt
); i
++)
496 print_gimple_stmt (file
, stmt
, 0, 0);
498 fprintf (file
, ")\n");
501 /* Dump conditions of a graphite basic block GBB on FILE. */
504 dump_gbb_conditions (FILE *file
, gimple_bb_p gbb
)
508 VEC (gimple
, heap
) *conditions
;
513 conditions
= GBB_CONDITIONS (gbb
);
514 if (VEC_empty (gimple
, conditions
))
517 fprintf (file
, "conditions bb_%d (", GBB_BB (gbb
)->index
);
519 for (i
= 0; VEC_iterate (gimple
, conditions
, i
, stmt
); i
++)
520 print_gimple_stmt (file
, stmt
, 0, 0);
522 fprintf (file
, ")\n");
525 /* Print to FILE all the data references of PBB. */
528 print_pdrs (FILE *file
, poly_bb_p pbb
)
533 for (i
= 0; VEC_iterate (poly_dr_p
, PBB_DRS (pbb
), i
, pdr
); i
++)
534 print_pdr (file
, pdr
);
537 /* Print to STDERR all the data references of PBB. */
540 debug_pdrs (poly_bb_p pbb
)
542 print_pdrs (stderr
, pbb
);
545 /* Print to FILE the domain and scattering function of PBB. */
548 print_pbb (FILE *file
, poly_bb_p pbb
)
550 fprintf (file
, "pbb (\n");
551 dump_gbb_conditions (file
, PBB_BLACK_BOX (pbb
));
552 dump_gbb_cases (file
, PBB_BLACK_BOX (pbb
));
553 print_pdrs (file
, pbb
);
554 print_pbb_domain (file
, pbb
);
555 print_scattering_function (file
, pbb
);
556 fprintf (file
, ")\n");
559 /* Print to FILE the parameters of SCOP. */
562 print_scop_params (FILE *file
, scop_p scop
)
567 fprintf (file
, "parameters (\n");
568 for (i
= 0; VEC_iterate (tree
, SESE_PARAMS (SCOP_REGION (scop
)), i
, t
); i
++)
570 fprintf (file
, "p_%d -> ", i
);
571 print_generic_expr (file
, t
, 0);
572 fprintf (file
, "\n");
574 fprintf (file
, ")\n");
577 /* Print to FILE the context of SCoP. */
579 print_scop_context (FILE *file
, scop_p scop
)
583 fprintf (file
, "context (\n");
584 fprintf (file
, "# eq");
586 for (i
= 0; i
< scop_nb_params (scop
); i
++)
587 fprintf (file
, " p%d", (int) i
);
589 fprintf (file
, " cst\n");
591 if (SCOP_CONTEXT (scop
))
592 ppl_print_powerset_matrix (file
, SCOP_CONTEXT (scop
));
594 fprintf (file
, ")\n");
597 /* Print to FILE the SCOP. */
600 print_scop (FILE *file
, scop_p scop
)
605 fprintf (file
, "scop (\n");
606 print_scop_params (file
, scop
);
607 print_scop_context (file
, scop
);
609 for (i
= 0; VEC_iterate (poly_bb_p
, SCOP_BBS (scop
), i
, pbb
); i
++)
610 print_pbb (file
, pbb
);
612 fprintf (file
, ")\n");
615 /* Print to STDERR the domain of PBB. */
618 debug_pbb_domain (poly_bb_p pbb
)
620 print_pbb_domain (stderr
, pbb
);
623 /* Print to FILE the domain and scattering function of PBB. */
626 debug_pbb (poly_bb_p pbb
)
628 print_pbb (stderr
, pbb
);
631 /* Print to STDERR the context of SCOP. */
634 debug_scop_context (scop_p scop
)
636 print_scop_context (stderr
, scop
);
639 /* Print to STDERR the SCOP. */
642 debug_scop (scop_p scop
)
644 print_scop (stderr
, scop
);
647 /* Print to STDERR the parameters of SCOP. */
650 debug_scop_params (scop_p scop
)
652 print_scop_params (stderr
, scop
);
656 /* The dimension in the transformed scattering polyhedron of PBB
657 containing the scattering iterator for the loop at depth LOOP_DEPTH. */
660 psct_scattering_dim_for_loop_depth (poly_bb_p pbb
, graphite_dim_t loop_depth
)
662 ppl_const_Constraint_System_t pcs
;
663 ppl_Constraint_System_const_iterator_t cit
, cend
;
664 ppl_const_Constraint_t cstr
;
665 ppl_Polyhedron_t ph
= PBB_TRANSFORMED_SCATTERING (pbb
);
666 ppl_dimension_type iter
= psct_iterator_dim (pbb
, loop_depth
);
667 ppl_Linear_Expression_t expr
;
668 ppl_Coefficient_t coef
;
673 ppl_new_Coefficient (&coef
);
674 ppl_Polyhedron_get_constraints (ph
, &pcs
);
675 ppl_new_Constraint_System_const_iterator (&cit
);
676 ppl_new_Constraint_System_const_iterator (&cend
);
678 for (ppl_Constraint_System_begin (pcs
, cit
),
679 ppl_Constraint_System_end (pcs
, cend
);
680 !ppl_Constraint_System_const_iterator_equal_test (cit
, cend
);
681 ppl_Constraint_System_const_iterator_increment (cit
))
683 ppl_Constraint_System_const_iterator_dereference (cit
, &cstr
);
684 ppl_new_Linear_Expression_from_Constraint (&expr
, cstr
);
685 ppl_Linear_Expression_coefficient (expr
, iter
, coef
);
686 ppl_Coefficient_to_mpz_t (coef
, val
);
688 if (value_zero_p (val
))
690 ppl_delete_Linear_Expression (expr
);
694 for (i
= 0; i
< pbb_nb_scattering_transform (pbb
); i
++)
696 ppl_dimension_type scatter
= psct_scattering_dim (pbb
, i
);
698 ppl_Linear_Expression_coefficient (expr
, scatter
, coef
);
699 ppl_Coefficient_to_mpz_t (coef
, val
);
701 if (value_notzero_p (val
))
704 ppl_delete_Linear_Expression (expr
);
705 ppl_delete_Coefficient (coef
);
706 ppl_delete_Constraint_System_const_iterator (cit
);
707 ppl_delete_Constraint_System_const_iterator (cend
);
717 /* Returns the number of iterations NITER of the loop around PBB at
721 pbb_number_of_iterations (poly_bb_p pbb
,
722 graphite_dim_t loop_depth
,
725 ppl_dimension_type loop_iter
= pbb_iterator_dim (pbb
, loop_depth
);
726 ppl_Linear_Expression_t le
;
727 ppl_Coefficient_t num
, denom
;
730 ppl_dimension_type dim
;
733 ppl_new_Coefficient (&num
);
734 ppl_new_Coefficient (&denom
);
735 ppl_Pointset_Powerset_C_Polyhedron_space_dimension (PBB_DOMAIN (pbb
), &dim
);
736 ppl_new_Linear_Expression_with_dimension (&le
, dim
);
737 ppl_set_coef (le
, loop_iter
, 1);
738 ppl_Pointset_Powerset_C_Polyhedron_maximize (PBB_DOMAIN (pbb
), le
,
739 num
, denom
, &maximum
);
743 ppl_Coefficient_to_mpz_t (num
, niter
);
744 ppl_Coefficient_to_mpz_t (denom
, dv
);
745 value_division (niter
, niter
, dv
);
748 value_set_si (niter
, -1);
751 ppl_delete_Linear_Expression (le
);
752 ppl_delete_Coefficient (num
);
753 ppl_delete_Coefficient (denom
);