Use SCEVAddRecExpr::isAffine.
[oota-llvm.git] / lib / Analysis / ScalarEvolutionExpander.cpp
1 //===- ScalarEvolutionExpander.cpp - Scalar Evolution Analysis --*- C++ -*-===//
2 //
3 //                     The LLVM Compiler Infrastructure
4 //
5 // This file is distributed under the University of Illinois Open Source
6 // License. See LICENSE.TXT for details.
7 //
8 //===----------------------------------------------------------------------===//
9 //
10 // This file contains the implementation of the scalar evolution expander,
11 // which is used to generate the code corresponding to a given scalar evolution
12 // expression.
13 //
14 //===----------------------------------------------------------------------===//
15
16 #include "llvm/Analysis/ScalarEvolutionExpander.h"
17 #include "llvm/Analysis/LoopInfo.h"
18 using namespace llvm;
19
20 /// InsertCastOfTo - Insert a cast of V to the specified type, doing what
21 /// we can to share the casts.
22 Value *SCEVExpander::InsertCastOfTo(Instruction::CastOps opcode, Value *V, 
23                                     const Type *Ty) {
24   // FIXME: keep track of the cast instruction.
25   if (Constant *C = dyn_cast<Constant>(V))
26     return ConstantExpr::getCast(opcode, C, Ty);
27   
28   if (Argument *A = dyn_cast<Argument>(V)) {
29     // Check to see if there is already a cast!
30     for (Value::use_iterator UI = A->use_begin(), E = A->use_end();
31          UI != E; ++UI) {
32       if ((*UI)->getType() == Ty)
33         if (CastInst *CI = dyn_cast<CastInst>(cast<Instruction>(*UI)))
34           if (CI->getOpcode() == opcode) {
35             // If the cast isn't the first instruction of the function, move it.
36             if (BasicBlock::iterator(CI) != 
37                 A->getParent()->getEntryBlock().begin()) {
38               CI->moveBefore(A->getParent()->getEntryBlock().begin());
39             }
40             return CI;
41           }
42     }
43     return CastInst::Create(opcode, V, Ty, V->getName(), 
44                             A->getParent()->getEntryBlock().begin());
45   }
46
47   Instruction *I = cast<Instruction>(V);
48
49   // Check to see if there is already a cast.  If there is, use it.
50   for (Value::use_iterator UI = I->use_begin(), E = I->use_end();
51        UI != E; ++UI) {
52     if ((*UI)->getType() == Ty)
53       if (CastInst *CI = dyn_cast<CastInst>(cast<Instruction>(*UI)))
54         if (CI->getOpcode() == opcode) {
55           BasicBlock::iterator It = I; ++It;
56           if (isa<InvokeInst>(I))
57             It = cast<InvokeInst>(I)->getNormalDest()->begin();
58           while (isa<PHINode>(It)) ++It;
59           if (It != BasicBlock::iterator(CI)) {
60             // Splice the cast immediately after the operand in question.
61             CI->moveBefore(It);
62           }
63           return CI;
64         }
65   }
66   BasicBlock::iterator IP = I; ++IP;
67   if (InvokeInst *II = dyn_cast<InvokeInst>(I))
68     IP = II->getNormalDest()->begin();
69   while (isa<PHINode>(IP)) ++IP;
70   return CastInst::Create(opcode, V, Ty, V->getName(), IP);
71 }
72
73 /// InsertBinop - Insert the specified binary operator, doing a small amount
74 /// of work to avoid inserting an obviously redundant operation.
75 Value *SCEVExpander::InsertBinop(Instruction::BinaryOps Opcode, Value *LHS,
76                                  Value *RHS, Instruction *InsertPt) {
77   // Fold a binop with constant operands.
78   if (Constant *CLHS = dyn_cast<Constant>(LHS))
79     if (Constant *CRHS = dyn_cast<Constant>(RHS))
80       return ConstantExpr::get(Opcode, CLHS, CRHS);
81
82   // Do a quick scan to see if we have this binop nearby.  If so, reuse it.
83   unsigned ScanLimit = 6;
84   BasicBlock::iterator BlockBegin = InsertPt->getParent()->begin();
85   if (InsertPt != BlockBegin) {
86     // Scanning starts from the last instruction before InsertPt.
87     BasicBlock::iterator IP = InsertPt;
88     --IP;
89     for (; ScanLimit; --IP, --ScanLimit) {
90       if (BinaryOperator *BinOp = dyn_cast<BinaryOperator>(IP))
91         if (BinOp->getOpcode() == Opcode && BinOp->getOperand(0) == LHS &&
92             BinOp->getOperand(1) == RHS)
93           return BinOp;
94       if (IP == BlockBegin) break;
95     }
96   }
97   
98   // If we haven't found this binop, insert it.
99   return BinaryOperator::Create(Opcode, LHS, RHS, "tmp", InsertPt);
100 }
101
102 Value *SCEVExpander::visitAddExpr(SCEVAddExpr *S) {
103   Value *V = expand(S->getOperand(S->getNumOperands()-1));
104
105   // Emit a bunch of add instructions
106   for (int i = S->getNumOperands()-2; i >= 0; --i)
107     V = InsertBinop(Instruction::Add, V, expand(S->getOperand(i)),
108                     InsertPt);
109   return V;
110 }
111     
112 Value *SCEVExpander::visitMulExpr(SCEVMulExpr *S) {
113   int FirstOp = 0;  // Set if we should emit a subtract.
114   if (SCEVConstant *SC = dyn_cast<SCEVConstant>(S->getOperand(0)))
115     if (SC->getValue()->isAllOnesValue())
116       FirstOp = 1;
117
118   int i = S->getNumOperands()-2;
119   Value *V = expand(S->getOperand(i+1));
120
121   // Emit a bunch of multiply instructions
122   for (; i >= FirstOp; --i)
123     V = InsertBinop(Instruction::Mul, V, expand(S->getOperand(i)),
124                     InsertPt);
125   // -1 * ...  --->  0 - ...
126   if (FirstOp == 1)
127     V = InsertBinop(Instruction::Sub, Constant::getNullValue(V->getType()), V,
128                     InsertPt);
129   return V;
130 }
131
132 Value *SCEVExpander::visitAddRecExpr(SCEVAddRecExpr *S) {
133   const Type *Ty = S->getType();
134   const Loop *L = S->getLoop();
135   // We cannot yet do fp recurrences, e.g. the xform of {X,+,F} --> X+{0,+,F}
136   assert(Ty->isInteger() && "Cannot expand fp recurrences yet!");
137
138   // {X,+,F} --> X + {0,+,F}
139   if (!S->getStart()->isZero()) {
140     Value *Start = expand(S->getStart());
141     std::vector<SCEVHandle> NewOps(S->op_begin(), S->op_end());
142     NewOps[0] = SE.getIntegerSCEV(0, Ty);
143     Value *Rest = expand(SE.getAddRecExpr(NewOps, L));
144
145     // FIXME: look for an existing add to use.
146     return InsertBinop(Instruction::Add, Rest, Start, InsertPt);
147   }
148
149   // {0,+,1} --> Insert a canonical induction variable into the loop!
150   if (S->isAffine() &&
151       S->getOperand(1) == SE.getIntegerSCEV(1, Ty)) {
152     // Create and insert the PHI node for the induction variable in the
153     // specified loop.
154     BasicBlock *Header = L->getHeader();
155     PHINode *PN = PHINode::Create(Ty, "indvar", Header->begin());
156     PN->addIncoming(Constant::getNullValue(Ty), L->getLoopPreheader());
157
158     pred_iterator HPI = pred_begin(Header);
159     assert(HPI != pred_end(Header) && "Loop with zero preds???");
160     if (!L->contains(*HPI)) ++HPI;
161     assert(HPI != pred_end(Header) && L->contains(*HPI) &&
162            "No backedge in loop?");
163
164     // Insert a unit add instruction right before the terminator corresponding
165     // to the back-edge.
166     Constant *One = ConstantInt::get(Ty, 1);
167     Instruction *Add = BinaryOperator::CreateAdd(PN, One, "indvar.next",
168                                                  (*HPI)->getTerminator());
169
170     pred_iterator PI = pred_begin(Header);
171     if (*PI == L->getLoopPreheader())
172       ++PI;
173     PN->addIncoming(Add, *PI);
174     return PN;
175   }
176
177   // Get the canonical induction variable I for this loop.
178   Value *I = getOrInsertCanonicalInductionVariable(L, Ty);
179
180   // If this is a simple linear addrec, emit it now as a special case.
181   if (S->isAffine()) {   // {0,+,F} --> i*F
182     Value *F = expand(S->getOperand(1));
183     
184     // IF the step is by one, just return the inserted IV.
185     if (ConstantInt *CI = dyn_cast<ConstantInt>(F))
186       if (CI->getValue() == 1)
187         return I;
188     
189     // If the insert point is directly inside of the loop, emit the multiply at
190     // the insert point.  Otherwise, L is a loop that is a parent of the insert
191     // point loop.  If we can, move the multiply to the outer most loop that it
192     // is safe to be in.
193     Instruction *MulInsertPt = InsertPt;
194     Loop *InsertPtLoop = LI.getLoopFor(MulInsertPt->getParent());
195     if (InsertPtLoop != L && InsertPtLoop &&
196         L->contains(InsertPtLoop->getHeader())) {
197       do {
198         // If we cannot hoist the multiply out of this loop, don't.
199         if (!InsertPtLoop->isLoopInvariant(F)) break;
200
201         BasicBlock *InsertPtLoopPH = InsertPtLoop->getLoopPreheader();
202
203         // If this loop hasn't got a preheader, we aren't able to hoist the
204         // multiply.
205         if (!InsertPtLoopPH)
206           break;
207
208         // Otherwise, move the insert point to the preheader.
209         MulInsertPt = InsertPtLoopPH->getTerminator();
210         InsertPtLoop = InsertPtLoop->getParentLoop();
211       } while (InsertPtLoop != L);
212     }
213     
214     return InsertBinop(Instruction::Mul, I, F, MulInsertPt);
215   }
216
217   // If this is a chain of recurrences, turn it into a closed form, using the
218   // folders, then expandCodeFor the closed form.  This allows the folders to
219   // simplify the expression without having to build a bunch of special code
220   // into this folder.
221   SCEVHandle IH = SE.getUnknown(I);   // Get I as a "symbolic" SCEV.
222
223   SCEVHandle V = S->evaluateAtIteration(IH, SE);
224   //cerr << "Evaluated: " << *this << "\n     to: " << *V << "\n";
225
226   return expand(V);
227 }
228
229 Value *SCEVExpander::visitTruncateExpr(SCEVTruncateExpr *S) {
230   Value *V = expand(S->getOperand());
231   return CastInst::createTruncOrBitCast(V, S->getType(), "tmp.", InsertPt);
232 }
233
234 Value *SCEVExpander::visitZeroExtendExpr(SCEVZeroExtendExpr *S) {
235   Value *V = expand(S->getOperand());
236   return CastInst::createZExtOrBitCast(V, S->getType(), "tmp.", InsertPt);
237 }
238
239 Value *SCEVExpander::visitSignExtendExpr(SCEVSignExtendExpr *S) {
240   Value *V = expand(S->getOperand());
241   return CastInst::createSExtOrBitCast(V, S->getType(), "tmp.", InsertPt);
242 }
243
244 Value *SCEVExpander::visitSMaxExpr(SCEVSMaxExpr *S) {
245   Value *LHS = expand(S->getOperand(0));
246   for (unsigned i = 1; i < S->getNumOperands(); ++i) {
247     Value *RHS = expand(S->getOperand(i));
248     Value *ICmp = new ICmpInst(ICmpInst::ICMP_SGT, LHS, RHS, "tmp", InsertPt);
249     LHS = SelectInst::Create(ICmp, LHS, RHS, "smax", InsertPt);
250   }
251   return LHS;
252 }
253
254 Value *SCEVExpander::visitUMaxExpr(SCEVUMaxExpr *S) {
255   Value *LHS = expand(S->getOperand(0));
256   for (unsigned i = 1; i < S->getNumOperands(); ++i) {
257     Value *RHS = expand(S->getOperand(i));
258     Value *ICmp = new ICmpInst(ICmpInst::ICMP_UGT, LHS, RHS, "tmp", InsertPt);
259     LHS = SelectInst::Create(ICmp, LHS, RHS, "umax", InsertPt);
260   }
261   return LHS;
262 }
263
264 Value *SCEVExpander::expandCodeFor(SCEVHandle SH, Instruction *IP) {
265   // Expand the code for this SCEV.
266   this->InsertPt = IP;
267   return expand(SH);
268 }
269
270 Value *SCEVExpander::expand(SCEV *S) {
271   // Check to see if we already expanded this.
272   std::map<SCEVHandle, Value*>::iterator I = InsertedExpressions.find(S);
273   if (I != InsertedExpressions.end())
274     return I->second;
275   
276   Value *V = visit(S);
277   InsertedExpressions[S] = V;
278   return V;
279 }