[InstCombine] fadd double (sitofp x), y check that the promotion is valid

[apilipenko]

Summary

Doing these transformations check that the result of integer addition is representable in the FP type.
(fadd double (sitofp x), fpcst) --> (sitofp (add int x, intcst))
(fadd double (sitofp x), (sitofp y)) --> (sitofp (add int x, y))

This is a fix for https://bugs.llvm.org//show_bug.cgi?id=27036

[boris.ulasevich]

test/Transforms/InstCombine/add-sitofp.ll

16

Why do we need lshr here? %n + 1 = ((%m >>> 24) && %b) + 1 takes range 1..256 - it fits quite well to 23-bit single precision mantissa, and optimisation can be applied here, isn't it?
Does it make sense for you? Do I miss something?

[spatel]

Before getting to any details about the patch, we need to address the question raised in PR27036: why are we doing this transform in InstCombine at all? The assumption is that an integer add is more canonical and/or cheaper than an FP add. Is that universally true?

[arsenm]

In D31182#706365, @spatel wrote:

Before getting to any details about the patch, we need to address the question raised in PR27036: why are we doing this transform in InstCombine at all? The assumption is that an integer add is more canonical and/or cheaper than an FP add. Is that universally true?

I would say no. For example on AMDGPU a f32 add is exactly as fast as an i32 add, but the i32 add has an additional carry output constraint to deal with

[scanon]

The assumption is that an integer add is more canonical and/or cheaper than an FP add. Is that universally true?

No. There exist GPU architectures where a 32b integer addition is a more expensive operation than a 32b FP add.

lib/Transforms/InstCombine/InstCombineAddSub.cpp

1401

s/mantissa/significand/ (https://en.wikipedia.org/wiki/Significand#Use_of_.22mantissa.22)

Also, semanticsPrecision returns the number of bits in the significand (which includes the implicit "integral" bit when present). There's no "plus one for the sign bit".

[scanon]

lib/Transforms/InstCombine/InstCombineAddSub.cpp

1429

Aren't we missing a check that RHSIntVal can be exactly in the floating-point type here? I.e. if (RHSIntVal->getType()->getIntegerBitWidth() <= MaxRepresentableBits), mirroring the check for LHS above?

[andrew.w.kaylor]

lib/Transforms/InstCombine/InstCombineAddSub.cpp

1404

This check seems more conservative than we'd want it to be. For instance, in the example in the comments we know that 'x & 1234' is can be represented as a 12-bit signed integer, but if x is i32 this code will return 32 for LHSIntVal->getType()->getIntegerBitWidth(). This directly relates to Boris' comment on the test case.

1410

The MaxRepresentableBits calculation and comparison won't be necessary if neither this condition nor the condition at line 1424 is true. I think the code should be organized to check these dyn_cast conditions first, particularly if you are going to do something more complicated with the number of bits required.

1413

This condition will block either transformation and so it could be checked at line 1395 and save us a bit of work sometimes.

test/Transforms/InstCombine/add-sitofp.ll

21

Can you add a negative test for the '(fadd double (sitofp x), (sitofp y)) --> (sitofp (add int x, y))' case also?

[andrew.w.kaylor]

lib/Transforms/InstCombine/InstCombineAddSub.cpp

1414

It looks like getFPToSI and getSIToFP will end up using the APFloat class. Wouldn't it be better to just use APFloat here and call APFloat::isInteger()?

[apilipenko]

In D31182#706365, @spatel wrote:

Before getting to any details about the patch, we need to address the question raised in PR27036: why are we doing this transform in InstCombine at all? The assumption is that an integer add is more canonical and/or cheaper than an FP add. Is that universally true?

It's difficult to tell. Even if FP addition is cheaper the second combine also saves us one sitofp conversion. It should be also taken into account.

These combines are buggy and there is a simple fix. Since I don't have a good understanding on the profitability of these combines I decided to go ahead with the fix rather then removing them for example.

[spatel]

In D31182#707598, @apilipenko wrote:

In D31182#706365, @spatel wrote:

Before getting to any details about the patch, we need to address the question raised in PR27036: why are we doing this transform in InstCombine at all? The assumption is that an integer add is more canonical and/or cheaper than an FP add. Is that universally true?

It's difficult to tell. Even if FP addition is cheaper the second combine also saves us one sitofp conversion. It should be also taken into account.

These combines are buggy and there is a simple fix. Since I don't have a good understanding on the profitability of these combines I decided to go ahead with the fix rather then removing them for example.

Agreed - solving the miscompiling is the most important thing. Since we can do that quickly with a small patch, let's go ahead. We can decide what is canonical form as a next step, but please add a 'TODO' comment about that.

This example leads into the same gray area we have with several instcombine folds: it's not clear if the value-track-ability is enough to justify the transform, whether a backend should be responsible for fixing something that is not optimal for all targets, if the removal of an IR instruction overrides either of those...

[apilipenko]

Address review comments.

lib/Transforms/InstCombine/InstCombineAddSub.cpp

1404

Given that we are uncertain if we need this transform at all (see the discussion above) I'd prefer not to make it too smart without a good reason. So, unless we have a clear reason why we want this code to be more aggressive I'd prefer to keep this fix as simple as it is now.

For more aggressive version of this transform we should rely on Float2Int pass and make it smarter if needed.

1413

The second transform looks for either RHSConv or LHSConv to have one use, so it doesn't always block the second transform.

1414

It sounds like a separate refactoring to the existing code.

1429

We should be fine as is. I assume that the LHS, RHS of the addition have the same types (there can be an assert) and there is also a check below that integer types are the same.

test/Transforms/InstCombine/add-sitofp.ll

16

I just duplicated the example above. I agree it is confusing in this context. I'll simplify the negative test case to make it more clear that the result of the operation is not guaranteed to fit into float type.

[andrew.w.kaylor]

lib/Transforms/InstCombine/InstCombineAddSub.cpp

1404

That seems reasonable, but can you add a comment indicating that it's being conservative and that the example transformation might not happen?

1413

I see. I missed that.

1414

Yes, that's true.

[apilipenko]

Ping.

Update, I didn't see that there was a request to add a comment. Will do and update the patch,

[scanon]

This addresses my concerns. Gentle ping to the other reviewers.

[spatel]

See inline for some nits.

lib/Transforms/InstCombine/InstCombineAddSub.cpp

1402–1403

Is that clang-formatted?

Please upload the version with the TODO/FIXME comments.

test/Transforms/InstCombine/add-sitofp.ll

17

"highes" -> "highest"; same typo is repeated below.

Also, I appreciate changing this test file to use FileCheck, but could you run utils/update_test_checks.py instead of hand-editing the CHECK lines? That makes the assertions stronger and easier to update (we're pretty sure we'll want to update these one way or another).

[apilipenko]

Add the comment about overly conservative legality check, add the test demonstrating it.

[spatel]

The code changes look good to me, but I'm not sure if the tests cover the case that Andy requested, so I'll let him comment.

[apilipenko]

@andrew.w.kaylor - ping, do you have any comments about the latest revision?

[andrew.w.kaylor]

I'm happy with this. Thanks for the improvements!

lib/Transforms/InstCombine/InstCombineAddSub.cpp

if (Value *V = dyn_castFNegVal(RHS)) {
       RI->copyFastMathFlags(&I);
       return RI;
     }
 
   // Check for (fadd double (sitofp x), y), see if we can merge this into an
   // integer add followed by a promotion.
   if (SIToFPInst *LHSConv = dyn_cast<SIToFPInst>(LHS)) {
     Value *LHSIntVal = LHSConv->getOperand(0);
+    Type *FPType = LHSConv->getType();
+
+    auto IsValidPromotion = [] (Type *FTy, Type *ITy) {
+      // Do we have enough bits in the significand to represent the result of the
+      // integer addition?
+      unsigned MaxRepresentableBits
+        = APFloat::semanticsPrecision(FTy->getFltSemantics());
+      return ITy->getIntegerBitWidth() <= MaxRepresentableBits;
+    };
 
     // (fadd double (sitofp x), fpcst) --> (sitofp (add int x, intcst))
     // ... if the constant fits in the integer value.  This is useful for things
     // like (double)(x & 1234) + 4.0 -> (double)((X & 1234)+4) which no longer
     // requires a constant pool load, and generally allows the add to be better
     // instcombined.
-    if (ConstantFP *CFP = dyn_cast<ConstantFP>(RHS)) {
+    if (ConstantFP *CFP = dyn_cast<ConstantFP>(RHS))
+      if (IsValidPromotion(FPType, LHSIntVal->getType())) {
         Constant *CI =
           ConstantExpr::getFPToSI(CFP, LHSIntVal->getType());
         if (LHSConv->hasOneUse() &&
             ConstantExpr::getSIToFP(CI, I.getType()) == CFP &&
             WillNotOverflowSignedAdd(LHSIntVal, CI, I)) {
           // Insert the new integer add.
           Value *NewAdd = Builder->CreateNSWAdd(LHSIntVal,
                                                 CI, "addconv");
           return new SIToFPInst(NewAdd, I.getType());
         }
       }
 
     // (fadd double (sitofp x), (sitofp y)) --> (sitofp (add int x, y))
     if (SIToFPInst *RHSConv = dyn_cast<SIToFPInst>(RHS)) {
       Value *RHSIntVal = RHSConv->getOperand(0);
-
+      // It's enough to check LHS types only because we require int types to
+      // be the same for this transform.
+      if (IsValidPromotion(FPType, LHSIntVal->getType())) {
         // Only do this if x/y have the same type, if at last one of them has a
         // single use (so we don't increase the number of int->fp conversions),
         // and if the integer add will not overflow.
         if (LHSIntVal->getType() == RHSIntVal->getType() &&
             (LHSConv->hasOneUse() || RHSConv->hasOneUse()) &&
             WillNotOverflowSignedAdd(LHSIntVal, RHSIntVal, I)) {
           // Insert the new integer add.
           Value *NewAdd = Builder->CreateNSWAdd(LHSIntVal,
                                                 RHSIntVal, "addconv");
           return new SIToFPInst(NewAdd, I.getType());
         }
       }
     }
+  }
 
   // select C, 0, B + select C, A, 0 -> select C, A, B
   {
     Value *A1, *B1, *C1, *A2, *B2, *C2;
     if (match(LHS, m_Select(m_Value(C1), m_Value(A1), m_Value(B1))) &&
         match(RHS, m_Select(m_Value(C2), m_Value(A2), m_Value(B2)))) {
       if (C1 == C2) {
         Constant *Z1=nullptr, *Z2=nullptr;

test/Transforms/InstCombine/add-sitofp.ll

-; RUN: opt < %s -instcombine -S | grep "add nuw nsw i32"
+; RUN: opt < %s -instcombine -S | FileCheck %s
 
 define double @x(i32 %a, i32 %b) nounwind {
+; CHECK-LABEL: @x
+; CHECK: add nuw nsw i32
   %m = lshr i32 %a, 24
   %n = and i32 %m, %b
   %o = sitofp i32 %n to double
   %p = fadd double %o, 1.0
   ret double %p
 }
+
+define double @test(i32 %a) nounwind {
+; CHECK-LABEL: @test
+; CHECK: add nuw nsw i32
+
+  ; Drop two highes bits to guarantee that %a + 1 doesn't overflow
+  %a_and = and i32 %a, 1073741823
+  %a_and_fp = sitofp i32 %a_and to double
+  %res = fadd double %a_and_fp, 1.0
+  ret double %res
+}
+
+define float @test_neg(i32 %a) nounwind {
+; CHECK-LABEL: @test_neg
+; CHECK: fadd float
+
+  ; Drop two highes bits to guarantee that %a + 1 doesn't overflow
+  %a_and = and i32 %a, 1073741823
+  %a_and_fp = sitofp i32 %a_and to float
+  %res = fadd float %a_and_fp, 1.0
+  ret float %res
+}
+
+define double @test_2(i32 %a, i32 %b) nounwind {
+; CHECK-LABEL: @test
+; CHECK: add nuw nsw i32
+
+  ; Drop two highes bits to guarantee that %a + %b doesn't overflow
+  %a_and = and i32 %a, 1073741823
+  %b_and = and i32 %b, 1073741823
+
+  %a_and_fp = sitofp i32 %a_and to double
+  %b_and_fp = sitofp i32 %b_and to double
+
+  %res = fadd double %a_and_fp, %b_and_fp
+  ret double %res
+}
+
+define float @test_2_neg(i32 %a, i32 %b) nounwind {
+; CHECK-LABEL: @test
+; CHECK: fadd float
+
+  ; Drop two highes bits to guarantee that %a + %b doesn't overflow
+  %a_and = and i32 %a, 1073741823
+  %b_and = and i32 %b, 1073741823
+
+  %a_and_fp = sitofp i32 %a_and to float
+  %b_and_fp = sitofp i32 %b_and to float
+
+  %res = fadd float %a_and_fp, %b_and_fp
+  ret float %res
+}
 No newline at end of file