On Mon, Aug 15, 2011 at 02:12:54AM -0500, Sebastian Pop wrote: > @@ -160,6 +182,14 @@ pbb_strip_mine_time_depth (poly_bb_p pbb, int > time_depth, int stride) > ppl_delete_Linear_Expression (expr); > ppl_Polyhedron_add_constraint (res, new_cstr); > ppl_delete_Constraint (new_cstr); > + > + { > + isl_dim *d = isl_map_get_dim (pbb->transformed); > + isl_constraint *c = isl_equality_alloc (d); > + > + c = isl_constraint_set_coefficient_si (c, isl_dim_out, strip + 1, 1); > + pbb->transformed = isl_map_add_constraint (pbb->transformed, c); > + } > } > }
I'm not sure what the point is of introducing an extra dimension fixed to zero, but it can be accomplished more easily using isl_map_fix_si. Also, it may clarify the code if you first construct the transformation and then apply it, like you do in the interchange. skimo