isl_aff_ceil: compute ceil(e/m) as floor((e + m - 1)/m)