Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] Annotation of /branches/charisee/src/compiler/mid-to-low/helper.sml
ViewVC logotype

Annotation of /branches/charisee/src/compiler/mid-to-low/helper.sml

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3195 - (view) (download)

1 : cchiw 2859 (*Helper functions
2 :     *)
3 :     structure Helper = struct
4 :     local
5 :    
6 :     structure DstIL = LowIL
7 :     structure DstTy = LowILTypes
8 :     structure DstOp = LowOps
9 :     structure Var = LowIL.Var
10 :     structure SrcIL = MidIL
11 :     structure SrcOp = MidOps
12 :     structure E = Ein
13 :     structure LowToS= LowToString
14 : cchiw 2870 structure FL=FloatLit
15 : jhr 2924 structure IMap = IntRedBlackMap
16 : cchiw 2859 in
17 :    
18 : cchiw 3194 val testing=1
19 : cchiw 2859 val bV= ref 0
20 :     fun err str=raise Fail(str)
21 :     val realTy=DstTy.TensorTy []
22 :     val intTy=DstTy.intTy
23 :     fun iTos e1=Int.toString e1
24 :     fun iToss es=String.concat(List.map iTos es)
25 :     fun testp n= (case testing
26 : jhr 2924 of 0 => 0
27 : cchiw 2859 | _ =>((print (String.concat n));1)
28 :     (*end case*))
29 :    
30 : jhr 2924 val empty = IMap.empty
31 :     fun lookup k d = IMap.find(d, k)
32 :     fun insert (k, v) d = IMap.insert(d, k, v)
33 :     fun find (v, mapp) = (case IMap.find(mapp, v)
34 :     of NONE => raise Fail(concat["Outside Bound(", Int.toString v, ")"])
35 :     | SOME s => s
36 :     (* end *))
37 : cchiw 2859
38 :     (*mapIndex:E.mu * dict-> int
39 :     * lookup
40 :     *)
41 : jhr 2924 fun mapIndex (e1, mapp)=(case e1
42 : cchiw 2859 of E.V e => find(e,mapp)
43 : jhr 2924 | E.C c => c
44 : cchiw 2859 (*end case*))
45 :    
46 :    
47 : cchiw 2870 (* intToReal:int->Var*LowIL.ASSGN list
48 : cchiw 2859 *)
49 : cchiw 2870 fun intToReal n=let
50 :     val a=DstIL.Var.new("real" ,realTy)
51 :     val b=DstIL.Var.new("cast" ,realTy)
52 :     val code=[DstIL.ASSGN (a,DstIL.LIT(Literal.Int(IntInf.fromInt n))),
53 :     DstIL.ASSGN (b,DstIL.OP(DstOp.IntToReal,[a]))]
54 : cchiw 2859 in
55 : cchiw 2870 (b,code)
56 : cchiw 2859 end
57 :    
58 : cchiw 2870
59 :    
60 : cchiw 2859 (* mkINt:int->Var*LowIL.ASSGN list
61 :     *)
62 :     fun mkInt n=let
63 :     val a=DstIL.Var.new("Int" ,intTy)
64 :     val code=DstIL.ASSGN (a,DstIL.LIT(Literal.Int(IntInf.fromInt n)))
65 :     val _ =testp[LowToS.toStringAll(intTy,code)]
66 :     in
67 :     (a,[code])
68 :     end
69 :    
70 :     (*assgnCons int list * Var list->Var*LowIL.ASSGN list
71 :     * cons elements on list
72 :     *)
73 :     fun assgnCons(pre,shape, args)=let
74 :     val ty=DstTy.TensorTy shape
75 : cchiw 2870 val a=DstIL.Var.new("cons"^"_" ,ty)
76 : cchiw 2859 val code=DstIL.ASSGN (a,DstIL.CONS(ty ,args))
77 :     val _ =testp[LowToS.toStringAll(ty,code)]
78 :     in
79 :     (a, [code])
80 :     end
81 :    
82 :     (*LowOps.Op.op * var list*string*LowIL.Ty
83 :     * -> Var*LowIL.ASSGN list
84 :     * Make lowIL assignment
85 :     *)
86 :     fun assgn(opss,args,pre,ty)=let
87 :     val a=DstIL.Var.new(pre ,ty)
88 :     val code=DstIL.ASSGN (a,DstIL.OP(opss,args))
89 :     val _ =testp[LowToS.toStringAll(ty,code)]
90 :     in
91 :     (a,[code])
92 :     end
93 :    
94 :     (*getTensorTy:E.params*E.tensor_id-> LowIL.Ty
95 :     * Integer, or Generic Tensor
96 :     *)
97 :     fun getTensorTy(params, id)=(case List.nth(params,id)
98 :     of E.TEN(3,[shape])=> DstTy.iVecTy(shape) (*FIX HERE*)
99 :     | E.TEN(_,shape)=> DstTy.TensorTy shape
100 :     |_=> err"NONE Tensor Param"
101 :     (*end case*))
102 :    
103 :     (* indexTensor:dict*string*E.params*Var list*E.tensor_id*E.alpha
104 :     * ->Var*LowIL.ASSGN list
105 :     * Index Tensor at specific indices to give a scalar result
106 :     *)
107 : cchiw 2923 fun indexTensor(_,(_,_,args,id, [],ty)) =
108 :     (List.nth(args,id),[])
109 : cchiw 2859 | indexTensor(mapp,(lhs,params,args,id,ix,ty))= let
110 :     val nU=List.nth(args,id)
111 :     val ixx=(List.map (fn (e1)=> mapIndex(e1,mapp)) ix)
112 :     val ix'=DstTy.indexTy ixx
113 :     val argTy=getTensorTy(params,id)
114 :     val opp=DstOp.IndexTensor(id,ix',argTy)
115 : cchiw 2870 val name=String.concat["Indx_",iToss ixx,"_"]
116 : cchiw 2859 in
117 :     assgn(opp,[nU],name,ty)
118 :     end
119 :    
120 :     (*projTensor:dict*(string* E.params*Var list*int*E.tensor_id*E.alpha)->Var*LowIL.ASSGN list
121 :     * projects tensor to a vector
122 :     *just used by EintoVecOps but made sense to keep it here
123 :     *)
124 :     fun projTensor(_,(_,_,args,_,id,[]))= (List.nth(args,id),[])
125 :     | projTensor(mapp,(lhs,params,args,vecIX,id,ix))= let
126 :     val nU=List.nth(args,id)
127 :     val ixx=(List.map (fn (e1)=> mapIndex(e1,mapp)) ix)
128 :     val ix'=DstTy.indexTy ixx
129 :     val argTy= getTensorTy(params,id)
130 :     val vecTy=DstTy.TensorTy [vecIX]
131 : cchiw 3195 val opp=DstOp.ProjectLast(id,vecIX,ix',argTy)
132 :     val name=String.concat["ProjLast_",iToss ixx,"_"]
133 : cchiw 2859 in
134 :     assgn(opp,[nU],name,vecTy)
135 :     end
136 :    
137 : cchiw 3195
138 :     fun projFirst(mapp,(lhs,params,args,vecIX,id,ix))= let
139 :     val nU=List.nth(args,id)
140 :     val argTy= getTensorTy(params,id)
141 :     val vecTy=DstTy.TensorTy [vecIX]
142 :     val ixx=(List.map (fn (e1)=> mapIndex(e1,mapp)) ix)
143 :     fun f cnt = let
144 :     val ix'=DstTy.indexTy ([cnt]@ixx)
145 :     val opp=DstOp.IndexTensor(id,ix',argTy)
146 :     val name=String.concat["IndexTensor_",iToss ixx,"_"]
147 :     in
148 :     assgn(opp,[nU],name,realTy)
149 :     end
150 :     val ops=List.tabulate( vecIX, fn e=> f e)
151 :     fun iter ([],vCs,Cs)=(vCs,Cs)
152 :     | iter((vB,B)::es,vCs,Cs)= iter(es,vCs@[vB],Cs@B)
153 :    
154 :     val (vCs,Cs)=iter(ops,[],[])
155 :     val (vD,D)=assgnCons("projFirstCons", [vecIX], vCs)
156 :     in
157 :     (vD, Cs@D)
158 :     end
159 :    
160 :    
161 :    
162 :     fun projFirst2(mapp,(lhs,params,args,vecIX,id,ix))= let
163 :     val nU=List.nth(args,id)
164 :     val argTy= getTensorTy(params,id)
165 :     val vecTy=DstTy.TensorTy [vecIX]
166 :     val ixx=(List.map (fn (e1)=> mapIndex(e1,mapp)) ix)
167 :     val ix'=DstTy.indexTy ixx
168 :     val opp=DstOp.ProjectFirst(id,vecIX,ix',argTy)
169 :     val name=String.concat["ProjFirst_",iToss ixx,"_"]
170 :     in assgn(opp,[nU],name,vecTy)
171 :     end
172 :    
173 :    
174 :    
175 : cchiw 2867 fun mkSqrt(nU,code)= let
176 :     val opp=DstOp.Sqrt
177 : cchiw 2870 val name=String.concat["_Sqrt_"]
178 : cchiw 2867 val (vA,A)=assgn(opp,[nU],name,realTy)
179 :     in
180 :     (vA,code@A)
181 :     end
182 :    
183 : cchiw 3138 fun mkCosine(nU,code)= let
184 :     val opp=DstOp.Cosine
185 :     val name=String.concat["_Cosine_"]
186 :     val (vA,A)=assgn(opp,[nU],name,realTy)
187 :     in
188 :     (vA,code@A)
189 :     end
190 : cchiw 2867
191 : cchiw 3138 fun mkArcCosine(nU,code)= let
192 :     val opp=DstOp.ArcCosine
193 :     val name=String.concat["_ArcCosine_"]
194 :     val (vA,A)=assgn(opp,[nU],name,realTy)
195 :     in
196 :     (vA,code@A)
197 :     end
198 :    
199 :     fun mkSine(nU,code)= let
200 :     val opp=DstOp.Sine
201 :     val name=String.concat["_Sine_"]
202 :     val (vA,A)=assgn(opp,[nU],name,realTy)
203 :     in
204 :     (vA,code@A)
205 :     end
206 :    
207 :     fun mkArcSine(nU,code)= let
208 :     val opp=DstOp.ArcSine
209 :     val name=String.concat["_ArcSine_"]
210 :     val (vA,A)=assgn(opp,[nU],name,realTy)
211 :     in
212 :     (vA,code@A)
213 :     end
214 :    
215 :    
216 : cchiw 2870 fun mkPowInt((nU,code),n)= let
217 :     val opp=DstOp.powInt
218 :     val name=String.concat["_PowInt_"]
219 :     val (r,rcode)=mkInt n
220 :     val (vA,A)=assgn(opp,[nU,r],name,realTy)
221 :     in
222 :     (vA,code@rcode@A)
223 :     end
224 : cchiw 2867
225 : cchiw 2870 fun mkPowRat((nU,code),rat)= let
226 :     val opp=DstOp.powRat(DstTy.R rat)
227 :     val name=String.concat["_PowRat_"]
228 :    
229 :     val (vA,A)=assgn(opp,[nU],name,realTy)
230 :     in
231 :     (vA,code@A)
232 :     end
233 :    
234 :    
235 :    
236 :    
237 :    
238 :    
239 : cchiw 2859 (* Some shortcuts. Arguements are Low-IL variables already indexed/projected
240 :     * string*Var list ->Var*LowIL.ASSGN list
241 :     *)
242 : cchiw 2870 fun mkAddSca(lhs,args)= assgn(DstOp.addSca,args,"addSca",realTy)
243 :     fun mkAddInt(lhs,args)= assgn(DstOp.addSca,args,"addInt",intTy)
244 :     fun mkAddPtr(lhs,args,ty)= assgn(DstOp.addSca,args,"addPtr",ty)
245 :     fun mkAddVec(lhs,vecIX,args)=assgn(DstOp.addVec vecIX,args,"addV",DstTy.TensorTy([vecIX]))
246 :     fun mkSubSca(lhs,args)= assgn(DstOp.subSca,args,"subSca",realTy)
247 :     fun mkProdSca(lhs,args)=assgn(DstOp.prodSca,args,"prodSca",realTy)
248 :     fun mkProdInt(lhs,args)=assgn(DstOp.prodSca,args,"prodInt",intTy)
249 :     fun mkProdVec(lhs,vecIX,args)=assgn(DstOp.prodVec vecIX,args,"prodV",DstTy.TensorTy([vecIX]))
250 :     fun mkDivSca(lhs,args)= assgn(DstOp.divSca,args,"divSca",realTy)
251 :     fun mkSumVec(lhs,vecIX,args)= assgn(DstOp.sumVec vecIX,args,"sumVec",realTy)
252 : cchiw 2859 fun mkDotVec(lhs,vecIX,args)=let
253 : cchiw 2870 val (vD, D)=mkProdVec("",vecIX,args)
254 :     val (vE, E)=mkSumVec("",vecIX,[vD])
255 : cchiw 2859 in (vE,D@E) end
256 :    
257 :    
258 :     (*mkMultiple:string*Var list*LowOps.Op *ListIL.Ty -> Var*LowIL.ASSGN list
259 :     *apply rator between each items on list1
260 :     *)
261 :     fun mkMultiple(lhs,list1,rator,ty)=let
262 :     fun add([],_,_) = err"no element in mkMultiple"
263 :     | add([e1],_,_) = (e1,[])
264 :     | add([e1,e2],code,_) = let
265 : cchiw 2870 val (vA,A)=assgn(rator,[e1,e2],"mult_2",ty)
266 : cchiw 2859 in (vA,code@A) end
267 :     | add(e1::e2::es,code,count) = let
268 : cchiw 2870 val (vA,A)=assgn(rator,[e1,e2],String.concat["mult_",iTos count],ty)
269 : cchiw 2859 in add(vA::es,code@A,count-1)
270 :     end
271 :     in
272 :     add(list1,[],List.length list1)
273 :     end
274 :    
275 :     (* deltaToInt:dict*E.mu*E.mu->int
276 :     * delta function
277 :     *)
278 :     fun deltaToInt(mapp,a,b)= let
279 :     val i=mapIndex(a,mapp)
280 :     val j=mapIndex(b,mapp)
281 :     in
282 :     if(i=j) then 1 else 0
283 :     end
284 :    
285 : cchiw 2870 fun evalDelta(mapp,a,b)= intToReal(deltaToInt(mapp,a,b))
286 : cchiw 2859
287 :     (*eval Epsilon-2d*)
288 :     fun evalEps2(mapp,a,b)=let
289 :     val i=mapIndex(E.V a,mapp)
290 :     val j=mapIndex(E.V b,mapp)
291 :     in
292 : cchiw 2870 if(i=j) then intToReal 0
293 : cchiw 2859 else
294 : cchiw 2870 if(j>i) then intToReal 1
295 :     else intToReal ~1
296 : cchiw 2859 end
297 :    
298 :     (*eval Epsilon-3d*)
299 :     fun evalEps3(mapp,a,b,c)=let
300 :     val i=mapIndex(E.V a,mapp)
301 :     val j=mapIndex(E.V b,mapp)
302 :     val k=mapIndex(E.V c,mapp)
303 :     in
304 : cchiw 2870 if(i=j orelse j=k orelse i=k) then intToReal 0
305 : cchiw 2859 else
306 :     if(j>i) then
307 : cchiw 2870 if(j>k andalso k>i) then intToReal ~1 else intToReal 1
308 :     else if(i>k andalso k>j) then intToReal 1 else intToReal ~1
309 : cchiw 2859
310 :     end
311 :    
312 :     end
313 :    
314 : jhr 2924 end

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0