1 //===- LoopIdiomRecognize.cpp - Loop idiom recognition --------------------===// 2 // 3 // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. 4 // See https://llvm.org/LICENSE.txt for license information. 5 // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception 6 // 7 //===----------------------------------------------------------------------===// 8 // 9 // This pass implements an idiom recognizer that transforms simple loops into a 10 // non-loop form. In cases that this kicks in, it can be a significant 11 // performance win. 12 // 13 // If compiling for code size we avoid idiom recognition if the resulting 14 // code could be larger than the code for the original loop. One way this could 15 // happen is if the loop is not removable after idiom recognition due to the 16 // presence of non-idiom instructions. The initial implementation of the 17 // heuristics applies to idioms in multi-block loops. 18 // 19 //===----------------------------------------------------------------------===// 20 // 21 // TODO List: 22 // 23 // Future loop memory idioms to recognize: 24 // memcmp, strlen, etc. 25 // Future floating point idioms to recognize in -ffast-math mode: 26 // fpowi 27 // Future integer operation idioms to recognize: 28 // ctpop 29 // 30 // Beware that isel's default lowering for ctpop is highly inefficient for 31 // i64 and larger types when i64 is legal and the value has few bits set. It 32 // would be good to enhance isel to emit a loop for ctpop in this case. 33 // 34 // This could recognize common matrix multiplies and dot product idioms and 35 // replace them with calls to BLAS (if linked in??). 36 // 37 //===----------------------------------------------------------------------===// 38 39 #include "llvm/Transforms/Scalar/LoopIdiomRecognize.h" 40 #include "llvm/ADT/APInt.h" 41 #include "llvm/ADT/ArrayRef.h" 42 #include "llvm/ADT/DenseMap.h" 43 #include "llvm/ADT/MapVector.h" 44 #include "llvm/ADT/SetVector.h" 45 #include "llvm/ADT/SmallPtrSet.h" 46 #include "llvm/ADT/SmallVector.h" 47 #include "llvm/ADT/Statistic.h" 48 #include "llvm/ADT/StringRef.h" 49 #include "llvm/Analysis/AliasAnalysis.h" 50 #include "llvm/Analysis/CmpInstAnalysis.h" 51 #include "llvm/Analysis/LoopAccessAnalysis.h" 52 #include "llvm/Analysis/LoopInfo.h" 53 #include "llvm/Analysis/LoopPass.h" 54 #include "llvm/Analysis/MemoryLocation.h" 55 #include "llvm/Analysis/MemorySSA.h" 56 #include "llvm/Analysis/MemorySSAUpdater.h" 57 #include "llvm/Analysis/MustExecute.h" 58 #include "llvm/Analysis/OptimizationRemarkEmitter.h" 59 #include "llvm/Analysis/ScalarEvolution.h" 60 #include "llvm/Analysis/ScalarEvolutionExpressions.h" 61 #include "llvm/Analysis/TargetLibraryInfo.h" 62 #include "llvm/Analysis/TargetTransformInfo.h" 63 #include "llvm/Analysis/ValueTracking.h" 64 #include "llvm/IR/Attributes.h" 65 #include "llvm/IR/BasicBlock.h" 66 #include "llvm/IR/Constant.h" 67 #include "llvm/IR/Constants.h" 68 #include "llvm/IR/DataLayout.h" 69 #include "llvm/IR/DebugLoc.h" 70 #include "llvm/IR/DerivedTypes.h" 71 #include "llvm/IR/Dominators.h" 72 #include "llvm/IR/GlobalValue.h" 73 #include "llvm/IR/GlobalVariable.h" 74 #include "llvm/IR/IRBuilder.h" 75 #include "llvm/IR/InstrTypes.h" 76 #include "llvm/IR/Instruction.h" 77 #include "llvm/IR/Instructions.h" 78 #include "llvm/IR/IntrinsicInst.h" 79 #include "llvm/IR/Intrinsics.h" 80 #include "llvm/IR/LLVMContext.h" 81 #include "llvm/IR/Module.h" 82 #include "llvm/IR/PassManager.h" 83 #include "llvm/IR/PatternMatch.h" 84 #include "llvm/IR/Type.h" 85 #include "llvm/IR/User.h" 86 #include "llvm/IR/Value.h" 87 #include "llvm/IR/ValueHandle.h" 88 #include "llvm/InitializePasses.h" 89 #include "llvm/Pass.h" 90 #include "llvm/Support/Casting.h" 91 #include "llvm/Support/CommandLine.h" 92 #include "llvm/Support/Debug.h" 93 #include "llvm/Support/InstructionCost.h" 94 #include "llvm/Support/raw_ostream.h" 95 #include "llvm/Transforms/Scalar.h" 96 #include "llvm/Transforms/Utils/BuildLibCalls.h" 97 #include "llvm/Transforms/Utils/Local.h" 98 #include "llvm/Transforms/Utils/LoopUtils.h" 99 #include "llvm/Transforms/Utils/ScalarEvolutionExpander.h" 100 #include <algorithm> 101 #include <cassert> 102 #include <cstdint> 103 #include <utility> 104 #include <vector> 105 106 using namespace llvm; 107 108 #define DEBUG_TYPE "loop-idiom" 109 110 STATISTIC(NumMemSet, "Number of memset's formed from loop stores"); 111 STATISTIC(NumMemCpy, "Number of memcpy's formed from loop load+stores"); 112 STATISTIC(NumMemMove, "Number of memmove's formed from loop load+stores"); 113 STATISTIC( 114 NumShiftUntilBitTest, 115 "Number of uncountable loops recognized as 'shift until bitttest' idiom"); 116 STATISTIC(NumShiftUntilZero, 117 "Number of uncountable loops recognized as 'shift until zero' idiom"); 118 119 bool DisableLIRP::All; 120 static cl::opt<bool, true> 121 DisableLIRPAll("disable-" DEBUG_TYPE "-all", 122 cl::desc("Options to disable Loop Idiom Recognize Pass."), 123 cl::location(DisableLIRP::All), cl::init(false), 124 cl::ReallyHidden); 125 126 bool DisableLIRP::Memset; 127 static cl::opt<bool, true> 128 DisableLIRPMemset("disable-" DEBUG_TYPE "-memset", 129 cl::desc("Proceed with loop idiom recognize pass, but do " 130 "not convert loop(s) to memset."), 131 cl::location(DisableLIRP::Memset), cl::init(false), 132 cl::ReallyHidden); 133 134 bool DisableLIRP::Memcpy; 135 static cl::opt<bool, true> 136 DisableLIRPMemcpy("disable-" DEBUG_TYPE "-memcpy", 137 cl::desc("Proceed with loop idiom recognize pass, but do " 138 "not convert loop(s) to memcpy."), 139 cl::location(DisableLIRP::Memcpy), cl::init(false), 140 cl::ReallyHidden); 141 142 static cl::opt<bool> UseLIRCodeSizeHeurs( 143 "use-lir-code-size-heurs", 144 cl::desc("Use loop idiom recognition code size heuristics when compiling" 145 "with -Os/-Oz"), 146 cl::init(true), cl::Hidden); 147 148 namespace { 149 150 class LoopIdiomRecognize { 151 Loop *CurLoop = nullptr; 152 AliasAnalysis *AA; 153 DominatorTree *DT; 154 LoopInfo *LI; 155 ScalarEvolution *SE; 156 TargetLibraryInfo *TLI; 157 const TargetTransformInfo *TTI; 158 const DataLayout *DL; 159 OptimizationRemarkEmitter &ORE; 160 bool ApplyCodeSizeHeuristics; 161 std::unique_ptr<MemorySSAUpdater> MSSAU; 162 163 public: 164 explicit LoopIdiomRecognize(AliasAnalysis *AA, DominatorTree *DT, 165 LoopInfo *LI, ScalarEvolution *SE, 166 TargetLibraryInfo *TLI, 167 const TargetTransformInfo *TTI, MemorySSA *MSSA, 168 const DataLayout *DL, 169 OptimizationRemarkEmitter &ORE) 170 : AA(AA), DT(DT), LI(LI), SE(SE), TLI(TLI), TTI(TTI), DL(DL), ORE(ORE) { 171 if (MSSA) 172 MSSAU = std::make_unique<MemorySSAUpdater>(MSSA); 173 } 174 175 bool runOnLoop(Loop *L); 176 177 private: 178 using StoreList = SmallVector<StoreInst *, 8>; 179 using StoreListMap = MapVector<Value *, StoreList>; 180 181 StoreListMap StoreRefsForMemset; 182 StoreListMap StoreRefsForMemsetPattern; 183 StoreList StoreRefsForMemcpy; 184 bool HasMemset; 185 bool HasMemsetPattern; 186 bool HasMemcpy; 187 188 /// Return code for isLegalStore() 189 enum LegalStoreKind { 190 None = 0, 191 Memset, 192 MemsetPattern, 193 Memcpy, 194 UnorderedAtomicMemcpy, 195 DontUse // Dummy retval never to be used. Allows catching errors in retval 196 // handling. 197 }; 198 199 /// \name Countable Loop Idiom Handling 200 /// @{ 201 202 bool runOnCountableLoop(); 203 bool runOnLoopBlock(BasicBlock *BB, const SCEV *BECount, 204 SmallVectorImpl<BasicBlock *> &ExitBlocks); 205 206 void collectStores(BasicBlock *BB); 207 LegalStoreKind isLegalStore(StoreInst *SI); 208 enum class ForMemset { No, Yes }; 209 bool processLoopStores(SmallVectorImpl<StoreInst *> &SL, const SCEV *BECount, 210 ForMemset For); 211 212 template <typename MemInst> 213 bool processLoopMemIntrinsic( 214 BasicBlock *BB, 215 bool (LoopIdiomRecognize::*Processor)(MemInst *, const SCEV *), 216 const SCEV *BECount); 217 bool processLoopMemCpy(MemCpyInst *MCI, const SCEV *BECount); 218 bool processLoopMemSet(MemSetInst *MSI, const SCEV *BECount); 219 220 bool processLoopStridedStore(Value *DestPtr, const SCEV *StoreSizeSCEV, 221 MaybeAlign StoreAlignment, Value *StoredVal, 222 Instruction *TheStore, 223 SmallPtrSetImpl<Instruction *> &Stores, 224 const SCEVAddRecExpr *Ev, const SCEV *BECount, 225 bool IsNegStride, bool IsLoopMemset = false); 226 bool processLoopStoreOfLoopLoad(StoreInst *SI, const SCEV *BECount); 227 bool processLoopStoreOfLoopLoad(Value *DestPtr, Value *SourcePtr, 228 const SCEV *StoreSize, MaybeAlign StoreAlign, 229 MaybeAlign LoadAlign, Instruction *TheStore, 230 Instruction *TheLoad, 231 const SCEVAddRecExpr *StoreEv, 232 const SCEVAddRecExpr *LoadEv, 233 const SCEV *BECount); 234 bool avoidLIRForMultiBlockLoop(bool IsMemset = false, 235 bool IsLoopMemset = false); 236 237 /// @} 238 /// \name Noncountable Loop Idiom Handling 239 /// @{ 240 241 bool runOnNoncountableLoop(); 242 243 bool recognizePopcount(); 244 void transformLoopToPopcount(BasicBlock *PreCondBB, Instruction *CntInst, 245 PHINode *CntPhi, Value *Var); 246 bool recognizeAndInsertFFS(); /// Find First Set: ctlz or cttz 247 void transformLoopToCountable(Intrinsic::ID IntrinID, BasicBlock *PreCondBB, 248 Instruction *CntInst, PHINode *CntPhi, 249 Value *Var, Instruction *DefX, 250 const DebugLoc &DL, bool ZeroCheck, 251 bool IsCntPhiUsedOutsideLoop); 252 253 bool recognizeShiftUntilBitTest(); 254 bool recognizeShiftUntilZero(); 255 256 /// @} 257 }; 258 259 class LoopIdiomRecognizeLegacyPass : public LoopPass { 260 public: 261 static char ID; 262 263 explicit LoopIdiomRecognizeLegacyPass() : LoopPass(ID) { 264 initializeLoopIdiomRecognizeLegacyPassPass( 265 *PassRegistry::getPassRegistry()); 266 } 267 268 bool runOnLoop(Loop *L, LPPassManager &LPM) override { 269 if (DisableLIRP::All) 270 return false; 271 272 if (skipLoop(L)) 273 return false; 274 275 AliasAnalysis *AA = &getAnalysis<AAResultsWrapperPass>().getAAResults(); 276 DominatorTree *DT = &getAnalysis<DominatorTreeWrapperPass>().getDomTree(); 277 LoopInfo *LI = &getAnalysis<LoopInfoWrapperPass>().getLoopInfo(); 278 ScalarEvolution *SE = &getAnalysis<ScalarEvolutionWrapperPass>().getSE(); 279 TargetLibraryInfo *TLI = 280 &getAnalysis<TargetLibraryInfoWrapperPass>().getTLI( 281 *L->getHeader()->getParent()); 282 const TargetTransformInfo *TTI = 283 &getAnalysis<TargetTransformInfoWrapperPass>().getTTI( 284 *L->getHeader()->getParent()); 285 const DataLayout *DL = &L->getHeader()->getModule()->getDataLayout(); 286 auto *MSSAAnalysis = getAnalysisIfAvailable<MemorySSAWrapperPass>(); 287 MemorySSA *MSSA = nullptr; 288 if (MSSAAnalysis) 289 MSSA = &MSSAAnalysis->getMSSA(); 290 291 // For the old PM, we can't use OptimizationRemarkEmitter as an analysis 292 // pass. Function analyses need to be preserved across loop transformations 293 // but ORE cannot be preserved (see comment before the pass definition). 294 OptimizationRemarkEmitter ORE(L->getHeader()->getParent()); 295 296 LoopIdiomRecognize LIR(AA, DT, LI, SE, TLI, TTI, MSSA, DL, ORE); 297 return LIR.runOnLoop(L); 298 } 299 300 /// This transformation requires natural loop information & requires that 301 /// loop preheaders be inserted into the CFG. 302 void getAnalysisUsage(AnalysisUsage &AU) const override { 303 AU.addRequired<TargetLibraryInfoWrapperPass>(); 304 AU.addRequired<TargetTransformInfoWrapperPass>(); 305 AU.addPreserved<MemorySSAWrapperPass>(); 306 getLoopAnalysisUsage(AU); 307 } 308 }; 309 310 } // end anonymous namespace 311 312 char LoopIdiomRecognizeLegacyPass::ID = 0; 313 314 PreservedAnalyses LoopIdiomRecognizePass::run(Loop &L, LoopAnalysisManager &AM, 315 LoopStandardAnalysisResults &AR, 316 LPMUpdater &) { 317 if (DisableLIRP::All) 318 return PreservedAnalyses::all(); 319 320 const auto *DL = &L.getHeader()->getModule()->getDataLayout(); 321 322 // For the new PM, we also can't use OptimizationRemarkEmitter as an analysis 323 // pass. Function analyses need to be preserved across loop transformations 324 // but ORE cannot be preserved (see comment before the pass definition). 325 OptimizationRemarkEmitter ORE(L.getHeader()->getParent()); 326 327 LoopIdiomRecognize LIR(&AR.AA, &AR.DT, &AR.LI, &AR.SE, &AR.TLI, &AR.TTI, 328 AR.MSSA, DL, ORE); 329 if (!LIR.runOnLoop(&L)) 330 return PreservedAnalyses::all(); 331 332 auto PA = getLoopPassPreservedAnalyses(); 333 if (AR.MSSA) 334 PA.preserve<MemorySSAAnalysis>(); 335 return PA; 336 } 337 338 INITIALIZE_PASS_BEGIN(LoopIdiomRecognizeLegacyPass, "loop-idiom", 339 "Recognize loop idioms", false, false) 340 INITIALIZE_PASS_DEPENDENCY(LoopPass) 341 INITIALIZE_PASS_DEPENDENCY(TargetLibraryInfoWrapperPass) 342 INITIALIZE_PASS_DEPENDENCY(TargetTransformInfoWrapperPass) 343 INITIALIZE_PASS_END(LoopIdiomRecognizeLegacyPass, "loop-idiom", 344 "Recognize loop idioms", false, false) 345 346 Pass *llvm::createLoopIdiomPass() { return new LoopIdiomRecognizeLegacyPass(); } 347 348 static void deleteDeadInstruction(Instruction *I) { 349 I->replaceAllUsesWith(UndefValue::get(I->getType())); 350 I->eraseFromParent(); 351 } 352 353 //===----------------------------------------------------------------------===// 354 // 355 // Implementation of LoopIdiomRecognize 356 // 357 //===----------------------------------------------------------------------===// 358 359 bool LoopIdiomRecognize::runOnLoop(Loop *L) { 360 CurLoop = L; 361 // If the loop could not be converted to canonical form, it must have an 362 // indirectbr in it, just give up. 363 if (!L->getLoopPreheader()) 364 return false; 365 366 // Disable loop idiom recognition if the function's name is a common idiom. 367 StringRef Name = L->getHeader()->getParent()->getName(); 368 if (Name == "memset" || Name == "memcpy") 369 return false; 370 371 // Determine if code size heuristics need to be applied. 372 ApplyCodeSizeHeuristics = 373 L->getHeader()->getParent()->hasOptSize() && UseLIRCodeSizeHeurs; 374 375 HasMemset = TLI->has(LibFunc_memset); 376 HasMemsetPattern = TLI->has(LibFunc_memset_pattern16); 377 HasMemcpy = TLI->has(LibFunc_memcpy); 378 379 if (HasMemset || HasMemsetPattern || HasMemcpy) 380 if (SE->hasLoopInvariantBackedgeTakenCount(L)) 381 return runOnCountableLoop(); 382 383 return runOnNoncountableLoop(); 384 } 385 386 bool LoopIdiomRecognize::runOnCountableLoop() { 387 const SCEV *BECount = SE->getBackedgeTakenCount(CurLoop); 388 assert(!isa<SCEVCouldNotCompute>(BECount) && 389 "runOnCountableLoop() called on a loop without a predictable" 390 "backedge-taken count"); 391 392 // If this loop executes exactly one time, then it should be peeled, not 393 // optimized by this pass. 394 if (const SCEVConstant *BECst = dyn_cast<SCEVConstant>(BECount)) 395 if (BECst->getAPInt() == 0) 396 return false; 397 398 SmallVector<BasicBlock *, 8> ExitBlocks; 399 CurLoop->getUniqueExitBlocks(ExitBlocks); 400 401 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Scanning: F[" 402 << CurLoop->getHeader()->getParent()->getName() 403 << "] Countable Loop %" << CurLoop->getHeader()->getName() 404 << "\n"); 405 406 // The following transforms hoist stores/memsets into the loop pre-header. 407 // Give up if the loop has instructions that may throw. 408 SimpleLoopSafetyInfo SafetyInfo; 409 SafetyInfo.computeLoopSafetyInfo(CurLoop); 410 if (SafetyInfo.anyBlockMayThrow()) 411 return false; 412 413 bool MadeChange = false; 414 415 // Scan all the blocks in the loop that are not in subloops. 416 for (auto *BB : CurLoop->getBlocks()) { 417 // Ignore blocks in subloops. 418 if (LI->getLoopFor(BB) != CurLoop) 419 continue; 420 421 MadeChange |= runOnLoopBlock(BB, BECount, ExitBlocks); 422 } 423 return MadeChange; 424 } 425 426 static APInt getStoreStride(const SCEVAddRecExpr *StoreEv) { 427 const SCEVConstant *ConstStride = cast<SCEVConstant>(StoreEv->getOperand(1)); 428 return ConstStride->getAPInt(); 429 } 430 431 /// getMemSetPatternValue - If a strided store of the specified value is safe to 432 /// turn into a memset_pattern16, return a ConstantArray of 16 bytes that should 433 /// be passed in. Otherwise, return null. 434 /// 435 /// Note that we don't ever attempt to use memset_pattern8 or 4, because these 436 /// just replicate their input array and then pass on to memset_pattern16. 437 static Constant *getMemSetPatternValue(Value *V, const DataLayout *DL) { 438 // FIXME: This could check for UndefValue because it can be merged into any 439 // other valid pattern. 440 441 // If the value isn't a constant, we can't promote it to being in a constant 442 // array. We could theoretically do a store to an alloca or something, but 443 // that doesn't seem worthwhile. 444 Constant *C = dyn_cast<Constant>(V); 445 if (!C) 446 return nullptr; 447 448 // Only handle simple values that are a power of two bytes in size. 449 uint64_t Size = DL->getTypeSizeInBits(V->getType()); 450 if (Size == 0 || (Size & 7) || (Size & (Size - 1))) 451 return nullptr; 452 453 // Don't care enough about darwin/ppc to implement this. 454 if (DL->isBigEndian()) 455 return nullptr; 456 457 // Convert to size in bytes. 458 Size /= 8; 459 460 // TODO: If CI is larger than 16-bytes, we can try slicing it in half to see 461 // if the top and bottom are the same (e.g. for vectors and large integers). 462 if (Size > 16) 463 return nullptr; 464 465 // If the constant is exactly 16 bytes, just use it. 466 if (Size == 16) 467 return C; 468 469 // Otherwise, we'll use an array of the constants. 470 unsigned ArraySize = 16 / Size; 471 ArrayType *AT = ArrayType::get(V->getType(), ArraySize); 472 return ConstantArray::get(AT, std::vector<Constant *>(ArraySize, C)); 473 } 474 475 LoopIdiomRecognize::LegalStoreKind 476 LoopIdiomRecognize::isLegalStore(StoreInst *SI) { 477 // Don't touch volatile stores. 478 if (SI->isVolatile()) 479 return LegalStoreKind::None; 480 // We only want simple or unordered-atomic stores. 481 if (!SI->isUnordered()) 482 return LegalStoreKind::None; 483 484 // Avoid merging nontemporal stores. 485 if (SI->getMetadata(LLVMContext::MD_nontemporal)) 486 return LegalStoreKind::None; 487 488 Value *StoredVal = SI->getValueOperand(); 489 Value *StorePtr = SI->getPointerOperand(); 490 491 // Don't convert stores of non-integral pointer types to memsets (which stores 492 // integers). 493 if (DL->isNonIntegralPointerType(StoredVal->getType()->getScalarType())) 494 return LegalStoreKind::None; 495 496 // Reject stores that are so large that they overflow an unsigned. 497 // When storing out scalable vectors we bail out for now, since the code 498 // below currently only works for constant strides. 499 TypeSize SizeInBits = DL->getTypeSizeInBits(StoredVal->getType()); 500 if (SizeInBits.isScalable() || (SizeInBits.getFixedSize() & 7) || 501 (SizeInBits.getFixedSize() >> 32) != 0) 502 return LegalStoreKind::None; 503 504 // See if the pointer expression is an AddRec like {base,+,1} on the current 505 // loop, which indicates a strided store. If we have something else, it's a 506 // random store we can't handle. 507 const SCEVAddRecExpr *StoreEv = 508 dyn_cast<SCEVAddRecExpr>(SE->getSCEV(StorePtr)); 509 if (!StoreEv || StoreEv->getLoop() != CurLoop || !StoreEv->isAffine()) 510 return LegalStoreKind::None; 511 512 // Check to see if we have a constant stride. 513 if (!isa<SCEVConstant>(StoreEv->getOperand(1))) 514 return LegalStoreKind::None; 515 516 // See if the store can be turned into a memset. 517 518 // If the stored value is a byte-wise value (like i32 -1), then it may be 519 // turned into a memset of i8 -1, assuming that all the consecutive bytes 520 // are stored. A store of i32 0x01020304 can never be turned into a memset, 521 // but it can be turned into memset_pattern if the target supports it. 522 Value *SplatValue = isBytewiseValue(StoredVal, *DL); 523 524 // Note: memset and memset_pattern on unordered-atomic is yet not supported 525 bool UnorderedAtomic = SI->isUnordered() && !SI->isSimple(); 526 527 // If we're allowed to form a memset, and the stored value would be 528 // acceptable for memset, use it. 529 if (!UnorderedAtomic && HasMemset && SplatValue && !DisableLIRP::Memset && 530 // Verify that the stored value is loop invariant. If not, we can't 531 // promote the memset. 532 CurLoop->isLoopInvariant(SplatValue)) { 533 // It looks like we can use SplatValue. 534 return LegalStoreKind::Memset; 535 } 536 if (!UnorderedAtomic && HasMemsetPattern && !DisableLIRP::Memset && 537 // Don't create memset_pattern16s with address spaces. 538 StorePtr->getType()->getPointerAddressSpace() == 0 && 539 getMemSetPatternValue(StoredVal, DL)) { 540 // It looks like we can use PatternValue! 541 return LegalStoreKind::MemsetPattern; 542 } 543 544 // Otherwise, see if the store can be turned into a memcpy. 545 if (HasMemcpy && !DisableLIRP::Memcpy) { 546 // Check to see if the stride matches the size of the store. If so, then we 547 // know that every byte is touched in the loop. 548 APInt Stride = getStoreStride(StoreEv); 549 unsigned StoreSize = DL->getTypeStoreSize(SI->getValueOperand()->getType()); 550 if (StoreSize != Stride && StoreSize != -Stride) 551 return LegalStoreKind::None; 552 553 // The store must be feeding a non-volatile load. 554 LoadInst *LI = dyn_cast<LoadInst>(SI->getValueOperand()); 555 556 // Only allow non-volatile loads 557 if (!LI || LI->isVolatile()) 558 return LegalStoreKind::None; 559 // Only allow simple or unordered-atomic loads 560 if (!LI->isUnordered()) 561 return LegalStoreKind::None; 562 563 // See if the pointer expression is an AddRec like {base,+,1} on the current 564 // loop, which indicates a strided load. If we have something else, it's a 565 // random load we can't handle. 566 const SCEVAddRecExpr *LoadEv = 567 dyn_cast<SCEVAddRecExpr>(SE->getSCEV(LI->getPointerOperand())); 568 if (!LoadEv || LoadEv->getLoop() != CurLoop || !LoadEv->isAffine()) 569 return LegalStoreKind::None; 570 571 // The store and load must share the same stride. 572 if (StoreEv->getOperand(1) != LoadEv->getOperand(1)) 573 return LegalStoreKind::None; 574 575 // Success. This store can be converted into a memcpy. 576 UnorderedAtomic = UnorderedAtomic || LI->isAtomic(); 577 return UnorderedAtomic ? LegalStoreKind::UnorderedAtomicMemcpy 578 : LegalStoreKind::Memcpy; 579 } 580 // This store can't be transformed into a memset/memcpy. 581 return LegalStoreKind::None; 582 } 583 584 void LoopIdiomRecognize::collectStores(BasicBlock *BB) { 585 StoreRefsForMemset.clear(); 586 StoreRefsForMemsetPattern.clear(); 587 StoreRefsForMemcpy.clear(); 588 for (Instruction &I : *BB) { 589 StoreInst *SI = dyn_cast<StoreInst>(&I); 590 if (!SI) 591 continue; 592 593 // Make sure this is a strided store with a constant stride. 594 switch (isLegalStore(SI)) { 595 case LegalStoreKind::None: 596 // Nothing to do 597 break; 598 case LegalStoreKind::Memset: { 599 // Find the base pointer. 600 Value *Ptr = getUnderlyingObject(SI->getPointerOperand()); 601 StoreRefsForMemset[Ptr].push_back(SI); 602 } break; 603 case LegalStoreKind::MemsetPattern: { 604 // Find the base pointer. 605 Value *Ptr = getUnderlyingObject(SI->getPointerOperand()); 606 StoreRefsForMemsetPattern[Ptr].push_back(SI); 607 } break; 608 case LegalStoreKind::Memcpy: 609 case LegalStoreKind::UnorderedAtomicMemcpy: 610 StoreRefsForMemcpy.push_back(SI); 611 break; 612 default: 613 assert(false && "unhandled return value"); 614 break; 615 } 616 } 617 } 618 619 /// runOnLoopBlock - Process the specified block, which lives in a counted loop 620 /// with the specified backedge count. This block is known to be in the current 621 /// loop and not in any subloops. 622 bool LoopIdiomRecognize::runOnLoopBlock( 623 BasicBlock *BB, const SCEV *BECount, 624 SmallVectorImpl<BasicBlock *> &ExitBlocks) { 625 // We can only promote stores in this block if they are unconditionally 626 // executed in the loop. For a block to be unconditionally executed, it has 627 // to dominate all the exit blocks of the loop. Verify this now. 628 for (BasicBlock *ExitBlock : ExitBlocks) 629 if (!DT->dominates(BB, ExitBlock)) 630 return false; 631 632 bool MadeChange = false; 633 // Look for store instructions, which may be optimized to memset/memcpy. 634 collectStores(BB); 635 636 // Look for a single store or sets of stores with a common base, which can be 637 // optimized into a memset (memset_pattern). The latter most commonly happens 638 // with structs and handunrolled loops. 639 for (auto &SL : StoreRefsForMemset) 640 MadeChange |= processLoopStores(SL.second, BECount, ForMemset::Yes); 641 642 for (auto &SL : StoreRefsForMemsetPattern) 643 MadeChange |= processLoopStores(SL.second, BECount, ForMemset::No); 644 645 // Optimize the store into a memcpy, if it feeds an similarly strided load. 646 for (auto &SI : StoreRefsForMemcpy) 647 MadeChange |= processLoopStoreOfLoopLoad(SI, BECount); 648 649 MadeChange |= processLoopMemIntrinsic<MemCpyInst>( 650 BB, &LoopIdiomRecognize::processLoopMemCpy, BECount); 651 MadeChange |= processLoopMemIntrinsic<MemSetInst>( 652 BB, &LoopIdiomRecognize::processLoopMemSet, BECount); 653 654 return MadeChange; 655 } 656 657 /// See if this store(s) can be promoted to a memset. 658 bool LoopIdiomRecognize::processLoopStores(SmallVectorImpl<StoreInst *> &SL, 659 const SCEV *BECount, ForMemset For) { 660 // Try to find consecutive stores that can be transformed into memsets. 661 SetVector<StoreInst *> Heads, Tails; 662 SmallDenseMap<StoreInst *, StoreInst *> ConsecutiveChain; 663 664 // Do a quadratic search on all of the given stores and find 665 // all of the pairs of stores that follow each other. 666 SmallVector<unsigned, 16> IndexQueue; 667 for (unsigned i = 0, e = SL.size(); i < e; ++i) { 668 assert(SL[i]->isSimple() && "Expected only non-volatile stores."); 669 670 Value *FirstStoredVal = SL[i]->getValueOperand(); 671 Value *FirstStorePtr = SL[i]->getPointerOperand(); 672 const SCEVAddRecExpr *FirstStoreEv = 673 cast<SCEVAddRecExpr>(SE->getSCEV(FirstStorePtr)); 674 APInt FirstStride = getStoreStride(FirstStoreEv); 675 unsigned FirstStoreSize = DL->getTypeStoreSize(SL[i]->getValueOperand()->getType()); 676 677 // See if we can optimize just this store in isolation. 678 if (FirstStride == FirstStoreSize || -FirstStride == FirstStoreSize) { 679 Heads.insert(SL[i]); 680 continue; 681 } 682 683 Value *FirstSplatValue = nullptr; 684 Constant *FirstPatternValue = nullptr; 685 686 if (For == ForMemset::Yes) 687 FirstSplatValue = isBytewiseValue(FirstStoredVal, *DL); 688 else 689 FirstPatternValue = getMemSetPatternValue(FirstStoredVal, DL); 690 691 assert((FirstSplatValue || FirstPatternValue) && 692 "Expected either splat value or pattern value."); 693 694 IndexQueue.clear(); 695 // If a store has multiple consecutive store candidates, search Stores 696 // array according to the sequence: from i+1 to e, then from i-1 to 0. 697 // This is because usually pairing with immediate succeeding or preceding 698 // candidate create the best chance to find memset opportunity. 699 unsigned j = 0; 700 for (j = i + 1; j < e; ++j) 701 IndexQueue.push_back(j); 702 for (j = i; j > 0; --j) 703 IndexQueue.push_back(j - 1); 704 705 for (auto &k : IndexQueue) { 706 assert(SL[k]->isSimple() && "Expected only non-volatile stores."); 707 Value *SecondStorePtr = SL[k]->getPointerOperand(); 708 const SCEVAddRecExpr *SecondStoreEv = 709 cast<SCEVAddRecExpr>(SE->getSCEV(SecondStorePtr)); 710 APInt SecondStride = getStoreStride(SecondStoreEv); 711 712 if (FirstStride != SecondStride) 713 continue; 714 715 Value *SecondStoredVal = SL[k]->getValueOperand(); 716 Value *SecondSplatValue = nullptr; 717 Constant *SecondPatternValue = nullptr; 718 719 if (For == ForMemset::Yes) 720 SecondSplatValue = isBytewiseValue(SecondStoredVal, *DL); 721 else 722 SecondPatternValue = getMemSetPatternValue(SecondStoredVal, DL); 723 724 assert((SecondSplatValue || SecondPatternValue) && 725 "Expected either splat value or pattern value."); 726 727 if (isConsecutiveAccess(SL[i], SL[k], *DL, *SE, false)) { 728 if (For == ForMemset::Yes) { 729 if (isa<UndefValue>(FirstSplatValue)) 730 FirstSplatValue = SecondSplatValue; 731 if (FirstSplatValue != SecondSplatValue) 732 continue; 733 } else { 734 if (isa<UndefValue>(FirstPatternValue)) 735 FirstPatternValue = SecondPatternValue; 736 if (FirstPatternValue != SecondPatternValue) 737 continue; 738 } 739 Tails.insert(SL[k]); 740 Heads.insert(SL[i]); 741 ConsecutiveChain[SL[i]] = SL[k]; 742 break; 743 } 744 } 745 } 746 747 // We may run into multiple chains that merge into a single chain. We mark the 748 // stores that we transformed so that we don't visit the same store twice. 749 SmallPtrSet<Value *, 16> TransformedStores; 750 bool Changed = false; 751 752 // For stores that start but don't end a link in the chain: 753 for (StoreInst *I : Heads) { 754 if (Tails.count(I)) 755 continue; 756 757 // We found a store instr that starts a chain. Now follow the chain and try 758 // to transform it. 759 SmallPtrSet<Instruction *, 8> AdjacentStores; 760 StoreInst *HeadStore = I; 761 unsigned StoreSize = 0; 762 763 // Collect the chain into a list. 764 while (Tails.count(I) || Heads.count(I)) { 765 if (TransformedStores.count(I)) 766 break; 767 AdjacentStores.insert(I); 768 769 StoreSize += DL->getTypeStoreSize(I->getValueOperand()->getType()); 770 // Move to the next value in the chain. 771 I = ConsecutiveChain[I]; 772 } 773 774 Value *StoredVal = HeadStore->getValueOperand(); 775 Value *StorePtr = HeadStore->getPointerOperand(); 776 const SCEVAddRecExpr *StoreEv = cast<SCEVAddRecExpr>(SE->getSCEV(StorePtr)); 777 APInt Stride = getStoreStride(StoreEv); 778 779 // Check to see if the stride matches the size of the stores. If so, then 780 // we know that every byte is touched in the loop. 781 if (StoreSize != Stride && StoreSize != -Stride) 782 continue; 783 784 bool IsNegStride = StoreSize == -Stride; 785 786 Type *IntIdxTy = DL->getIndexType(StorePtr->getType()); 787 const SCEV *StoreSizeSCEV = SE->getConstant(IntIdxTy, StoreSize); 788 if (processLoopStridedStore(StorePtr, StoreSizeSCEV, 789 MaybeAlign(HeadStore->getAlign()), StoredVal, 790 HeadStore, AdjacentStores, StoreEv, BECount, 791 IsNegStride)) { 792 TransformedStores.insert(AdjacentStores.begin(), AdjacentStores.end()); 793 Changed = true; 794 } 795 } 796 797 return Changed; 798 } 799 800 /// processLoopMemIntrinsic - Template function for calling different processor 801 /// functions based on mem instrinsic type. 802 template <typename MemInst> 803 bool LoopIdiomRecognize::processLoopMemIntrinsic( 804 BasicBlock *BB, 805 bool (LoopIdiomRecognize::*Processor)(MemInst *, const SCEV *), 806 const SCEV *BECount) { 807 bool MadeChange = false; 808 for (BasicBlock::iterator I = BB->begin(), E = BB->end(); I != E;) { 809 Instruction *Inst = &*I++; 810 // Look for memory instructions, which may be optimized to a larger one. 811 if (MemInst *MI = dyn_cast<MemInst>(Inst)) { 812 WeakTrackingVH InstPtr(&*I); 813 if (!(this->*Processor)(MI, BECount)) 814 continue; 815 MadeChange = true; 816 817 // If processing the instruction invalidated our iterator, start over from 818 // the top of the block. 819 if (!InstPtr) 820 I = BB->begin(); 821 } 822 } 823 return MadeChange; 824 } 825 826 /// processLoopMemCpy - See if this memcpy can be promoted to a large memcpy 827 bool LoopIdiomRecognize::processLoopMemCpy(MemCpyInst *MCI, 828 const SCEV *BECount) { 829 // We can only handle non-volatile memcpys with a constant size. 830 if (MCI->isVolatile() || !isa<ConstantInt>(MCI->getLength())) 831 return false; 832 833 // If we're not allowed to hack on memcpy, we fail. 834 if ((!HasMemcpy && !isa<MemCpyInlineInst>(MCI)) || DisableLIRP::Memcpy) 835 return false; 836 837 Value *Dest = MCI->getDest(); 838 Value *Source = MCI->getSource(); 839 if (!Dest || !Source) 840 return false; 841 842 // See if the load and store pointer expressions are AddRec like {base,+,1} on 843 // the current loop, which indicates a strided load and store. If we have 844 // something else, it's a random load or store we can't handle. 845 const SCEVAddRecExpr *StoreEv = dyn_cast<SCEVAddRecExpr>(SE->getSCEV(Dest)); 846 if (!StoreEv || StoreEv->getLoop() != CurLoop || !StoreEv->isAffine()) 847 return false; 848 const SCEVAddRecExpr *LoadEv = dyn_cast<SCEVAddRecExpr>(SE->getSCEV(Source)); 849 if (!LoadEv || LoadEv->getLoop() != CurLoop || !LoadEv->isAffine()) 850 return false; 851 852 // Reject memcpys that are so large that they overflow an unsigned. 853 uint64_t SizeInBytes = cast<ConstantInt>(MCI->getLength())->getZExtValue(); 854 if ((SizeInBytes >> 32) != 0) 855 return false; 856 857 // Check if the stride matches the size of the memcpy. If so, then we know 858 // that every byte is touched in the loop. 859 const SCEVConstant *ConstStoreStride = 860 dyn_cast<SCEVConstant>(StoreEv->getOperand(1)); 861 const SCEVConstant *ConstLoadStride = 862 dyn_cast<SCEVConstant>(LoadEv->getOperand(1)); 863 if (!ConstStoreStride || !ConstLoadStride) 864 return false; 865 866 APInt StoreStrideValue = ConstStoreStride->getAPInt(); 867 APInt LoadStrideValue = ConstLoadStride->getAPInt(); 868 // Huge stride value - give up 869 if (StoreStrideValue.getBitWidth() > 64 || LoadStrideValue.getBitWidth() > 64) 870 return false; 871 872 if (SizeInBytes != StoreStrideValue && SizeInBytes != -StoreStrideValue) { 873 ORE.emit([&]() { 874 return OptimizationRemarkMissed(DEBUG_TYPE, "SizeStrideUnequal", MCI) 875 << ore::NV("Inst", "memcpy") << " in " 876 << ore::NV("Function", MCI->getFunction()) 877 << " function will not be hoisted: " 878 << ore::NV("Reason", "memcpy size is not equal to stride"); 879 }); 880 return false; 881 } 882 883 int64_t StoreStrideInt = StoreStrideValue.getSExtValue(); 884 int64_t LoadStrideInt = LoadStrideValue.getSExtValue(); 885 // Check if the load stride matches the store stride. 886 if (StoreStrideInt != LoadStrideInt) 887 return false; 888 889 return processLoopStoreOfLoopLoad( 890 Dest, Source, SE->getConstant(Dest->getType(), SizeInBytes), 891 MCI->getDestAlign(), MCI->getSourceAlign(), MCI, MCI, StoreEv, LoadEv, 892 BECount); 893 } 894 895 /// processLoopMemSet - See if this memset can be promoted to a large memset. 896 bool LoopIdiomRecognize::processLoopMemSet(MemSetInst *MSI, 897 const SCEV *BECount) { 898 // We can only handle non-volatile memsets. 899 if (MSI->isVolatile()) 900 return false; 901 902 // If we're not allowed to hack on memset, we fail. 903 if (!HasMemset || DisableLIRP::Memset) 904 return false; 905 906 Value *Pointer = MSI->getDest(); 907 908 // See if the pointer expression is an AddRec like {base,+,1} on the current 909 // loop, which indicates a strided store. If we have something else, it's a 910 // random store we can't handle. 911 const SCEVAddRecExpr *Ev = dyn_cast<SCEVAddRecExpr>(SE->getSCEV(Pointer)); 912 if (!Ev || Ev->getLoop() != CurLoop) 913 return false; 914 if (!Ev->isAffine()) { 915 LLVM_DEBUG(dbgs() << " Pointer is not affine, abort\n"); 916 return false; 917 } 918 919 const SCEV *PointerStrideSCEV = Ev->getOperand(1); 920 const SCEV *MemsetSizeSCEV = SE->getSCEV(MSI->getLength()); 921 if (!PointerStrideSCEV || !MemsetSizeSCEV) 922 return false; 923 924 bool IsNegStride = false; 925 const bool IsConstantSize = isa<ConstantInt>(MSI->getLength()); 926 927 if (IsConstantSize) { 928 // Memset size is constant. 929 // Check if the pointer stride matches the memset size. If so, then 930 // we know that every byte is touched in the loop. 931 LLVM_DEBUG(dbgs() << " memset size is constant\n"); 932 uint64_t SizeInBytes = cast<ConstantInt>(MSI->getLength())->getZExtValue(); 933 const SCEVConstant *ConstStride = dyn_cast<SCEVConstant>(Ev->getOperand(1)); 934 if (!ConstStride) 935 return false; 936 937 APInt Stride = ConstStride->getAPInt(); 938 if (SizeInBytes != Stride && SizeInBytes != -Stride) 939 return false; 940 941 IsNegStride = SizeInBytes == -Stride; 942 } else { 943 // Memset size is non-constant. 944 // Check if the pointer stride matches the memset size. 945 // To be conservative, the pass would not promote pointers that aren't in 946 // address space zero. Also, the pass only handles memset length and stride 947 // that are invariant for the top level loop. 948 LLVM_DEBUG(dbgs() << " memset size is non-constant\n"); 949 if (Pointer->getType()->getPointerAddressSpace() != 0) { 950 LLVM_DEBUG(dbgs() << " pointer is not in address space zero, " 951 << "abort\n"); 952 return false; 953 } 954 if (!SE->isLoopInvariant(MemsetSizeSCEV, CurLoop)) { 955 LLVM_DEBUG(dbgs() << " memset size is not a loop-invariant, " 956 << "abort\n"); 957 return false; 958 } 959 960 // Compare positive direction PointerStrideSCEV with MemsetSizeSCEV 961 IsNegStride = PointerStrideSCEV->isNonConstantNegative(); 962 const SCEV *PositiveStrideSCEV = 963 IsNegStride ? SE->getNegativeSCEV(PointerStrideSCEV) 964 : PointerStrideSCEV; 965 LLVM_DEBUG(dbgs() << " MemsetSizeSCEV: " << *MemsetSizeSCEV << "\n" 966 << " PositiveStrideSCEV: " << *PositiveStrideSCEV 967 << "\n"); 968 969 if (PositiveStrideSCEV != MemsetSizeSCEV) { 970 // If an expression is covered by the loop guard, compare again and 971 // proceed with optimization if equal. 972 const SCEV *FoldedPositiveStride = 973 SE->applyLoopGuards(PositiveStrideSCEV, CurLoop); 974 const SCEV *FoldedMemsetSize = 975 SE->applyLoopGuards(MemsetSizeSCEV, CurLoop); 976 977 LLVM_DEBUG(dbgs() << " Try to fold SCEV based on loop guard\n" 978 << " FoldedMemsetSize: " << *FoldedMemsetSize << "\n" 979 << " FoldedPositiveStride: " << *FoldedPositiveStride 980 << "\n"); 981 982 if (FoldedPositiveStride != FoldedMemsetSize) { 983 LLVM_DEBUG(dbgs() << " SCEV don't match, abort\n"); 984 return false; 985 } 986 } 987 } 988 989 // Verify that the memset value is loop invariant. If not, we can't promote 990 // the memset. 991 Value *SplatValue = MSI->getValue(); 992 if (!SplatValue || !CurLoop->isLoopInvariant(SplatValue)) 993 return false; 994 995 SmallPtrSet<Instruction *, 1> MSIs; 996 MSIs.insert(MSI); 997 return processLoopStridedStore(Pointer, SE->getSCEV(MSI->getLength()), 998 MaybeAlign(MSI->getDestAlignment()), 999 SplatValue, MSI, MSIs, Ev, BECount, 1000 IsNegStride, /*IsLoopMemset=*/true); 1001 } 1002 1003 /// mayLoopAccessLocation - Return true if the specified loop might access the 1004 /// specified pointer location, which is a loop-strided access. The 'Access' 1005 /// argument specifies what the verboten forms of access are (read or write). 1006 static bool 1007 mayLoopAccessLocation(Value *Ptr, ModRefInfo Access, Loop *L, 1008 const SCEV *BECount, const SCEV *StoreSizeSCEV, 1009 AliasAnalysis &AA, 1010 SmallPtrSetImpl<Instruction *> &IgnoredInsts) { 1011 // Get the location that may be stored across the loop. Since the access is 1012 // strided positively through memory, we say that the modified location starts 1013 // at the pointer and has infinite size. 1014 LocationSize AccessSize = LocationSize::afterPointer(); 1015 1016 // If the loop iterates a fixed number of times, we can refine the access size 1017 // to be exactly the size of the memset, which is (BECount+1)*StoreSize 1018 const SCEVConstant *BECst = dyn_cast<SCEVConstant>(BECount); 1019 const SCEVConstant *ConstSize = dyn_cast<SCEVConstant>(StoreSizeSCEV); 1020 if (BECst && ConstSize) 1021 AccessSize = LocationSize::precise((BECst->getValue()->getZExtValue() + 1) * 1022 ConstSize->getValue()->getZExtValue()); 1023 1024 // TODO: For this to be really effective, we have to dive into the pointer 1025 // operand in the store. Store to &A[i] of 100 will always return may alias 1026 // with store of &A[100], we need to StoreLoc to be "A" with size of 100, 1027 // which will then no-alias a store to &A[100]. 1028 MemoryLocation StoreLoc(Ptr, AccessSize); 1029 1030 for (BasicBlock *B : L->blocks()) 1031 for (Instruction &I : *B) 1032 if (!IgnoredInsts.contains(&I) && 1033 isModOrRefSet( 1034 intersectModRef(AA.getModRefInfo(&I, StoreLoc), Access))) 1035 return true; 1036 return false; 1037 } 1038 1039 // If we have a negative stride, Start refers to the end of the memory location 1040 // we're trying to memset. Therefore, we need to recompute the base pointer, 1041 // which is just Start - BECount*Size. 1042 static const SCEV *getStartForNegStride(const SCEV *Start, const SCEV *BECount, 1043 Type *IntPtr, const SCEV *StoreSizeSCEV, 1044 ScalarEvolution *SE) { 1045 const SCEV *Index = SE->getTruncateOrZeroExtend(BECount, IntPtr); 1046 if (!StoreSizeSCEV->isOne()) { 1047 // index = back edge count * store size 1048 Index = SE->getMulExpr(Index, 1049 SE->getTruncateOrZeroExtend(StoreSizeSCEV, IntPtr), 1050 SCEV::FlagNUW); 1051 } 1052 // base pointer = start - index * store size 1053 return SE->getMinusSCEV(Start, Index); 1054 } 1055 1056 /// Compute trip count from the backedge taken count. 1057 static const SCEV *getTripCount(const SCEV *BECount, Type *IntPtr, 1058 Loop *CurLoop, const DataLayout *DL, 1059 ScalarEvolution *SE) { 1060 const SCEV *TripCountS = nullptr; 1061 // The # stored bytes is (BECount+1). Expand the trip count out to 1062 // pointer size if it isn't already. 1063 // 1064 // If we're going to need to zero extend the BE count, check if we can add 1065 // one to it prior to zero extending without overflow. Provided this is safe, 1066 // it allows better simplification of the +1. 1067 if (DL->getTypeSizeInBits(BECount->getType()) < 1068 DL->getTypeSizeInBits(IntPtr) && 1069 SE->isLoopEntryGuardedByCond( 1070 CurLoop, ICmpInst::ICMP_NE, BECount, 1071 SE->getNegativeSCEV(SE->getOne(BECount->getType())))) { 1072 TripCountS = SE->getZeroExtendExpr( 1073 SE->getAddExpr(BECount, SE->getOne(BECount->getType()), SCEV::FlagNUW), 1074 IntPtr); 1075 } else { 1076 TripCountS = SE->getAddExpr(SE->getTruncateOrZeroExtend(BECount, IntPtr), 1077 SE->getOne(IntPtr), SCEV::FlagNUW); 1078 } 1079 1080 return TripCountS; 1081 } 1082 1083 /// Compute the number of bytes as a SCEV from the backedge taken count. 1084 /// 1085 /// This also maps the SCEV into the provided type and tries to handle the 1086 /// computation in a way that will fold cleanly. 1087 static const SCEV *getNumBytes(const SCEV *BECount, Type *IntPtr, 1088 const SCEV *StoreSizeSCEV, Loop *CurLoop, 1089 const DataLayout *DL, ScalarEvolution *SE) { 1090 const SCEV *TripCountSCEV = getTripCount(BECount, IntPtr, CurLoop, DL, SE); 1091 1092 return SE->getMulExpr(TripCountSCEV, 1093 SE->getTruncateOrZeroExtend(StoreSizeSCEV, IntPtr), 1094 SCEV::FlagNUW); 1095 } 1096 1097 /// processLoopStridedStore - We see a strided store of some value. If we can 1098 /// transform this into a memset or memset_pattern in the loop preheader, do so. 1099 bool LoopIdiomRecognize::processLoopStridedStore( 1100 Value *DestPtr, const SCEV *StoreSizeSCEV, MaybeAlign StoreAlignment, 1101 Value *StoredVal, Instruction *TheStore, 1102 SmallPtrSetImpl<Instruction *> &Stores, const SCEVAddRecExpr *Ev, 1103 const SCEV *BECount, bool IsNegStride, bool IsLoopMemset) { 1104 Value *SplatValue = isBytewiseValue(StoredVal, *DL); 1105 Constant *PatternValue = nullptr; 1106 1107 if (!SplatValue) 1108 PatternValue = getMemSetPatternValue(StoredVal, DL); 1109 1110 assert((SplatValue || PatternValue) && 1111 "Expected either splat value or pattern value."); 1112 1113 // The trip count of the loop and the base pointer of the addrec SCEV is 1114 // guaranteed to be loop invariant, which means that it should dominate the 1115 // header. This allows us to insert code for it in the preheader. 1116 unsigned DestAS = DestPtr->getType()->getPointerAddressSpace(); 1117 BasicBlock *Preheader = CurLoop->getLoopPreheader(); 1118 IRBuilder<> Builder(Preheader->getTerminator()); 1119 SCEVExpander Expander(*SE, *DL, "loop-idiom"); 1120 SCEVExpanderCleaner ExpCleaner(Expander); 1121 1122 Type *DestInt8PtrTy = Builder.getInt8PtrTy(DestAS); 1123 Type *IntIdxTy = DL->getIndexType(DestPtr->getType()); 1124 1125 bool Changed = false; 1126 const SCEV *Start = Ev->getStart(); 1127 // Handle negative strided loops. 1128 if (IsNegStride) 1129 Start = getStartForNegStride(Start, BECount, IntIdxTy, StoreSizeSCEV, SE); 1130 1131 // TODO: ideally we should still be able to generate memset if SCEV expander 1132 // is taught to generate the dependencies at the latest point. 1133 if (!isSafeToExpand(Start, *SE)) 1134 return Changed; 1135 1136 // Okay, we have a strided store "p[i]" of a splattable value. We can turn 1137 // this into a memset in the loop preheader now if we want. However, this 1138 // would be unsafe to do if there is anything else in the loop that may read 1139 // or write to the aliased location. Check for any overlap by generating the 1140 // base pointer and checking the region. 1141 Value *BasePtr = 1142 Expander.expandCodeFor(Start, DestInt8PtrTy, Preheader->getTerminator()); 1143 1144 // From here on out, conservatively report to the pass manager that we've 1145 // changed the IR, even if we later clean up these added instructions. There 1146 // may be structural differences e.g. in the order of use lists not accounted 1147 // for in just a textual dump of the IR. This is written as a variable, even 1148 // though statically all the places this dominates could be replaced with 1149 // 'true', with the hope that anyone trying to be clever / "more precise" with 1150 // the return value will read this comment, and leave them alone. 1151 Changed = true; 1152 1153 if (mayLoopAccessLocation(BasePtr, ModRefInfo::ModRef, CurLoop, BECount, 1154 StoreSizeSCEV, *AA, Stores)) 1155 return Changed; 1156 1157 if (avoidLIRForMultiBlockLoop(/*IsMemset=*/true, IsLoopMemset)) 1158 return Changed; 1159 1160 // Okay, everything looks good, insert the memset. 1161 1162 const SCEV *NumBytesS = 1163 getNumBytes(BECount, IntIdxTy, StoreSizeSCEV, CurLoop, DL, SE); 1164 1165 // TODO: ideally we should still be able to generate memset if SCEV expander 1166 // is taught to generate the dependencies at the latest point. 1167 if (!isSafeToExpand(NumBytesS, *SE)) 1168 return Changed; 1169 1170 Value *NumBytes = 1171 Expander.expandCodeFor(NumBytesS, IntIdxTy, Preheader->getTerminator()); 1172 1173 CallInst *NewCall; 1174 if (SplatValue) { 1175 AAMDNodes AATags = TheStore->getAAMetadata(); 1176 if (auto CI = dyn_cast<ConstantInt>(NumBytes)) 1177 AATags = AATags.extendTo(CI->getZExtValue()); 1178 else 1179 AATags = AATags.extendTo(-1); 1180 1181 NewCall = Builder.CreateMemSet( 1182 BasePtr, SplatValue, NumBytes, MaybeAlign(StoreAlignment), 1183 /*isVolatile=*/false, AATags.TBAA, AATags.Scope, AATags.NoAlias); 1184 } else { 1185 // Everything is emitted in default address space 1186 Type *Int8PtrTy = DestInt8PtrTy; 1187 1188 Module *M = TheStore->getModule(); 1189 StringRef FuncName = "memset_pattern16"; 1190 FunctionCallee MSP = M->getOrInsertFunction(FuncName, Builder.getVoidTy(), 1191 Int8PtrTy, Int8PtrTy, IntIdxTy); 1192 inferLibFuncAttributes(M, FuncName, *TLI); 1193 1194 // Otherwise we should form a memset_pattern16. PatternValue is known to be 1195 // an constant array of 16-bytes. Plop the value into a mergable global. 1196 GlobalVariable *GV = new GlobalVariable(*M, PatternValue->getType(), true, 1197 GlobalValue::PrivateLinkage, 1198 PatternValue, ".memset_pattern"); 1199 GV->setUnnamedAddr(GlobalValue::UnnamedAddr::Global); // Ok to merge these. 1200 GV->setAlignment(Align(16)); 1201 Value *PatternPtr = ConstantExpr::getBitCast(GV, Int8PtrTy); 1202 NewCall = Builder.CreateCall(MSP, {BasePtr, PatternPtr, NumBytes}); 1203 } 1204 NewCall->setDebugLoc(TheStore->getDebugLoc()); 1205 1206 if (MSSAU) { 1207 MemoryAccess *NewMemAcc = MSSAU->createMemoryAccessInBB( 1208 NewCall, nullptr, NewCall->getParent(), MemorySSA::BeforeTerminator); 1209 MSSAU->insertDef(cast<MemoryDef>(NewMemAcc), true); 1210 } 1211 1212 LLVM_DEBUG(dbgs() << " Formed memset: " << *NewCall << "\n" 1213 << " from store to: " << *Ev << " at: " << *TheStore 1214 << "\n"); 1215 1216 ORE.emit([&]() { 1217 OptimizationRemark R(DEBUG_TYPE, "ProcessLoopStridedStore", 1218 NewCall->getDebugLoc(), Preheader); 1219 R << "Transformed loop-strided store in " 1220 << ore::NV("Function", TheStore->getFunction()) 1221 << " function into a call to " 1222 << ore::NV("NewFunction", NewCall->getCalledFunction()) 1223 << "() intrinsic"; 1224 if (!Stores.empty()) 1225 R << ore::setExtraArgs(); 1226 for (auto *I : Stores) { 1227 R << ore::NV("FromBlock", I->getParent()->getName()) 1228 << ore::NV("ToBlock", Preheader->getName()); 1229 } 1230 return R; 1231 }); 1232 1233 // Okay, the memset has been formed. Zap the original store and anything that 1234 // feeds into it. 1235 for (auto *I : Stores) { 1236 if (MSSAU) 1237 MSSAU->removeMemoryAccess(I, true); 1238 deleteDeadInstruction(I); 1239 } 1240 if (MSSAU && VerifyMemorySSA) 1241 MSSAU->getMemorySSA()->verifyMemorySSA(); 1242 ++NumMemSet; 1243 ExpCleaner.markResultUsed(); 1244 return true; 1245 } 1246 1247 /// If the stored value is a strided load in the same loop with the same stride 1248 /// this may be transformable into a memcpy. This kicks in for stuff like 1249 /// for (i) A[i] = B[i]; 1250 bool LoopIdiomRecognize::processLoopStoreOfLoopLoad(StoreInst *SI, 1251 const SCEV *BECount) { 1252 assert(SI->isUnordered() && "Expected only non-volatile non-ordered stores."); 1253 1254 Value *StorePtr = SI->getPointerOperand(); 1255 const SCEVAddRecExpr *StoreEv = cast<SCEVAddRecExpr>(SE->getSCEV(StorePtr)); 1256 unsigned StoreSize = DL->getTypeStoreSize(SI->getValueOperand()->getType()); 1257 1258 // The store must be feeding a non-volatile load. 1259 LoadInst *LI = cast<LoadInst>(SI->getValueOperand()); 1260 assert(LI->isUnordered() && "Expected only non-volatile non-ordered loads."); 1261 1262 // See if the pointer expression is an AddRec like {base,+,1} on the current 1263 // loop, which indicates a strided load. If we have something else, it's a 1264 // random load we can't handle. 1265 Value *LoadPtr = LI->getPointerOperand(); 1266 const SCEVAddRecExpr *LoadEv = cast<SCEVAddRecExpr>(SE->getSCEV(LoadPtr)); 1267 1268 const SCEV *StoreSizeSCEV = SE->getConstant(StorePtr->getType(), StoreSize); 1269 return processLoopStoreOfLoopLoad(StorePtr, LoadPtr, StoreSizeSCEV, 1270 SI->getAlign(), LI->getAlign(), SI, LI, 1271 StoreEv, LoadEv, BECount); 1272 } 1273 1274 class MemmoveVerifier { 1275 public: 1276 explicit MemmoveVerifier(const Value &LoadBasePtr, const Value &StoreBasePtr, 1277 const DataLayout &DL) 1278 : DL(DL), BP1(llvm::GetPointerBaseWithConstantOffset( 1279 LoadBasePtr.stripPointerCasts(), LoadOff, DL)), 1280 BP2(llvm::GetPointerBaseWithConstantOffset( 1281 StoreBasePtr.stripPointerCasts(), StoreOff, DL)), 1282 IsSameObject(BP1 == BP2) {} 1283 1284 bool loadAndStoreMayFormMemmove(unsigned StoreSize, bool IsNegStride, 1285 const Instruction &TheLoad, 1286 bool IsMemCpy) const { 1287 if (IsMemCpy) { 1288 // Ensure that LoadBasePtr is after StoreBasePtr or before StoreBasePtr 1289 // for negative stride. 1290 if ((!IsNegStride && LoadOff <= StoreOff) || 1291 (IsNegStride && LoadOff >= StoreOff)) 1292 return false; 1293 } else { 1294 // Ensure that LoadBasePtr is after StoreBasePtr or before StoreBasePtr 1295 // for negative stride. LoadBasePtr shouldn't overlap with StoreBasePtr. 1296 int64_t LoadSize = 1297 DL.getTypeSizeInBits(TheLoad.getType()).getFixedSize() / 8; 1298 if (BP1 != BP2 || LoadSize != int64_t(StoreSize)) 1299 return false; 1300 if ((!IsNegStride && LoadOff < StoreOff + int64_t(StoreSize)) || 1301 (IsNegStride && LoadOff + LoadSize > StoreOff)) 1302 return false; 1303 } 1304 return true; 1305 } 1306 1307 private: 1308 const DataLayout &DL; 1309 int64_t LoadOff = 0; 1310 int64_t StoreOff = 0; 1311 const Value *BP1; 1312 const Value *BP2; 1313 1314 public: 1315 const bool IsSameObject; 1316 }; 1317 1318 bool LoopIdiomRecognize::processLoopStoreOfLoopLoad( 1319 Value *DestPtr, Value *SourcePtr, const SCEV *StoreSizeSCEV, 1320 MaybeAlign StoreAlign, MaybeAlign LoadAlign, Instruction *TheStore, 1321 Instruction *TheLoad, const SCEVAddRecExpr *StoreEv, 1322 const SCEVAddRecExpr *LoadEv, const SCEV *BECount) { 1323 1324 // FIXME: until llvm.memcpy.inline supports dynamic sizes, we need to 1325 // conservatively bail here, since otherwise we may have to transform 1326 // llvm.memcpy.inline into llvm.memcpy which is illegal. 1327 if (isa<MemCpyInlineInst>(TheStore)) 1328 return false; 1329 1330 // The trip count of the loop and the base pointer of the addrec SCEV is 1331 // guaranteed to be loop invariant, which means that it should dominate the 1332 // header. This allows us to insert code for it in the preheader. 1333 BasicBlock *Preheader = CurLoop->getLoopPreheader(); 1334 IRBuilder<> Builder(Preheader->getTerminator()); 1335 SCEVExpander Expander(*SE, *DL, "loop-idiom"); 1336 1337 SCEVExpanderCleaner ExpCleaner(Expander); 1338 1339 bool Changed = false; 1340 const SCEV *StrStart = StoreEv->getStart(); 1341 unsigned StrAS = DestPtr->getType()->getPointerAddressSpace(); 1342 Type *IntIdxTy = Builder.getIntNTy(DL->getIndexSizeInBits(StrAS)); 1343 1344 APInt Stride = getStoreStride(StoreEv); 1345 const SCEVConstant *ConstStoreSize = dyn_cast<SCEVConstant>(StoreSizeSCEV); 1346 1347 // TODO: Deal with non-constant size; Currently expect constant store size 1348 assert(ConstStoreSize && "store size is expected to be a constant"); 1349 1350 int64_t StoreSize = ConstStoreSize->getValue()->getZExtValue(); 1351 bool IsNegStride = StoreSize == -Stride; 1352 1353 // Handle negative strided loops. 1354 if (IsNegStride) 1355 StrStart = 1356 getStartForNegStride(StrStart, BECount, IntIdxTy, StoreSizeSCEV, SE); 1357 1358 // Okay, we have a strided store "p[i]" of a loaded value. We can turn 1359 // this into a memcpy in the loop preheader now if we want. However, this 1360 // would be unsafe to do if there is anything else in the loop that may read 1361 // or write the memory region we're storing to. This includes the load that 1362 // feeds the stores. Check for an alias by generating the base address and 1363 // checking everything. 1364 Value *StoreBasePtr = Expander.expandCodeFor( 1365 StrStart, Builder.getInt8PtrTy(StrAS), Preheader->getTerminator()); 1366 1367 // From here on out, conservatively report to the pass manager that we've 1368 // changed the IR, even if we later clean up these added instructions. There 1369 // may be structural differences e.g. in the order of use lists not accounted 1370 // for in just a textual dump of the IR. This is written as a variable, even 1371 // though statically all the places this dominates could be replaced with 1372 // 'true', with the hope that anyone trying to be clever / "more precise" with 1373 // the return value will read this comment, and leave them alone. 1374 Changed = true; 1375 1376 SmallPtrSet<Instruction *, 2> IgnoredInsts; 1377 IgnoredInsts.insert(TheStore); 1378 1379 bool IsMemCpy = isa<MemCpyInst>(TheStore); 1380 const StringRef InstRemark = IsMemCpy ? "memcpy" : "load and store"; 1381 1382 bool LoopAccessStore = 1383 mayLoopAccessLocation(StoreBasePtr, ModRefInfo::ModRef, CurLoop, BECount, 1384 StoreSizeSCEV, *AA, IgnoredInsts); 1385 if (LoopAccessStore) { 1386 // For memmove case it's not enough to guarantee that loop doesn't access 1387 // TheStore and TheLoad. Additionally we need to make sure that TheStore is 1388 // the only user of TheLoad. 1389 if (!TheLoad->hasOneUse()) 1390 return Changed; 1391 IgnoredInsts.insert(TheLoad); 1392 if (mayLoopAccessLocation(StoreBasePtr, ModRefInfo::ModRef, CurLoop, 1393 BECount, StoreSizeSCEV, *AA, IgnoredInsts)) { 1394 ORE.emit([&]() { 1395 return OptimizationRemarkMissed(DEBUG_TYPE, "LoopMayAccessStore", 1396 TheStore) 1397 << ore::NV("Inst", InstRemark) << " in " 1398 << ore::NV("Function", TheStore->getFunction()) 1399 << " function will not be hoisted: " 1400 << ore::NV("Reason", "The loop may access store location"); 1401 }); 1402 return Changed; 1403 } 1404 IgnoredInsts.erase(TheLoad); 1405 } 1406 1407 const SCEV *LdStart = LoadEv->getStart(); 1408 unsigned LdAS = SourcePtr->getType()->getPointerAddressSpace(); 1409 1410 // Handle negative strided loops. 1411 if (IsNegStride) 1412 LdStart = 1413 getStartForNegStride(LdStart, BECount, IntIdxTy, StoreSizeSCEV, SE); 1414 1415 // For a memcpy, we have to make sure that the input array is not being 1416 // mutated by the loop. 1417 Value *LoadBasePtr = Expander.expandCodeFor( 1418 LdStart, Builder.getInt8PtrTy(LdAS), Preheader->getTerminator()); 1419 1420 // If the store is a memcpy instruction, we must check if it will write to 1421 // the load memory locations. So remove it from the ignored stores. 1422 if (IsMemCpy) 1423 IgnoredInsts.erase(TheStore); 1424 MemmoveVerifier Verifier(*LoadBasePtr, *StoreBasePtr, *DL); 1425 if (mayLoopAccessLocation(LoadBasePtr, ModRefInfo::Mod, CurLoop, BECount, 1426 StoreSizeSCEV, *AA, IgnoredInsts)) { 1427 if (!IsMemCpy) { 1428 ORE.emit([&]() { 1429 return OptimizationRemarkMissed(DEBUG_TYPE, "LoopMayAccessLoad", 1430 TheLoad) 1431 << ore::NV("Inst", InstRemark) << " in " 1432 << ore::NV("Function", TheStore->getFunction()) 1433 << " function will not be hoisted: " 1434 << ore::NV("Reason", "The loop may access load location"); 1435 }); 1436 return Changed; 1437 } 1438 // At this point loop may access load only for memcpy in same underlying 1439 // object. If that's not the case bail out. 1440 if (!Verifier.IsSameObject) 1441 return Changed; 1442 } 1443 1444 bool UseMemMove = IsMemCpy ? Verifier.IsSameObject : LoopAccessStore; 1445 if (UseMemMove) 1446 if (!Verifier.loadAndStoreMayFormMemmove(StoreSize, IsNegStride, *TheLoad, 1447 IsMemCpy)) 1448 return Changed; 1449 1450 if (avoidLIRForMultiBlockLoop()) 1451 return Changed; 1452 1453 // Okay, everything is safe, we can transform this! 1454 1455 const SCEV *NumBytesS = 1456 getNumBytes(BECount, IntIdxTy, StoreSizeSCEV, CurLoop, DL, SE); 1457 1458 Value *NumBytes = 1459 Expander.expandCodeFor(NumBytesS, IntIdxTy, Preheader->getTerminator()); 1460 1461 AAMDNodes AATags = TheLoad->getAAMetadata(); 1462 AAMDNodes StoreAATags = TheStore->getAAMetadata(); 1463 AATags = AATags.merge(StoreAATags); 1464 if (auto CI = dyn_cast<ConstantInt>(NumBytes)) 1465 AATags = AATags.extendTo(CI->getZExtValue()); 1466 else 1467 AATags = AATags.extendTo(-1); 1468 1469 CallInst *NewCall = nullptr; 1470 // Check whether to generate an unordered atomic memcpy: 1471 // If the load or store are atomic, then they must necessarily be unordered 1472 // by previous checks. 1473 if (!TheStore->isAtomic() && !TheLoad->isAtomic()) { 1474 if (UseMemMove) 1475 NewCall = Builder.CreateMemMove( 1476 StoreBasePtr, StoreAlign, LoadBasePtr, LoadAlign, NumBytes, 1477 /*isVolatile=*/false, AATags.TBAA, AATags.Scope, AATags.NoAlias); 1478 else 1479 NewCall = 1480 Builder.CreateMemCpy(StoreBasePtr, StoreAlign, LoadBasePtr, LoadAlign, 1481 NumBytes, /*isVolatile=*/false, AATags.TBAA, 1482 AATags.TBAAStruct, AATags.Scope, AATags.NoAlias); 1483 } else { 1484 // For now don't support unordered atomic memmove. 1485 if (UseMemMove) 1486 return Changed; 1487 // We cannot allow unaligned ops for unordered load/store, so reject 1488 // anything where the alignment isn't at least the element size. 1489 assert((StoreAlign.hasValue() && LoadAlign.hasValue()) && 1490 "Expect unordered load/store to have align."); 1491 if (StoreAlign.getValue() < StoreSize || LoadAlign.getValue() < StoreSize) 1492 return Changed; 1493 1494 // If the element.atomic memcpy is not lowered into explicit 1495 // loads/stores later, then it will be lowered into an element-size 1496 // specific lib call. If the lib call doesn't exist for our store size, then 1497 // we shouldn't generate the memcpy. 1498 if (StoreSize > TTI->getAtomicMemIntrinsicMaxElementSize()) 1499 return Changed; 1500 1501 // Create the call. 1502 // Note that unordered atomic loads/stores are *required* by the spec to 1503 // have an alignment but non-atomic loads/stores may not. 1504 NewCall = Builder.CreateElementUnorderedAtomicMemCpy( 1505 StoreBasePtr, StoreAlign.getValue(), LoadBasePtr, LoadAlign.getValue(), 1506 NumBytes, StoreSize, AATags.TBAA, AATags.TBAAStruct, AATags.Scope, 1507 AATags.NoAlias); 1508 } 1509 NewCall->setDebugLoc(TheStore->getDebugLoc()); 1510 1511 if (MSSAU) { 1512 MemoryAccess *NewMemAcc = MSSAU->createMemoryAccessInBB( 1513 NewCall, nullptr, NewCall->getParent(), MemorySSA::BeforeTerminator); 1514 MSSAU->insertDef(cast<MemoryDef>(NewMemAcc), true); 1515 } 1516 1517 LLVM_DEBUG(dbgs() << " Formed new call: " << *NewCall << "\n" 1518 << " from load ptr=" << *LoadEv << " at: " << *TheLoad 1519 << "\n" 1520 << " from store ptr=" << *StoreEv << " at: " << *TheStore 1521 << "\n"); 1522 1523 ORE.emit([&]() { 1524 return OptimizationRemark(DEBUG_TYPE, "ProcessLoopStoreOfLoopLoad", 1525 NewCall->getDebugLoc(), Preheader) 1526 << "Formed a call to " 1527 << ore::NV("NewFunction", NewCall->getCalledFunction()) 1528 << "() intrinsic from " << ore::NV("Inst", InstRemark) 1529 << " instruction in " << ore::NV("Function", TheStore->getFunction()) 1530 << " function" 1531 << ore::setExtraArgs() 1532 << ore::NV("FromBlock", TheStore->getParent()->getName()) 1533 << ore::NV("ToBlock", Preheader->getName()); 1534 }); 1535 1536 // Okay, a new call to memcpy/memmove has been formed. Zap the original store 1537 // and anything that feeds into it. 1538 if (MSSAU) 1539 MSSAU->removeMemoryAccess(TheStore, true); 1540 deleteDeadInstruction(TheStore); 1541 if (MSSAU && VerifyMemorySSA) 1542 MSSAU->getMemorySSA()->verifyMemorySSA(); 1543 if (UseMemMove) 1544 ++NumMemMove; 1545 else 1546 ++NumMemCpy; 1547 ExpCleaner.markResultUsed(); 1548 return true; 1549 } 1550 1551 // When compiling for codesize we avoid idiom recognition for a multi-block loop 1552 // unless it is a loop_memset idiom or a memset/memcpy idiom in a nested loop. 1553 // 1554 bool LoopIdiomRecognize::avoidLIRForMultiBlockLoop(bool IsMemset, 1555 bool IsLoopMemset) { 1556 if (ApplyCodeSizeHeuristics && CurLoop->getNumBlocks() > 1) { 1557 if (CurLoop->isOutermost() && (!IsMemset || !IsLoopMemset)) { 1558 LLVM_DEBUG(dbgs() << " " << CurLoop->getHeader()->getParent()->getName() 1559 << " : LIR " << (IsMemset ? "Memset" : "Memcpy") 1560 << " avoided: multi-block top-level loop\n"); 1561 return true; 1562 } 1563 } 1564 1565 return false; 1566 } 1567 1568 bool LoopIdiomRecognize::runOnNoncountableLoop() { 1569 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Scanning: F[" 1570 << CurLoop->getHeader()->getParent()->getName() 1571 << "] Noncountable Loop %" 1572 << CurLoop->getHeader()->getName() << "\n"); 1573 1574 return recognizePopcount() || recognizeAndInsertFFS() || 1575 recognizeShiftUntilBitTest() || recognizeShiftUntilZero(); 1576 } 1577 1578 /// Check if the given conditional branch is based on the comparison between 1579 /// a variable and zero, and if the variable is non-zero or zero (JmpOnZero is 1580 /// true), the control yields to the loop entry. If the branch matches the 1581 /// behavior, the variable involved in the comparison is returned. This function 1582 /// will be called to see if the precondition and postcondition of the loop are 1583 /// in desirable form. 1584 static Value *matchCondition(BranchInst *BI, BasicBlock *LoopEntry, 1585 bool JmpOnZero = false) { 1586 if (!BI || !BI->isConditional()) 1587 return nullptr; 1588 1589 ICmpInst *Cond = dyn_cast<ICmpInst>(BI->getCondition()); 1590 if (!Cond) 1591 return nullptr; 1592 1593 ConstantInt *CmpZero = dyn_cast<ConstantInt>(Cond->getOperand(1)); 1594 if (!CmpZero || !CmpZero->isZero()) 1595 return nullptr; 1596 1597 BasicBlock *TrueSucc = BI->getSuccessor(0); 1598 BasicBlock *FalseSucc = BI->getSuccessor(1); 1599 if (JmpOnZero) 1600 std::swap(TrueSucc, FalseSucc); 1601 1602 ICmpInst::Predicate Pred = Cond->getPredicate(); 1603 if ((Pred == ICmpInst::ICMP_NE && TrueSucc == LoopEntry) || 1604 (Pred == ICmpInst::ICMP_EQ && FalseSucc == LoopEntry)) 1605 return Cond->getOperand(0); 1606 1607 return nullptr; 1608 } 1609 1610 // Check if the recurrence variable `VarX` is in the right form to create 1611 // the idiom. Returns the value coerced to a PHINode if so. 1612 static PHINode *getRecurrenceVar(Value *VarX, Instruction *DefX, 1613 BasicBlock *LoopEntry) { 1614 auto *PhiX = dyn_cast<PHINode>(VarX); 1615 if (PhiX && PhiX->getParent() == LoopEntry && 1616 (PhiX->getOperand(0) == DefX || PhiX->getOperand(1) == DefX)) 1617 return PhiX; 1618 return nullptr; 1619 } 1620 1621 /// Return true iff the idiom is detected in the loop. 1622 /// 1623 /// Additionally: 1624 /// 1) \p CntInst is set to the instruction counting the population bit. 1625 /// 2) \p CntPhi is set to the corresponding phi node. 1626 /// 3) \p Var is set to the value whose population bits are being counted. 1627 /// 1628 /// The core idiom we are trying to detect is: 1629 /// \code 1630 /// if (x0 != 0) 1631 /// goto loop-exit // the precondition of the loop 1632 /// cnt0 = init-val; 1633 /// do { 1634 /// x1 = phi (x0, x2); 1635 /// cnt1 = phi(cnt0, cnt2); 1636 /// 1637 /// cnt2 = cnt1 + 1; 1638 /// ... 1639 /// x2 = x1 & (x1 - 1); 1640 /// ... 1641 /// } while(x != 0); 1642 /// 1643 /// loop-exit: 1644 /// \endcode 1645 static bool detectPopcountIdiom(Loop *CurLoop, BasicBlock *PreCondBB, 1646 Instruction *&CntInst, PHINode *&CntPhi, 1647 Value *&Var) { 1648 // step 1: Check to see if the look-back branch match this pattern: 1649 // "if (a!=0) goto loop-entry". 1650 BasicBlock *LoopEntry; 1651 Instruction *DefX2, *CountInst; 1652 Value *VarX1, *VarX0; 1653 PHINode *PhiX, *CountPhi; 1654 1655 DefX2 = CountInst = nullptr; 1656 VarX1 = VarX0 = nullptr; 1657 PhiX = CountPhi = nullptr; 1658 LoopEntry = *(CurLoop->block_begin()); 1659 1660 // step 1: Check if the loop-back branch is in desirable form. 1661 { 1662 if (Value *T = matchCondition( 1663 dyn_cast<BranchInst>(LoopEntry->getTerminator()), LoopEntry)) 1664 DefX2 = dyn_cast<Instruction>(T); 1665 else 1666 return false; 1667 } 1668 1669 // step 2: detect instructions corresponding to "x2 = x1 & (x1 - 1)" 1670 { 1671 if (!DefX2 || DefX2->getOpcode() != Instruction::And) 1672 return false; 1673 1674 BinaryOperator *SubOneOp; 1675 1676 if ((SubOneOp = dyn_cast<BinaryOperator>(DefX2->getOperand(0)))) 1677 VarX1 = DefX2->getOperand(1); 1678 else { 1679 VarX1 = DefX2->getOperand(0); 1680 SubOneOp = dyn_cast<BinaryOperator>(DefX2->getOperand(1)); 1681 } 1682 if (!SubOneOp || SubOneOp->getOperand(0) != VarX1) 1683 return false; 1684 1685 ConstantInt *Dec = dyn_cast<ConstantInt>(SubOneOp->getOperand(1)); 1686 if (!Dec || 1687 !((SubOneOp->getOpcode() == Instruction::Sub && Dec->isOne()) || 1688 (SubOneOp->getOpcode() == Instruction::Add && 1689 Dec->isMinusOne()))) { 1690 return false; 1691 } 1692 } 1693 1694 // step 3: Check the recurrence of variable X 1695 PhiX = getRecurrenceVar(VarX1, DefX2, LoopEntry); 1696 if (!PhiX) 1697 return false; 1698 1699 // step 4: Find the instruction which count the population: cnt2 = cnt1 + 1 1700 { 1701 CountInst = nullptr; 1702 for (Instruction &Inst : llvm::make_range( 1703 LoopEntry->getFirstNonPHI()->getIterator(), LoopEntry->end())) { 1704 if (Inst.getOpcode() != Instruction::Add) 1705 continue; 1706 1707 ConstantInt *Inc = dyn_cast<ConstantInt>(Inst.getOperand(1)); 1708 if (!Inc || !Inc->isOne()) 1709 continue; 1710 1711 PHINode *Phi = getRecurrenceVar(Inst.getOperand(0), &Inst, LoopEntry); 1712 if (!Phi) 1713 continue; 1714 1715 // Check if the result of the instruction is live of the loop. 1716 bool LiveOutLoop = false; 1717 for (User *U : Inst.users()) { 1718 if ((cast<Instruction>(U))->getParent() != LoopEntry) { 1719 LiveOutLoop = true; 1720 break; 1721 } 1722 } 1723 1724 if (LiveOutLoop) { 1725 CountInst = &Inst; 1726 CountPhi = Phi; 1727 break; 1728 } 1729 } 1730 1731 if (!CountInst) 1732 return false; 1733 } 1734 1735 // step 5: check if the precondition is in this form: 1736 // "if (x != 0) goto loop-head ; else goto somewhere-we-don't-care;" 1737 { 1738 auto *PreCondBr = dyn_cast<BranchInst>(PreCondBB->getTerminator()); 1739 Value *T = matchCondition(PreCondBr, CurLoop->getLoopPreheader()); 1740 if (T != PhiX->getOperand(0) && T != PhiX->getOperand(1)) 1741 return false; 1742 1743 CntInst = CountInst; 1744 CntPhi = CountPhi; 1745 Var = T; 1746 } 1747 1748 return true; 1749 } 1750 1751 /// Return true if the idiom is detected in the loop. 1752 /// 1753 /// Additionally: 1754 /// 1) \p CntInst is set to the instruction Counting Leading Zeros (CTLZ) 1755 /// or nullptr if there is no such. 1756 /// 2) \p CntPhi is set to the corresponding phi node 1757 /// or nullptr if there is no such. 1758 /// 3) \p Var is set to the value whose CTLZ could be used. 1759 /// 4) \p DefX is set to the instruction calculating Loop exit condition. 1760 /// 1761 /// The core idiom we are trying to detect is: 1762 /// \code 1763 /// if (x0 == 0) 1764 /// goto loop-exit // the precondition of the loop 1765 /// cnt0 = init-val; 1766 /// do { 1767 /// x = phi (x0, x.next); //PhiX 1768 /// cnt = phi(cnt0, cnt.next); 1769 /// 1770 /// cnt.next = cnt + 1; 1771 /// ... 1772 /// x.next = x >> 1; // DefX 1773 /// ... 1774 /// } while(x.next != 0); 1775 /// 1776 /// loop-exit: 1777 /// \endcode 1778 static bool detectShiftUntilZeroIdiom(Loop *CurLoop, const DataLayout &DL, 1779 Intrinsic::ID &IntrinID, Value *&InitX, 1780 Instruction *&CntInst, PHINode *&CntPhi, 1781 Instruction *&DefX) { 1782 BasicBlock *LoopEntry; 1783 Value *VarX = nullptr; 1784 1785 DefX = nullptr; 1786 CntInst = nullptr; 1787 CntPhi = nullptr; 1788 LoopEntry = *(CurLoop->block_begin()); 1789 1790 // step 1: Check if the loop-back branch is in desirable form. 1791 if (Value *T = matchCondition( 1792 dyn_cast<BranchInst>(LoopEntry->getTerminator()), LoopEntry)) 1793 DefX = dyn_cast<Instruction>(T); 1794 else 1795 return false; 1796 1797 // step 2: detect instructions corresponding to "x.next = x >> 1 or x << 1" 1798 if (!DefX || !DefX->isShift()) 1799 return false; 1800 IntrinID = DefX->getOpcode() == Instruction::Shl ? Intrinsic::cttz : 1801 Intrinsic::ctlz; 1802 ConstantInt *Shft = dyn_cast<ConstantInt>(DefX->getOperand(1)); 1803 if (!Shft || !Shft->isOne()) 1804 return false; 1805 VarX = DefX->getOperand(0); 1806 1807 // step 3: Check the recurrence of variable X 1808 PHINode *PhiX = getRecurrenceVar(VarX, DefX, LoopEntry); 1809 if (!PhiX) 1810 return false; 1811 1812 InitX = PhiX->getIncomingValueForBlock(CurLoop->getLoopPreheader()); 1813 1814 // Make sure the initial value can't be negative otherwise the ashr in the 1815 // loop might never reach zero which would make the loop infinite. 1816 if (DefX->getOpcode() == Instruction::AShr && !isKnownNonNegative(InitX, DL)) 1817 return false; 1818 1819 // step 4: Find the instruction which count the CTLZ: cnt.next = cnt + 1 1820 // or cnt.next = cnt + -1. 1821 // TODO: We can skip the step. If loop trip count is known (CTLZ), 1822 // then all uses of "cnt.next" could be optimized to the trip count 1823 // plus "cnt0". Currently it is not optimized. 1824 // This step could be used to detect POPCNT instruction: 1825 // cnt.next = cnt + (x.next & 1) 1826 for (Instruction &Inst : llvm::make_range( 1827 LoopEntry->getFirstNonPHI()->getIterator(), LoopEntry->end())) { 1828 if (Inst.getOpcode() != Instruction::Add) 1829 continue; 1830 1831 ConstantInt *Inc = dyn_cast<ConstantInt>(Inst.getOperand(1)); 1832 if (!Inc || (!Inc->isOne() && !Inc->isMinusOne())) 1833 continue; 1834 1835 PHINode *Phi = getRecurrenceVar(Inst.getOperand(0), &Inst, LoopEntry); 1836 if (!Phi) 1837 continue; 1838 1839 CntInst = &Inst; 1840 CntPhi = Phi; 1841 break; 1842 } 1843 if (!CntInst) 1844 return false; 1845 1846 return true; 1847 } 1848 1849 /// Recognize CTLZ or CTTZ idiom in a non-countable loop and convert the loop 1850 /// to countable (with CTLZ / CTTZ trip count). If CTLZ / CTTZ inserted as a new 1851 /// trip count returns true; otherwise, returns false. 1852 bool LoopIdiomRecognize::recognizeAndInsertFFS() { 1853 // Give up if the loop has multiple blocks or multiple backedges. 1854 if (CurLoop->getNumBackEdges() != 1 || CurLoop->getNumBlocks() != 1) 1855 return false; 1856 1857 Intrinsic::ID IntrinID; 1858 Value *InitX; 1859 Instruction *DefX = nullptr; 1860 PHINode *CntPhi = nullptr; 1861 Instruction *CntInst = nullptr; 1862 // Help decide if transformation is profitable. For ShiftUntilZero idiom, 1863 // this is always 6. 1864 size_t IdiomCanonicalSize = 6; 1865 1866 if (!detectShiftUntilZeroIdiom(CurLoop, *DL, IntrinID, InitX, 1867 CntInst, CntPhi, DefX)) 1868 return false; 1869 1870 bool IsCntPhiUsedOutsideLoop = false; 1871 for (User *U : CntPhi->users()) 1872 if (!CurLoop->contains(cast<Instruction>(U))) { 1873 IsCntPhiUsedOutsideLoop = true; 1874 break; 1875 } 1876 bool IsCntInstUsedOutsideLoop = false; 1877 for (User *U : CntInst->users()) 1878 if (!CurLoop->contains(cast<Instruction>(U))) { 1879 IsCntInstUsedOutsideLoop = true; 1880 break; 1881 } 1882 // If both CntInst and CntPhi are used outside the loop the profitability 1883 // is questionable. 1884 if (IsCntInstUsedOutsideLoop && IsCntPhiUsedOutsideLoop) 1885 return false; 1886 1887 // For some CPUs result of CTLZ(X) intrinsic is undefined 1888 // when X is 0. If we can not guarantee X != 0, we need to check this 1889 // when expand. 1890 bool ZeroCheck = false; 1891 // It is safe to assume Preheader exist as it was checked in 1892 // parent function RunOnLoop. 1893 BasicBlock *PH = CurLoop->getLoopPreheader(); 1894 1895 // If we are using the count instruction outside the loop, make sure we 1896 // have a zero check as a precondition. Without the check the loop would run 1897 // one iteration for before any check of the input value. This means 0 and 1 1898 // would have identical behavior in the original loop and thus 1899 if (!IsCntPhiUsedOutsideLoop) { 1900 auto *PreCondBB = PH->getSinglePredecessor(); 1901 if (!PreCondBB) 1902 return false; 1903 auto *PreCondBI = dyn_cast<BranchInst>(PreCondBB->getTerminator()); 1904 if (!PreCondBI) 1905 return false; 1906 if (matchCondition(PreCondBI, PH) != InitX) 1907 return false; 1908 ZeroCheck = true; 1909 } 1910 1911 // Check if CTLZ / CTTZ intrinsic is profitable. Assume it is always 1912 // profitable if we delete the loop. 1913 1914 // the loop has only 6 instructions: 1915 // %n.addr.0 = phi [ %n, %entry ], [ %shr, %while.cond ] 1916 // %i.0 = phi [ %i0, %entry ], [ %inc, %while.cond ] 1917 // %shr = ashr %n.addr.0, 1 1918 // %tobool = icmp eq %shr, 0 1919 // %inc = add nsw %i.0, 1 1920 // br i1 %tobool 1921 1922 const Value *Args[] = {InitX, 1923 ConstantInt::getBool(InitX->getContext(), ZeroCheck)}; 1924 1925 // @llvm.dbg doesn't count as they have no semantic effect. 1926 auto InstWithoutDebugIt = CurLoop->getHeader()->instructionsWithoutDebug(); 1927 uint32_t HeaderSize = 1928 std::distance(InstWithoutDebugIt.begin(), InstWithoutDebugIt.end()); 1929 1930 IntrinsicCostAttributes Attrs(IntrinID, InitX->getType(), Args); 1931 InstructionCost Cost = 1932 TTI->getIntrinsicInstrCost(Attrs, TargetTransformInfo::TCK_SizeAndLatency); 1933 if (HeaderSize != IdiomCanonicalSize && 1934 Cost > TargetTransformInfo::TCC_Basic) 1935 return false; 1936 1937 transformLoopToCountable(IntrinID, PH, CntInst, CntPhi, InitX, DefX, 1938 DefX->getDebugLoc(), ZeroCheck, 1939 IsCntPhiUsedOutsideLoop); 1940 return true; 1941 } 1942 1943 /// Recognizes a population count idiom in a non-countable loop. 1944 /// 1945 /// If detected, transforms the relevant code to issue the popcount intrinsic 1946 /// function call, and returns true; otherwise, returns false. 1947 bool LoopIdiomRecognize::recognizePopcount() { 1948 if (TTI->getPopcntSupport(32) != TargetTransformInfo::PSK_FastHardware) 1949 return false; 1950 1951 // Counting population are usually conducted by few arithmetic instructions. 1952 // Such instructions can be easily "absorbed" by vacant slots in a 1953 // non-compact loop. Therefore, recognizing popcount idiom only makes sense 1954 // in a compact loop. 1955 1956 // Give up if the loop has multiple blocks or multiple backedges. 1957 if (CurLoop->getNumBackEdges() != 1 || CurLoop->getNumBlocks() != 1) 1958 return false; 1959 1960 BasicBlock *LoopBody = *(CurLoop->block_begin()); 1961 if (LoopBody->size() >= 20) { 1962 // The loop is too big, bail out. 1963 return false; 1964 } 1965 1966 // It should have a preheader containing nothing but an unconditional branch. 1967 BasicBlock *PH = CurLoop->getLoopPreheader(); 1968 if (!PH || &PH->front() != PH->getTerminator()) 1969 return false; 1970 auto *EntryBI = dyn_cast<BranchInst>(PH->getTerminator()); 1971 if (!EntryBI || EntryBI->isConditional()) 1972 return false; 1973 1974 // It should have a precondition block where the generated popcount intrinsic 1975 // function can be inserted. 1976 auto *PreCondBB = PH->getSinglePredecessor(); 1977 if (!PreCondBB) 1978 return false; 1979 auto *PreCondBI = dyn_cast<BranchInst>(PreCondBB->getTerminator()); 1980 if (!PreCondBI || PreCondBI->isUnconditional()) 1981 return false; 1982 1983 Instruction *CntInst; 1984 PHINode *CntPhi; 1985 Value *Val; 1986 if (!detectPopcountIdiom(CurLoop, PreCondBB, CntInst, CntPhi, Val)) 1987 return false; 1988 1989 transformLoopToPopcount(PreCondBB, CntInst, CntPhi, Val); 1990 return true; 1991 } 1992 1993 static CallInst *createPopcntIntrinsic(IRBuilder<> &IRBuilder, Value *Val, 1994 const DebugLoc &DL) { 1995 Value *Ops[] = {Val}; 1996 Type *Tys[] = {Val->getType()}; 1997 1998 Module *M = IRBuilder.GetInsertBlock()->getParent()->getParent(); 1999 Function *Func = Intrinsic::getDeclaration(M, Intrinsic::ctpop, Tys); 2000 CallInst *CI = IRBuilder.CreateCall(Func, Ops); 2001 CI->setDebugLoc(DL); 2002 2003 return CI; 2004 } 2005 2006 static CallInst *createFFSIntrinsic(IRBuilder<> &IRBuilder, Value *Val, 2007 const DebugLoc &DL, bool ZeroCheck, 2008 Intrinsic::ID IID) { 2009 Value *Ops[] = {Val, IRBuilder.getInt1(ZeroCheck)}; 2010 Type *Tys[] = {Val->getType()}; 2011 2012 Module *M = IRBuilder.GetInsertBlock()->getParent()->getParent(); 2013 Function *Func = Intrinsic::getDeclaration(M, IID, Tys); 2014 CallInst *CI = IRBuilder.CreateCall(Func, Ops); 2015 CI->setDebugLoc(DL); 2016 2017 return CI; 2018 } 2019 2020 /// Transform the following loop (Using CTLZ, CTTZ is similar): 2021 /// loop: 2022 /// CntPhi = PHI [Cnt0, CntInst] 2023 /// PhiX = PHI [InitX, DefX] 2024 /// CntInst = CntPhi + 1 2025 /// DefX = PhiX >> 1 2026 /// LOOP_BODY 2027 /// Br: loop if (DefX != 0) 2028 /// Use(CntPhi) or Use(CntInst) 2029 /// 2030 /// Into: 2031 /// If CntPhi used outside the loop: 2032 /// CountPrev = BitWidth(InitX) - CTLZ(InitX >> 1) 2033 /// Count = CountPrev + 1 2034 /// else 2035 /// Count = BitWidth(InitX) - CTLZ(InitX) 2036 /// loop: 2037 /// CntPhi = PHI [Cnt0, CntInst] 2038 /// PhiX = PHI [InitX, DefX] 2039 /// PhiCount = PHI [Count, Dec] 2040 /// CntInst = CntPhi + 1 2041 /// DefX = PhiX >> 1 2042 /// Dec = PhiCount - 1 2043 /// LOOP_BODY 2044 /// Br: loop if (Dec != 0) 2045 /// Use(CountPrev + Cnt0) // Use(CntPhi) 2046 /// or 2047 /// Use(Count + Cnt0) // Use(CntInst) 2048 /// 2049 /// If LOOP_BODY is empty the loop will be deleted. 2050 /// If CntInst and DefX are not used in LOOP_BODY they will be removed. 2051 void LoopIdiomRecognize::transformLoopToCountable( 2052 Intrinsic::ID IntrinID, BasicBlock *Preheader, Instruction *CntInst, 2053 PHINode *CntPhi, Value *InitX, Instruction *DefX, const DebugLoc &DL, 2054 bool ZeroCheck, bool IsCntPhiUsedOutsideLoop) { 2055 BranchInst *PreheaderBr = cast<BranchInst>(Preheader->getTerminator()); 2056 2057 // Step 1: Insert the CTLZ/CTTZ instruction at the end of the preheader block 2058 IRBuilder<> Builder(PreheaderBr); 2059 Builder.SetCurrentDebugLocation(DL); 2060 2061 // If there are no uses of CntPhi crate: 2062 // Count = BitWidth - CTLZ(InitX); 2063 // NewCount = Count; 2064 // If there are uses of CntPhi create: 2065 // NewCount = BitWidth - CTLZ(InitX >> 1); 2066 // Count = NewCount + 1; 2067 Value *InitXNext; 2068 if (IsCntPhiUsedOutsideLoop) { 2069 if (DefX->getOpcode() == Instruction::AShr) 2070 InitXNext = Builder.CreateAShr(InitX, 1); 2071 else if (DefX->getOpcode() == Instruction::LShr) 2072 InitXNext = Builder.CreateLShr(InitX, 1); 2073 else if (DefX->getOpcode() == Instruction::Shl) // cttz 2074 InitXNext = Builder.CreateShl(InitX, 1); 2075 else 2076 llvm_unreachable("Unexpected opcode!"); 2077 } else 2078 InitXNext = InitX; 2079 Value *Count = 2080 createFFSIntrinsic(Builder, InitXNext, DL, ZeroCheck, IntrinID); 2081 Type *CountTy = Count->getType(); 2082 Count = Builder.CreateSub( 2083 ConstantInt::get(CountTy, CountTy->getIntegerBitWidth()), Count); 2084 Value *NewCount = Count; 2085 if (IsCntPhiUsedOutsideLoop) 2086 Count = Builder.CreateAdd(Count, ConstantInt::get(CountTy, 1)); 2087 2088 NewCount = Builder.CreateZExtOrTrunc(NewCount, CntInst->getType()); 2089 2090 Value *CntInitVal = CntPhi->getIncomingValueForBlock(Preheader); 2091 if (cast<ConstantInt>(CntInst->getOperand(1))->isOne()) { 2092 // If the counter was being incremented in the loop, add NewCount to the 2093 // counter's initial value, but only if the initial value is not zero. 2094 ConstantInt *InitConst = dyn_cast<ConstantInt>(CntInitVal); 2095 if (!InitConst || !InitConst->isZero()) 2096 NewCount = Builder.CreateAdd(NewCount, CntInitVal); 2097 } else { 2098 // If the count was being decremented in the loop, subtract NewCount from 2099 // the counter's initial value. 2100 NewCount = Builder.CreateSub(CntInitVal, NewCount); 2101 } 2102 2103 // Step 2: Insert new IV and loop condition: 2104 // loop: 2105 // ... 2106 // PhiCount = PHI [Count, Dec] 2107 // ... 2108 // Dec = PhiCount - 1 2109 // ... 2110 // Br: loop if (Dec != 0) 2111 BasicBlock *Body = *(CurLoop->block_begin()); 2112 auto *LbBr = cast<BranchInst>(Body->getTerminator()); 2113 ICmpInst *LbCond = cast<ICmpInst>(LbBr->getCondition()); 2114 2115 PHINode *TcPhi = PHINode::Create(CountTy, 2, "tcphi", &Body->front()); 2116 2117 Builder.SetInsertPoint(LbCond); 2118 Instruction *TcDec = cast<Instruction>(Builder.CreateSub( 2119 TcPhi, ConstantInt::get(CountTy, 1), "tcdec", false, true)); 2120 2121 TcPhi->addIncoming(Count, Preheader); 2122 TcPhi->addIncoming(TcDec, Body); 2123 2124 CmpInst::Predicate Pred = 2125 (LbBr->getSuccessor(0) == Body) ? CmpInst::ICMP_NE : CmpInst::ICMP_EQ; 2126 LbCond->setPredicate(Pred); 2127 LbCond->setOperand(0, TcDec); 2128 LbCond->setOperand(1, ConstantInt::get(CountTy, 0)); 2129 2130 // Step 3: All the references to the original counter outside 2131 // the loop are replaced with the NewCount 2132 if (IsCntPhiUsedOutsideLoop) 2133 CntPhi->replaceUsesOutsideBlock(NewCount, Body); 2134 else 2135 CntInst->replaceUsesOutsideBlock(NewCount, Body); 2136 2137 // step 4: Forget the "non-computable" trip-count SCEV associated with the 2138 // loop. The loop would otherwise not be deleted even if it becomes empty. 2139 SE->forgetLoop(CurLoop); 2140 } 2141 2142 void LoopIdiomRecognize::transformLoopToPopcount(BasicBlock *PreCondBB, 2143 Instruction *CntInst, 2144 PHINode *CntPhi, Value *Var) { 2145 BasicBlock *PreHead = CurLoop->getLoopPreheader(); 2146 auto *PreCondBr = cast<BranchInst>(PreCondBB->getTerminator()); 2147 const DebugLoc &DL = CntInst->getDebugLoc(); 2148 2149 // Assuming before transformation, the loop is following: 2150 // if (x) // the precondition 2151 // do { cnt++; x &= x - 1; } while(x); 2152 2153 // Step 1: Insert the ctpop instruction at the end of the precondition block 2154 IRBuilder<> Builder(PreCondBr); 2155 Value *PopCnt, *PopCntZext, *NewCount, *TripCnt; 2156 { 2157 PopCnt = createPopcntIntrinsic(Builder, Var, DL); 2158 NewCount = PopCntZext = 2159 Builder.CreateZExtOrTrunc(PopCnt, cast<IntegerType>(CntPhi->getType())); 2160 2161 if (NewCount != PopCnt) 2162 (cast<Instruction>(NewCount))->setDebugLoc(DL); 2163 2164 // TripCnt is exactly the number of iterations the loop has 2165 TripCnt = NewCount; 2166 2167 // If the population counter's initial value is not zero, insert Add Inst. 2168 Value *CntInitVal = CntPhi->getIncomingValueForBlock(PreHead); 2169 ConstantInt *InitConst = dyn_cast<ConstantInt>(CntInitVal); 2170 if (!InitConst || !InitConst->isZero()) { 2171 NewCount = Builder.CreateAdd(NewCount, CntInitVal); 2172 (cast<Instruction>(NewCount))->setDebugLoc(DL); 2173 } 2174 } 2175 2176 // Step 2: Replace the precondition from "if (x == 0) goto loop-exit" to 2177 // "if (NewCount == 0) loop-exit". Without this change, the intrinsic 2178 // function would be partial dead code, and downstream passes will drag 2179 // it back from the precondition block to the preheader. 2180 { 2181 ICmpInst *PreCond = cast<ICmpInst>(PreCondBr->getCondition()); 2182 2183 Value *Opnd0 = PopCntZext; 2184 Value *Opnd1 = ConstantInt::get(PopCntZext->getType(), 0); 2185 if (PreCond->getOperand(0) != Var) 2186 std::swap(Opnd0, Opnd1); 2187 2188 ICmpInst *NewPreCond = cast<ICmpInst>( 2189 Builder.CreateICmp(PreCond->getPredicate(), Opnd0, Opnd1)); 2190 PreCondBr->setCondition(NewPreCond); 2191 2192 RecursivelyDeleteTriviallyDeadInstructions(PreCond, TLI); 2193 } 2194 2195 // Step 3: Note that the population count is exactly the trip count of the 2196 // loop in question, which enable us to convert the loop from noncountable 2197 // loop into a countable one. The benefit is twofold: 2198 // 2199 // - If the loop only counts population, the entire loop becomes dead after 2200 // the transformation. It is a lot easier to prove a countable loop dead 2201 // than to prove a noncountable one. (In some C dialects, an infinite loop 2202 // isn't dead even if it computes nothing useful. In general, DCE needs 2203 // to prove a noncountable loop finite before safely delete it.) 2204 // 2205 // - If the loop also performs something else, it remains alive. 2206 // Since it is transformed to countable form, it can be aggressively 2207 // optimized by some optimizations which are in general not applicable 2208 // to a noncountable loop. 2209 // 2210 // After this step, this loop (conceptually) would look like following: 2211 // newcnt = __builtin_ctpop(x); 2212 // t = newcnt; 2213 // if (x) 2214 // do { cnt++; x &= x-1; t--) } while (t > 0); 2215 BasicBlock *Body = *(CurLoop->block_begin()); 2216 { 2217 auto *LbBr = cast<BranchInst>(Body->getTerminator()); 2218 ICmpInst *LbCond = cast<ICmpInst>(LbBr->getCondition()); 2219 Type *Ty = TripCnt->getType(); 2220 2221 PHINode *TcPhi = PHINode::Create(Ty, 2, "tcphi", &Body->front()); 2222 2223 Builder.SetInsertPoint(LbCond); 2224 Instruction *TcDec = cast<Instruction>( 2225 Builder.CreateSub(TcPhi, ConstantInt::get(Ty, 1), 2226 "tcdec", false, true)); 2227 2228 TcPhi->addIncoming(TripCnt, PreHead); 2229 TcPhi->addIncoming(TcDec, Body); 2230 2231 CmpInst::Predicate Pred = 2232 (LbBr->getSuccessor(0) == Body) ? CmpInst::ICMP_UGT : CmpInst::ICMP_SLE; 2233 LbCond->setPredicate(Pred); 2234 LbCond->setOperand(0, TcDec); 2235 LbCond->setOperand(1, ConstantInt::get(Ty, 0)); 2236 } 2237 2238 // Step 4: All the references to the original population counter outside 2239 // the loop are replaced with the NewCount -- the value returned from 2240 // __builtin_ctpop(). 2241 CntInst->replaceUsesOutsideBlock(NewCount, Body); 2242 2243 // step 5: Forget the "non-computable" trip-count SCEV associated with the 2244 // loop. The loop would otherwise not be deleted even if it becomes empty. 2245 SE->forgetLoop(CurLoop); 2246 } 2247 2248 /// Match loop-invariant value. 2249 template <typename SubPattern_t> struct match_LoopInvariant { 2250 SubPattern_t SubPattern; 2251 const Loop *L; 2252 2253 match_LoopInvariant(const SubPattern_t &SP, const Loop *L) 2254 : SubPattern(SP), L(L) {} 2255 2256 template <typename ITy> bool match(ITy *V) { 2257 return L->isLoopInvariant(V) && SubPattern.match(V); 2258 } 2259 }; 2260 2261 /// Matches if the value is loop-invariant. 2262 template <typename Ty> 2263 inline match_LoopInvariant<Ty> m_LoopInvariant(const Ty &M, const Loop *L) { 2264 return match_LoopInvariant<Ty>(M, L); 2265 } 2266 2267 /// Return true if the idiom is detected in the loop. 2268 /// 2269 /// The core idiom we are trying to detect is: 2270 /// \code 2271 /// entry: 2272 /// <...> 2273 /// %bitmask = shl i32 1, %bitpos 2274 /// br label %loop 2275 /// 2276 /// loop: 2277 /// %x.curr = phi i32 [ %x, %entry ], [ %x.next, %loop ] 2278 /// %x.curr.bitmasked = and i32 %x.curr, %bitmask 2279 /// %x.curr.isbitunset = icmp eq i32 %x.curr.bitmasked, 0 2280 /// %x.next = shl i32 %x.curr, 1 2281 /// <...> 2282 /// br i1 %x.curr.isbitunset, label %loop, label %end 2283 /// 2284 /// end: 2285 /// %x.curr.res = phi i32 [ %x.curr, %loop ] <...> 2286 /// %x.next.res = phi i32 [ %x.next, %loop ] <...> 2287 /// <...> 2288 /// \endcode 2289 static bool detectShiftUntilBitTestIdiom(Loop *CurLoop, Value *&BaseX, 2290 Value *&BitMask, Value *&BitPos, 2291 Value *&CurrX, Instruction *&NextX) { 2292 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2293 " Performing shift-until-bittest idiom detection.\n"); 2294 2295 // Give up if the loop has multiple blocks or multiple backedges. 2296 if (CurLoop->getNumBlocks() != 1 || CurLoop->getNumBackEdges() != 1) { 2297 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad block/backedge count.\n"); 2298 return false; 2299 } 2300 2301 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2302 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2303 assert(LoopPreheaderBB && "There is always a loop preheader."); 2304 2305 using namespace PatternMatch; 2306 2307 // Step 1: Check if the loop backedge is in desirable form. 2308 2309 ICmpInst::Predicate Pred; 2310 Value *CmpLHS, *CmpRHS; 2311 BasicBlock *TrueBB, *FalseBB; 2312 if (!match(LoopHeaderBB->getTerminator(), 2313 m_Br(m_ICmp(Pred, m_Value(CmpLHS), m_Value(CmpRHS)), 2314 m_BasicBlock(TrueBB), m_BasicBlock(FalseBB)))) { 2315 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge structure.\n"); 2316 return false; 2317 } 2318 2319 // Step 2: Check if the backedge's condition is in desirable form. 2320 2321 auto MatchVariableBitMask = [&]() { 2322 return ICmpInst::isEquality(Pred) && match(CmpRHS, m_Zero()) && 2323 match(CmpLHS, 2324 m_c_And(m_Value(CurrX), 2325 m_CombineAnd( 2326 m_Value(BitMask), 2327 m_LoopInvariant(m_Shl(m_One(), m_Value(BitPos)), 2328 CurLoop)))); 2329 }; 2330 auto MatchConstantBitMask = [&]() { 2331 return ICmpInst::isEquality(Pred) && match(CmpRHS, m_Zero()) && 2332 match(CmpLHS, m_And(m_Value(CurrX), 2333 m_CombineAnd(m_Value(BitMask), m_Power2()))) && 2334 (BitPos = ConstantExpr::getExactLogBase2(cast<Constant>(BitMask))); 2335 }; 2336 auto MatchDecomposableConstantBitMask = [&]() { 2337 APInt Mask; 2338 return llvm::decomposeBitTestICmp(CmpLHS, CmpRHS, Pred, CurrX, Mask) && 2339 ICmpInst::isEquality(Pred) && Mask.isPowerOf2() && 2340 (BitMask = ConstantInt::get(CurrX->getType(), Mask)) && 2341 (BitPos = ConstantInt::get(CurrX->getType(), Mask.logBase2())); 2342 }; 2343 2344 if (!MatchVariableBitMask() && !MatchConstantBitMask() && 2345 !MatchDecomposableConstantBitMask()) { 2346 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge comparison.\n"); 2347 return false; 2348 } 2349 2350 // Step 3: Check if the recurrence is in desirable form. 2351 auto *CurrXPN = dyn_cast<PHINode>(CurrX); 2352 if (!CurrXPN || CurrXPN->getParent() != LoopHeaderBB) { 2353 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Not an expected PHI node.\n"); 2354 return false; 2355 } 2356 2357 BaseX = CurrXPN->getIncomingValueForBlock(LoopPreheaderBB); 2358 NextX = 2359 dyn_cast<Instruction>(CurrXPN->getIncomingValueForBlock(LoopHeaderBB)); 2360 2361 assert(CurLoop->isLoopInvariant(BaseX) && 2362 "Expected BaseX to be avaliable in the preheader!"); 2363 2364 if (!NextX || !match(NextX, m_Shl(m_Specific(CurrX), m_One()))) { 2365 // FIXME: support right-shift? 2366 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad recurrence.\n"); 2367 return false; 2368 } 2369 2370 // Step 4: Check if the backedge's destinations are in desirable form. 2371 2372 assert(ICmpInst::isEquality(Pred) && 2373 "Should only get equality predicates here."); 2374 2375 // cmp-br is commutative, so canonicalize to a single variant. 2376 if (Pred != ICmpInst::Predicate::ICMP_EQ) { 2377 Pred = ICmpInst::getInversePredicate(Pred); 2378 std::swap(TrueBB, FalseBB); 2379 } 2380 2381 // We expect to exit loop when comparison yields false, 2382 // so when it yields true we should branch back to loop header. 2383 if (TrueBB != LoopHeaderBB) { 2384 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge flow.\n"); 2385 return false; 2386 } 2387 2388 // Okay, idiom checks out. 2389 return true; 2390 } 2391 2392 /// Look for the following loop: 2393 /// \code 2394 /// entry: 2395 /// <...> 2396 /// %bitmask = shl i32 1, %bitpos 2397 /// br label %loop 2398 /// 2399 /// loop: 2400 /// %x.curr = phi i32 [ %x, %entry ], [ %x.next, %loop ] 2401 /// %x.curr.bitmasked = and i32 %x.curr, %bitmask 2402 /// %x.curr.isbitunset = icmp eq i32 %x.curr.bitmasked, 0 2403 /// %x.next = shl i32 %x.curr, 1 2404 /// <...> 2405 /// br i1 %x.curr.isbitunset, label %loop, label %end 2406 /// 2407 /// end: 2408 /// %x.curr.res = phi i32 [ %x.curr, %loop ] <...> 2409 /// %x.next.res = phi i32 [ %x.next, %loop ] <...> 2410 /// <...> 2411 /// \endcode 2412 /// 2413 /// And transform it into: 2414 /// \code 2415 /// entry: 2416 /// %bitmask = shl i32 1, %bitpos 2417 /// %lowbitmask = add i32 %bitmask, -1 2418 /// %mask = or i32 %lowbitmask, %bitmask 2419 /// %x.masked = and i32 %x, %mask 2420 /// %x.masked.numleadingzeros = call i32 @llvm.ctlz.i32(i32 %x.masked, 2421 /// i1 true) 2422 /// %x.masked.numactivebits = sub i32 32, %x.masked.numleadingzeros 2423 /// %x.masked.leadingonepos = add i32 %x.masked.numactivebits, -1 2424 /// %backedgetakencount = sub i32 %bitpos, %x.masked.leadingonepos 2425 /// %tripcount = add i32 %backedgetakencount, 1 2426 /// %x.curr = shl i32 %x, %backedgetakencount 2427 /// %x.next = shl i32 %x, %tripcount 2428 /// br label %loop 2429 /// 2430 /// loop: 2431 /// %loop.iv = phi i32 [ 0, %entry ], [ %loop.iv.next, %loop ] 2432 /// %loop.iv.next = add nuw i32 %loop.iv, 1 2433 /// %loop.ivcheck = icmp eq i32 %loop.iv.next, %tripcount 2434 /// <...> 2435 /// br i1 %loop.ivcheck, label %end, label %loop 2436 /// 2437 /// end: 2438 /// %x.curr.res = phi i32 [ %x.curr, %loop ] <...> 2439 /// %x.next.res = phi i32 [ %x.next, %loop ] <...> 2440 /// <...> 2441 /// \endcode 2442 bool LoopIdiomRecognize::recognizeShiftUntilBitTest() { 2443 bool MadeChange = false; 2444 2445 Value *X, *BitMask, *BitPos, *XCurr; 2446 Instruction *XNext; 2447 if (!detectShiftUntilBitTestIdiom(CurLoop, X, BitMask, BitPos, XCurr, 2448 XNext)) { 2449 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2450 " shift-until-bittest idiom detection failed.\n"); 2451 return MadeChange; 2452 } 2453 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-bittest idiom detected!\n"); 2454 2455 // Ok, it is the idiom we were looking for, we *could* transform this loop, 2456 // but is it profitable to transform? 2457 2458 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2459 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2460 assert(LoopPreheaderBB && "There is always a loop preheader."); 2461 2462 BasicBlock *SuccessorBB = CurLoop->getExitBlock(); 2463 assert(SuccessorBB && "There is only a single successor."); 2464 2465 IRBuilder<> Builder(LoopPreheaderBB->getTerminator()); 2466 Builder.SetCurrentDebugLocation(cast<Instruction>(XCurr)->getDebugLoc()); 2467 2468 Intrinsic::ID IntrID = Intrinsic::ctlz; 2469 Type *Ty = X->getType(); 2470 unsigned Bitwidth = Ty->getScalarSizeInBits(); 2471 2472 TargetTransformInfo::TargetCostKind CostKind = 2473 TargetTransformInfo::TCK_SizeAndLatency; 2474 2475 // The rewrite is considered to be unprofitable iff and only iff the 2476 // intrinsic/shift we'll use are not cheap. Note that we are okay with *just* 2477 // making the loop countable, even if nothing else changes. 2478 IntrinsicCostAttributes Attrs( 2479 IntrID, Ty, {UndefValue::get(Ty), /*is_zero_undef=*/Builder.getTrue()}); 2480 InstructionCost Cost = TTI->getIntrinsicInstrCost(Attrs, CostKind); 2481 if (Cost > TargetTransformInfo::TCC_Basic) { 2482 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2483 " Intrinsic is too costly, not beneficial\n"); 2484 return MadeChange; 2485 } 2486 if (TTI->getArithmeticInstrCost(Instruction::Shl, Ty, CostKind) > 2487 TargetTransformInfo::TCC_Basic) { 2488 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Shift is too costly, not beneficial\n"); 2489 return MadeChange; 2490 } 2491 2492 // Ok, transform appears worthwhile. 2493 MadeChange = true; 2494 2495 // Step 1: Compute the loop trip count. 2496 2497 Value *LowBitMask = Builder.CreateAdd(BitMask, Constant::getAllOnesValue(Ty), 2498 BitPos->getName() + ".lowbitmask"); 2499 Value *Mask = 2500 Builder.CreateOr(LowBitMask, BitMask, BitPos->getName() + ".mask"); 2501 Value *XMasked = Builder.CreateAnd(X, Mask, X->getName() + ".masked"); 2502 CallInst *XMaskedNumLeadingZeros = Builder.CreateIntrinsic( 2503 IntrID, Ty, {XMasked, /*is_zero_undef=*/Builder.getTrue()}, 2504 /*FMFSource=*/nullptr, XMasked->getName() + ".numleadingzeros"); 2505 Value *XMaskedNumActiveBits = Builder.CreateSub( 2506 ConstantInt::get(Ty, Ty->getScalarSizeInBits()), XMaskedNumLeadingZeros, 2507 XMasked->getName() + ".numactivebits", /*HasNUW=*/true, 2508 /*HasNSW=*/Bitwidth != 2); 2509 Value *XMaskedLeadingOnePos = 2510 Builder.CreateAdd(XMaskedNumActiveBits, Constant::getAllOnesValue(Ty), 2511 XMasked->getName() + ".leadingonepos", /*HasNUW=*/false, 2512 /*HasNSW=*/Bitwidth > 2); 2513 2514 Value *LoopBackedgeTakenCount = Builder.CreateSub( 2515 BitPos, XMaskedLeadingOnePos, CurLoop->getName() + ".backedgetakencount", 2516 /*HasNUW=*/true, /*HasNSW=*/true); 2517 // We know loop's backedge-taken count, but what's loop's trip count? 2518 // Note that while NUW is always safe, while NSW is only for bitwidths != 2. 2519 Value *LoopTripCount = 2520 Builder.CreateAdd(LoopBackedgeTakenCount, ConstantInt::get(Ty, 1), 2521 CurLoop->getName() + ".tripcount", /*HasNUW=*/true, 2522 /*HasNSW=*/Bitwidth != 2); 2523 2524 // Step 2: Compute the recurrence's final value without a loop. 2525 2526 // NewX is always safe to compute, because `LoopBackedgeTakenCount` 2527 // will always be smaller than `bitwidth(X)`, i.e. we never get poison. 2528 Value *NewX = Builder.CreateShl(X, LoopBackedgeTakenCount); 2529 NewX->takeName(XCurr); 2530 if (auto *I = dyn_cast<Instruction>(NewX)) 2531 I->copyIRFlags(XNext, /*IncludeWrapFlags=*/true); 2532 2533 Value *NewXNext; 2534 // Rewriting XNext is more complicated, however, because `X << LoopTripCount` 2535 // will be poison iff `LoopTripCount == bitwidth(X)` (which will happen 2536 // iff `BitPos` is `bitwidth(x) - 1` and `X` is `1`). So unless we know 2537 // that isn't the case, we'll need to emit an alternative, safe IR. 2538 if (XNext->hasNoSignedWrap() || XNext->hasNoUnsignedWrap() || 2539 PatternMatch::match( 2540 BitPos, PatternMatch::m_SpecificInt_ICMP( 2541 ICmpInst::ICMP_NE, APInt(Ty->getScalarSizeInBits(), 2542 Ty->getScalarSizeInBits() - 1)))) 2543 NewXNext = Builder.CreateShl(X, LoopTripCount); 2544 else { 2545 // Otherwise, just additionally shift by one. It's the smallest solution, 2546 // alternatively, we could check that NewX is INT_MIN (or BitPos is ) 2547 // and select 0 instead. 2548 NewXNext = Builder.CreateShl(NewX, ConstantInt::get(Ty, 1)); 2549 } 2550 2551 NewXNext->takeName(XNext); 2552 if (auto *I = dyn_cast<Instruction>(NewXNext)) 2553 I->copyIRFlags(XNext, /*IncludeWrapFlags=*/true); 2554 2555 // Step 3: Adjust the successor basic block to recieve the computed 2556 // recurrence's final value instead of the recurrence itself. 2557 2558 XCurr->replaceUsesOutsideBlock(NewX, LoopHeaderBB); 2559 XNext->replaceUsesOutsideBlock(NewXNext, LoopHeaderBB); 2560 2561 // Step 4: Rewrite the loop into a countable form, with canonical IV. 2562 2563 // The new canonical induction variable. 2564 Builder.SetInsertPoint(&LoopHeaderBB->front()); 2565 auto *IV = Builder.CreatePHI(Ty, 2, CurLoop->getName() + ".iv"); 2566 2567 // The induction itself. 2568 // Note that while NUW is always safe, while NSW is only for bitwidths != 2. 2569 Builder.SetInsertPoint(LoopHeaderBB->getTerminator()); 2570 auto *IVNext = 2571 Builder.CreateAdd(IV, ConstantInt::get(Ty, 1), IV->getName() + ".next", 2572 /*HasNUW=*/true, /*HasNSW=*/Bitwidth != 2); 2573 2574 // The loop trip count check. 2575 auto *IVCheck = Builder.CreateICmpEQ(IVNext, LoopTripCount, 2576 CurLoop->getName() + ".ivcheck"); 2577 Builder.CreateCondBr(IVCheck, SuccessorBB, LoopHeaderBB); 2578 LoopHeaderBB->getTerminator()->eraseFromParent(); 2579 2580 // Populate the IV PHI. 2581 IV->addIncoming(ConstantInt::get(Ty, 0), LoopPreheaderBB); 2582 IV->addIncoming(IVNext, LoopHeaderBB); 2583 2584 // Step 5: Forget the "non-computable" trip-count SCEV associated with the 2585 // loop. The loop would otherwise not be deleted even if it becomes empty. 2586 2587 SE->forgetLoop(CurLoop); 2588 2589 // Other passes will take care of actually deleting the loop if possible. 2590 2591 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-bittest idiom optimized!\n"); 2592 2593 ++NumShiftUntilBitTest; 2594 return MadeChange; 2595 } 2596 2597 /// Return true if the idiom is detected in the loop. 2598 /// 2599 /// The core idiom we are trying to detect is: 2600 /// \code 2601 /// entry: 2602 /// <...> 2603 /// %start = <...> 2604 /// %extraoffset = <...> 2605 /// <...> 2606 /// br label %for.cond 2607 /// 2608 /// loop: 2609 /// %iv = phi i8 [ %start, %entry ], [ %iv.next, %for.cond ] 2610 /// %nbits = add nsw i8 %iv, %extraoffset 2611 /// %val.shifted = {{l,a}shr,shl} i8 %val, %nbits 2612 /// %val.shifted.iszero = icmp eq i8 %val.shifted, 0 2613 /// %iv.next = add i8 %iv, 1 2614 /// <...> 2615 /// br i1 %val.shifted.iszero, label %end, label %loop 2616 /// 2617 /// end: 2618 /// %iv.res = phi i8 [ %iv, %loop ] <...> 2619 /// %nbits.res = phi i8 [ %nbits, %loop ] <...> 2620 /// %val.shifted.res = phi i8 [ %val.shifted, %loop ] <...> 2621 /// %val.shifted.iszero.res = phi i1 [ %val.shifted.iszero, %loop ] <...> 2622 /// %iv.next.res = phi i8 [ %iv.next, %loop ] <...> 2623 /// <...> 2624 /// \endcode 2625 static bool detectShiftUntilZeroIdiom(Loop *CurLoop, ScalarEvolution *SE, 2626 Instruction *&ValShiftedIsZero, 2627 Intrinsic::ID &IntrinID, Instruction *&IV, 2628 Value *&Start, Value *&Val, 2629 const SCEV *&ExtraOffsetExpr, 2630 bool &InvertedCond) { 2631 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2632 " Performing shift-until-zero idiom detection.\n"); 2633 2634 // Give up if the loop has multiple blocks or multiple backedges. 2635 if (CurLoop->getNumBlocks() != 1 || CurLoop->getNumBackEdges() != 1) { 2636 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad block/backedge count.\n"); 2637 return false; 2638 } 2639 2640 Instruction *ValShifted, *NBits, *IVNext; 2641 Value *ExtraOffset; 2642 2643 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2644 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2645 assert(LoopPreheaderBB && "There is always a loop preheader."); 2646 2647 using namespace PatternMatch; 2648 2649 // Step 1: Check if the loop backedge, condition is in desirable form. 2650 2651 ICmpInst::Predicate Pred; 2652 BasicBlock *TrueBB, *FalseBB; 2653 if (!match(LoopHeaderBB->getTerminator(), 2654 m_Br(m_Instruction(ValShiftedIsZero), m_BasicBlock(TrueBB), 2655 m_BasicBlock(FalseBB))) || 2656 !match(ValShiftedIsZero, 2657 m_ICmp(Pred, m_Instruction(ValShifted), m_Zero())) || 2658 !ICmpInst::isEquality(Pred)) { 2659 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge structure.\n"); 2660 return false; 2661 } 2662 2663 // Step 2: Check if the comparison's operand is in desirable form. 2664 // FIXME: Val could be a one-input PHI node, which we should look past. 2665 if (!match(ValShifted, m_Shift(m_LoopInvariant(m_Value(Val), CurLoop), 2666 m_Instruction(NBits)))) { 2667 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad comparisons value computation.\n"); 2668 return false; 2669 } 2670 IntrinID = ValShifted->getOpcode() == Instruction::Shl ? Intrinsic::cttz 2671 : Intrinsic::ctlz; 2672 2673 // Step 3: Check if the shift amount is in desirable form. 2674 2675 if (match(NBits, m_c_Add(m_Instruction(IV), 2676 m_LoopInvariant(m_Value(ExtraOffset), CurLoop))) && 2677 (NBits->hasNoSignedWrap() || NBits->hasNoUnsignedWrap())) 2678 ExtraOffsetExpr = SE->getNegativeSCEV(SE->getSCEV(ExtraOffset)); 2679 else if (match(NBits, 2680 m_Sub(m_Instruction(IV), 2681 m_LoopInvariant(m_Value(ExtraOffset), CurLoop))) && 2682 NBits->hasNoSignedWrap()) 2683 ExtraOffsetExpr = SE->getSCEV(ExtraOffset); 2684 else { 2685 IV = NBits; 2686 ExtraOffsetExpr = SE->getZero(NBits->getType()); 2687 } 2688 2689 // Step 4: Check if the recurrence is in desirable form. 2690 auto *IVPN = dyn_cast<PHINode>(IV); 2691 if (!IVPN || IVPN->getParent() != LoopHeaderBB) { 2692 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Not an expected PHI node.\n"); 2693 return false; 2694 } 2695 2696 Start = IVPN->getIncomingValueForBlock(LoopPreheaderBB); 2697 IVNext = dyn_cast<Instruction>(IVPN->getIncomingValueForBlock(LoopHeaderBB)); 2698 2699 if (!IVNext || !match(IVNext, m_Add(m_Specific(IVPN), m_One()))) { 2700 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad recurrence.\n"); 2701 return false; 2702 } 2703 2704 // Step 4: Check if the backedge's destinations are in desirable form. 2705 2706 assert(ICmpInst::isEquality(Pred) && 2707 "Should only get equality predicates here."); 2708 2709 // cmp-br is commutative, so canonicalize to a single variant. 2710 InvertedCond = Pred != ICmpInst::Predicate::ICMP_EQ; 2711 if (InvertedCond) { 2712 Pred = ICmpInst::getInversePredicate(Pred); 2713 std::swap(TrueBB, FalseBB); 2714 } 2715 2716 // We expect to exit loop when comparison yields true, 2717 // so when it yields false we should branch back to loop header. 2718 if (FalseBB != LoopHeaderBB) { 2719 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge flow.\n"); 2720 return false; 2721 } 2722 2723 // The new, countable, loop will certainly only run a known number of 2724 // iterations, It won't be infinite. But the old loop might be infinite 2725 // under certain conditions. For logical shifts, the value will become zero 2726 // after at most bitwidth(%Val) loop iterations. However, for arithmetic 2727 // right-shift, iff the sign bit was set, the value will never become zero, 2728 // and the loop may never finish. 2729 if (ValShifted->getOpcode() == Instruction::AShr && 2730 !isMustProgress(CurLoop) && !SE->isKnownNonNegative(SE->getSCEV(Val))) { 2731 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Can not prove the loop is finite.\n"); 2732 return false; 2733 } 2734 2735 // Okay, idiom checks out. 2736 return true; 2737 } 2738 2739 /// Look for the following loop: 2740 /// \code 2741 /// entry: 2742 /// <...> 2743 /// %start = <...> 2744 /// %extraoffset = <...> 2745 /// <...> 2746 /// br label %for.cond 2747 /// 2748 /// loop: 2749 /// %iv = phi i8 [ %start, %entry ], [ %iv.next, %for.cond ] 2750 /// %nbits = add nsw i8 %iv, %extraoffset 2751 /// %val.shifted = {{l,a}shr,shl} i8 %val, %nbits 2752 /// %val.shifted.iszero = icmp eq i8 %val.shifted, 0 2753 /// %iv.next = add i8 %iv, 1 2754 /// <...> 2755 /// br i1 %val.shifted.iszero, label %end, label %loop 2756 /// 2757 /// end: 2758 /// %iv.res = phi i8 [ %iv, %loop ] <...> 2759 /// %nbits.res = phi i8 [ %nbits, %loop ] <...> 2760 /// %val.shifted.res = phi i8 [ %val.shifted, %loop ] <...> 2761 /// %val.shifted.iszero.res = phi i1 [ %val.shifted.iszero, %loop ] <...> 2762 /// %iv.next.res = phi i8 [ %iv.next, %loop ] <...> 2763 /// <...> 2764 /// \endcode 2765 /// 2766 /// And transform it into: 2767 /// \code 2768 /// entry: 2769 /// <...> 2770 /// %start = <...> 2771 /// %extraoffset = <...> 2772 /// <...> 2773 /// %val.numleadingzeros = call i8 @llvm.ct{l,t}z.i8(i8 %val, i1 0) 2774 /// %val.numactivebits = sub i8 8, %val.numleadingzeros 2775 /// %extraoffset.neg = sub i8 0, %extraoffset 2776 /// %tmp = add i8 %val.numactivebits, %extraoffset.neg 2777 /// %iv.final = call i8 @llvm.smax.i8(i8 %tmp, i8 %start) 2778 /// %loop.tripcount = sub i8 %iv.final, %start 2779 /// br label %loop 2780 /// 2781 /// loop: 2782 /// %loop.iv = phi i8 [ 0, %entry ], [ %loop.iv.next, %loop ] 2783 /// %loop.iv.next = add i8 %loop.iv, 1 2784 /// %loop.ivcheck = icmp eq i8 %loop.iv.next, %loop.tripcount 2785 /// %iv = add i8 %loop.iv, %start 2786 /// <...> 2787 /// br i1 %loop.ivcheck, label %end, label %loop 2788 /// 2789 /// end: 2790 /// %iv.res = phi i8 [ %iv.final, %loop ] <...> 2791 /// <...> 2792 /// \endcode 2793 bool LoopIdiomRecognize::recognizeShiftUntilZero() { 2794 bool MadeChange = false; 2795 2796 Instruction *ValShiftedIsZero; 2797 Intrinsic::ID IntrID; 2798 Instruction *IV; 2799 Value *Start, *Val; 2800 const SCEV *ExtraOffsetExpr; 2801 bool InvertedCond; 2802 if (!detectShiftUntilZeroIdiom(CurLoop, SE, ValShiftedIsZero, IntrID, IV, 2803 Start, Val, ExtraOffsetExpr, InvertedCond)) { 2804 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2805 " shift-until-zero idiom detection failed.\n"); 2806 return MadeChange; 2807 } 2808 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-zero idiom detected!\n"); 2809 2810 // Ok, it is the idiom we were looking for, we *could* transform this loop, 2811 // but is it profitable to transform? 2812 2813 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2814 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2815 assert(LoopPreheaderBB && "There is always a loop preheader."); 2816 2817 BasicBlock *SuccessorBB = CurLoop->getExitBlock(); 2818 assert(SuccessorBB && "There is only a single successor."); 2819 2820 IRBuilder<> Builder(LoopPreheaderBB->getTerminator()); 2821 Builder.SetCurrentDebugLocation(IV->getDebugLoc()); 2822 2823 Type *Ty = Val->getType(); 2824 unsigned Bitwidth = Ty->getScalarSizeInBits(); 2825 2826 TargetTransformInfo::TargetCostKind CostKind = 2827 TargetTransformInfo::TCK_SizeAndLatency; 2828 2829 // The rewrite is considered to be unprofitable iff and only iff the 2830 // intrinsic we'll use are not cheap. Note that we are okay with *just* 2831 // making the loop countable, even if nothing else changes. 2832 IntrinsicCostAttributes Attrs( 2833 IntrID, Ty, {UndefValue::get(Ty), /*is_zero_undef=*/Builder.getFalse()}); 2834 InstructionCost Cost = TTI->getIntrinsicInstrCost(Attrs, CostKind); 2835 if (Cost > TargetTransformInfo::TCC_Basic) { 2836 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2837 " Intrinsic is too costly, not beneficial\n"); 2838 return MadeChange; 2839 } 2840 2841 // Ok, transform appears worthwhile. 2842 MadeChange = true; 2843 2844 bool OffsetIsZero = false; 2845 if (auto *ExtraOffsetExprC = dyn_cast<SCEVConstant>(ExtraOffsetExpr)) 2846 OffsetIsZero = ExtraOffsetExprC->isZero(); 2847 2848 // Step 1: Compute the loop's final IV value / trip count. 2849 2850 CallInst *ValNumLeadingZeros = Builder.CreateIntrinsic( 2851 IntrID, Ty, {Val, /*is_zero_undef=*/Builder.getFalse()}, 2852 /*FMFSource=*/nullptr, Val->getName() + ".numleadingzeros"); 2853 Value *ValNumActiveBits = Builder.CreateSub( 2854 ConstantInt::get(Ty, Ty->getScalarSizeInBits()), ValNumLeadingZeros, 2855 Val->getName() + ".numactivebits", /*HasNUW=*/true, 2856 /*HasNSW=*/Bitwidth != 2); 2857 2858 SCEVExpander Expander(*SE, *DL, "loop-idiom"); 2859 Expander.setInsertPoint(&*Builder.GetInsertPoint()); 2860 Value *ExtraOffset = Expander.expandCodeFor(ExtraOffsetExpr); 2861 2862 Value *ValNumActiveBitsOffset = Builder.CreateAdd( 2863 ValNumActiveBits, ExtraOffset, ValNumActiveBits->getName() + ".offset", 2864 /*HasNUW=*/OffsetIsZero, /*HasNSW=*/true); 2865 Value *IVFinal = Builder.CreateIntrinsic(Intrinsic::smax, {Ty}, 2866 {ValNumActiveBitsOffset, Start}, 2867 /*FMFSource=*/nullptr, "iv.final"); 2868 2869 auto *LoopBackedgeTakenCount = cast<Instruction>(Builder.CreateSub( 2870 IVFinal, Start, CurLoop->getName() + ".backedgetakencount", 2871 /*HasNUW=*/OffsetIsZero, /*HasNSW=*/true)); 2872 // FIXME: or when the offset was `add nuw` 2873 2874 // We know loop's backedge-taken count, but what's loop's trip count? 2875 Value *LoopTripCount = 2876 Builder.CreateAdd(LoopBackedgeTakenCount, ConstantInt::get(Ty, 1), 2877 CurLoop->getName() + ".tripcount", /*HasNUW=*/true, 2878 /*HasNSW=*/Bitwidth != 2); 2879 2880 // Step 2: Adjust the successor basic block to recieve the original 2881 // induction variable's final value instead of the orig. IV itself. 2882 2883 IV->replaceUsesOutsideBlock(IVFinal, LoopHeaderBB); 2884 2885 // Step 3: Rewrite the loop into a countable form, with canonical IV. 2886 2887 // The new canonical induction variable. 2888 Builder.SetInsertPoint(&LoopHeaderBB->front()); 2889 auto *CIV = Builder.CreatePHI(Ty, 2, CurLoop->getName() + ".iv"); 2890 2891 // The induction itself. 2892 Builder.SetInsertPoint(LoopHeaderBB->getFirstNonPHI()); 2893 auto *CIVNext = 2894 Builder.CreateAdd(CIV, ConstantInt::get(Ty, 1), CIV->getName() + ".next", 2895 /*HasNUW=*/true, /*HasNSW=*/Bitwidth != 2); 2896 2897 // The loop trip count check. 2898 auto *CIVCheck = Builder.CreateICmpEQ(CIVNext, LoopTripCount, 2899 CurLoop->getName() + ".ivcheck"); 2900 auto *NewIVCheck = CIVCheck; 2901 if (InvertedCond) { 2902 NewIVCheck = Builder.CreateNot(CIVCheck); 2903 NewIVCheck->takeName(ValShiftedIsZero); 2904 } 2905 2906 // The original IV, but rebased to be an offset to the CIV. 2907 auto *IVDePHId = Builder.CreateAdd(CIV, Start, "", /*HasNUW=*/false, 2908 /*HasNSW=*/true); // FIXME: what about NUW? 2909 IVDePHId->takeName(IV); 2910 2911 // The loop terminator. 2912 Builder.SetInsertPoint(LoopHeaderBB->getTerminator()); 2913 Builder.CreateCondBr(CIVCheck, SuccessorBB, LoopHeaderBB); 2914 LoopHeaderBB->getTerminator()->eraseFromParent(); 2915 2916 // Populate the IV PHI. 2917 CIV->addIncoming(ConstantInt::get(Ty, 0), LoopPreheaderBB); 2918 CIV->addIncoming(CIVNext, LoopHeaderBB); 2919 2920 // Step 4: Forget the "non-computable" trip-count SCEV associated with the 2921 // loop. The loop would otherwise not be deleted even if it becomes empty. 2922 2923 SE->forgetLoop(CurLoop); 2924 2925 // Step 5: Try to cleanup the loop's body somewhat. 2926 IV->replaceAllUsesWith(IVDePHId); 2927 IV->eraseFromParent(); 2928 2929 ValShiftedIsZero->replaceAllUsesWith(NewIVCheck); 2930 ValShiftedIsZero->eraseFromParent(); 2931 2932 // Other passes will take care of actually deleting the loop if possible. 2933 2934 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-zero idiom optimized!\n"); 2935 2936 ++NumShiftUntilZero; 2937 return MadeChange; 2938 } 2939