gcc/graphite-poly.c

   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>.
   5
   6 This file is part of GCC.
   7
   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)
  11 any later version.
  12
  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.
  17
  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/>.  */
  21 #include "config.h"
  22 #include "system.h"
  23 #include "coretypes.h"
  24 #include "tm.h"
  25 #include "ggc.h"
  26 #include "tree.h"
  27 #include "rtl.h"
  28 #include "output.h"
  29 #include "basic-block.h"
  30 #include "diagnostic.h"
  31 #include "tree-flow.h"
  32 #include "toplev.h"
  33 #include "tree-dump.h"
  34 #include "timevar.h"
  35 #include "cfgloop.h"
  36 #include "tree-chrec.h"
  37 #include "tree-data-ref.h"
  38 #include "tree-scalar-evolution.h"
  39 #include "tree-pass.h"
  40 #include "domwalk.h"
  41 #include "value-prof.h"
  42 #include "pointer-set.h"
  43 #include "gimple.h"
  44 #include "params.h"
  45
  46 #ifdef HAVE_cloog
  47 #include "cloog/cloog.h"
  48 #include "ppl_c.h"
  49 #include "sese.h"
  50 #include "graphite-ppl.h"
  51 #include "graphite.h"
  52 #include "graphite-poly.h"
  53 #include "graphite-dependences.h"
  54
  55 /* Return the maximal loop depth in SCOP.  */
  56
  57 int
  58 scop_max_loop_depth (scop_p scop)
  59 {
  60   int i;
  61   poly_bb_p pbb;
  62   int max_nb_loops = 0;
  63
  64   for (i = 0; VEC_iterate (poly_bb_p, SCOP_BBS (scop), i, pbb); i++)
  65     {
  66       int nb_loops = pbb_dim_iter_domain (pbb);
  67       if (max_nb_loops < nb_loops)
  68         max_nb_loops = nb_loops;
  69     }
  70
  71   return max_nb_loops;
  72 }
  73
  74 /* Extend the scattering matrix of PBB to MAX_SCATTERING scattering
  75    dimensions.  */
  76
  77 static void
  78 extend_scattering (poly_bb_p pbb, int max_scattering)
  79 {
  80   ppl_dimension_type nb_old_dims, nb_new_dims;
  81   int nb_added_dims, i;
  82   ppl_Coefficient_t coef;
  83   Value one;
  84
  85   nb_added_dims = max_scattering - pbb_nb_scattering_transform (pbb);
  86   value_init (one);
  87   value_set_si (one, 1);
  88   ppl_new_Coefficient (&coef);
  89   ppl_assign_Coefficient_from_mpz_t (coef, one);
  90
  91   gcc_assert (nb_added_dims >= 0);
  92
  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;
  96
  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;
 100
 101   /* Add identity matrix for the added dimensions.  */
 102   for (i = max_scattering - nb_added_dims; i < max_scattering; i++)
 103     {
 104       ppl_Constraint_t cstr;
 105       ppl_Linear_Expression_t expr;
 106
 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);
 113     }
 114
 115   ppl_delete_Coefficient (coef);
 116   value_clear (one);
 117 }
 118
 119 /* All scattering matrices in SCOP will have the same number of scattering
 120    dimensions.  */
 121
 122 int
 123 unify_scattering_dimensions (scop_p scop)
 124 {
 125   int i;
 126   poly_bb_p pbb;
 127   graphite_dim_t max_scattering = 0;
 128
 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);
 131
 132   for (i = 0; VEC_iterate (poly_bb_p, SCOP_BBS (scop), i, pbb); i++)
 133     extend_scattering (pbb, max_scattering);
 134
 135   return max_scattering;
 136 }
 137
 138 /* Prints to FILE the scattering function of PBB.  */
 139
 140 void
 141 print_scattering_function (FILE *file, poly_bb_p pbb)
 142 {
 143   graphite_dim_t i;
 144
 145   if (!PBB_TRANSFORMED (pbb))
 146     return;
 147
 148   fprintf (file, "scattering bb_%d (\n", GBB_BB (PBB_BLACK_BOX (pbb))->index);
 149   fprintf (file, "#  eq");
 150
 151   for (i = 0; i < pbb_nb_scattering_transform (pbb); i++)
 152     fprintf (file, "     s%d", (int) i);
 153
 154   for (i = 0; i < pbb_nb_local_vars (pbb); i++)
 155     fprintf (file, "    lv%d", (int) i);
 156
 157   for (i = 0; i < pbb_dim_iter_domain (pbb); i++)
 158     fprintf (file, "     i%d", (int) i);
 159
 160   for (i = 0; i < pbb_nb_params (pbb); i++)
 161     fprintf (file, "     p%d", (int) i);
 162
 163   fprintf (file, "    cst\n");
 164
 165   ppl_print_polyhedron_matrix (file, PBB_TRANSFORMED_SCATTERING (pbb));
 166
 167   fprintf (file, ")\n");
 168 }
 169
 170 /* Prints to FILE the iteration domain of PBB.  */
 171
 172 void
 173 print_iteration_domain (FILE *file, poly_bb_p pbb)
 174 {
 175   print_pbb_domain (file, pbb);
 176 }
 177
 178 /* Prints to FILE the scattering functions of every PBB of SCOP.  */
 179
 180 void
 181 print_scattering_functions (FILE *file, scop_p scop)
 182 {
 183   int i;
 184   poly_bb_p pbb;
 185
 186   for (i = 0; VEC_iterate (poly_bb_p, SCOP_BBS (scop), i, pbb); i++)
 187     print_scattering_function (file, pbb);
 188 }
 189
 190 /* Prints to FILE the iteration domains of every PBB of SCOP.  */
 191
 192 void
 193 print_iteration_domains (FILE *file, scop_p scop)
 194 {
 195   int i;
 196   poly_bb_p pbb;
 197
 198   for (i = 0; VEC_iterate (poly_bb_p, SCOP_BBS (scop), i, pbb); i++)
 199     print_iteration_domain (file, pbb);
 200 }
 201
 202 /* Prints to STDERR the scattering function of PBB.  */
 203
 204 void
 205 debug_scattering_function (poly_bb_p pbb)
 206 {
 207   print_scattering_function (stderr, pbb);
 208 }
 209
 210 /* Prints to STDERR the iteration domain of PBB.  */
 211
 212 void
 213 debug_iteration_domain (poly_bb_p pbb)
 214 {
 215   print_iteration_domain (stderr, pbb);
 216 }
 217
 218 /* Prints to STDERR the scattering functions of every PBB of SCOP.  */
 219
 220 void
 221 debug_scattering_functions (scop_p scop)
 222 {
 223   print_scattering_functions (stderr, scop);
 224 }
 225
 226 /* Prints to STDERR the iteration domains of every PBB of SCOP.  */
 227
 228 void
 229 debug_iteration_domains (scop_p scop)
 230 {
 231   print_iteration_domains (stderr, scop);
 232 }
 233
 234 /* Apply graphite transformations to all the basic blocks of SCOP.  */
 235
 236 bool
 237 apply_poly_transforms (scop_p scop)
 238 {
 239   bool transform_done = false;
 240
 241   /* Generate code even if we did not apply any real transformation.
 242      This also allows to check the performance for the identity
 243      transformation: GIMPLE -> GRAPHITE -> GIMPLE
 244      Keep in mind that CLooG optimizes in control, so the loop structure
 245      may change, even if we only use -fgraphite-identity.  */
 246   if (flag_graphite_identity)
 247     transform_done = true;
 248
 249   if (flag_loop_parallelize_all)
 250     transform_done = true;
 251
 252   if (flag_loop_block)
 253     gcc_unreachable (); /* Not yet supported.  */
 254
 255   if (flag_loop_strip_mine)
 256     transform_done |= scop_do_strip_mine (scop);
 257
 258   if (flag_loop_interchange)
 259     transform_done |= scop_do_interchange (scop);
 260
 261   return transform_done;
 262 }
 263
 264 /* Returns true when it PDR1 is a duplicate of PDR2: same PBB, and
 265    their ACCESSES, TYPE, and NB_SUBSCRIPTS are the same.  */
 266
 267 static inline bool
 268 can_collapse_pdrs (poly_dr_p pdr1, poly_dr_p pdr2)
 269 {
 270   bool res;
 271   ppl_Pointset_Powerset_C_Polyhedron_t af1, af2, diff;
 272
 273   if (PDR_PBB (pdr1) != PDR_PBB (pdr2)
 274       || PDR_NB_SUBSCRIPTS (pdr1) != PDR_NB_SUBSCRIPTS (pdr2)
 275       || PDR_TYPE (pdr1) != PDR_TYPE (pdr2))
 276     return false;
 277
 278   af1 = PDR_ACCESSES (pdr1);
 279   af2 = PDR_ACCESSES (pdr2);
 280   ppl_new_Pointset_Powerset_C_Polyhedron_from_Pointset_Powerset_C_Polyhedron
 281     (&diff, af1);
 282   ppl_Pointset_Powerset_C_Polyhedron_difference_assign (diff, af2);
 283
 284   res = ppl_Pointset_Powerset_C_Polyhedron_is_empty (diff);
 285   ppl_delete_Pointset_Powerset_C_Polyhedron (diff);
 286   return res;
 287 }
 288
 289 /* Removes duplicated data references in PBB.  */
 290
 291 void
 292 pbb_remove_duplicate_pdrs (poly_bb_p pbb)
 293 {
 294   int i, j;
 295   poly_dr_p pdr1, pdr2;
 296   unsigned n = VEC_length (poly_dr_p, PBB_DRS (pbb));
 297   VEC (poly_dr_p, heap) *collapsed = VEC_alloc (poly_dr_p, heap, n);
 298
 299   for (i = 0; VEC_iterate (poly_dr_p, PBB_DRS (pbb), i, pdr1); i++)
 300     for (j = 0; VEC_iterate (poly_dr_p, collapsed, j, pdr2); j++)
 301       if (!can_collapse_pdrs (pdr1, pdr2))
 302         VEC_quick_push (poly_dr_p, collapsed, pdr1);
 303 }
 304
 305 /* Create a new polyhedral data reference and add it to PBB.  It is
 306    defined by its ACCESSES, its TYPE, and the number of subscripts
 307    NB_SUBSCRIPTS.  */
 308
 309 void
 310 new_poly_dr (poly_bb_p pbb,
 311              ppl_Pointset_Powerset_C_Polyhedron_t accesses,
 312              enum poly_dr_type type, void *cdr, graphite_dim_t nb_subscripts)
 313 {
 314   static int id = 0;
 315   poly_dr_p pdr = XNEW (struct poly_dr);
 316
 317   PDR_ID (pdr) = id++;
 318   PDR_NB_REFS (pdr) = 1;
 319   PDR_PBB (pdr) = pbb;
 320   PDR_ACCESSES (pdr) = accesses;
 321   PDR_TYPE (pdr) = type;
 322   PDR_CDR (pdr) = cdr;
 323   PDR_NB_SUBSCRIPTS (pdr) = nb_subscripts;
 324   VEC_safe_push (poly_dr_p, heap, PBB_DRS (pbb), pdr);
 325 }
 326
 327 /* Free polyhedral data reference PDR.  */
 328
 329 void
 330 free_poly_dr (poly_dr_p pdr)
 331 {
 332   ppl_delete_Pointset_Powerset_C_Polyhedron (PDR_ACCESSES (pdr));
 333   XDELETE (pdr);
 334 }
 335
 336 /* Create a new polyhedral black box.  */
 337
 338 void
 339 new_poly_bb (scop_p scop, void *black_box)
 340 {
 341   poly_bb_p pbb = XNEW (struct poly_bb);
 342
 343   PBB_DOMAIN (pbb) = NULL;
 344   PBB_SCOP (pbb) = scop;
 345   pbb_set_black_box (pbb, black_box);
 346   PBB_TRANSFORMED (pbb) = NULL;
 347   PBB_SAVED (pbb) = NULL;
 348   PBB_ORIGINAL (pbb) = NULL;
 349   PBB_DRS (pbb) = VEC_alloc (poly_dr_p, heap, 3);
 350   VEC_safe_push (poly_bb_p, heap, SCOP_BBS (scop), pbb);
 351 }
 352
 353 /* Free polyhedral black box.  */
 354
 355 void
 356 free_poly_bb (poly_bb_p pbb)
 357 {
 358   int i;
 359   poly_dr_p pdr;
 360
 361   ppl_delete_Pointset_Powerset_C_Polyhedron (PBB_DOMAIN (pbb));
 362
 363   if (PBB_TRANSFORMED (pbb))
 364     poly_scattering_free (PBB_TRANSFORMED (pbb));
 365
 366   if (PBB_SAVED (pbb))
 367     poly_scattering_free (PBB_SAVED (pbb));
 368
 369   if (PBB_ORIGINAL (pbb))
 370     poly_scattering_free (PBB_ORIGINAL (pbb));
 371
 372   if (PBB_DRS (pbb))
 373     for (i = 0; VEC_iterate (poly_dr_p, PBB_DRS (pbb), i, pdr); i++)
 374       free_poly_dr (pdr);
 375
 376   VEC_free (poly_dr_p, heap, PBB_DRS (pbb));
 377   XDELETE (pbb);
 378 }
 379
 380 static void
 381 print_pdr_access_layout (FILE *file, poly_dr_p pdr)
 382 {
 383   graphite_dim_t i;
 384
 385   fprintf (file, "#  eq");
 386
 387   for (i = 0; i < pdr_dim_iter_domain (pdr); i++)
 388     fprintf (file, "     i%d", (int) i);
 389
 390   for (i = 0; i < pdr_nb_params (pdr); i++)
 391     fprintf (file, "     p%d", (int) i);
 392
 393   fprintf (file, "  alias");
 394
 395   for (i = 0; i < PDR_NB_SUBSCRIPTS (pdr); i++)
 396     fprintf (file, "   sub%d", (int) i);
 397
 398   fprintf (file, "    cst\n");
 399 }
 400
 401 /* Prints to FILE the polyhedral data reference PDR.  */
 402
 403 void
 404 print_pdr (FILE *file, poly_dr_p pdr)
 405 {
 406   fprintf (file, "pdr_%d (", PDR_ID (pdr));
 407
 408   switch (PDR_TYPE (pdr))
 409     {
 410     case PDR_READ:
 411       fprintf (file, "read \n");
 412       break;
 413
 414     case PDR_WRITE:
 415       fprintf (file, "write \n");
 416       break;
 417
 418     case PDR_MAY_WRITE:
 419       fprintf (file, "may_write \n");
 420       break;
 421
 422     default:
 423       gcc_unreachable ();
 424     }
 425
 426   dump_data_reference (file, (data_reference_p) PDR_CDR (pdr));
 427
 428   fprintf (file, "data accesses (\n");
 429   print_pdr_access_layout (file, pdr);
 430   ppl_print_powerset_matrix (file, PDR_ACCESSES (pdr));
 431   fprintf (file, ")\n");
 432
 433   fprintf (file, ")\n");
 434 }
 435
 436 /* Prints to STDERR the polyhedral data reference PDR.  */
 437
 438 void
 439 debug_pdr (poly_dr_p pdr)
 440 {
 441   print_pdr (stderr, pdr);
 442 }
 443
 444 /* Creates a new SCOP containing REGION.  */
 445
 446 scop_p
 447 new_scop (void *region)
 448 {
 449   scop_p scop = XNEW (struct scop);
 450
 451   SCOP_DEP_GRAPH (scop) = NULL;
 452   SCOP_CONTEXT (scop) = NULL;
 453   scop_set_region (scop, region);
 454   SCOP_BBS (scop) = VEC_alloc (poly_bb_p, heap, 3);
 455   SCOP_ORIGINAL_PDDRS (scop) = htab_create (10, hash_poly_ddr_p,
 456                                             eq_poly_ddr_p, free_poly_ddr);
 457   return scop;
 458 }
 459
 460 /* Deletes SCOP.  */
 461
 462 void
 463 free_scop (scop_p scop)
 464 {
 465   int i;
 466   poly_bb_p pbb;
 467
 468   for (i = 0; VEC_iterate (poly_bb_p, SCOP_BBS (scop), i, pbb); i++)
 469     free_poly_bb (pbb);
 470
 471   VEC_free (poly_bb_p, heap, SCOP_BBS (scop));
 472
 473   if (SCOP_CONTEXT (scop))
 474     ppl_delete_Pointset_Powerset_C_Polyhedron (SCOP_CONTEXT (scop));
 475
 476   htab_delete (SCOP_ORIGINAL_PDDRS (scop));
 477   XDELETE (scop);
 478 }
 479
 480 /* Print to FILE the domain of PBB.  */
 481
 482 void
 483 print_pbb_domain (FILE *file, poly_bb_p pbb)
 484 {
 485   graphite_dim_t i;
 486   gimple_bb_p gbb = PBB_BLACK_BOX (pbb);
 487
 488   if (!PBB_DOMAIN (pbb))
 489     return;
 490
 491   fprintf (file, "domains bb_%d (\n", GBB_BB (gbb)->index);
 492   fprintf (file, "#  eq");
 493
 494   for (i = 0; i < pbb_dim_iter_domain (pbb); i++)
 495     fprintf (file, "     i%d", (int) i);
 496
 497   for (i = 0; i < pbb_nb_params (pbb); i++)
 498     fprintf (file, "     p%d", (int) i);
 499
 500   fprintf (file, "    cst\n");
 501
 502   if (PBB_DOMAIN (pbb))
 503     ppl_print_powerset_matrix (file, PBB_DOMAIN (pbb));
 504
 505   fprintf (file, ")\n");
 506 }
 507
 508 /* Dump the cases of a graphite basic block GBB on FILE.  */
 509
 510 static void
 511 dump_gbb_cases (FILE *file, gimple_bb_p gbb)
 512 {
 513   int i;
 514   gimple stmt;
 515   VEC (gimple, heap) *cases;
 516
 517   if (!gbb)
 518     return;
 519
 520   cases = GBB_CONDITION_CASES (gbb);
 521   if (VEC_empty (gimple, cases))
 522     return;
 523
 524   fprintf (file, "cases bb_%d (", GBB_BB (gbb)->index);
 525
 526   for (i = 0; VEC_iterate (gimple, cases, i, stmt); i++)
 527     print_gimple_stmt (file, stmt, 0, 0);
 528
 529   fprintf (file, ")\n");
 530 }
 531
 532 /* Dump conditions of a graphite basic block GBB on FILE.  */
 533
 534 static void
 535 dump_gbb_conditions (FILE *file, gimple_bb_p gbb)
 536 {
 537   int i;
 538   gimple stmt;
 539   VEC (gimple, heap) *conditions;
 540
 541   if (!gbb)
 542     return;
 543
 544   conditions = GBB_CONDITIONS (gbb);
 545   if (VEC_empty (gimple, conditions))
 546     return;
 547
 548   fprintf (file, "conditions bb_%d (", GBB_BB (gbb)->index);
 549
 550   for (i = 0; VEC_iterate (gimple, conditions, i, stmt); i++)
 551     print_gimple_stmt (file, stmt, 0, 0);
 552
 553   fprintf (file, ")\n");
 554 }
 555
 556 /* Print to FILE all the data references of PBB.  */
 557
 558 void
 559 print_pdrs (FILE *file, poly_bb_p pbb)
 560 {
 561   int i;
 562   poly_dr_p pdr;
 563
 564   for (i = 0; VEC_iterate (poly_dr_p, PBB_DRS (pbb), i, pdr); i++)
 565     print_pdr (file, pdr);
 566 }
 567
 568 /* Print to STDERR all the data references of PBB.  */
 569
 570 void
 571 debug_pdrs (poly_bb_p pbb)
 572 {
 573   print_pdrs (stderr, pbb);
 574 }
 575
 576 /* Print to FILE the domain and scattering function of PBB.  */
 577
 578 void
 579 print_pbb (FILE *file, poly_bb_p pbb)
 580 {
 581   fprintf (file, "pbb_%d (\n", GBB_BB (PBB_BLACK_BOX (pbb))->index);
 582   dump_gbb_conditions (file, PBB_BLACK_BOX (pbb));
 583   dump_gbb_cases (file, PBB_BLACK_BOX (pbb));
 584   print_pdrs (file, pbb);
 585   print_pbb_domain (file, pbb);
 586   print_scattering_function (file, pbb);
 587   fprintf (file, ")\n");
 588 }
 589
 590 /* Print to FILE the parameters of SCOP.  */
 591
 592 void
 593 print_scop_params (FILE *file, scop_p scop)
 594 {
 595   int i;
 596   tree t;
 597
 598   fprintf (file, "parameters (\n");
 599   for (i = 0; VEC_iterate (tree, SESE_PARAMS (SCOP_REGION (scop)), i, t); i++)
 600     {
 601       fprintf (file, "p_%d -> ", i);
 602       print_generic_expr (file, t, 0);
 603       fprintf (file, "\n");
 604     }
 605   fprintf (file, ")\n");
 606 }
 607
 608 /* Print to FILE the context of SCoP.  */
 609 void
 610 print_scop_context (FILE *file, scop_p scop)
 611 {
 612   graphite_dim_t i;
 613
 614   fprintf (file, "context (\n");
 615   fprintf (file, "#  eq");
 616
 617   for (i = 0; i < scop_nb_params (scop); i++)
 618     fprintf (file, "     p%d", (int) i);
 619
 620   fprintf (file, "    cst\n");
 621
 622   if (SCOP_CONTEXT (scop))
 623     ppl_print_powerset_matrix (file, SCOP_CONTEXT (scop));
 624
 625   fprintf (file, ")\n");
 626 }
 627
 628 /* Print to FILE the SCOP.  */
 629
 630 void
 631 print_scop (FILE *file, scop_p scop)
 632 {
 633   int i;
 634   poly_bb_p pbb;
 635
 636   fprintf (file, "scop (\n");
 637   print_scop_params (file, scop);
 638   print_scop_context (file, scop);
 639
 640   for (i = 0; VEC_iterate (poly_bb_p, SCOP_BBS (scop), i, pbb); i++)
 641     print_pbb (file, pbb);
 642
 643   fprintf (file, ")\n");
 644 }
 645
 646 /* Print to STDERR the domain of PBB.  */
 647
 648 void
 649 debug_pbb_domain (poly_bb_p pbb)
 650 {
 651   print_pbb_domain (stderr, pbb);
 652 }
 653
 654 /* Print to FILE the domain and scattering function of PBB.  */
 655
 656 void
 657 debug_pbb (poly_bb_p pbb)
 658 {
 659   print_pbb (stderr, pbb);
 660 }
 661
 662 /* Print to STDERR the context of SCOP.  */
 663
 664 void
 665 debug_scop_context (scop_p scop)
 666 {
 667   print_scop_context (stderr, scop);
 668 }
 669
 670 /* Print to STDERR the SCOP.  */
 671
 672 void
 673 debug_scop (scop_p scop)
 674 {
 675   print_scop (stderr, scop);
 676 }
 677
 678 /* Print to STDERR the parameters of SCOP.  */
 679
 680 void
 681 debug_scop_params (scop_p scop)
 682 {
 683   print_scop_params (stderr, scop);
 684 }
 685
 686
 687 /* The dimension in the transformed scattering polyhedron of PBB
 688    containing the scattering iterator for the loop at depth LOOP_DEPTH.  */
 689
 690 ppl_dimension_type
 691 psct_scattering_dim_for_loop_depth (poly_bb_p pbb, graphite_dim_t loop_depth)
 692 {
 693   ppl_const_Constraint_System_t pcs;
 694   ppl_Constraint_System_const_iterator_t cit, cend;
 695   ppl_const_Constraint_t cstr;
 696   ppl_Polyhedron_t ph = PBB_TRANSFORMED_SCATTERING (pbb);
 697   ppl_dimension_type iter = psct_iterator_dim (pbb, loop_depth);
 698   ppl_Linear_Expression_t expr;
 699   ppl_Coefficient_t coef;
 700   Value val;
 701   graphite_dim_t i;
 702
 703   value_init (val);
 704   ppl_new_Coefficient (&coef);
 705   ppl_Polyhedron_get_constraints (ph, &pcs);
 706   ppl_new_Constraint_System_const_iterator (&cit);
 707   ppl_new_Constraint_System_const_iterator (&cend);
 708
 709   for (ppl_Constraint_System_begin (pcs, cit),
 710          ppl_Constraint_System_end (pcs, cend);
 711        !ppl_Constraint_System_const_iterator_equal_test (cit, cend);
 712        ppl_Constraint_System_const_iterator_increment (cit))
 713     {
 714       ppl_Constraint_System_const_iterator_dereference (cit, &cstr);
 715       ppl_new_Linear_Expression_from_Constraint (&expr, cstr);
 716       ppl_Linear_Expression_coefficient (expr, iter, coef);
 717       ppl_Coefficient_to_mpz_t (coef, val);
 718
 719       if (value_zero_p (val))
 720         {
 721           ppl_delete_Linear_Expression (expr);
 722           continue;
 723         }
 724
 725       for (i = 0; i < pbb_nb_scattering_transform (pbb); i++)
 726         {
 727           ppl_dimension_type scatter = psct_scattering_dim (pbb, i);
 728
 729           ppl_Linear_Expression_coefficient (expr, scatter, coef);
 730           ppl_Coefficient_to_mpz_t (coef, val);
 731
 732           if (value_notzero_p (val))
 733             {
 734               value_clear (val);
 735               ppl_delete_Linear_Expression (expr);
 736               ppl_delete_Coefficient (coef);
 737               ppl_delete_Constraint_System_const_iterator (cit);
 738               ppl_delete_Constraint_System_const_iterator (cend);
 739
 740               return scatter;
 741             }
 742         }
 743     }
 744
 745   gcc_unreachable ();
 746 }
 747
 748 /* Returns the number of iterations NITER of the loop around PBB at
 749    depth LOOP_DEPTH.  */
 750
 751 void
 752 pbb_number_of_iterations (poly_bb_p pbb,
 753                           graphite_dim_t loop_depth,
 754                           Value niter)
 755 {
 756   ppl_Linear_Expression_t le;
 757   ppl_dimension_type dim;
 758
 759   ppl_Pointset_Powerset_C_Polyhedron_space_dimension (PBB_DOMAIN (pbb), &dim);
 760   ppl_new_Linear_Expression_with_dimension (&le, dim);
 761   ppl_set_coef (le, pbb_iterator_dim (pbb, loop_depth), 1);
 762   value_set_si (niter, -1);
 763   ppl_max_for_le_pointset (PBB_DOMAIN (pbb), le, niter);
 764   ppl_delete_Linear_Expression (le);
 765 }
 766
 767 /* Returns the number of iterations NITER of the loop around PBB at
 768    time(scattering) dimension TIME_DEPTH.  */
 769
 770 void
 771 pbb_number_of_iterations_at_time (poly_bb_p pbb,
 772                                   graphite_dim_t time_depth,
 773                                   Value niter)
 774 {
 775   ppl_Pointset_Powerset_C_Polyhedron_t ext_domain, sctr;
 776   ppl_Linear_Expression_t le;
 777   ppl_dimension_type dim;
 778
 779   value_set_si (niter, -1);
 780
 781   /* Takes together domain and scattering polyhedrons, and composes
 782      them into the bigger polyhedron that has the following format:
 783      t0..t_{n-1} | l0..l_{nlcl-1} | i0..i_{niter-1} | g0..g_{nparm-1}.
 784      t0..t_{n-1} are time dimensions (scattering dimensions)
 785      l0..l_{nclc-1} are local variables in scatterin function
 786      i0..i_{niter-1} are original iteration variables
 787      g0..g_{nparam-1} are global parameters.  */
 788
 789   ppl_new_Pointset_Powerset_C_Polyhedron_from_Pointset_Powerset_C_Polyhedron
 790     (&ext_domain, PBB_DOMAIN (pbb));
 791   ppl_insert_dimensions_pointset (ext_domain, 0,
 792                                   pbb_nb_scattering_transform (pbb)
 793                                   + pbb_nb_local_vars (pbb));
 794   ppl_new_Pointset_Powerset_C_Polyhedron_from_C_Polyhedron (&sctr,
 795       PBB_TRANSFORMED_SCATTERING (pbb));
 796   ppl_Pointset_Powerset_C_Polyhedron_intersection_assign (sctr, ext_domain);
 797
 798   ppl_Pointset_Powerset_C_Polyhedron_space_dimension (sctr, &dim);
 799   ppl_new_Linear_Expression_with_dimension (&le, dim);
 800   ppl_set_coef (le, time_depth, 1);
 801   ppl_max_for_le_pointset (sctr, le, niter);
 802
 803   ppl_delete_Linear_Expression (le);
 804   ppl_delete_Pointset_Powerset_C_Polyhedron (sctr);
 805   ppl_delete_Pointset_Powerset_C_Polyhedron (ext_domain);
 806 }
 807
 808 #endif
 809