llvm dev - Jul 2016 - Regarding ScalarEvolution's loop backedge computation

Hi Pankaj,

Chawla, Pankaj via llvm-dev wrote:
 > It looks like ScalarEvolution bails out of loop backedge computation if
 > it cannot prove the IV stride as either positive or negative (based on
 > loop control condition). I think this logic can be refined for signed 
IVs.
 >
 > Consider this simple loop-
 >
 > void foo(int *A, int n, int s) {
 >
 > int i;
 >
 > for(i=0; i<n; i += s) {
 >
 > A[i]++;
 >
 > }
 >
 > }
 >
 > The IV of this loop has this SCEV form-
 >
 > {0,+,%s}<nsw><%for.body>

This looks valid -- we already do things like this for

for (i = A; i != B; i += 5)
   ...

and compute the backedge taken count as "(B - A) / 5" (roughly :) )
since if (B - A) is not divisible by 5 then we have UB due to
overflow.  We just have to be careful around cases like:

   for(i = 0; i < 60; i += s) {
     may_exit();
   }

"s" can be (say) -3 and the loop can take 160 backedges and then
"exit(0)", avoiding the undefined behavior due to underflow. 
"s" can
also be zero, in which case the loop can potentially take an infinite
number of backedges.

However, in the example you gave (written in LLVM's canonical rotated
form):

   if (0 < N) {
     i = 0;
     do {
       a[i]++;
       i += s; // NSW add
     } while (i < N);
   }

For any s <= 0 we have undefined behavior, so it is sound to assume s
 > 0.

Do you want to take a crack at fixing this?  I'm traveling till the
10th of July, but I can review your change once I'm back.

-- Sanjoy

> On Jul 5, 2016, at 6:01 PM, Chawla, Pankaj via llvm-dev <llvm-dev at
lists.llvm.org> wrote:
> 
> Hi Sanjoy,
> 
> The following trivial change in howManyLessThans() seems to resolve the
problem with the original loop-
> 
>   // Avoid negative or zero stride values
> -  if (!isKnownPositive(Stride))
> +  if (!NoWrap && !isKnownPositive(Stride))
>     return getCouldNotCompute();
> 
> However, I was experimenting with a few variants of the loop I posted and
they seem to have different issues which may require more involved fixes. I am
listing them here-
> 
> 1) Changing the loop control condition from '<' to
'<='.
> 
> The canonical form of this loop is something like this-
> 
> if (0 < N) {
>  i = 0;
>  do {
>    a[i]++;
>    i += s; // NSW add
>  } while (! (i > N));  // sgt compare
> }
> 
> The 'sgt' compare is inverted to 'sle' for analysis.
ScalarEvolution isn't really expecting 'sle' in canonicalized loops
so it reverts to brute force exit count computation using
computeExitCountExhaustively() which doesn't work. This looks like a
canonicalization issue.
> 
> 2) Variants with '>' and '>='. For example-
> 
>  for(i=n; i>=0; i-=s) {
>    A[i]++;
>  }
> 
> In this case the SCEV form of IV does not have 'nsw' flag- 
> {%n,+,(-1 * %s)}<%for.body>
> 
> 
> For now, I can submit a patch which fixes the issue with the original loop.
> Please let me know how to proceed.
If you’re asking about how to submit a patch to LLVM, this may help:
http://llvm.org/docs/DeveloperPolicy.html#making-and-submitting-a-patch
<http://llvm.org/docs/DeveloperPolicy.html#making-and-submitting-a-patch>
(And possibly this: http://llvm.org/docs/Phabricator.html
<http://llvm.org/docs/Phabricator.html> )

— 
Mehdi

> 
> Thanks,
> Pankaj
> 
> -----Original Message-----
> From: Chawla, Pankaj 
> Sent: Thursday, June 30, 2016 12:13 PM
> To: 'Sanjoy Das'
> Cc: llvm-dev at lists.llvm.org
> Subject: RE: [llvm-dev] Regarding ScalarEvolution's loop backedge
computation
> 
> Hi Sanjoy,
> 
> Thank you for the clarification!
> I will give it a try and put up the changes for review.
> 
> -Pankaj
> 
> -----Original Message-----
> From: Sanjoy Das [mailto:sanjoy at playingwithpointers.com] 
> Sent: Wednesday, June 29, 2016 5:02 PM
> To: Chawla, Pankaj
> Cc: llvm-dev at lists.llvm.org
> Subject: Re: [llvm-dev] Regarding ScalarEvolution's loop backedge
computation
> 
> Hi Pankaj,
> 
> Chawla, Pankaj via llvm-dev wrote:
>> It looks like ScalarEvolution bails out of loop backedge computation if
> it cannot prove the IV stride as either positive or negative (based on 
> loop control condition). I think this logic can be refined for signed IVs.
>> 
>> Consider this simple loop-
>> 
>> void foo(int *A, int n, int s) {
>> 
>> int i;
>> 
>> for(i=0; i<n; i += s) {
>> 
>> A[i]++;
>> 
>> }
>> 
>> }
>> 
>> The IV of this loop has this SCEV form-  >  >
{0,+,%s}<nsw><%for.body>
> 
> This looks valid -- we already do things like this for
> 
> for (i = A; i != B; i += 5)
>   ...
> 
> and compute the backedge taken count as "(B - A) / 5" (roughly :)
) since if (B - A) is not divisible by 5 then we have UB due to overflow.  We
just have to be careful around cases like:
> 
>   for(i = 0; i < 60; i += s) {
>     may_exit();
>   }
> 
> "s" can be (say) -3 and the loop can take 160 backedges and then
"exit(0)", avoiding the undefined behavior due to underflow. 
"s" can also be zero, in which case the loop can potentially take an
infinite number of backedges.
> 
> However, in the example you gave (written in LLVM's canonical rotated
> form):
> 
>   if (0 < N) {
>     i = 0;
>     do {
>       a[i]++;
>       i += s; // NSW add
>     } while (i < N);
>   }
> 
> For any s <= 0 we have undefined behavior, so it is sound to assume s 
> 0.
> 
> Do you want to take a crack at fixing this?  I'm traveling till the
10th of July, but I can review your change once I'm back.
> 
> -- Sanjoy
> _______________________________________________
> LLVM Developers mailing list
> llvm-dev at lists.llvm.org
> http://lists.llvm.org/cgi-bin/mailman/listinfo/llvm-dev
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://lists.llvm.org/pipermail/llvm-dev/attachments/20160706/1ffefb26/attachment.html>