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 (unsigned i = 0, e = ExitBlocks.size(); i != e; ++i) 629 if (!DT->dominates(BB, ExitBlocks[i])) 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 (SetVector<StoreInst *>::iterator it = Heads.begin(), e = Heads.end(); 754 it != e; ++it) { 755 if (Tails.count(*it)) 756 continue; 757 758 // We found a store instr that starts a chain. Now follow the chain and try 759 // to transform it. 760 SmallPtrSet<Instruction *, 8> AdjacentStores; 761 StoreInst *I = *it; 762 763 StoreInst *HeadStore = I; 764 unsigned StoreSize = 0; 765 766 // Collect the chain into a list. 767 while (Tails.count(I) || Heads.count(I)) { 768 if (TransformedStores.count(I)) 769 break; 770 AdjacentStores.insert(I); 771 772 StoreSize += DL->getTypeStoreSize(I->getValueOperand()->getType()); 773 // Move to the next value in the chain. 774 I = ConsecutiveChain[I]; 775 } 776 777 Value *StoredVal = HeadStore->getValueOperand(); 778 Value *StorePtr = HeadStore->getPointerOperand(); 779 const SCEVAddRecExpr *StoreEv = cast<SCEVAddRecExpr>(SE->getSCEV(StorePtr)); 780 APInt Stride = getStoreStride(StoreEv); 781 782 // Check to see if the stride matches the size of the stores. If so, then 783 // we know that every byte is touched in the loop. 784 if (StoreSize != Stride && StoreSize != -Stride) 785 continue; 786 787 bool IsNegStride = StoreSize == -Stride; 788 789 const SCEV *StoreSizeSCEV = SE->getConstant(BECount->getType(), StoreSize); 790 if (processLoopStridedStore(StorePtr, StoreSizeSCEV, 791 MaybeAlign(HeadStore->getAlignment()), 792 StoredVal, HeadStore, AdjacentStores, StoreEv, 793 BECount, IsNegStride)) { 794 TransformedStores.insert(AdjacentStores.begin(), AdjacentStores.end()); 795 Changed = true; 796 } 797 } 798 799 return Changed; 800 } 801 802 /// processLoopMemIntrinsic - Template function for calling different processor 803 /// functions based on mem instrinsic type. 804 template <typename MemInst> 805 bool LoopIdiomRecognize::processLoopMemIntrinsic( 806 BasicBlock *BB, 807 bool (LoopIdiomRecognize::*Processor)(MemInst *, const SCEV *), 808 const SCEV *BECount) { 809 bool MadeChange = false; 810 for (BasicBlock::iterator I = BB->begin(), E = BB->end(); I != E;) { 811 Instruction *Inst = &*I++; 812 // Look for memory instructions, which may be optimized to a larger one. 813 if (MemInst *MI = dyn_cast<MemInst>(Inst)) { 814 WeakTrackingVH InstPtr(&*I); 815 if (!(this->*Processor)(MI, BECount)) 816 continue; 817 MadeChange = true; 818 819 // If processing the instruction invalidated our iterator, start over from 820 // the top of the block. 821 if (!InstPtr) 822 I = BB->begin(); 823 } 824 } 825 return MadeChange; 826 } 827 828 /// processLoopMemCpy - See if this memcpy can be promoted to a large memcpy 829 bool LoopIdiomRecognize::processLoopMemCpy(MemCpyInst *MCI, 830 const SCEV *BECount) { 831 // We can only handle non-volatile memcpys with a constant size. 832 if (MCI->isVolatile() || !isa<ConstantInt>(MCI->getLength())) 833 return false; 834 835 // If we're not allowed to hack on memcpy, we fail. 836 if ((!HasMemcpy && !isa<MemCpyInlineInst>(MCI)) || DisableLIRP::Memcpy) 837 return false; 838 839 Value *Dest = MCI->getDest(); 840 Value *Source = MCI->getSource(); 841 if (!Dest || !Source) 842 return false; 843 844 // See if the load and store pointer expressions are AddRec like {base,+,1} on 845 // the current loop, which indicates a strided load and store. If we have 846 // something else, it's a random load or store we can't handle. 847 const SCEVAddRecExpr *StoreEv = dyn_cast<SCEVAddRecExpr>(SE->getSCEV(Dest)); 848 if (!StoreEv || StoreEv->getLoop() != CurLoop || !StoreEv->isAffine()) 849 return false; 850 const SCEVAddRecExpr *LoadEv = dyn_cast<SCEVAddRecExpr>(SE->getSCEV(Source)); 851 if (!LoadEv || LoadEv->getLoop() != CurLoop || !LoadEv->isAffine()) 852 return false; 853 854 // Reject memcpys that are so large that they overflow an unsigned. 855 uint64_t SizeInBytes = cast<ConstantInt>(MCI->getLength())->getZExtValue(); 856 if ((SizeInBytes >> 32) != 0) 857 return false; 858 859 // Check if the stride matches the size of the memcpy. If so, then we know 860 // that every byte is touched in the loop. 861 const SCEVConstant *ConstStoreStride = 862 dyn_cast<SCEVConstant>(StoreEv->getOperand(1)); 863 const SCEVConstant *ConstLoadStride = 864 dyn_cast<SCEVConstant>(LoadEv->getOperand(1)); 865 if (!ConstStoreStride || !ConstLoadStride) 866 return false; 867 868 APInt StoreStrideValue = ConstStoreStride->getAPInt(); 869 APInt LoadStrideValue = ConstLoadStride->getAPInt(); 870 // Huge stride value - give up 871 if (StoreStrideValue.getBitWidth() > 64 || LoadStrideValue.getBitWidth() > 64) 872 return false; 873 874 if (SizeInBytes != StoreStrideValue && SizeInBytes != -StoreStrideValue) { 875 ORE.emit([&]() { 876 return OptimizationRemarkMissed(DEBUG_TYPE, "SizeStrideUnequal", MCI) 877 << ore::NV("Inst", "memcpy") << " in " 878 << ore::NV("Function", MCI->getFunction()) 879 << " function will not be hoisted: " 880 << ore::NV("Reason", "memcpy size is not equal to stride"); 881 }); 882 return false; 883 } 884 885 int64_t StoreStrideInt = StoreStrideValue.getSExtValue(); 886 int64_t LoadStrideInt = LoadStrideValue.getSExtValue(); 887 // Check if the load stride matches the store stride. 888 if (StoreStrideInt != LoadStrideInt) 889 return false; 890 891 return processLoopStoreOfLoopLoad( 892 Dest, Source, SE->getConstant(Dest->getType(), SizeInBytes), 893 MCI->getDestAlign(), MCI->getSourceAlign(), MCI, MCI, StoreEv, LoadEv, 894 BECount); 895 } 896 897 /// processLoopMemSet - See if this memset can be promoted to a large memset. 898 bool LoopIdiomRecognize::processLoopMemSet(MemSetInst *MSI, 899 const SCEV *BECount) { 900 // We can only handle non-volatile memsets. 901 if (MSI->isVolatile()) 902 return false; 903 904 // If we're not allowed to hack on memset, we fail. 905 if (!HasMemset || DisableLIRP::Memset) 906 return false; 907 908 Value *Pointer = MSI->getDest(); 909 910 // See if the pointer expression is an AddRec like {base,+,1} on the current 911 // loop, which indicates a strided store. If we have something else, it's a 912 // random store we can't handle. 913 const SCEVAddRecExpr *Ev = dyn_cast<SCEVAddRecExpr>(SE->getSCEV(Pointer)); 914 if (!Ev || Ev->getLoop() != CurLoop) 915 return false; 916 if (!Ev->isAffine()) { 917 LLVM_DEBUG(dbgs() << " Pointer is not affine, abort\n"); 918 return false; 919 } 920 921 const SCEV *PointerStrideSCEV = Ev->getOperand(1); 922 const SCEV *MemsetSizeSCEV = SE->getSCEV(MSI->getLength()); 923 if (!PointerStrideSCEV || !MemsetSizeSCEV) 924 return false; 925 926 bool IsNegStride = false; 927 const bool IsConstantSize = isa<ConstantInt>(MSI->getLength()); 928 929 if (IsConstantSize) { 930 // Memset size is constant. 931 // Check if the pointer stride matches the memset size. If so, then 932 // we know that every byte is touched in the loop. 933 LLVM_DEBUG(dbgs() << " memset size is constant\n"); 934 uint64_t SizeInBytes = cast<ConstantInt>(MSI->getLength())->getZExtValue(); 935 const SCEVConstant *ConstStride = dyn_cast<SCEVConstant>(Ev->getOperand(1)); 936 if (!ConstStride) 937 return false; 938 939 APInt Stride = ConstStride->getAPInt(); 940 if (SizeInBytes != Stride && SizeInBytes != -Stride) 941 return false; 942 943 IsNegStride = SizeInBytes == -Stride; 944 } else { 945 // Memset size is non-constant. 946 // Check if the pointer stride matches the memset size. 947 // To be conservative, the pass would not promote pointers that aren't in 948 // address space zero. Also, the pass only handles memset length and stride 949 // that are invariant for the top level loop. 950 LLVM_DEBUG(dbgs() << " memset size is non-constant\n"); 951 if (Pointer->getType()->getPointerAddressSpace() != 0) { 952 LLVM_DEBUG(dbgs() << " pointer is not in address space zero, " 953 << "abort\n"); 954 return false; 955 } 956 if (!SE->isLoopInvariant(MemsetSizeSCEV, CurLoop)) { 957 LLVM_DEBUG(dbgs() << " memset size is not a loop-invariant, " 958 << "abort\n"); 959 return false; 960 } 961 962 // Compare positive direction PointerStrideSCEV with MemsetSizeSCEV 963 IsNegStride = PointerStrideSCEV->isNonConstantNegative(); 964 const SCEV *PositiveStrideSCEV = 965 IsNegStride ? SE->getNegativeSCEV(PointerStrideSCEV) 966 : PointerStrideSCEV; 967 LLVM_DEBUG(dbgs() << " MemsetSizeSCEV: " << *MemsetSizeSCEV << "\n" 968 << " PositiveStrideSCEV: " << *PositiveStrideSCEV 969 << "\n"); 970 971 if (PositiveStrideSCEV != MemsetSizeSCEV) { 972 // TODO: folding can be done to the SCEVs 973 // The folding is to fold expressions that is covered by the loop guard 974 // at loop entry. After the folding, compare again and proceed 975 // optimization if equal. 976 LLVM_DEBUG(dbgs() << " SCEV don't match, abort\n"); 977 return false; 978 } 979 } 980 981 // Verify that the memset value is loop invariant. If not, we can't promote 982 // the memset. 983 Value *SplatValue = MSI->getValue(); 984 if (!SplatValue || !CurLoop->isLoopInvariant(SplatValue)) 985 return false; 986 987 SmallPtrSet<Instruction *, 1> MSIs; 988 MSIs.insert(MSI); 989 return processLoopStridedStore(Pointer, SE->getSCEV(MSI->getLength()), 990 MaybeAlign(MSI->getDestAlignment()), 991 SplatValue, MSI, MSIs, Ev, BECount, 992 IsNegStride, /*IsLoopMemset=*/true); 993 } 994 995 /// mayLoopAccessLocation - Return true if the specified loop might access the 996 /// specified pointer location, which is a loop-strided access. The 'Access' 997 /// argument specifies what the verboten forms of access are (read or write). 998 static bool 999 mayLoopAccessLocation(Value *Ptr, ModRefInfo Access, Loop *L, 1000 const SCEV *BECount, const SCEV *StoreSizeSCEV, 1001 AliasAnalysis &AA, 1002 SmallPtrSetImpl<Instruction *> &IgnoredInsts) { 1003 // Get the location that may be stored across the loop. Since the access is 1004 // strided positively through memory, we say that the modified location starts 1005 // at the pointer and has infinite size. 1006 LocationSize AccessSize = LocationSize::afterPointer(); 1007 1008 // If the loop iterates a fixed number of times, we can refine the access size 1009 // to be exactly the size of the memset, which is (BECount+1)*StoreSize 1010 const SCEVConstant *BECst = dyn_cast<SCEVConstant>(BECount); 1011 const SCEVConstant *ConstSize = dyn_cast<SCEVConstant>(StoreSizeSCEV); 1012 if (BECst && ConstSize) 1013 AccessSize = LocationSize::precise((BECst->getValue()->getZExtValue() + 1) * 1014 ConstSize->getValue()->getZExtValue()); 1015 1016 // TODO: For this to be really effective, we have to dive into the pointer 1017 // operand in the store. Store to &A[i] of 100 will always return may alias 1018 // with store of &A[100], we need to StoreLoc to be "A" with size of 100, 1019 // which will then no-alias a store to &A[100]. 1020 MemoryLocation StoreLoc(Ptr, AccessSize); 1021 1022 for (Loop::block_iterator BI = L->block_begin(), E = L->block_end(); BI != E; 1023 ++BI) 1024 for (Instruction &I : **BI) 1025 if (IgnoredInsts.count(&I) == 0 && 1026 isModOrRefSet( 1027 intersectModRef(AA.getModRefInfo(&I, StoreLoc), Access))) 1028 return true; 1029 return false; 1030 } 1031 1032 // If we have a negative stride, Start refers to the end of the memory location 1033 // we're trying to memset. Therefore, we need to recompute the base pointer, 1034 // which is just Start - BECount*Size. 1035 static const SCEV *getStartForNegStride(const SCEV *Start, const SCEV *BECount, 1036 Type *IntPtr, const SCEV *StoreSizeSCEV, 1037 ScalarEvolution *SE) { 1038 const SCEV *Index = SE->getTruncateOrZeroExtend(BECount, IntPtr); 1039 if (!StoreSizeSCEV->isOne()) { 1040 // index = back edge count * store size 1041 Index = SE->getMulExpr(Index, 1042 SE->getTruncateOrZeroExtend(StoreSizeSCEV, IntPtr), 1043 SCEV::FlagNUW); 1044 } 1045 // base pointer = start - index * store size 1046 return SE->getMinusSCEV(Start, Index); 1047 } 1048 1049 /// Compute trip count from the backedge taken count. 1050 static const SCEV *getTripCount(const SCEV *BECount, Type *IntPtr, 1051 Loop *CurLoop, const DataLayout *DL, 1052 ScalarEvolution *SE) { 1053 const SCEV *TripCountS = nullptr; 1054 // The # stored bytes is (BECount+1). Expand the trip count out to 1055 // pointer size if it isn't already. 1056 // 1057 // If we're going to need to zero extend the BE count, check if we can add 1058 // one to it prior to zero extending without overflow. Provided this is safe, 1059 // it allows better simplification of the +1. 1060 if (DL->getTypeSizeInBits(BECount->getType()) < 1061 DL->getTypeSizeInBits(IntPtr) && 1062 SE->isLoopEntryGuardedByCond( 1063 CurLoop, ICmpInst::ICMP_NE, BECount, 1064 SE->getNegativeSCEV(SE->getOne(BECount->getType())))) { 1065 TripCountS = SE->getZeroExtendExpr( 1066 SE->getAddExpr(BECount, SE->getOne(BECount->getType()), SCEV::FlagNUW), 1067 IntPtr); 1068 } else { 1069 TripCountS = SE->getAddExpr(SE->getTruncateOrZeroExtend(BECount, IntPtr), 1070 SE->getOne(IntPtr), SCEV::FlagNUW); 1071 } 1072 1073 return TripCountS; 1074 } 1075 1076 /// Compute the number of bytes as a SCEV from the backedge taken count. 1077 /// 1078 /// This also maps the SCEV into the provided type and tries to handle the 1079 /// computation in a way that will fold cleanly. 1080 static const SCEV *getNumBytes(const SCEV *BECount, Type *IntPtr, 1081 const SCEV *StoreSizeSCEV, Loop *CurLoop, 1082 const DataLayout *DL, ScalarEvolution *SE) { 1083 const SCEV *TripCountSCEV = getTripCount(BECount, IntPtr, CurLoop, DL, SE); 1084 1085 return SE->getMulExpr(TripCountSCEV, 1086 SE->getTruncateOrZeroExtend(StoreSizeSCEV, IntPtr), 1087 SCEV::FlagNUW); 1088 } 1089 1090 /// processLoopStridedStore - We see a strided store of some value. If we can 1091 /// transform this into a memset or memset_pattern in the loop preheader, do so. 1092 bool LoopIdiomRecognize::processLoopStridedStore( 1093 Value *DestPtr, const SCEV *StoreSizeSCEV, MaybeAlign StoreAlignment, 1094 Value *StoredVal, Instruction *TheStore, 1095 SmallPtrSetImpl<Instruction *> &Stores, const SCEVAddRecExpr *Ev, 1096 const SCEV *BECount, bool IsNegStride, bool IsLoopMemset) { 1097 Value *SplatValue = isBytewiseValue(StoredVal, *DL); 1098 Constant *PatternValue = nullptr; 1099 1100 if (!SplatValue) 1101 PatternValue = getMemSetPatternValue(StoredVal, DL); 1102 1103 assert((SplatValue || PatternValue) && 1104 "Expected either splat value or pattern value."); 1105 1106 // The trip count of the loop and the base pointer of the addrec SCEV is 1107 // guaranteed to be loop invariant, which means that it should dominate the 1108 // header. This allows us to insert code for it in the preheader. 1109 unsigned DestAS = DestPtr->getType()->getPointerAddressSpace(); 1110 BasicBlock *Preheader = CurLoop->getLoopPreheader(); 1111 IRBuilder<> Builder(Preheader->getTerminator()); 1112 SCEVExpander Expander(*SE, *DL, "loop-idiom"); 1113 SCEVExpanderCleaner ExpCleaner(Expander, *DT); 1114 1115 Type *DestInt8PtrTy = Builder.getInt8PtrTy(DestAS); 1116 Type *IntIdxTy = DL->getIndexType(DestPtr->getType()); 1117 1118 bool Changed = false; 1119 const SCEV *Start = Ev->getStart(); 1120 // Handle negative strided loops. 1121 if (IsNegStride) 1122 Start = getStartForNegStride(Start, BECount, IntIdxTy, StoreSizeSCEV, SE); 1123 1124 // TODO: ideally we should still be able to generate memset if SCEV expander 1125 // is taught to generate the dependencies at the latest point. 1126 if (!isSafeToExpand(Start, *SE)) 1127 return Changed; 1128 1129 // Okay, we have a strided store "p[i]" of a splattable value. We can turn 1130 // this into a memset in the loop preheader now if we want. However, this 1131 // would be unsafe to do if there is anything else in the loop that may read 1132 // or write to the aliased location. Check for any overlap by generating the 1133 // base pointer and checking the region. 1134 Value *BasePtr = 1135 Expander.expandCodeFor(Start, DestInt8PtrTy, Preheader->getTerminator()); 1136 1137 // From here on out, conservatively report to the pass manager that we've 1138 // changed the IR, even if we later clean up these added instructions. There 1139 // may be structural differences e.g. in the order of use lists not accounted 1140 // for in just a textual dump of the IR. This is written as a variable, even 1141 // though statically all the places this dominates could be replaced with 1142 // 'true', with the hope that anyone trying to be clever / "more precise" with 1143 // the return value will read this comment, and leave them alone. 1144 Changed = true; 1145 1146 if (mayLoopAccessLocation(BasePtr, ModRefInfo::ModRef, CurLoop, BECount, 1147 StoreSizeSCEV, *AA, Stores)) 1148 return Changed; 1149 1150 if (avoidLIRForMultiBlockLoop(/*IsMemset=*/true, IsLoopMemset)) 1151 return Changed; 1152 1153 // Okay, everything looks good, insert the memset. 1154 1155 const SCEV *NumBytesS = 1156 getNumBytes(BECount, IntIdxTy, StoreSizeSCEV, CurLoop, DL, SE); 1157 1158 // TODO: ideally we should still be able to generate memset if SCEV expander 1159 // is taught to generate the dependencies at the latest point. 1160 if (!isSafeToExpand(NumBytesS, *SE)) 1161 return Changed; 1162 1163 Value *NumBytes = 1164 Expander.expandCodeFor(NumBytesS, IntIdxTy, Preheader->getTerminator()); 1165 1166 CallInst *NewCall; 1167 if (SplatValue) { 1168 NewCall = Builder.CreateMemSet(BasePtr, SplatValue, NumBytes, 1169 MaybeAlign(StoreAlignment)); 1170 } else { 1171 // Everything is emitted in default address space 1172 Type *Int8PtrTy = DestInt8PtrTy; 1173 1174 Module *M = TheStore->getModule(); 1175 StringRef FuncName = "memset_pattern16"; 1176 FunctionCallee MSP = M->getOrInsertFunction(FuncName, Builder.getVoidTy(), 1177 Int8PtrTy, Int8PtrTy, IntIdxTy); 1178 inferLibFuncAttributes(M, FuncName, *TLI); 1179 1180 // Otherwise we should form a memset_pattern16. PatternValue is known to be 1181 // an constant array of 16-bytes. Plop the value into a mergable global. 1182 GlobalVariable *GV = new GlobalVariable(*M, PatternValue->getType(), true, 1183 GlobalValue::PrivateLinkage, 1184 PatternValue, ".memset_pattern"); 1185 GV->setUnnamedAddr(GlobalValue::UnnamedAddr::Global); // Ok to merge these. 1186 GV->setAlignment(Align(16)); 1187 Value *PatternPtr = ConstantExpr::getBitCast(GV, Int8PtrTy); 1188 NewCall = Builder.CreateCall(MSP, {BasePtr, PatternPtr, NumBytes}); 1189 } 1190 NewCall->setDebugLoc(TheStore->getDebugLoc()); 1191 1192 if (MSSAU) { 1193 MemoryAccess *NewMemAcc = MSSAU->createMemoryAccessInBB( 1194 NewCall, nullptr, NewCall->getParent(), MemorySSA::BeforeTerminator); 1195 MSSAU->insertDef(cast<MemoryDef>(NewMemAcc), true); 1196 } 1197 1198 LLVM_DEBUG(dbgs() << " Formed memset: " << *NewCall << "\n" 1199 << " from store to: " << *Ev << " at: " << *TheStore 1200 << "\n"); 1201 1202 ORE.emit([&]() { 1203 OptimizationRemark R(DEBUG_TYPE, "ProcessLoopStridedStore", 1204 NewCall->getDebugLoc(), Preheader); 1205 R << "Transformed loop-strided store in " 1206 << ore::NV("Function", TheStore->getFunction()) 1207 << " function into a call to " 1208 << ore::NV("NewFunction", NewCall->getCalledFunction()) 1209 << "() intrinsic"; 1210 if (!Stores.empty()) 1211 R << ore::setExtraArgs(); 1212 for (auto *I : Stores) { 1213 R << ore::NV("FromBlock", I->getParent()->getName()) 1214 << ore::NV("ToBlock", Preheader->getName()); 1215 } 1216 return R; 1217 }); 1218 1219 // Okay, the memset has been formed. Zap the original store and anything that 1220 // feeds into it. 1221 for (auto *I : Stores) { 1222 if (MSSAU) 1223 MSSAU->removeMemoryAccess(I, true); 1224 deleteDeadInstruction(I); 1225 } 1226 if (MSSAU && VerifyMemorySSA) 1227 MSSAU->getMemorySSA()->verifyMemorySSA(); 1228 ++NumMemSet; 1229 ExpCleaner.markResultUsed(); 1230 return true; 1231 } 1232 1233 /// If the stored value is a strided load in the same loop with the same stride 1234 /// this may be transformable into a memcpy. This kicks in for stuff like 1235 /// for (i) A[i] = B[i]; 1236 bool LoopIdiomRecognize::processLoopStoreOfLoopLoad(StoreInst *SI, 1237 const SCEV *BECount) { 1238 assert(SI->isUnordered() && "Expected only non-volatile non-ordered stores."); 1239 1240 Value *StorePtr = SI->getPointerOperand(); 1241 const SCEVAddRecExpr *StoreEv = cast<SCEVAddRecExpr>(SE->getSCEV(StorePtr)); 1242 unsigned StoreSize = DL->getTypeStoreSize(SI->getValueOperand()->getType()); 1243 1244 // The store must be feeding a non-volatile load. 1245 LoadInst *LI = cast<LoadInst>(SI->getValueOperand()); 1246 assert(LI->isUnordered() && "Expected only non-volatile non-ordered loads."); 1247 1248 // See if the pointer expression is an AddRec like {base,+,1} on the current 1249 // loop, which indicates a strided load. If we have something else, it's a 1250 // random load we can't handle. 1251 Value *LoadPtr = LI->getPointerOperand(); 1252 const SCEVAddRecExpr *LoadEv = cast<SCEVAddRecExpr>(SE->getSCEV(LoadPtr)); 1253 1254 const SCEV *StoreSizeSCEV = SE->getConstant(StorePtr->getType(), StoreSize); 1255 return processLoopStoreOfLoopLoad(StorePtr, LoadPtr, StoreSizeSCEV, 1256 SI->getAlign(), LI->getAlign(), SI, LI, 1257 StoreEv, LoadEv, BECount); 1258 } 1259 1260 class MemmoveVerifier { 1261 public: 1262 explicit MemmoveVerifier(const Value &LoadBasePtr, const Value &StoreBasePtr, 1263 const DataLayout &DL) 1264 : DL(DL), LoadOff(0), StoreOff(0), 1265 BP1(llvm::GetPointerBaseWithConstantOffset( 1266 LoadBasePtr.stripPointerCasts(), LoadOff, DL)), 1267 BP2(llvm::GetPointerBaseWithConstantOffset( 1268 StoreBasePtr.stripPointerCasts(), StoreOff, DL)), 1269 IsSameObject(BP1 == BP2) {} 1270 1271 bool loadAndStoreMayFormMemmove(unsigned StoreSize, bool IsNegStride, 1272 const Instruction &TheLoad, 1273 bool IsMemCpy) const { 1274 if (IsMemCpy) { 1275 // Ensure that LoadBasePtr is after StoreBasePtr or before StoreBasePtr 1276 // for negative stride. 1277 if ((!IsNegStride && LoadOff <= StoreOff) || 1278 (IsNegStride && LoadOff >= StoreOff)) 1279 return false; 1280 } else { 1281 // Ensure that LoadBasePtr is after StoreBasePtr or before StoreBasePtr 1282 // for negative stride. LoadBasePtr shouldn't overlap with StoreBasePtr. 1283 int64_t LoadSize = 1284 DL.getTypeSizeInBits(TheLoad.getType()).getFixedSize() / 8; 1285 if (BP1 != BP2 || LoadSize != int64_t(StoreSize)) 1286 return false; 1287 if ((!IsNegStride && LoadOff < StoreOff + int64_t(StoreSize)) || 1288 (IsNegStride && LoadOff + LoadSize > StoreOff)) 1289 return false; 1290 } 1291 return true; 1292 } 1293 1294 private: 1295 const DataLayout &DL; 1296 int64_t LoadOff; 1297 int64_t StoreOff; 1298 const Value *BP1; 1299 const Value *BP2; 1300 1301 public: 1302 const bool IsSameObject; 1303 }; 1304 1305 bool LoopIdiomRecognize::processLoopStoreOfLoopLoad( 1306 Value *DestPtr, Value *SourcePtr, const SCEV *StoreSizeSCEV, 1307 MaybeAlign StoreAlign, MaybeAlign LoadAlign, Instruction *TheStore, 1308 Instruction *TheLoad, const SCEVAddRecExpr *StoreEv, 1309 const SCEVAddRecExpr *LoadEv, const SCEV *BECount) { 1310 1311 // FIXME: until llvm.memcpy.inline supports dynamic sizes, we need to 1312 // conservatively bail here, since otherwise we may have to transform 1313 // llvm.memcpy.inline into llvm.memcpy which is illegal. 1314 if (isa<MemCpyInlineInst>(TheStore)) 1315 return false; 1316 1317 // The trip count of the loop and the base pointer of the addrec SCEV is 1318 // guaranteed to be loop invariant, which means that it should dominate the 1319 // header. This allows us to insert code for it in the preheader. 1320 BasicBlock *Preheader = CurLoop->getLoopPreheader(); 1321 IRBuilder<> Builder(Preheader->getTerminator()); 1322 SCEVExpander Expander(*SE, *DL, "loop-idiom"); 1323 1324 SCEVExpanderCleaner ExpCleaner(Expander, *DT); 1325 1326 bool Changed = false; 1327 const SCEV *StrStart = StoreEv->getStart(); 1328 unsigned StrAS = DestPtr->getType()->getPointerAddressSpace(); 1329 Type *IntIdxTy = Builder.getIntNTy(DL->getIndexSizeInBits(StrAS)); 1330 1331 APInt Stride = getStoreStride(StoreEv); 1332 const SCEVConstant *ConstStoreSize = dyn_cast<SCEVConstant>(StoreSizeSCEV); 1333 1334 // TODO: Deal with non-constant size; Currently expect constant store size 1335 assert(ConstStoreSize && "store size is expected to be a constant"); 1336 1337 int64_t StoreSize = ConstStoreSize->getValue()->getZExtValue(); 1338 bool IsNegStride = StoreSize == -Stride; 1339 1340 // Handle negative strided loops. 1341 if (IsNegStride) 1342 StrStart = 1343 getStartForNegStride(StrStart, BECount, IntIdxTy, StoreSizeSCEV, SE); 1344 1345 // Okay, we have a strided store "p[i]" of a loaded value. We can turn 1346 // this into a memcpy in the loop preheader now if we want. However, this 1347 // would be unsafe to do if there is anything else in the loop that may read 1348 // or write the memory region we're storing to. This includes the load that 1349 // feeds the stores. Check for an alias by generating the base address and 1350 // checking everything. 1351 Value *StoreBasePtr = Expander.expandCodeFor( 1352 StrStart, Builder.getInt8PtrTy(StrAS), Preheader->getTerminator()); 1353 1354 // From here on out, conservatively report to the pass manager that we've 1355 // changed the IR, even if we later clean up these added instructions. There 1356 // may be structural differences e.g. in the order of use lists not accounted 1357 // for in just a textual dump of the IR. This is written as a variable, even 1358 // though statically all the places this dominates could be replaced with 1359 // 'true', with the hope that anyone trying to be clever / "more precise" with 1360 // the return value will read this comment, and leave them alone. 1361 Changed = true; 1362 1363 SmallPtrSet<Instruction *, 2> IgnoredInsts; 1364 IgnoredInsts.insert(TheStore); 1365 1366 bool IsMemCpy = isa<MemCpyInst>(TheStore); 1367 const StringRef InstRemark = IsMemCpy ? "memcpy" : "load and store"; 1368 1369 bool LoopAccessStore = 1370 mayLoopAccessLocation(StoreBasePtr, ModRefInfo::ModRef, CurLoop, BECount, 1371 StoreSizeSCEV, *AA, IgnoredInsts); 1372 if (LoopAccessStore) { 1373 // For memmove case it's not enough to guarantee that loop doesn't access 1374 // TheStore and TheLoad. Additionally we need to make sure that TheStore is 1375 // the only user of TheLoad. 1376 if (!TheLoad->hasOneUse()) 1377 return Changed; 1378 IgnoredInsts.insert(TheLoad); 1379 if (mayLoopAccessLocation(StoreBasePtr, ModRefInfo::ModRef, CurLoop, 1380 BECount, StoreSizeSCEV, *AA, IgnoredInsts)) { 1381 ORE.emit([&]() { 1382 return OptimizationRemarkMissed(DEBUG_TYPE, "LoopMayAccessStore", 1383 TheStore) 1384 << ore::NV("Inst", InstRemark) << " in " 1385 << ore::NV("Function", TheStore->getFunction()) 1386 << " function will not be hoisted: " 1387 << ore::NV("Reason", "The loop may access store location"); 1388 }); 1389 return Changed; 1390 } 1391 IgnoredInsts.erase(TheLoad); 1392 } 1393 1394 const SCEV *LdStart = LoadEv->getStart(); 1395 unsigned LdAS = SourcePtr->getType()->getPointerAddressSpace(); 1396 1397 // Handle negative strided loops. 1398 if (IsNegStride) 1399 LdStart = 1400 getStartForNegStride(LdStart, BECount, IntIdxTy, StoreSizeSCEV, SE); 1401 1402 // For a memcpy, we have to make sure that the input array is not being 1403 // mutated by the loop. 1404 Value *LoadBasePtr = Expander.expandCodeFor( 1405 LdStart, Builder.getInt8PtrTy(LdAS), Preheader->getTerminator()); 1406 1407 // If the store is a memcpy instruction, we must check if it will write to 1408 // the load memory locations. So remove it from the ignored stores. 1409 if (IsMemCpy) 1410 IgnoredInsts.erase(TheStore); 1411 MemmoveVerifier Verifier(*LoadBasePtr, *StoreBasePtr, *DL); 1412 if (mayLoopAccessLocation(LoadBasePtr, ModRefInfo::Mod, CurLoop, BECount, 1413 StoreSizeSCEV, *AA, IgnoredInsts)) { 1414 if (!IsMemCpy) { 1415 ORE.emit([&]() { 1416 return OptimizationRemarkMissed(DEBUG_TYPE, "LoopMayAccessLoad", 1417 TheLoad) 1418 << ore::NV("Inst", InstRemark) << " in " 1419 << ore::NV("Function", TheStore->getFunction()) 1420 << " function will not be hoisted: " 1421 << ore::NV("Reason", "The loop may access load location"); 1422 }); 1423 return Changed; 1424 } 1425 // At this point loop may access load only for memcpy in same underlying 1426 // object. If that's not the case bail out. 1427 if (!Verifier.IsSameObject) 1428 return Changed; 1429 } 1430 1431 bool UseMemMove = IsMemCpy ? Verifier.IsSameObject : LoopAccessStore; 1432 if (UseMemMove) 1433 if (!Verifier.loadAndStoreMayFormMemmove(StoreSize, IsNegStride, *TheLoad, 1434 IsMemCpy)) 1435 return Changed; 1436 1437 if (avoidLIRForMultiBlockLoop()) 1438 return Changed; 1439 1440 // Okay, everything is safe, we can transform this! 1441 1442 const SCEV *NumBytesS = 1443 getNumBytes(BECount, IntIdxTy, StoreSizeSCEV, CurLoop, DL, SE); 1444 1445 Value *NumBytes = 1446 Expander.expandCodeFor(NumBytesS, IntIdxTy, Preheader->getTerminator()); 1447 1448 CallInst *NewCall = nullptr; 1449 // Check whether to generate an unordered atomic memcpy: 1450 // If the load or store are atomic, then they must necessarily be unordered 1451 // by previous checks. 1452 if (!TheStore->isAtomic() && !TheLoad->isAtomic()) { 1453 if (UseMemMove) 1454 NewCall = Builder.CreateMemMove(StoreBasePtr, StoreAlign, LoadBasePtr, 1455 LoadAlign, NumBytes); 1456 else 1457 NewCall = Builder.CreateMemCpy(StoreBasePtr, StoreAlign, LoadBasePtr, 1458 LoadAlign, NumBytes); 1459 } else { 1460 // For now don't support unordered atomic memmove. 1461 if (UseMemMove) 1462 return Changed; 1463 // We cannot allow unaligned ops for unordered load/store, so reject 1464 // anything where the alignment isn't at least the element size. 1465 assert((StoreAlign.hasValue() && LoadAlign.hasValue()) && 1466 "Expect unordered load/store to have align."); 1467 if (StoreAlign.getValue() < StoreSize || LoadAlign.getValue() < StoreSize) 1468 return Changed; 1469 1470 // If the element.atomic memcpy is not lowered into explicit 1471 // loads/stores later, then it will be lowered into an element-size 1472 // specific lib call. If the lib call doesn't exist for our store size, then 1473 // we shouldn't generate the memcpy. 1474 if (StoreSize > TTI->getAtomicMemIntrinsicMaxElementSize()) 1475 return Changed; 1476 1477 // Create the call. 1478 // Note that unordered atomic loads/stores are *required* by the spec to 1479 // have an alignment but non-atomic loads/stores may not. 1480 NewCall = Builder.CreateElementUnorderedAtomicMemCpy( 1481 StoreBasePtr, StoreAlign.getValue(), LoadBasePtr, LoadAlign.getValue(), 1482 NumBytes, StoreSize); 1483 } 1484 NewCall->setDebugLoc(TheStore->getDebugLoc()); 1485 1486 if (MSSAU) { 1487 MemoryAccess *NewMemAcc = MSSAU->createMemoryAccessInBB( 1488 NewCall, nullptr, NewCall->getParent(), MemorySSA::BeforeTerminator); 1489 MSSAU->insertDef(cast<MemoryDef>(NewMemAcc), true); 1490 } 1491 1492 LLVM_DEBUG(dbgs() << " Formed new call: " << *NewCall << "\n" 1493 << " from load ptr=" << *LoadEv << " at: " << *TheLoad 1494 << "\n" 1495 << " from store ptr=" << *StoreEv << " at: " << *TheStore 1496 << "\n"); 1497 1498 ORE.emit([&]() { 1499 return OptimizationRemark(DEBUG_TYPE, "ProcessLoopStoreOfLoopLoad", 1500 NewCall->getDebugLoc(), Preheader) 1501 << "Formed a call to " 1502 << ore::NV("NewFunction", NewCall->getCalledFunction()) 1503 << "() intrinsic from " << ore::NV("Inst", InstRemark) 1504 << " instruction in " << ore::NV("Function", TheStore->getFunction()) 1505 << " function" 1506 << ore::setExtraArgs() 1507 << ore::NV("FromBlock", TheStore->getParent()->getName()) 1508 << ore::NV("ToBlock", Preheader->getName()); 1509 }); 1510 1511 // Okay, a new call to memcpy/memmove has been formed. Zap the original store 1512 // and anything that feeds into it. 1513 if (MSSAU) 1514 MSSAU->removeMemoryAccess(TheStore, true); 1515 deleteDeadInstruction(TheStore); 1516 if (MSSAU && VerifyMemorySSA) 1517 MSSAU->getMemorySSA()->verifyMemorySSA(); 1518 if (UseMemMove) 1519 ++NumMemMove; 1520 else 1521 ++NumMemCpy; 1522 ExpCleaner.markResultUsed(); 1523 return true; 1524 } 1525 1526 // When compiling for codesize we avoid idiom recognition for a multi-block loop 1527 // unless it is a loop_memset idiom or a memset/memcpy idiom in a nested loop. 1528 // 1529 bool LoopIdiomRecognize::avoidLIRForMultiBlockLoop(bool IsMemset, 1530 bool IsLoopMemset) { 1531 if (ApplyCodeSizeHeuristics && CurLoop->getNumBlocks() > 1) { 1532 if (CurLoop->isOutermost() && (!IsMemset || !IsLoopMemset)) { 1533 LLVM_DEBUG(dbgs() << " " << CurLoop->getHeader()->getParent()->getName() 1534 << " : LIR " << (IsMemset ? "Memset" : "Memcpy") 1535 << " avoided: multi-block top-level loop\n"); 1536 return true; 1537 } 1538 } 1539 1540 return false; 1541 } 1542 1543 bool LoopIdiomRecognize::runOnNoncountableLoop() { 1544 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Scanning: F[" 1545 << CurLoop->getHeader()->getParent()->getName() 1546 << "] Noncountable Loop %" 1547 << CurLoop->getHeader()->getName() << "\n"); 1548 1549 return recognizePopcount() || recognizeAndInsertFFS() || 1550 recognizeShiftUntilBitTest() || recognizeShiftUntilZero(); 1551 } 1552 1553 /// Check if the given conditional branch is based on the comparison between 1554 /// a variable and zero, and if the variable is non-zero or zero (JmpOnZero is 1555 /// true), the control yields to the loop entry. If the branch matches the 1556 /// behavior, the variable involved in the comparison is returned. This function 1557 /// will be called to see if the precondition and postcondition of the loop are 1558 /// in desirable form. 1559 static Value *matchCondition(BranchInst *BI, BasicBlock *LoopEntry, 1560 bool JmpOnZero = false) { 1561 if (!BI || !BI->isConditional()) 1562 return nullptr; 1563 1564 ICmpInst *Cond = dyn_cast<ICmpInst>(BI->getCondition()); 1565 if (!Cond) 1566 return nullptr; 1567 1568 ConstantInt *CmpZero = dyn_cast<ConstantInt>(Cond->getOperand(1)); 1569 if (!CmpZero || !CmpZero->isZero()) 1570 return nullptr; 1571 1572 BasicBlock *TrueSucc = BI->getSuccessor(0); 1573 BasicBlock *FalseSucc = BI->getSuccessor(1); 1574 if (JmpOnZero) 1575 std::swap(TrueSucc, FalseSucc); 1576 1577 ICmpInst::Predicate Pred = Cond->getPredicate(); 1578 if ((Pred == ICmpInst::ICMP_NE && TrueSucc == LoopEntry) || 1579 (Pred == ICmpInst::ICMP_EQ && FalseSucc == LoopEntry)) 1580 return Cond->getOperand(0); 1581 1582 return nullptr; 1583 } 1584 1585 // Check if the recurrence variable `VarX` is in the right form to create 1586 // the idiom. Returns the value coerced to a PHINode if so. 1587 static PHINode *getRecurrenceVar(Value *VarX, Instruction *DefX, 1588 BasicBlock *LoopEntry) { 1589 auto *PhiX = dyn_cast<PHINode>(VarX); 1590 if (PhiX && PhiX->getParent() == LoopEntry && 1591 (PhiX->getOperand(0) == DefX || PhiX->getOperand(1) == DefX)) 1592 return PhiX; 1593 return nullptr; 1594 } 1595 1596 /// Return true iff the idiom is detected in the loop. 1597 /// 1598 /// Additionally: 1599 /// 1) \p CntInst is set to the instruction counting the population bit. 1600 /// 2) \p CntPhi is set to the corresponding phi node. 1601 /// 3) \p Var is set to the value whose population bits are being counted. 1602 /// 1603 /// The core idiom we are trying to detect is: 1604 /// \code 1605 /// if (x0 != 0) 1606 /// goto loop-exit // the precondition of the loop 1607 /// cnt0 = init-val; 1608 /// do { 1609 /// x1 = phi (x0, x2); 1610 /// cnt1 = phi(cnt0, cnt2); 1611 /// 1612 /// cnt2 = cnt1 + 1; 1613 /// ... 1614 /// x2 = x1 & (x1 - 1); 1615 /// ... 1616 /// } while(x != 0); 1617 /// 1618 /// loop-exit: 1619 /// \endcode 1620 static bool detectPopcountIdiom(Loop *CurLoop, BasicBlock *PreCondBB, 1621 Instruction *&CntInst, PHINode *&CntPhi, 1622 Value *&Var) { 1623 // step 1: Check to see if the look-back branch match this pattern: 1624 // "if (a!=0) goto loop-entry". 1625 BasicBlock *LoopEntry; 1626 Instruction *DefX2, *CountInst; 1627 Value *VarX1, *VarX0; 1628 PHINode *PhiX, *CountPhi; 1629 1630 DefX2 = CountInst = nullptr; 1631 VarX1 = VarX0 = nullptr; 1632 PhiX = CountPhi = nullptr; 1633 LoopEntry = *(CurLoop->block_begin()); 1634 1635 // step 1: Check if the loop-back branch is in desirable form. 1636 { 1637 if (Value *T = matchCondition( 1638 dyn_cast<BranchInst>(LoopEntry->getTerminator()), LoopEntry)) 1639 DefX2 = dyn_cast<Instruction>(T); 1640 else 1641 return false; 1642 } 1643 1644 // step 2: detect instructions corresponding to "x2 = x1 & (x1 - 1)" 1645 { 1646 if (!DefX2 || DefX2->getOpcode() != Instruction::And) 1647 return false; 1648 1649 BinaryOperator *SubOneOp; 1650 1651 if ((SubOneOp = dyn_cast<BinaryOperator>(DefX2->getOperand(0)))) 1652 VarX1 = DefX2->getOperand(1); 1653 else { 1654 VarX1 = DefX2->getOperand(0); 1655 SubOneOp = dyn_cast<BinaryOperator>(DefX2->getOperand(1)); 1656 } 1657 if (!SubOneOp || SubOneOp->getOperand(0) != VarX1) 1658 return false; 1659 1660 ConstantInt *Dec = dyn_cast<ConstantInt>(SubOneOp->getOperand(1)); 1661 if (!Dec || 1662 !((SubOneOp->getOpcode() == Instruction::Sub && Dec->isOne()) || 1663 (SubOneOp->getOpcode() == Instruction::Add && 1664 Dec->isMinusOne()))) { 1665 return false; 1666 } 1667 } 1668 1669 // step 3: Check the recurrence of variable X 1670 PhiX = getRecurrenceVar(VarX1, DefX2, LoopEntry); 1671 if (!PhiX) 1672 return false; 1673 1674 // step 4: Find the instruction which count the population: cnt2 = cnt1 + 1 1675 { 1676 CountInst = nullptr; 1677 for (BasicBlock::iterator Iter = LoopEntry->getFirstNonPHI()->getIterator(), 1678 IterE = LoopEntry->end(); 1679 Iter != IterE; Iter++) { 1680 Instruction *Inst = &*Iter; 1681 if (Inst->getOpcode() != Instruction::Add) 1682 continue; 1683 1684 ConstantInt *Inc = dyn_cast<ConstantInt>(Inst->getOperand(1)); 1685 if (!Inc || !Inc->isOne()) 1686 continue; 1687 1688 PHINode *Phi = getRecurrenceVar(Inst->getOperand(0), Inst, LoopEntry); 1689 if (!Phi) 1690 continue; 1691 1692 // Check if the result of the instruction is live of the loop. 1693 bool LiveOutLoop = false; 1694 for (User *U : Inst->users()) { 1695 if ((cast<Instruction>(U))->getParent() != LoopEntry) { 1696 LiveOutLoop = true; 1697 break; 1698 } 1699 } 1700 1701 if (LiveOutLoop) { 1702 CountInst = Inst; 1703 CountPhi = Phi; 1704 break; 1705 } 1706 } 1707 1708 if (!CountInst) 1709 return false; 1710 } 1711 1712 // step 5: check if the precondition is in this form: 1713 // "if (x != 0) goto loop-head ; else goto somewhere-we-don't-care;" 1714 { 1715 auto *PreCondBr = dyn_cast<BranchInst>(PreCondBB->getTerminator()); 1716 Value *T = matchCondition(PreCondBr, CurLoop->getLoopPreheader()); 1717 if (T != PhiX->getOperand(0) && T != PhiX->getOperand(1)) 1718 return false; 1719 1720 CntInst = CountInst; 1721 CntPhi = CountPhi; 1722 Var = T; 1723 } 1724 1725 return true; 1726 } 1727 1728 /// Return true if the idiom is detected in the loop. 1729 /// 1730 /// Additionally: 1731 /// 1) \p CntInst is set to the instruction Counting Leading Zeros (CTLZ) 1732 /// or nullptr if there is no such. 1733 /// 2) \p CntPhi is set to the corresponding phi node 1734 /// or nullptr if there is no such. 1735 /// 3) \p Var is set to the value whose CTLZ could be used. 1736 /// 4) \p DefX is set to the instruction calculating Loop exit condition. 1737 /// 1738 /// The core idiom we are trying to detect is: 1739 /// \code 1740 /// if (x0 == 0) 1741 /// goto loop-exit // the precondition of the loop 1742 /// cnt0 = init-val; 1743 /// do { 1744 /// x = phi (x0, x.next); //PhiX 1745 /// cnt = phi(cnt0, cnt.next); 1746 /// 1747 /// cnt.next = cnt + 1; 1748 /// ... 1749 /// x.next = x >> 1; // DefX 1750 /// ... 1751 /// } while(x.next != 0); 1752 /// 1753 /// loop-exit: 1754 /// \endcode 1755 static bool detectShiftUntilZeroIdiom(Loop *CurLoop, const DataLayout &DL, 1756 Intrinsic::ID &IntrinID, Value *&InitX, 1757 Instruction *&CntInst, PHINode *&CntPhi, 1758 Instruction *&DefX) { 1759 BasicBlock *LoopEntry; 1760 Value *VarX = nullptr; 1761 1762 DefX = nullptr; 1763 CntInst = nullptr; 1764 CntPhi = nullptr; 1765 LoopEntry = *(CurLoop->block_begin()); 1766 1767 // step 1: Check if the loop-back branch is in desirable form. 1768 if (Value *T = matchCondition( 1769 dyn_cast<BranchInst>(LoopEntry->getTerminator()), LoopEntry)) 1770 DefX = dyn_cast<Instruction>(T); 1771 else 1772 return false; 1773 1774 // step 2: detect instructions corresponding to "x.next = x >> 1 or x << 1" 1775 if (!DefX || !DefX->isShift()) 1776 return false; 1777 IntrinID = DefX->getOpcode() == Instruction::Shl ? Intrinsic::cttz : 1778 Intrinsic::ctlz; 1779 ConstantInt *Shft = dyn_cast<ConstantInt>(DefX->getOperand(1)); 1780 if (!Shft || !Shft->isOne()) 1781 return false; 1782 VarX = DefX->getOperand(0); 1783 1784 // step 3: Check the recurrence of variable X 1785 PHINode *PhiX = getRecurrenceVar(VarX, DefX, LoopEntry); 1786 if (!PhiX) 1787 return false; 1788 1789 InitX = PhiX->getIncomingValueForBlock(CurLoop->getLoopPreheader()); 1790 1791 // Make sure the initial value can't be negative otherwise the ashr in the 1792 // loop might never reach zero which would make the loop infinite. 1793 if (DefX->getOpcode() == Instruction::AShr && !isKnownNonNegative(InitX, DL)) 1794 return false; 1795 1796 // step 4: Find the instruction which count the CTLZ: cnt.next = cnt + 1 1797 // or cnt.next = cnt + -1. 1798 // TODO: We can skip the step. If loop trip count is known (CTLZ), 1799 // then all uses of "cnt.next" could be optimized to the trip count 1800 // plus "cnt0". Currently it is not optimized. 1801 // This step could be used to detect POPCNT instruction: 1802 // cnt.next = cnt + (x.next & 1) 1803 for (BasicBlock::iterator Iter = LoopEntry->getFirstNonPHI()->getIterator(), 1804 IterE = LoopEntry->end(); 1805 Iter != IterE; Iter++) { 1806 Instruction *Inst = &*Iter; 1807 if (Inst->getOpcode() != Instruction::Add) 1808 continue; 1809 1810 ConstantInt *Inc = dyn_cast<ConstantInt>(Inst->getOperand(1)); 1811 if (!Inc || (!Inc->isOne() && !Inc->isMinusOne())) 1812 continue; 1813 1814 PHINode *Phi = getRecurrenceVar(Inst->getOperand(0), Inst, LoopEntry); 1815 if (!Phi) 1816 continue; 1817 1818 CntInst = Inst; 1819 CntPhi = Phi; 1820 break; 1821 } 1822 if (!CntInst) 1823 return false; 1824 1825 return true; 1826 } 1827 1828 /// Recognize CTLZ or CTTZ idiom in a non-countable loop and convert the loop 1829 /// to countable (with CTLZ / CTTZ trip count). If CTLZ / CTTZ inserted as a new 1830 /// trip count returns true; otherwise, returns false. 1831 bool LoopIdiomRecognize::recognizeAndInsertFFS() { 1832 // Give up if the loop has multiple blocks or multiple backedges. 1833 if (CurLoop->getNumBackEdges() != 1 || CurLoop->getNumBlocks() != 1) 1834 return false; 1835 1836 Intrinsic::ID IntrinID; 1837 Value *InitX; 1838 Instruction *DefX = nullptr; 1839 PHINode *CntPhi = nullptr; 1840 Instruction *CntInst = nullptr; 1841 // Help decide if transformation is profitable. For ShiftUntilZero idiom, 1842 // this is always 6. 1843 size_t IdiomCanonicalSize = 6; 1844 1845 if (!detectShiftUntilZeroIdiom(CurLoop, *DL, IntrinID, InitX, 1846 CntInst, CntPhi, DefX)) 1847 return false; 1848 1849 bool IsCntPhiUsedOutsideLoop = false; 1850 for (User *U : CntPhi->users()) 1851 if (!CurLoop->contains(cast<Instruction>(U))) { 1852 IsCntPhiUsedOutsideLoop = true; 1853 break; 1854 } 1855 bool IsCntInstUsedOutsideLoop = false; 1856 for (User *U : CntInst->users()) 1857 if (!CurLoop->contains(cast<Instruction>(U))) { 1858 IsCntInstUsedOutsideLoop = true; 1859 break; 1860 } 1861 // If both CntInst and CntPhi are used outside the loop the profitability 1862 // is questionable. 1863 if (IsCntInstUsedOutsideLoop && IsCntPhiUsedOutsideLoop) 1864 return false; 1865 1866 // For some CPUs result of CTLZ(X) intrinsic is undefined 1867 // when X is 0. If we can not guarantee X != 0, we need to check this 1868 // when expand. 1869 bool ZeroCheck = false; 1870 // It is safe to assume Preheader exist as it was checked in 1871 // parent function RunOnLoop. 1872 BasicBlock *PH = CurLoop->getLoopPreheader(); 1873 1874 // If we are using the count instruction outside the loop, make sure we 1875 // have a zero check as a precondition. Without the check the loop would run 1876 // one iteration for before any check of the input value. This means 0 and 1 1877 // would have identical behavior in the original loop and thus 1878 if (!IsCntPhiUsedOutsideLoop) { 1879 auto *PreCondBB = PH->getSinglePredecessor(); 1880 if (!PreCondBB) 1881 return false; 1882 auto *PreCondBI = dyn_cast<BranchInst>(PreCondBB->getTerminator()); 1883 if (!PreCondBI) 1884 return false; 1885 if (matchCondition(PreCondBI, PH) != InitX) 1886 return false; 1887 ZeroCheck = true; 1888 } 1889 1890 // Check if CTLZ / CTTZ intrinsic is profitable. Assume it is always 1891 // profitable if we delete the loop. 1892 1893 // the loop has only 6 instructions: 1894 // %n.addr.0 = phi [ %n, %entry ], [ %shr, %while.cond ] 1895 // %i.0 = phi [ %i0, %entry ], [ %inc, %while.cond ] 1896 // %shr = ashr %n.addr.0, 1 1897 // %tobool = icmp eq %shr, 0 1898 // %inc = add nsw %i.0, 1 1899 // br i1 %tobool 1900 1901 const Value *Args[] = {InitX, 1902 ConstantInt::getBool(InitX->getContext(), ZeroCheck)}; 1903 1904 // @llvm.dbg doesn't count as they have no semantic effect. 1905 auto InstWithoutDebugIt = CurLoop->getHeader()->instructionsWithoutDebug(); 1906 uint32_t HeaderSize = 1907 std::distance(InstWithoutDebugIt.begin(), InstWithoutDebugIt.end()); 1908 1909 IntrinsicCostAttributes Attrs(IntrinID, InitX->getType(), Args); 1910 InstructionCost Cost = 1911 TTI->getIntrinsicInstrCost(Attrs, TargetTransformInfo::TCK_SizeAndLatency); 1912 if (HeaderSize != IdiomCanonicalSize && 1913 Cost > TargetTransformInfo::TCC_Basic) 1914 return false; 1915 1916 transformLoopToCountable(IntrinID, PH, CntInst, CntPhi, InitX, DefX, 1917 DefX->getDebugLoc(), ZeroCheck, 1918 IsCntPhiUsedOutsideLoop); 1919 return true; 1920 } 1921 1922 /// Recognizes a population count idiom in a non-countable loop. 1923 /// 1924 /// If detected, transforms the relevant code to issue the popcount intrinsic 1925 /// function call, and returns true; otherwise, returns false. 1926 bool LoopIdiomRecognize::recognizePopcount() { 1927 if (TTI->getPopcntSupport(32) != TargetTransformInfo::PSK_FastHardware) 1928 return false; 1929 1930 // Counting population are usually conducted by few arithmetic instructions. 1931 // Such instructions can be easily "absorbed" by vacant slots in a 1932 // non-compact loop. Therefore, recognizing popcount idiom only makes sense 1933 // in a compact loop. 1934 1935 // Give up if the loop has multiple blocks or multiple backedges. 1936 if (CurLoop->getNumBackEdges() != 1 || CurLoop->getNumBlocks() != 1) 1937 return false; 1938 1939 BasicBlock *LoopBody = *(CurLoop->block_begin()); 1940 if (LoopBody->size() >= 20) { 1941 // The loop is too big, bail out. 1942 return false; 1943 } 1944 1945 // It should have a preheader containing nothing but an unconditional branch. 1946 BasicBlock *PH = CurLoop->getLoopPreheader(); 1947 if (!PH || &PH->front() != PH->getTerminator()) 1948 return false; 1949 auto *EntryBI = dyn_cast<BranchInst>(PH->getTerminator()); 1950 if (!EntryBI || EntryBI->isConditional()) 1951 return false; 1952 1953 // It should have a precondition block where the generated popcount intrinsic 1954 // function can be inserted. 1955 auto *PreCondBB = PH->getSinglePredecessor(); 1956 if (!PreCondBB) 1957 return false; 1958 auto *PreCondBI = dyn_cast<BranchInst>(PreCondBB->getTerminator()); 1959 if (!PreCondBI || PreCondBI->isUnconditional()) 1960 return false; 1961 1962 Instruction *CntInst; 1963 PHINode *CntPhi; 1964 Value *Val; 1965 if (!detectPopcountIdiom(CurLoop, PreCondBB, CntInst, CntPhi, Val)) 1966 return false; 1967 1968 transformLoopToPopcount(PreCondBB, CntInst, CntPhi, Val); 1969 return true; 1970 } 1971 1972 static CallInst *createPopcntIntrinsic(IRBuilder<> &IRBuilder, Value *Val, 1973 const DebugLoc &DL) { 1974 Value *Ops[] = {Val}; 1975 Type *Tys[] = {Val->getType()}; 1976 1977 Module *M = IRBuilder.GetInsertBlock()->getParent()->getParent(); 1978 Function *Func = Intrinsic::getDeclaration(M, Intrinsic::ctpop, Tys); 1979 CallInst *CI = IRBuilder.CreateCall(Func, Ops); 1980 CI->setDebugLoc(DL); 1981 1982 return CI; 1983 } 1984 1985 static CallInst *createFFSIntrinsic(IRBuilder<> &IRBuilder, Value *Val, 1986 const DebugLoc &DL, bool ZeroCheck, 1987 Intrinsic::ID IID) { 1988 Value *Ops[] = {Val, IRBuilder.getInt1(ZeroCheck)}; 1989 Type *Tys[] = {Val->getType()}; 1990 1991 Module *M = IRBuilder.GetInsertBlock()->getParent()->getParent(); 1992 Function *Func = Intrinsic::getDeclaration(M, IID, Tys); 1993 CallInst *CI = IRBuilder.CreateCall(Func, Ops); 1994 CI->setDebugLoc(DL); 1995 1996 return CI; 1997 } 1998 1999 /// Transform the following loop (Using CTLZ, CTTZ is similar): 2000 /// loop: 2001 /// CntPhi = PHI [Cnt0, CntInst] 2002 /// PhiX = PHI [InitX, DefX] 2003 /// CntInst = CntPhi + 1 2004 /// DefX = PhiX >> 1 2005 /// LOOP_BODY 2006 /// Br: loop if (DefX != 0) 2007 /// Use(CntPhi) or Use(CntInst) 2008 /// 2009 /// Into: 2010 /// If CntPhi used outside the loop: 2011 /// CountPrev = BitWidth(InitX) - CTLZ(InitX >> 1) 2012 /// Count = CountPrev + 1 2013 /// else 2014 /// Count = BitWidth(InitX) - CTLZ(InitX) 2015 /// loop: 2016 /// CntPhi = PHI [Cnt0, CntInst] 2017 /// PhiX = PHI [InitX, DefX] 2018 /// PhiCount = PHI [Count, Dec] 2019 /// CntInst = CntPhi + 1 2020 /// DefX = PhiX >> 1 2021 /// Dec = PhiCount - 1 2022 /// LOOP_BODY 2023 /// Br: loop if (Dec != 0) 2024 /// Use(CountPrev + Cnt0) // Use(CntPhi) 2025 /// or 2026 /// Use(Count + Cnt0) // Use(CntInst) 2027 /// 2028 /// If LOOP_BODY is empty the loop will be deleted. 2029 /// If CntInst and DefX are not used in LOOP_BODY they will be removed. 2030 void LoopIdiomRecognize::transformLoopToCountable( 2031 Intrinsic::ID IntrinID, BasicBlock *Preheader, Instruction *CntInst, 2032 PHINode *CntPhi, Value *InitX, Instruction *DefX, const DebugLoc &DL, 2033 bool ZeroCheck, bool IsCntPhiUsedOutsideLoop) { 2034 BranchInst *PreheaderBr = cast<BranchInst>(Preheader->getTerminator()); 2035 2036 // Step 1: Insert the CTLZ/CTTZ instruction at the end of the preheader block 2037 IRBuilder<> Builder(PreheaderBr); 2038 Builder.SetCurrentDebugLocation(DL); 2039 2040 // If there are no uses of CntPhi crate: 2041 // Count = BitWidth - CTLZ(InitX); 2042 // NewCount = Count; 2043 // If there are uses of CntPhi create: 2044 // NewCount = BitWidth - CTLZ(InitX >> 1); 2045 // Count = NewCount + 1; 2046 Value *InitXNext; 2047 if (IsCntPhiUsedOutsideLoop) { 2048 if (DefX->getOpcode() == Instruction::AShr) 2049 InitXNext = Builder.CreateAShr(InitX, 1); 2050 else if (DefX->getOpcode() == Instruction::LShr) 2051 InitXNext = Builder.CreateLShr(InitX, 1); 2052 else if (DefX->getOpcode() == Instruction::Shl) // cttz 2053 InitXNext = Builder.CreateShl(InitX, 1); 2054 else 2055 llvm_unreachable("Unexpected opcode!"); 2056 } else 2057 InitXNext = InitX; 2058 Value *Count = 2059 createFFSIntrinsic(Builder, InitXNext, DL, ZeroCheck, IntrinID); 2060 Type *CountTy = Count->getType(); 2061 Count = Builder.CreateSub( 2062 ConstantInt::get(CountTy, CountTy->getIntegerBitWidth()), Count); 2063 Value *NewCount = Count; 2064 if (IsCntPhiUsedOutsideLoop) 2065 Count = Builder.CreateAdd(Count, ConstantInt::get(CountTy, 1)); 2066 2067 NewCount = Builder.CreateZExtOrTrunc(NewCount, CntInst->getType()); 2068 2069 Value *CntInitVal = CntPhi->getIncomingValueForBlock(Preheader); 2070 if (cast<ConstantInt>(CntInst->getOperand(1))->isOne()) { 2071 // If the counter was being incremented in the loop, add NewCount to the 2072 // counter's initial value, but only if the initial value is not zero. 2073 ConstantInt *InitConst = dyn_cast<ConstantInt>(CntInitVal); 2074 if (!InitConst || !InitConst->isZero()) 2075 NewCount = Builder.CreateAdd(NewCount, CntInitVal); 2076 } else { 2077 // If the count was being decremented in the loop, subtract NewCount from 2078 // the counter's initial value. 2079 NewCount = Builder.CreateSub(CntInitVal, NewCount); 2080 } 2081 2082 // Step 2: Insert new IV and loop condition: 2083 // loop: 2084 // ... 2085 // PhiCount = PHI [Count, Dec] 2086 // ... 2087 // Dec = PhiCount - 1 2088 // ... 2089 // Br: loop if (Dec != 0) 2090 BasicBlock *Body = *(CurLoop->block_begin()); 2091 auto *LbBr = cast<BranchInst>(Body->getTerminator()); 2092 ICmpInst *LbCond = cast<ICmpInst>(LbBr->getCondition()); 2093 2094 PHINode *TcPhi = PHINode::Create(CountTy, 2, "tcphi", &Body->front()); 2095 2096 Builder.SetInsertPoint(LbCond); 2097 Instruction *TcDec = cast<Instruction>(Builder.CreateSub( 2098 TcPhi, ConstantInt::get(CountTy, 1), "tcdec", false, true)); 2099 2100 TcPhi->addIncoming(Count, Preheader); 2101 TcPhi->addIncoming(TcDec, Body); 2102 2103 CmpInst::Predicate Pred = 2104 (LbBr->getSuccessor(0) == Body) ? CmpInst::ICMP_NE : CmpInst::ICMP_EQ; 2105 LbCond->setPredicate(Pred); 2106 LbCond->setOperand(0, TcDec); 2107 LbCond->setOperand(1, ConstantInt::get(CountTy, 0)); 2108 2109 // Step 3: All the references to the original counter outside 2110 // the loop are replaced with the NewCount 2111 if (IsCntPhiUsedOutsideLoop) 2112 CntPhi->replaceUsesOutsideBlock(NewCount, Body); 2113 else 2114 CntInst->replaceUsesOutsideBlock(NewCount, Body); 2115 2116 // step 4: Forget the "non-computable" trip-count SCEV associated with the 2117 // loop. The loop would otherwise not be deleted even if it becomes empty. 2118 SE->forgetLoop(CurLoop); 2119 } 2120 2121 void LoopIdiomRecognize::transformLoopToPopcount(BasicBlock *PreCondBB, 2122 Instruction *CntInst, 2123 PHINode *CntPhi, Value *Var) { 2124 BasicBlock *PreHead = CurLoop->getLoopPreheader(); 2125 auto *PreCondBr = cast<BranchInst>(PreCondBB->getTerminator()); 2126 const DebugLoc &DL = CntInst->getDebugLoc(); 2127 2128 // Assuming before transformation, the loop is following: 2129 // if (x) // the precondition 2130 // do { cnt++; x &= x - 1; } while(x); 2131 2132 // Step 1: Insert the ctpop instruction at the end of the precondition block 2133 IRBuilder<> Builder(PreCondBr); 2134 Value *PopCnt, *PopCntZext, *NewCount, *TripCnt; 2135 { 2136 PopCnt = createPopcntIntrinsic(Builder, Var, DL); 2137 NewCount = PopCntZext = 2138 Builder.CreateZExtOrTrunc(PopCnt, cast<IntegerType>(CntPhi->getType())); 2139 2140 if (NewCount != PopCnt) 2141 (cast<Instruction>(NewCount))->setDebugLoc(DL); 2142 2143 // TripCnt is exactly the number of iterations the loop has 2144 TripCnt = NewCount; 2145 2146 // If the population counter's initial value is not zero, insert Add Inst. 2147 Value *CntInitVal = CntPhi->getIncomingValueForBlock(PreHead); 2148 ConstantInt *InitConst = dyn_cast<ConstantInt>(CntInitVal); 2149 if (!InitConst || !InitConst->isZero()) { 2150 NewCount = Builder.CreateAdd(NewCount, CntInitVal); 2151 (cast<Instruction>(NewCount))->setDebugLoc(DL); 2152 } 2153 } 2154 2155 // Step 2: Replace the precondition from "if (x == 0) goto loop-exit" to 2156 // "if (NewCount == 0) loop-exit". Without this change, the intrinsic 2157 // function would be partial dead code, and downstream passes will drag 2158 // it back from the precondition block to the preheader. 2159 { 2160 ICmpInst *PreCond = cast<ICmpInst>(PreCondBr->getCondition()); 2161 2162 Value *Opnd0 = PopCntZext; 2163 Value *Opnd1 = ConstantInt::get(PopCntZext->getType(), 0); 2164 if (PreCond->getOperand(0) != Var) 2165 std::swap(Opnd0, Opnd1); 2166 2167 ICmpInst *NewPreCond = cast<ICmpInst>( 2168 Builder.CreateICmp(PreCond->getPredicate(), Opnd0, Opnd1)); 2169 PreCondBr->setCondition(NewPreCond); 2170 2171 RecursivelyDeleteTriviallyDeadInstructions(PreCond, TLI); 2172 } 2173 2174 // Step 3: Note that the population count is exactly the trip count of the 2175 // loop in question, which enable us to convert the loop from noncountable 2176 // loop into a countable one. The benefit is twofold: 2177 // 2178 // - If the loop only counts population, the entire loop becomes dead after 2179 // the transformation. It is a lot easier to prove a countable loop dead 2180 // than to prove a noncountable one. (In some C dialects, an infinite loop 2181 // isn't dead even if it computes nothing useful. In general, DCE needs 2182 // to prove a noncountable loop finite before safely delete it.) 2183 // 2184 // - If the loop also performs something else, it remains alive. 2185 // Since it is transformed to countable form, it can be aggressively 2186 // optimized by some optimizations which are in general not applicable 2187 // to a noncountable loop. 2188 // 2189 // After this step, this loop (conceptually) would look like following: 2190 // newcnt = __builtin_ctpop(x); 2191 // t = newcnt; 2192 // if (x) 2193 // do { cnt++; x &= x-1; t--) } while (t > 0); 2194 BasicBlock *Body = *(CurLoop->block_begin()); 2195 { 2196 auto *LbBr = cast<BranchInst>(Body->getTerminator()); 2197 ICmpInst *LbCond = cast<ICmpInst>(LbBr->getCondition()); 2198 Type *Ty = TripCnt->getType(); 2199 2200 PHINode *TcPhi = PHINode::Create(Ty, 2, "tcphi", &Body->front()); 2201 2202 Builder.SetInsertPoint(LbCond); 2203 Instruction *TcDec = cast<Instruction>( 2204 Builder.CreateSub(TcPhi, ConstantInt::get(Ty, 1), 2205 "tcdec", false, true)); 2206 2207 TcPhi->addIncoming(TripCnt, PreHead); 2208 TcPhi->addIncoming(TcDec, Body); 2209 2210 CmpInst::Predicate Pred = 2211 (LbBr->getSuccessor(0) == Body) ? CmpInst::ICMP_UGT : CmpInst::ICMP_SLE; 2212 LbCond->setPredicate(Pred); 2213 LbCond->setOperand(0, TcDec); 2214 LbCond->setOperand(1, ConstantInt::get(Ty, 0)); 2215 } 2216 2217 // Step 4: All the references to the original population counter outside 2218 // the loop are replaced with the NewCount -- the value returned from 2219 // __builtin_ctpop(). 2220 CntInst->replaceUsesOutsideBlock(NewCount, Body); 2221 2222 // step 5: Forget the "non-computable" trip-count SCEV associated with the 2223 // loop. The loop would otherwise not be deleted even if it becomes empty. 2224 SE->forgetLoop(CurLoop); 2225 } 2226 2227 /// Match loop-invariant value. 2228 template <typename SubPattern_t> struct match_LoopInvariant { 2229 SubPattern_t SubPattern; 2230 const Loop *L; 2231 2232 match_LoopInvariant(const SubPattern_t &SP, const Loop *L) 2233 : SubPattern(SP), L(L) {} 2234 2235 template <typename ITy> bool match(ITy *V) { 2236 return L->isLoopInvariant(V) && SubPattern.match(V); 2237 } 2238 }; 2239 2240 /// Matches if the value is loop-invariant. 2241 template <typename Ty> 2242 inline match_LoopInvariant<Ty> m_LoopInvariant(const Ty &M, const Loop *L) { 2243 return match_LoopInvariant<Ty>(M, L); 2244 } 2245 2246 /// Return true if the idiom is detected in the loop. 2247 /// 2248 /// The core idiom we are trying to detect is: 2249 /// \code 2250 /// entry: 2251 /// <...> 2252 /// %bitmask = shl i32 1, %bitpos 2253 /// br label %loop 2254 /// 2255 /// loop: 2256 /// %x.curr = phi i32 [ %x, %entry ], [ %x.next, %loop ] 2257 /// %x.curr.bitmasked = and i32 %x.curr, %bitmask 2258 /// %x.curr.isbitunset = icmp eq i32 %x.curr.bitmasked, 0 2259 /// %x.next = shl i32 %x.curr, 1 2260 /// <...> 2261 /// br i1 %x.curr.isbitunset, label %loop, label %end 2262 /// 2263 /// end: 2264 /// %x.curr.res = phi i32 [ %x.curr, %loop ] <...> 2265 /// %x.next.res = phi i32 [ %x.next, %loop ] <...> 2266 /// <...> 2267 /// \endcode 2268 static bool detectShiftUntilBitTestIdiom(Loop *CurLoop, Value *&BaseX, 2269 Value *&BitMask, Value *&BitPos, 2270 Value *&CurrX, Instruction *&NextX) { 2271 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2272 " Performing shift-until-bittest idiom detection.\n"); 2273 2274 // Give up if the loop has multiple blocks or multiple backedges. 2275 if (CurLoop->getNumBlocks() != 1 || CurLoop->getNumBackEdges() != 1) { 2276 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad block/backedge count.\n"); 2277 return false; 2278 } 2279 2280 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2281 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2282 assert(LoopPreheaderBB && "There is always a loop preheader."); 2283 2284 using namespace PatternMatch; 2285 2286 // Step 1: Check if the loop backedge is in desirable form. 2287 2288 ICmpInst::Predicate Pred; 2289 Value *CmpLHS, *CmpRHS; 2290 BasicBlock *TrueBB, *FalseBB; 2291 if (!match(LoopHeaderBB->getTerminator(), 2292 m_Br(m_ICmp(Pred, m_Value(CmpLHS), m_Value(CmpRHS)), 2293 m_BasicBlock(TrueBB), m_BasicBlock(FalseBB)))) { 2294 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge structure.\n"); 2295 return false; 2296 } 2297 2298 // Step 2: Check if the backedge's condition is in desirable form. 2299 2300 auto MatchVariableBitMask = [&]() { 2301 return ICmpInst::isEquality(Pred) && match(CmpRHS, m_Zero()) && 2302 match(CmpLHS, 2303 m_c_And(m_Value(CurrX), 2304 m_CombineAnd( 2305 m_Value(BitMask), 2306 m_LoopInvariant(m_Shl(m_One(), m_Value(BitPos)), 2307 CurLoop)))); 2308 }; 2309 auto MatchConstantBitMask = [&]() { 2310 return ICmpInst::isEquality(Pred) && match(CmpRHS, m_Zero()) && 2311 match(CmpLHS, m_And(m_Value(CurrX), 2312 m_CombineAnd(m_Value(BitMask), m_Power2()))) && 2313 (BitPos = ConstantExpr::getExactLogBase2(cast<Constant>(BitMask))); 2314 }; 2315 auto MatchDecomposableConstantBitMask = [&]() { 2316 APInt Mask; 2317 return llvm::decomposeBitTestICmp(CmpLHS, CmpRHS, Pred, CurrX, Mask) && 2318 ICmpInst::isEquality(Pred) && Mask.isPowerOf2() && 2319 (BitMask = ConstantInt::get(CurrX->getType(), Mask)) && 2320 (BitPos = ConstantInt::get(CurrX->getType(), Mask.logBase2())); 2321 }; 2322 2323 if (!MatchVariableBitMask() && !MatchConstantBitMask() && 2324 !MatchDecomposableConstantBitMask()) { 2325 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge comparison.\n"); 2326 return false; 2327 } 2328 2329 // Step 3: Check if the recurrence is in desirable form. 2330 auto *CurrXPN = dyn_cast<PHINode>(CurrX); 2331 if (!CurrXPN || CurrXPN->getParent() != LoopHeaderBB) { 2332 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Not an expected PHI node.\n"); 2333 return false; 2334 } 2335 2336 BaseX = CurrXPN->getIncomingValueForBlock(LoopPreheaderBB); 2337 NextX = 2338 dyn_cast<Instruction>(CurrXPN->getIncomingValueForBlock(LoopHeaderBB)); 2339 2340 assert(CurLoop->isLoopInvariant(BaseX) && 2341 "Expected BaseX to be avaliable in the preheader!"); 2342 2343 if (!NextX || !match(NextX, m_Shl(m_Specific(CurrX), m_One()))) { 2344 // FIXME: support right-shift? 2345 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad recurrence.\n"); 2346 return false; 2347 } 2348 2349 // Step 4: Check if the backedge's destinations are in desirable form. 2350 2351 assert(ICmpInst::isEquality(Pred) && 2352 "Should only get equality predicates here."); 2353 2354 // cmp-br is commutative, so canonicalize to a single variant. 2355 if (Pred != ICmpInst::Predicate::ICMP_EQ) { 2356 Pred = ICmpInst::getInversePredicate(Pred); 2357 std::swap(TrueBB, FalseBB); 2358 } 2359 2360 // We expect to exit loop when comparison yields false, 2361 // so when it yields true we should branch back to loop header. 2362 if (TrueBB != LoopHeaderBB) { 2363 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge flow.\n"); 2364 return false; 2365 } 2366 2367 // Okay, idiom checks out. 2368 return true; 2369 } 2370 2371 /// Look for the following loop: 2372 /// \code 2373 /// entry: 2374 /// <...> 2375 /// %bitmask = shl i32 1, %bitpos 2376 /// br label %loop 2377 /// 2378 /// loop: 2379 /// %x.curr = phi i32 [ %x, %entry ], [ %x.next, %loop ] 2380 /// %x.curr.bitmasked = and i32 %x.curr, %bitmask 2381 /// %x.curr.isbitunset = icmp eq i32 %x.curr.bitmasked, 0 2382 /// %x.next = shl i32 %x.curr, 1 2383 /// <...> 2384 /// br i1 %x.curr.isbitunset, label %loop, label %end 2385 /// 2386 /// end: 2387 /// %x.curr.res = phi i32 [ %x.curr, %loop ] <...> 2388 /// %x.next.res = phi i32 [ %x.next, %loop ] <...> 2389 /// <...> 2390 /// \endcode 2391 /// 2392 /// And transform it into: 2393 /// \code 2394 /// entry: 2395 /// %bitmask = shl i32 1, %bitpos 2396 /// %lowbitmask = add i32 %bitmask, -1 2397 /// %mask = or i32 %lowbitmask, %bitmask 2398 /// %x.masked = and i32 %x, %mask 2399 /// %x.masked.numleadingzeros = call i32 @llvm.ctlz.i32(i32 %x.masked, 2400 /// i1 true) 2401 /// %x.masked.numactivebits = sub i32 32, %x.masked.numleadingzeros 2402 /// %x.masked.leadingonepos = add i32 %x.masked.numactivebits, -1 2403 /// %backedgetakencount = sub i32 %bitpos, %x.masked.leadingonepos 2404 /// %tripcount = add i32 %backedgetakencount, 1 2405 /// %x.curr = shl i32 %x, %backedgetakencount 2406 /// %x.next = shl i32 %x, %tripcount 2407 /// br label %loop 2408 /// 2409 /// loop: 2410 /// %loop.iv = phi i32 [ 0, %entry ], [ %loop.iv.next, %loop ] 2411 /// %loop.iv.next = add nuw i32 %loop.iv, 1 2412 /// %loop.ivcheck = icmp eq i32 %loop.iv.next, %tripcount 2413 /// <...> 2414 /// br i1 %loop.ivcheck, label %end, label %loop 2415 /// 2416 /// end: 2417 /// %x.curr.res = phi i32 [ %x.curr, %loop ] <...> 2418 /// %x.next.res = phi i32 [ %x.next, %loop ] <...> 2419 /// <...> 2420 /// \endcode 2421 bool LoopIdiomRecognize::recognizeShiftUntilBitTest() { 2422 bool MadeChange = false; 2423 2424 Value *X, *BitMask, *BitPos, *XCurr; 2425 Instruction *XNext; 2426 if (!detectShiftUntilBitTestIdiom(CurLoop, X, BitMask, BitPos, XCurr, 2427 XNext)) { 2428 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2429 " shift-until-bittest idiom detection failed.\n"); 2430 return MadeChange; 2431 } 2432 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-bittest idiom detected!\n"); 2433 2434 // Ok, it is the idiom we were looking for, we *could* transform this loop, 2435 // but is it profitable to transform? 2436 2437 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2438 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2439 assert(LoopPreheaderBB && "There is always a loop preheader."); 2440 2441 BasicBlock *SuccessorBB = CurLoop->getExitBlock(); 2442 assert(SuccessorBB && "There is only a single successor."); 2443 2444 IRBuilder<> Builder(LoopPreheaderBB->getTerminator()); 2445 Builder.SetCurrentDebugLocation(cast<Instruction>(XCurr)->getDebugLoc()); 2446 2447 Intrinsic::ID IntrID = Intrinsic::ctlz; 2448 Type *Ty = X->getType(); 2449 unsigned Bitwidth = Ty->getScalarSizeInBits(); 2450 2451 TargetTransformInfo::TargetCostKind CostKind = 2452 TargetTransformInfo::TCK_SizeAndLatency; 2453 2454 // The rewrite is considered to be unprofitable iff and only iff the 2455 // intrinsic/shift we'll use are not cheap. Note that we are okay with *just* 2456 // making the loop countable, even if nothing else changes. 2457 IntrinsicCostAttributes Attrs( 2458 IntrID, Ty, {UndefValue::get(Ty), /*is_zero_undef=*/Builder.getTrue()}); 2459 InstructionCost Cost = TTI->getIntrinsicInstrCost(Attrs, CostKind); 2460 if (Cost > TargetTransformInfo::TCC_Basic) { 2461 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2462 " Intrinsic is too costly, not beneficial\n"); 2463 return MadeChange; 2464 } 2465 if (TTI->getArithmeticInstrCost(Instruction::Shl, Ty, CostKind) > 2466 TargetTransformInfo::TCC_Basic) { 2467 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Shift is too costly, not beneficial\n"); 2468 return MadeChange; 2469 } 2470 2471 // Ok, transform appears worthwhile. 2472 MadeChange = true; 2473 2474 // Step 1: Compute the loop trip count. 2475 2476 Value *LowBitMask = Builder.CreateAdd(BitMask, Constant::getAllOnesValue(Ty), 2477 BitPos->getName() + ".lowbitmask"); 2478 Value *Mask = 2479 Builder.CreateOr(LowBitMask, BitMask, BitPos->getName() + ".mask"); 2480 Value *XMasked = Builder.CreateAnd(X, Mask, X->getName() + ".masked"); 2481 CallInst *XMaskedNumLeadingZeros = Builder.CreateIntrinsic( 2482 IntrID, Ty, {XMasked, /*is_zero_undef=*/Builder.getTrue()}, 2483 /*FMFSource=*/nullptr, XMasked->getName() + ".numleadingzeros"); 2484 Value *XMaskedNumActiveBits = Builder.CreateSub( 2485 ConstantInt::get(Ty, Ty->getScalarSizeInBits()), XMaskedNumLeadingZeros, 2486 XMasked->getName() + ".numactivebits", /*HasNUW=*/true, 2487 /*HasNSW=*/Bitwidth != 2); 2488 Value *XMaskedLeadingOnePos = 2489 Builder.CreateAdd(XMaskedNumActiveBits, Constant::getAllOnesValue(Ty), 2490 XMasked->getName() + ".leadingonepos", /*HasNUW=*/false, 2491 /*HasNSW=*/Bitwidth > 2); 2492 2493 Value *LoopBackedgeTakenCount = Builder.CreateSub( 2494 BitPos, XMaskedLeadingOnePos, CurLoop->getName() + ".backedgetakencount", 2495 /*HasNUW=*/true, /*HasNSW=*/true); 2496 // We know loop's backedge-taken count, but what's loop's trip count? 2497 // Note that while NUW is always safe, while NSW is only for bitwidths != 2. 2498 Value *LoopTripCount = 2499 Builder.CreateAdd(LoopBackedgeTakenCount, ConstantInt::get(Ty, 1), 2500 CurLoop->getName() + ".tripcount", /*HasNUW=*/true, 2501 /*HasNSW=*/Bitwidth != 2); 2502 2503 // Step 2: Compute the recurrence's final value without a loop. 2504 2505 // NewX is always safe to compute, because `LoopBackedgeTakenCount` 2506 // will always be smaller than `bitwidth(X)`, i.e. we never get poison. 2507 Value *NewX = Builder.CreateShl(X, LoopBackedgeTakenCount); 2508 NewX->takeName(XCurr); 2509 if (auto *I = dyn_cast<Instruction>(NewX)) 2510 I->copyIRFlags(XNext, /*IncludeWrapFlags=*/true); 2511 2512 Value *NewXNext; 2513 // Rewriting XNext is more complicated, however, because `X << LoopTripCount` 2514 // will be poison iff `LoopTripCount == bitwidth(X)` (which will happen 2515 // iff `BitPos` is `bitwidth(x) - 1` and `X` is `1`). So unless we know 2516 // that isn't the case, we'll need to emit an alternative, safe IR. 2517 if (XNext->hasNoSignedWrap() || XNext->hasNoUnsignedWrap() || 2518 PatternMatch::match( 2519 BitPos, PatternMatch::m_SpecificInt_ICMP( 2520 ICmpInst::ICMP_NE, APInt(Ty->getScalarSizeInBits(), 2521 Ty->getScalarSizeInBits() - 1)))) 2522 NewXNext = Builder.CreateShl(X, LoopTripCount); 2523 else { 2524 // Otherwise, just additionally shift by one. It's the smallest solution, 2525 // alternatively, we could check that NewX is INT_MIN (or BitPos is ) 2526 // and select 0 instead. 2527 NewXNext = Builder.CreateShl(NewX, ConstantInt::get(Ty, 1)); 2528 } 2529 2530 NewXNext->takeName(XNext); 2531 if (auto *I = dyn_cast<Instruction>(NewXNext)) 2532 I->copyIRFlags(XNext, /*IncludeWrapFlags=*/true); 2533 2534 // Step 3: Adjust the successor basic block to recieve the computed 2535 // recurrence's final value instead of the recurrence itself. 2536 2537 XCurr->replaceUsesOutsideBlock(NewX, LoopHeaderBB); 2538 XNext->replaceUsesOutsideBlock(NewXNext, LoopHeaderBB); 2539 2540 // Step 4: Rewrite the loop into a countable form, with canonical IV. 2541 2542 // The new canonical induction variable. 2543 Builder.SetInsertPoint(&LoopHeaderBB->front()); 2544 auto *IV = Builder.CreatePHI(Ty, 2, CurLoop->getName() + ".iv"); 2545 2546 // The induction itself. 2547 // Note that while NUW is always safe, while NSW is only for bitwidths != 2. 2548 Builder.SetInsertPoint(LoopHeaderBB->getTerminator()); 2549 auto *IVNext = 2550 Builder.CreateAdd(IV, ConstantInt::get(Ty, 1), IV->getName() + ".next", 2551 /*HasNUW=*/true, /*HasNSW=*/Bitwidth != 2); 2552 2553 // The loop trip count check. 2554 auto *IVCheck = Builder.CreateICmpEQ(IVNext, LoopTripCount, 2555 CurLoop->getName() + ".ivcheck"); 2556 Builder.CreateCondBr(IVCheck, SuccessorBB, LoopHeaderBB); 2557 LoopHeaderBB->getTerminator()->eraseFromParent(); 2558 2559 // Populate the IV PHI. 2560 IV->addIncoming(ConstantInt::get(Ty, 0), LoopPreheaderBB); 2561 IV->addIncoming(IVNext, LoopHeaderBB); 2562 2563 // Step 5: Forget the "non-computable" trip-count SCEV associated with the 2564 // loop. The loop would otherwise not be deleted even if it becomes empty. 2565 2566 SE->forgetLoop(CurLoop); 2567 2568 // Other passes will take care of actually deleting the loop if possible. 2569 2570 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-bittest idiom optimized!\n"); 2571 2572 ++NumShiftUntilBitTest; 2573 return MadeChange; 2574 } 2575 2576 /// Return true if the idiom is detected in the loop. 2577 /// 2578 /// The core idiom we are trying to detect is: 2579 /// \code 2580 /// entry: 2581 /// <...> 2582 /// %start = <...> 2583 /// %extraoffset = <...> 2584 /// <...> 2585 /// br label %for.cond 2586 /// 2587 /// loop: 2588 /// %iv = phi i8 [ %start, %entry ], [ %iv.next, %for.cond ] 2589 /// %nbits = add nsw i8 %iv, %extraoffset 2590 /// %val.shifted = {{l,a}shr,shl} i8 %val, %nbits 2591 /// %val.shifted.iszero = icmp eq i8 %val.shifted, 0 2592 /// %iv.next = add i8 %iv, 1 2593 /// <...> 2594 /// br i1 %val.shifted.iszero, label %end, label %loop 2595 /// 2596 /// end: 2597 /// %iv.res = phi i8 [ %iv, %loop ] <...> 2598 /// %nbits.res = phi i8 [ %nbits, %loop ] <...> 2599 /// %val.shifted.res = phi i8 [ %val.shifted, %loop ] <...> 2600 /// %val.shifted.iszero.res = phi i1 [ %val.shifted.iszero, %loop ] <...> 2601 /// %iv.next.res = phi i8 [ %iv.next, %loop ] <...> 2602 /// <...> 2603 /// \endcode 2604 static bool detectShiftUntilZeroIdiom(Loop *CurLoop, ScalarEvolution *SE, 2605 Instruction *&ValShiftedIsZero, 2606 Intrinsic::ID &IntrinID, Instruction *&IV, 2607 Value *&Start, Value *&Val, 2608 const SCEV *&ExtraOffsetExpr, 2609 bool &InvertedCond) { 2610 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2611 " Performing shift-until-zero idiom detection.\n"); 2612 2613 // Give up if the loop has multiple blocks or multiple backedges. 2614 if (CurLoop->getNumBlocks() != 1 || CurLoop->getNumBackEdges() != 1) { 2615 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad block/backedge count.\n"); 2616 return false; 2617 } 2618 2619 Instruction *ValShifted, *NBits, *IVNext; 2620 Value *ExtraOffset; 2621 2622 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2623 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2624 assert(LoopPreheaderBB && "There is always a loop preheader."); 2625 2626 using namespace PatternMatch; 2627 2628 // Step 1: Check if the loop backedge, condition is in desirable form. 2629 2630 ICmpInst::Predicate Pred; 2631 BasicBlock *TrueBB, *FalseBB; 2632 if (!match(LoopHeaderBB->getTerminator(), 2633 m_Br(m_Instruction(ValShiftedIsZero), m_BasicBlock(TrueBB), 2634 m_BasicBlock(FalseBB))) || 2635 !match(ValShiftedIsZero, 2636 m_ICmp(Pred, m_Instruction(ValShifted), m_Zero())) || 2637 !ICmpInst::isEquality(Pred)) { 2638 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge structure.\n"); 2639 return false; 2640 } 2641 2642 // Step 2: Check if the comparison's operand is in desirable form. 2643 // FIXME: Val could be a one-input PHI node, which we should look past. 2644 if (!match(ValShifted, m_Shift(m_LoopInvariant(m_Value(Val), CurLoop), 2645 m_Instruction(NBits)))) { 2646 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad comparisons value computation.\n"); 2647 return false; 2648 } 2649 IntrinID = ValShifted->getOpcode() == Instruction::Shl ? Intrinsic::cttz 2650 : Intrinsic::ctlz; 2651 2652 // Step 3: Check if the shift amount is in desirable form. 2653 2654 if (match(NBits, m_c_Add(m_Instruction(IV), 2655 m_LoopInvariant(m_Value(ExtraOffset), CurLoop))) && 2656 (NBits->hasNoSignedWrap() || NBits->hasNoUnsignedWrap())) 2657 ExtraOffsetExpr = SE->getNegativeSCEV(SE->getSCEV(ExtraOffset)); 2658 else if (match(NBits, 2659 m_Sub(m_Instruction(IV), 2660 m_LoopInvariant(m_Value(ExtraOffset), CurLoop))) && 2661 NBits->hasNoSignedWrap()) 2662 ExtraOffsetExpr = SE->getSCEV(ExtraOffset); 2663 else { 2664 IV = NBits; 2665 ExtraOffsetExpr = SE->getZero(NBits->getType()); 2666 } 2667 2668 // Step 4: Check if the recurrence is in desirable form. 2669 auto *IVPN = dyn_cast<PHINode>(IV); 2670 if (!IVPN || IVPN->getParent() != LoopHeaderBB) { 2671 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Not an expected PHI node.\n"); 2672 return false; 2673 } 2674 2675 Start = IVPN->getIncomingValueForBlock(LoopPreheaderBB); 2676 IVNext = dyn_cast<Instruction>(IVPN->getIncomingValueForBlock(LoopHeaderBB)); 2677 2678 if (!IVNext || !match(IVNext, m_Add(m_Specific(IVPN), m_One()))) { 2679 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad recurrence.\n"); 2680 return false; 2681 } 2682 2683 // Step 4: Check if the backedge's destinations are in desirable form. 2684 2685 assert(ICmpInst::isEquality(Pred) && 2686 "Should only get equality predicates here."); 2687 2688 // cmp-br is commutative, so canonicalize to a single variant. 2689 InvertedCond = Pred != ICmpInst::Predicate::ICMP_EQ; 2690 if (InvertedCond) { 2691 Pred = ICmpInst::getInversePredicate(Pred); 2692 std::swap(TrueBB, FalseBB); 2693 } 2694 2695 // We expect to exit loop when comparison yields true, 2696 // so when it yields false we should branch back to loop header. 2697 if (FalseBB != LoopHeaderBB) { 2698 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Bad backedge flow.\n"); 2699 return false; 2700 } 2701 2702 // The new, countable, loop will certainly only run a known number of 2703 // iterations, It won't be infinite. But the old loop might be infinite 2704 // under certain conditions. For logical shifts, the value will become zero 2705 // after at most bitwidth(%Val) loop iterations. However, for arithmetic 2706 // right-shift, iff the sign bit was set, the value will never become zero, 2707 // and the loop may never finish. 2708 if (ValShifted->getOpcode() == Instruction::AShr && 2709 !isMustProgress(CurLoop) && !SE->isKnownNonNegative(SE->getSCEV(Val))) { 2710 LLVM_DEBUG(dbgs() << DEBUG_TYPE " Can not prove the loop is finite.\n"); 2711 return false; 2712 } 2713 2714 // Okay, idiom checks out. 2715 return true; 2716 } 2717 2718 /// Look for the following loop: 2719 /// \code 2720 /// entry: 2721 /// <...> 2722 /// %start = <...> 2723 /// %extraoffset = <...> 2724 /// <...> 2725 /// br label %for.cond 2726 /// 2727 /// loop: 2728 /// %iv = phi i8 [ %start, %entry ], [ %iv.next, %for.cond ] 2729 /// %nbits = add nsw i8 %iv, %extraoffset 2730 /// %val.shifted = {{l,a}shr,shl} i8 %val, %nbits 2731 /// %val.shifted.iszero = icmp eq i8 %val.shifted, 0 2732 /// %iv.next = add i8 %iv, 1 2733 /// <...> 2734 /// br i1 %val.shifted.iszero, label %end, label %loop 2735 /// 2736 /// end: 2737 /// %iv.res = phi i8 [ %iv, %loop ] <...> 2738 /// %nbits.res = phi i8 [ %nbits, %loop ] <...> 2739 /// %val.shifted.res = phi i8 [ %val.shifted, %loop ] <...> 2740 /// %val.shifted.iszero.res = phi i1 [ %val.shifted.iszero, %loop ] <...> 2741 /// %iv.next.res = phi i8 [ %iv.next, %loop ] <...> 2742 /// <...> 2743 /// \endcode 2744 /// 2745 /// And transform it into: 2746 /// \code 2747 /// entry: 2748 /// <...> 2749 /// %start = <...> 2750 /// %extraoffset = <...> 2751 /// <...> 2752 /// %val.numleadingzeros = call i8 @llvm.ct{l,t}z.i8(i8 %val, i1 0) 2753 /// %val.numactivebits = sub i8 8, %val.numleadingzeros 2754 /// %extraoffset.neg = sub i8 0, %extraoffset 2755 /// %tmp = add i8 %val.numactivebits, %extraoffset.neg 2756 /// %iv.final = call i8 @llvm.smax.i8(i8 %tmp, i8 %start) 2757 /// %loop.tripcount = sub i8 %iv.final, %start 2758 /// br label %loop 2759 /// 2760 /// loop: 2761 /// %loop.iv = phi i8 [ 0, %entry ], [ %loop.iv.next, %loop ] 2762 /// %loop.iv.next = add i8 %loop.iv, 1 2763 /// %loop.ivcheck = icmp eq i8 %loop.iv.next, %loop.tripcount 2764 /// %iv = add i8 %loop.iv, %start 2765 /// <...> 2766 /// br i1 %loop.ivcheck, label %end, label %loop 2767 /// 2768 /// end: 2769 /// %iv.res = phi i8 [ %iv.final, %loop ] <...> 2770 /// <...> 2771 /// \endcode 2772 bool LoopIdiomRecognize::recognizeShiftUntilZero() { 2773 bool MadeChange = false; 2774 2775 Instruction *ValShiftedIsZero; 2776 Intrinsic::ID IntrID; 2777 Instruction *IV; 2778 Value *Start, *Val; 2779 const SCEV *ExtraOffsetExpr; 2780 bool InvertedCond; 2781 if (!detectShiftUntilZeroIdiom(CurLoop, SE, ValShiftedIsZero, IntrID, IV, 2782 Start, Val, ExtraOffsetExpr, InvertedCond)) { 2783 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2784 " shift-until-zero idiom detection failed.\n"); 2785 return MadeChange; 2786 } 2787 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-zero idiom detected!\n"); 2788 2789 // Ok, it is the idiom we were looking for, we *could* transform this loop, 2790 // but is it profitable to transform? 2791 2792 BasicBlock *LoopHeaderBB = CurLoop->getHeader(); 2793 BasicBlock *LoopPreheaderBB = CurLoop->getLoopPreheader(); 2794 assert(LoopPreheaderBB && "There is always a loop preheader."); 2795 2796 BasicBlock *SuccessorBB = CurLoop->getExitBlock(); 2797 assert(SuccessorBB && "There is only a single successor."); 2798 2799 IRBuilder<> Builder(LoopPreheaderBB->getTerminator()); 2800 Builder.SetCurrentDebugLocation(IV->getDebugLoc()); 2801 2802 Type *Ty = Val->getType(); 2803 unsigned Bitwidth = Ty->getScalarSizeInBits(); 2804 2805 TargetTransformInfo::TargetCostKind CostKind = 2806 TargetTransformInfo::TCK_SizeAndLatency; 2807 2808 // The rewrite is considered to be unprofitable iff and only iff the 2809 // intrinsic we'll use are not cheap. Note that we are okay with *just* 2810 // making the loop countable, even if nothing else changes. 2811 IntrinsicCostAttributes Attrs( 2812 IntrID, Ty, {UndefValue::get(Ty), /*is_zero_undef=*/Builder.getFalse()}); 2813 InstructionCost Cost = TTI->getIntrinsicInstrCost(Attrs, CostKind); 2814 if (Cost > TargetTransformInfo::TCC_Basic) { 2815 LLVM_DEBUG(dbgs() << DEBUG_TYPE 2816 " Intrinsic is too costly, not beneficial\n"); 2817 return MadeChange; 2818 } 2819 2820 // Ok, transform appears worthwhile. 2821 MadeChange = true; 2822 2823 bool OffsetIsZero = false; 2824 if (auto *ExtraOffsetExprC = dyn_cast<SCEVConstant>(ExtraOffsetExpr)) 2825 OffsetIsZero = ExtraOffsetExprC->isZero(); 2826 2827 // Step 1: Compute the loop's final IV value / trip count. 2828 2829 CallInst *ValNumLeadingZeros = Builder.CreateIntrinsic( 2830 IntrID, Ty, {Val, /*is_zero_undef=*/Builder.getFalse()}, 2831 /*FMFSource=*/nullptr, Val->getName() + ".numleadingzeros"); 2832 Value *ValNumActiveBits = Builder.CreateSub( 2833 ConstantInt::get(Ty, Ty->getScalarSizeInBits()), ValNumLeadingZeros, 2834 Val->getName() + ".numactivebits", /*HasNUW=*/true, 2835 /*HasNSW=*/Bitwidth != 2); 2836 2837 SCEVExpander Expander(*SE, *DL, "loop-idiom"); 2838 Expander.setInsertPoint(&*Builder.GetInsertPoint()); 2839 Value *ExtraOffset = Expander.expandCodeFor(ExtraOffsetExpr); 2840 2841 Value *ValNumActiveBitsOffset = Builder.CreateAdd( 2842 ValNumActiveBits, ExtraOffset, ValNumActiveBits->getName() + ".offset", 2843 /*HasNUW=*/OffsetIsZero, /*HasNSW=*/true); 2844 Value *IVFinal = Builder.CreateIntrinsic(Intrinsic::smax, {Ty}, 2845 {ValNumActiveBitsOffset, Start}, 2846 /*FMFSource=*/nullptr, "iv.final"); 2847 2848 auto *LoopBackedgeTakenCount = cast<Instruction>(Builder.CreateSub( 2849 IVFinal, Start, CurLoop->getName() + ".backedgetakencount", 2850 /*HasNUW=*/OffsetIsZero, /*HasNSW=*/true)); 2851 // FIXME: or when the offset was `add nuw` 2852 2853 // We know loop's backedge-taken count, but what's loop's trip count? 2854 Value *LoopTripCount = 2855 Builder.CreateAdd(LoopBackedgeTakenCount, ConstantInt::get(Ty, 1), 2856 CurLoop->getName() + ".tripcount", /*HasNUW=*/true, 2857 /*HasNSW=*/Bitwidth != 2); 2858 2859 // Step 2: Adjust the successor basic block to recieve the original 2860 // induction variable's final value instead of the orig. IV itself. 2861 2862 IV->replaceUsesOutsideBlock(IVFinal, LoopHeaderBB); 2863 2864 // Step 3: Rewrite the loop into a countable form, with canonical IV. 2865 2866 // The new canonical induction variable. 2867 Builder.SetInsertPoint(&LoopHeaderBB->front()); 2868 auto *CIV = Builder.CreatePHI(Ty, 2, CurLoop->getName() + ".iv"); 2869 2870 // The induction itself. 2871 Builder.SetInsertPoint(LoopHeaderBB->getFirstNonPHI()); 2872 auto *CIVNext = 2873 Builder.CreateAdd(CIV, ConstantInt::get(Ty, 1), CIV->getName() + ".next", 2874 /*HasNUW=*/true, /*HasNSW=*/Bitwidth != 2); 2875 2876 // The loop trip count check. 2877 auto *CIVCheck = Builder.CreateICmpEQ(CIVNext, LoopTripCount, 2878 CurLoop->getName() + ".ivcheck"); 2879 auto *NewIVCheck = CIVCheck; 2880 if (InvertedCond) { 2881 NewIVCheck = Builder.CreateNot(CIVCheck); 2882 NewIVCheck->takeName(ValShiftedIsZero); 2883 } 2884 2885 // The original IV, but rebased to be an offset to the CIV. 2886 auto *IVDePHId = Builder.CreateAdd(CIV, Start, "", /*HasNUW=*/false, 2887 /*HasNSW=*/true); // FIXME: what about NUW? 2888 IVDePHId->takeName(IV); 2889 2890 // The loop terminator. 2891 Builder.SetInsertPoint(LoopHeaderBB->getTerminator()); 2892 Builder.CreateCondBr(CIVCheck, SuccessorBB, LoopHeaderBB); 2893 LoopHeaderBB->getTerminator()->eraseFromParent(); 2894 2895 // Populate the IV PHI. 2896 CIV->addIncoming(ConstantInt::get(Ty, 0), LoopPreheaderBB); 2897 CIV->addIncoming(CIVNext, LoopHeaderBB); 2898 2899 // Step 4: Forget the "non-computable" trip-count SCEV associated with the 2900 // loop. The loop would otherwise not be deleted even if it becomes empty. 2901 2902 SE->forgetLoop(CurLoop); 2903 2904 // Step 5: Try to cleanup the loop's body somewhat. 2905 IV->replaceAllUsesWith(IVDePHId); 2906 IV->eraseFromParent(); 2907 2908 ValShiftedIsZero->replaceAllUsesWith(NewIVCheck); 2909 ValShiftedIsZero->eraseFromParent(); 2910 2911 // Other passes will take care of actually deleting the loop if possible. 2912 2913 LLVM_DEBUG(dbgs() << DEBUG_TYPE " shift-until-zero idiom optimized!\n"); 2914 2915 ++NumShiftUntilZero; 2916 return MadeChange; 2917 } 2918