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

SCM Repository

[diderot] Diff of /branches/charisee_dev/src/compiler/high-to-mid/handleEin.sml
ViewVC logotype

Diff of /branches/charisee_dev/src/compiler/high-to-mid/handleEin.sml

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

branches/charisee/src/compiler/high-to-mid/handleEin.sml revision 2870, Wed Feb 25 21:47:43 2015 UTC branches/charisee_dev/src/compiler/high-to-mid/handleEin.sml revision 3503, Thu Dec 17 23:13:57 2015 UTC
# Line 1  Line 1 
1  (* Expands probe ein  (* Expands probe ein
2   *   *
3   * COPYRIGHT (c) 2012 The Diderot Project (http://diderot-language.cs.uchicago.edu)   * This code is part of the Diderot Project (http://diderot-language.cs.uchicago.edu)
4     *
5     * COPYRIGHT (c) 2015 The University of Chicago
6   * All rights reserved.   * All rights reserved.
7   *)   *)
8    
# Line 10  Line 12 
12    
13      structure E = Ein      structure E = Ein
14      structure DstIL = MidIL      structure DstIL = MidIL
     structure DstOp = MidOps  
15      structure P=Printer      structure P=Printer
     structure T=TransformEin  
16      structure MidToS=MidToString      structure MidToS=MidToString
17      structure DstV = DstIL.Var      structure L=LiftEin
18      in      in
19    
20        val liftflag=true
21        val splitflag=true
22        val splitlimit=NONE
23        val testinitial=0
24      val testing=0      val testing=0
     fun setEin(params,index,body)=Ein.EIN{params=params, index=index, body=body}  
     fun assignEinApp(y,params,index,body,args)= (y,DstIL.EINAPP(setEin(params,index,body),args))  
     fun assignEinApp(y,params,index,body,args)= (y,DstIL.EINAPP(setEin(params,index,body),args))  
       fun iTos i =Int.toString i  
25      fun testp n=(case testing      fun testp n=(case testing
26          of 0=> 1          of 0=> 1
27          | _ =>(print(String.concat n);1)          | _ =>(print(String.concat n);1)
28          (*end case*))          (*end case*))
29      val einappzero=DstIL.EINAPP(setEin([],[],E.Const 0),[])      fun test0 n=(case testinitial
30      fun setEinZero y=  (y,einappzero)          of 0=>testp n
31            | _ => (print(String.concat n);1))
32      (* isZero:ein_exp->int      fun toStringBind e= MidToS.toStringBind e
33      *  Checks to see if the entire body is 0. If so returns 1.      fun useCount (HighIL.V{useCnt, ...}) = !useCnt
34        fun expandEinOp (srcy,bind00)=let
35            (*val _= raise Fail"development branch- do not use"*)
36    
37            val star="************"
38            val _= ("\n Start"^star)
39            val _ =test0[star,"\n Original EinApp",star,"\n\n","start get test",toStringBind bind00]
40            val (y,einapp as DstIL.EINAPP(Ein.EIN{params, index, body},args))=bind00
41    
42    
43            (* ************* Sweep for 0's'*********** *)
44            val _ =testp["Pre Sweep\n",P.printbody body,"\n"]
45            val bodysweep=sweepEin.sweep body
46            val ein1=Ein.EIN{params=params, index=index, body=bodysweep}
47            val _ =testp["Post Sweep\n",P.printbody bodysweep,"\n"]
48    
49    
50            (* ************** distribute and push Summation*********** *)
51            val bind3 = (y,DstIL.EINAPP(SummationEin.main ein1,args))
52            val _ =testp["\n\n summation \n\n =>", toStringBind bind3]
53    
54            (* **************** lift - split ************* *)
55            val newbies5=  (case (splitflag,liftflag,splitlimit)
56                of (true,false,NONE) => Split.splitEinApp bind3
57                (* **************** split phase ************* *)
58                | (true,false,SOME limit) =>let
59                    val _ = (String.concat["\nSPLit with limit",Int.toString(limit)])
60                    in Split.limitSplit(bind3,limit) end
61                | (true,true,SOME limit) =>  let
62                    (* **************** lift and split Limit ************* *)
63                    val (einapp1,e2a,fieldset)=L.testLift  bind3
64                    val _ = (String.concat["\nSPLit with limit",Int.toString(limit),"\n **** lift Length:",Int.toString(length(e2a))])
65    
66                   (* val t2=List.map (fn t=>Split.splitEinApp t) e2*)
67                    val t2=List.map (fn t =>Split.limitSplit2(t,limit,fieldset)) e2a
68                    val e2= List.foldr op@ [] t2
69                    val newbies4=Split.limitSplit2(einapp1,limit,fieldset)
70                    val _ = print(String.concat["\n split made pieces(",Int.toString(length(e2)-length(e2a)),") orig (", Int.toString(length(newbies4)),")"])
71                    in (e2@newbies4) end
72                | (true,true, NONE) => let
73                    (* **************** lift and split All ************* *)
74                    val (einapp1,e2a,_)=L.testLift  bind3
75                    val _ =  print(String.concat["\n **** lift Length"^(Int.toString(length(e2a)))])
76                    val e2= List.foldr op@ [] (List.map (fn t=>Split.splitEinApp t) e2a)
77                    val newbies4=Split.splitEinApp einapp1
78                    val _ = print(String.concat["\n split made pieces(",Int.toString(length(e2)-length(e2a)),") orig (", Int.toString(length(newbies4)),")"])
79                    in (e2@newbies4) end
80    
81    (*
82                | (false,true,_) => let
83                    val (einapp1,e2)=L.testLift  bind3
84                    val _ = ("\n **** lift Length"^(Int.toString(length(e2))))
85                    in (e2@[einapp1]) end
86                |  (false,false,_)=>Split.splitEinApp bind3
87      *)      *)
     fun isZero e=let  
         val _=testp["\n Started Zero",P.printbody e]  
         fun iterList []= 1  
             | iterList(e1::list)=(case (isZero e1)  
             of 0=>0  
             | _ => iterList list  
             (*end case*))  
         fun iterProdList []= 0  
         | iterProdList [e1]=isZero e1  
         | iterProdList(e1::list)=(case (isZero e1)  
             of 1=>1  
             | _ => iterList list  
88              (*end case*))              (*end case*))
89          val z=(case e          val _ = "************** TASH"
90              of E.Tensor _         => 0          (**************** printing result ***************)
91              | E.Add e             => iterList e          val _ =testp["\n\n Returning \n\n =>",
92              | E.Sub(e1,e2)         => iterList [e1,e2]          String.concatWith",\n\t"(List.map toStringBind newbies5)]
93              | E.Div(e1,e2)         => iterProdList [e1,e2]          val n =length(newbies5)
94              | E.Sum(_ ,e)          => isZero e          val _ = if (n > 10) then print(String.concat["\n **** SPLIT Length",Int.toString(n)]) else print ""
95              | E.Prod e             => iterProdList e           (* ************** ProbeEIN *********** *)
96              | E.Neg e              => isZero e          (*gives probeEIN the set. Could split into a sepearate step*)
97              | E.Probe(E.Conv _,_)  => 0          val fieldset= einVarSet.EinVarSet.empty
98              | E.Sqrt e             => isZero e          val varset= einVarSet.VarSet.empty
99              | E.Const 0            => 1          val isField = ref 0
100              | E.Field _            => 1          val isFieldR = ref 0
101              | E.Partial _          => 1          val isReplaced = ref 0
102              | E.Apply _            => 1          fun iter([],_,_)=[]
103              | E.Lift _             => 1           | iter(e1::es,fieldset,varset)= let
104              | E.Conv _             => 1              val (e2,fieldset,varset,isFieldt,isReplacedt) = ProbeEin.expandEinOp(e1,fieldset,varset)
105              | E.PowInt(e ,_)            => isZero e              val rtn=[e2]@(iter(es,fieldset,varset))
106              | E.PowReal(e,_)             => isZero e              val isFieldRt=if ((isFieldt+isReplacedt)=2) then 1 else 0
107              | _                    => 0              in (isField:=(!isField)+isFieldt;isReplaced:=(!isReplaced )+isReplacedt;
108              (*end case*))              isFieldR:=(!isFieldR)+isFieldRt;rtn)
         val _=testp["\n Is it Zero",P.printbody e,"--", iTos z]  
         in  
             z  
109          end          end
110            val code=iter(newbies5, fieldset,varset)
111            val flatcode= List.foldr op@ [] code
112    
113            val n= List.length(code)
114            val m= List.length(flatcode)
115            val _ =testp[HighIL.Var.name srcy ," Pieces",Int.toString (n),"\n"]
116    
117              val _ = if (n>5)
118              then  print(String.concat["\n",star,":Name:",HighIL.Var.name srcy ,
119                ":VarCount(",Int.toString(useCount srcy),")","\n\tPieces now: ",Int.toString (n),"flatcode",Int.toString (m),
120                " Fields:",Int.toString(!isField)," isReplaced:",Int.toString(!isReplaced),
121                " Replaced_Fields:",Int.toString(!isFieldR)])
122                else print ""
123    
124                val flatcode= List.foldr op@ [] code
125    
126    
127            (*val _=List.map (fn(_,DstIL.EINAPP(e,_))=>checkEin.checkEIN e | _=> 1) flatcode*)
128            val _=(star^"end")
129            in
130                List.map (fn (y,rator)=> DstIL.ASSGN(y,rator)) flatcode
131            end
132    
     fun sweep e= (case e  
         of E.Tensor _           => e  
         | E.Add es             => E.Add(List.map sweep es)  
         | E.Sub(e1,e2)         => E.Sub(sweep e1,sweep e2)  
         | E.Div(e1,e2)         => E.Div(sweep e1,sweep e2)  
         | E.Sum(c ,e)          => E.Sum(c, sweep e)  
         | E.Prod es             => E.Prod(List.map sweep es)  
         | E.Neg e              => E.Neg(sweep e)  
         | E.Probe(E.Conv _,_)  => e  
         | E.Sqrt e             => E.Sqrt(sweep e)  
         | E.Const _            => e  
         | E.ConstR _           => e  
         | E.Delta _               => e  
         | E.Epsilon _             => e  
         | E.Eps2 _                => e  
         | E.Field _            => E.Const 0  
         | E.Partial _          => E.Const 0  
         | E.Apply _            => E.Const 0  
         | E.Lift _             => E.Const 0  
         | E.Conv _             => E.Const 0  
         | E.PowInt(e ,n)       => E.PowInt(sweep e ,n)  
         | E.PowReal(e,n)       => E.PowReal(sweep e ,n)  
         | E.Probe _            =>  
 raise Fail (String.concat["\n Incorrect probe, substition was not made. Is the Field in an if statement? :",P.printbody e])  
     (*end case*))  
133    end; (* local *)    end; (* local *)
134    
135  end (* local *)  end (* local *)

Legend:
Removed from v.2870  
changed lines
  Added in v.3503

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