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 3138 - (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 :     val testing=0
19 :     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 :     val opp=DstOp.ProjectTensor(id,vecIX,ix',argTy)
132 : cchiw 2870 val name=String.concat["Proj_",iToss ixx,"_"]
133 : cchiw 2859 in
134 :     assgn(opp,[nU],name,vecTy)
135 :     end
136 :    
137 : cchiw 2867 fun mkSqrt(nU,code)= let
138 :     val opp=DstOp.Sqrt
139 : cchiw 2870 val name=String.concat["_Sqrt_"]
140 : cchiw 2867 val (vA,A)=assgn(opp,[nU],name,realTy)
141 :     in
142 :     (vA,code@A)
143 :     end
144 :    
145 : cchiw 3138 fun mkCosine(nU,code)= let
146 :     val opp=DstOp.Cosine
147 :     val name=String.concat["_Cosine_"]
148 :     val (vA,A)=assgn(opp,[nU],name,realTy)
149 :     in
150 :     (vA,code@A)
151 :     end
152 : cchiw 2867
153 : cchiw 3138 fun mkArcCosine(nU,code)= let
154 :     val opp=DstOp.ArcCosine
155 :     val name=String.concat["_ArcCosine_"]
156 :     val (vA,A)=assgn(opp,[nU],name,realTy)
157 :     in
158 :     (vA,code@A)
159 :     end
160 :    
161 :     fun mkSine(nU,code)= let
162 :     val opp=DstOp.Sine
163 :     val name=String.concat["_Sine_"]
164 :     val (vA,A)=assgn(opp,[nU],name,realTy)
165 :     in
166 :     (vA,code@A)
167 :     end
168 :    
169 :     fun mkArcSine(nU,code)= let
170 :     val opp=DstOp.ArcSine
171 :     val name=String.concat["_ArcSine_"]
172 :     val (vA,A)=assgn(opp,[nU],name,realTy)
173 :     in
174 :     (vA,code@A)
175 :     end
176 :    
177 :    
178 : cchiw 2870 fun mkPowInt((nU,code),n)= let
179 :     val opp=DstOp.powInt
180 :     val name=String.concat["_PowInt_"]
181 :     val (r,rcode)=mkInt n
182 :     val (vA,A)=assgn(opp,[nU,r],name,realTy)
183 :     in
184 :     (vA,code@rcode@A)
185 :     end
186 : cchiw 2867
187 : cchiw 2870 fun mkPowRat((nU,code),rat)= let
188 :     val opp=DstOp.powRat(DstTy.R rat)
189 :     val name=String.concat["_PowRat_"]
190 :    
191 :     val (vA,A)=assgn(opp,[nU],name,realTy)
192 :     in
193 :     (vA,code@A)
194 :     end
195 :    
196 :    
197 :    
198 :    
199 :    
200 :    
201 : cchiw 2859 (* Some shortcuts. Arguements are Low-IL variables already indexed/projected
202 :     * string*Var list ->Var*LowIL.ASSGN list
203 :     *)
204 : cchiw 2870 fun mkAddSca(lhs,args)= assgn(DstOp.addSca,args,"addSca",realTy)
205 :     fun mkAddInt(lhs,args)= assgn(DstOp.addSca,args,"addInt",intTy)
206 :     fun mkAddPtr(lhs,args,ty)= assgn(DstOp.addSca,args,"addPtr",ty)
207 :     fun mkAddVec(lhs,vecIX,args)=assgn(DstOp.addVec vecIX,args,"addV",DstTy.TensorTy([vecIX]))
208 :     fun mkSubSca(lhs,args)= assgn(DstOp.subSca,args,"subSca",realTy)
209 :     fun mkProdSca(lhs,args)=assgn(DstOp.prodSca,args,"prodSca",realTy)
210 :     fun mkProdInt(lhs,args)=assgn(DstOp.prodSca,args,"prodInt",intTy)
211 :     fun mkProdVec(lhs,vecIX,args)=assgn(DstOp.prodVec vecIX,args,"prodV",DstTy.TensorTy([vecIX]))
212 :     fun mkDivSca(lhs,args)= assgn(DstOp.divSca,args,"divSca",realTy)
213 :     fun mkSumVec(lhs,vecIX,args)= assgn(DstOp.sumVec vecIX,args,"sumVec",realTy)
214 : cchiw 2859 fun mkDotVec(lhs,vecIX,args)=let
215 : cchiw 2870 val (vD, D)=mkProdVec("",vecIX,args)
216 :     val (vE, E)=mkSumVec("",vecIX,[vD])
217 : cchiw 2859 in (vE,D@E) end
218 :    
219 :    
220 :     (*mkMultiple:string*Var list*LowOps.Op *ListIL.Ty -> Var*LowIL.ASSGN list
221 :     *apply rator between each items on list1
222 :     *)
223 :     fun mkMultiple(lhs,list1,rator,ty)=let
224 :     fun add([],_,_) = err"no element in mkMultiple"
225 :     | add([e1],_,_) = (e1,[])
226 :     | add([e1,e2],code,_) = let
227 : cchiw 2870 val (vA,A)=assgn(rator,[e1,e2],"mult_2",ty)
228 : cchiw 2859 in (vA,code@A) end
229 :     | add(e1::e2::es,code,count) = let
230 : cchiw 2870 val (vA,A)=assgn(rator,[e1,e2],String.concat["mult_",iTos count],ty)
231 : cchiw 2859 in add(vA::es,code@A,count-1)
232 :     end
233 :     in
234 :     add(list1,[],List.length list1)
235 :     end
236 :    
237 :     (* deltaToInt:dict*E.mu*E.mu->int
238 :     * delta function
239 :     *)
240 :     fun deltaToInt(mapp,a,b)= let
241 :     val i=mapIndex(a,mapp)
242 :     val j=mapIndex(b,mapp)
243 :     in
244 :     if(i=j) then 1 else 0
245 :     end
246 :    
247 : cchiw 2870 fun evalDelta(mapp,a,b)= intToReal(deltaToInt(mapp,a,b))
248 : cchiw 2859
249 :     (*eval Epsilon-2d*)
250 :     fun evalEps2(mapp,a,b)=let
251 :     val i=mapIndex(E.V a,mapp)
252 :     val j=mapIndex(E.V b,mapp)
253 :     in
254 : cchiw 2870 if(i=j) then intToReal 0
255 : cchiw 2859 else
256 : cchiw 2870 if(j>i) then intToReal 1
257 :     else intToReal ~1
258 : cchiw 2859 end
259 :    
260 :     (*eval Epsilon-3d*)
261 :     fun evalEps3(mapp,a,b,c)=let
262 :     val i=mapIndex(E.V a,mapp)
263 :     val j=mapIndex(E.V b,mapp)
264 :     val k=mapIndex(E.V c,mapp)
265 :     in
266 : cchiw 2870 if(i=j orelse j=k orelse i=k) then intToReal 0
267 : cchiw 2859 else
268 :     if(j>i) then
269 : cchiw 2870 if(j>k andalso k>i) then intToReal ~1 else intToReal 1
270 :     else if(i>k andalso k>j) then intToReal 1 else intToReal ~1
271 : cchiw 2859
272 :     end
273 :    
274 :     end
275 :    
276 : jhr 2924 end

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