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

SCM Repository

[diderot] Diff of /branches/vis15/src/compiler/simplify/simplify.sml
ViewVC logotype

Diff of /branches/vis15/src/compiler/simplify/simplify.sml

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

revision 3465, Sun Nov 29 20:04:16 2015 UTC revision 4113, Thu Jun 30 13:52:19 2016 UTC
# Line 5  Line 5 
5   * COPYRIGHT (c) 2015 The University of Chicago   * COPYRIGHT (c) 2015 The University of Chicago
6   * All rights reserved.   * All rights reserved.
7   *   *
8   * Simplify the AST representation.   * Simplify the AST representation.  This phase involves the following transformations:
9     *
10     *      - types are simplified by removing meta variables (which will have been resolved)
11     *
12     *      - expressions are simplified to involve a single operation on variables
13     *
14     *      - global reductions are converted to MapReduce statements
15     *
16     *      - other comprehensions and reductions are converted to foreach loops
17     *
18     *      - unreachable code is pruned
19     *
20     *      - negation of literal integers and reals are constant folded
21   *)   *)
22    
23  structure Simplify : sig  structure Simplify : sig
# Line 31  Line 43 
43              | Ty.T_String => STy.T_String              | Ty.T_String => STy.T_String
44              | Ty.T_Sequence(ty, NONE) => STy.T_Sequence(cvtTy ty, NONE)              | Ty.T_Sequence(ty, NONE) => STy.T_Sequence(cvtTy ty, NONE)
45              | Ty.T_Sequence(ty, SOME dim) => STy.T_Sequence(cvtTy ty, SOME(TU.monoDim dim))              | Ty.T_Sequence(ty, SOME dim) => STy.T_Sequence(cvtTy ty, SOME(TU.monoDim dim))
46              | Ty.T_Named id => STy.T_Named id              | Ty.T_Strand id => STy.T_Strand id
47              | Ty.T_Kernel n => STy.T_Kernel(TU.monoDiff n)              | Ty.T_Kernel n => STy.T_Kernel(TU.monoDiff n)
48              | Ty.T_Tensor shape => STy.T_Tensor(TU.monoShape shape)              | Ty.T_Tensor shape => STy.T_Tensor(TU.monoShape shape)
49              | Ty.T_Image{dim, shape} => STy.T_Image{              | Ty.T_Image{dim, shape} => STy.T_Image{
# Line 47  Line 59 
59              | Ty.T_Error => raise Fail "unexpected T_Error in Simplify"              | Ty.T_Error => raise Fail "unexpected T_Error in Simplify"
60            (* end case *))            (* end case *))
61    
62        fun apiTypeOf x = let
63              fun cvtTy STy.T_Bool = APITypes.BoolTy
64                | cvtTy STy.T_Int = APITypes.IntTy
65                | cvtTy STy.T_String = APITypes.StringTy
66                | cvtTy (STy.T_Sequence(ty, len)) = APITypes.SeqTy(cvtTy ty, len)
67                | cvtTy (STy.T_Tensor shape) = APITypes.TensorTy shape
68                | cvtTy (STy.T_Image{dim, shape}) = APITypes.ImageTy(dim, shape)
69                | cvtTy ty = raise Fail "bogus API type"
70              in
71                cvtTy (SimpleVar.typeOf x)
72              end
73    
74      fun newTemp ty = SimpleVar.new ("_t", SimpleVar.LocalVar, ty)      fun newTemp ty = SimpleVar.new ("_t", SimpleVar.LocalVar, ty)
75    
76    (* a property to map AST variables to SimpleAST variables *)    (* a property to map AST variables to SimpleAST variables *)
# Line 59  Line 83 
83      fun cvtVars xs = List.map cvtVar xs      fun cvtVars xs = List.map cvtVar xs
84    
85    (* make a block out of a list of statements that are in reverse order *)    (* make a block out of a list of statements that are in reverse order *)
86      fun mkBlock stms = S.Block(List.rev stms)      fun mkBlock stms = S.Block{props = PropList.newHolder(), code = List.rev stms}
   
     fun inputImage (errStrm, nrrd, dim, shape) = (  
           case ImageInfo.fromNrrd(NrrdInfo.getInfo(errStrm, nrrd), dim, shape)  
            of NONE => raise Fail(concat["nrrd file \"", nrrd, "\" does not have expected type"])  
             | SOME info => S.Proxy(nrrd, info)  
           (* end case *))  
   
     datatype 'a ctl_flow_info  
       = EXIT                    (* stm sequence always exits; no pruning so far *)  
       | PRUNE of 'a             (* stm sequence always exits at last stm in argument, which  
                                  * is either a block or stm list *)  
       | CONT                    (* stm sequence falls through *)  
       | EDIT of 'a              (* pruned code that has non-exiting paths *)  
   
     fun pruneUnreachableCode (blk as S.Block stms) = let  
           fun isExit S.S_Die = true  
             | isExit S.S_Stabilize = true  
             | isExit (S.S_Return _) = true  
             | isExit _ = false  
           fun pruneStms [] = CONT  
             | pruneStms [S.S_IfThenElse(x, blk1, blk2)] = (  
                 case pruneIf(x, blk1, blk2)  
                  of EXIT => EXIT  
                   | PRUNE stm => PRUNE[stm]  
                   | CONT => CONT  
                   | EDIT stm => EDIT[stm]  
                 (* end case *))  
             | pruneStms [stm] = if isExit stm then EXIT else CONT  
             | pruneStms ((stm as S.S_IfThenElse(x, blk1, blk2))::stms) = (  
                 case pruneIf(x, blk1, blk2)  
                  of EXIT => PRUNE[stm]  
                   | PRUNE stm => PRUNE[stm]  
                   | CONT => (case pruneStms stms  
                        of PRUNE stms => PRUNE(stm::stms)  
                         | EDIT stms => EDIT(stm::stms)  
                         | EXIT => EXIT (* different instances of ctl_flow_info *)  
                         | CONT => CONT  
                       (* end case *))  
                   | EDIT stm => (case pruneStms stms  
                        of PRUNE stms => PRUNE(stm::stms)  
                         | EDIT stms => EDIT(stm::stms)  
                         | _ => EDIT(stm::stms)  
                       (* end case *))  
                 (* end case *))  
             | pruneStms (stm::stms) = if isExit stm  
                 then PRUNE[stm]  
                 else (case pruneStms stms  
                    of PRUNE stms => PRUNE(stm::stms)  
                     | EDIT stms => EDIT(stm::stms)  
                     | info => info  
                   (* end case *))  
           and pruneIf (x, blk1, blk2) = (case (pruneBlk blk1, pruneBlk blk2)  
                  of (EXIT,       EXIT      ) => EXIT  
                   | (CONT,       CONT      ) => CONT  
                   | (CONT,       EXIT      ) => CONT  
                   | (EXIT,       CONT      ) => CONT  
                   | (CONT,       EDIT blk2 ) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (EDIT blk1,  CONT      ) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (CONT,       PRUNE blk2) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (PRUNE blk1, CONT      ) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (EXIT,       EDIT blk2 ) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (EDIT blk1,  EXIT      ) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (EDIT blk1,  EDIT blk2 ) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (EDIT blk1,  PRUNE blk2) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (PRUNE blk1, EDIT blk2 ) => EDIT(S.S_IfThenElse(x, blk1, blk2))  
                   | (EXIT,       PRUNE blk2) => PRUNE(S.S_IfThenElse(x, blk1, blk2))  
                   | (PRUNE blk1, EXIT      ) => PRUNE(S.S_IfThenElse(x, blk1, blk2))  
                   | (PRUNE blk1, PRUNE blk2) => PRUNE(S.S_IfThenElse(x, blk1, blk2))  
                 (* end case *))  
           and pruneBlk (S.Block stms) = (case pruneStms stms  
                  of PRUNE stms => PRUNE(S.Block stms)  
                   | EDIT stms => EDIT(S.Block stms)  
                   | EXIT => EXIT (* different instances of ctl_flow_info *)  
                   | CONT => CONT  
                 (* end case *))  
           in  
             case pruneBlk blk  
              of PRUNE blk => blk  
               | EDIT blk => blk  
               | _=> blk  
             (* end case *)  
           end  
87    
88    (* simplify a statement into a single statement (i.e., a block if it expands    (* simplify a statement into a single statement (i.e., a block if it expands
89     * into more than one new statement).     * into more than one new statement).
90     *)     *)
91      fun simplifyBlock errStrm stm = mkBlock (simplifyStmt (errStrm, stm, []))      fun simplifyBlock (errStrm, stm) = mkBlock (simplifyStmt (errStrm, stm, []))
92    
93    (* simplify the statement stm where stms is a reverse-order list of preceeding simplified    (* simplify the statement stm where stms is a reverse-order list of preceeding simplified
94     * statements.  This function returns a reverse-order list of simplified statements.     * statements.  This function returns a reverse-order list of simplified statements.
# Line 173  Line 115 
115                  end                  end
116              | AST.S_IfThenElse(e, s1, s2) => let              | AST.S_IfThenElse(e, s1, s2) => let
117                  val (stms, x) = simplifyExpToVar (errStrm, e, stms)                  val (stms, x) = simplifyExpToVar (errStrm, e, stms)
118                  val s1 = simplifyBlock errStrm s1                  val s1 = simplifyBlock (errStrm, s1)
119                  val s2 = simplifyBlock errStrm s2                  val s2 = simplifyBlock (errStrm, s2)
120                  in                  in
121                    S.S_IfThenElse(x, s1, s2) :: stms                    S.S_IfThenElse(x, s1, s2) :: stms
122                  end                  end
123              | AST.S_Foreach((x, e), body) => let              | AST.S_Foreach((x, e), body) => let
124                  val (stms, xs') = simplifyExpToVar (errStrm, e, stms)                  val (stms, xs') = simplifyExpToVar (errStrm, e, stms)
125                  val body' = simplifyBlock errStrm body                  val body' = simplifyBlock (errStrm, body)
126                  in                  in
127                    S.S_Foreach(cvtVar x, xs', body') :: stms                    S.S_Foreach(cvtVar x, xs', body') :: stms
128                  end                  end
# Line 239  Line 181 
181                      | _ => (stms, S.E_Var(cvtVar x))                      | _ => (stms, S.E_Var(cvtVar x))
182                    (* end case *))                    (* end case *))
183                | AST.E_Lit lit => (stms, S.E_Lit lit)                | AST.E_Lit lit => (stms, S.E_Lit lit)
184                  | AST.E_Kernel h => (stms, S.E_Kernel h)
185                | AST.E_Select(e, (fld, _)) => let                | AST.E_Select(e, (fld, _)) => let
186                    val (stms, x) = simplifyExpToVar (errStrm, e, stms)                    val (stms, x) = simplifyExpToVar (errStrm, e, stms)
187                    in                    in
# Line 311  Line 254 
254                    in                    in
255                      (foreachStm :: initStm :: stms, S.E_Var acc)                      (foreachStm :: initStm :: stms, S.E_Var acc)
256                    end                    end
257                | AST.E_ParallelMap(e, x, xs, ty) => raise Fail "FIXME"                | AST.E_ParallelMap(e, x, xs, ty) => raise Fail "FIXME: ParallelMap"
258                | AST.E_Tensor(es, ty) => let                | AST.E_Tensor(es, ty) => let
259                    val (stms, xs) = simplifyExpsToVars (errStrm, es, stms)                    val (stms, xs) = simplifyExpsToVars (errStrm, es, stms)
260                    in                    in
# Line 324  Line 267 
267                    end                    end
268                | AST.E_Slice(e, indices, ty) => let (* tensor slicing *)                | AST.E_Slice(e, indices, ty) => let (* tensor slicing *)
269                    val (stms, x) = simplifyExpToVar (errStrm, e, stms)                    val (stms, x) = simplifyExpToVar (errStrm, e, stms)
270                    fun f ([], ys, stms) = (stms, List.rev ys)                    fun f NONE = NONE
271                      | f (NONE::es, ys, stms) = f (es, NONE::ys, stms)                      | f (SOME(AST.E_Lit(Literal.Int i))) = SOME(Int.fromLarge i)
272                      | f (SOME e::es, ys, stms) = let                      | f _ = raise Fail "expected integer literal in slice"
273                          val (stms, y) = simplifyExpToVar (errStrm, e, stms)                    val indices = List.map f indices
                         in  
                           f (es, SOME y::ys, stms)  
                         end  
                   val (stms, indices) = f (indices, [], stms)  
274                    in                    in
275                      (stms, S.E_Slice(x, indices, cvtTy ty))                      (stms, S.E_Slice(x, indices, cvtTy ty))
276                    end                    end
# Line 352  Line 291 
291                | AST.E_LoadNrrd(_, nrrd, ty) => (case cvtTy ty                | AST.E_LoadNrrd(_, nrrd, ty) => (case cvtTy ty
292                     of ty as SimpleTypes.T_Sequence(_, NONE) => (stms, S.E_LoadSeq(ty, nrrd))                     of ty as SimpleTypes.T_Sequence(_, NONE) => (stms, S.E_LoadSeq(ty, nrrd))
293                      | ty as SimpleTypes.T_Image{dim, shape} => (                      | ty as SimpleTypes.T_Image{dim, shape} => (
294                          case ImageInfo.fromNrrd(NrrdInfo.getInfo(errStrm, nrrd), dim, shape)                          case NrrdInfo.getInfo (errStrm, nrrd)
295                           of NONE => raise Fail(concat[                           of SOME info => (case ImageInfo.fromNrrd(info, dim, shape)
296                                   of NONE => (
297                                        Error.error (errStrm, [
298                                  "nrrd file \"", nrrd, "\" does not have expected type"                                  "nrrd file \"", nrrd, "\" does not have expected type"
299                                ])                                        ]);
300                                        (stms, S.E_LoadImage(ty, nrrd, ImageInfo.mkInfo(dim, shape))))
301                            | SOME info => (stms, S.E_LoadImage(ty, nrrd, info))                            | SOME info => (stms, S.E_LoadImage(ty, nrrd, info))
302                          (* end case *))                          (* end case *))
303                              | NONE => (
304                                  Error.warning (errStrm, [
305                                      "nrrd file \"", nrrd, "\" does not exist"
306                                    ]);
307                                  (stms, S.E_LoadImage(ty, nrrd, ImageInfo.mkInfo(dim, shape))))
308                            (* end case *))
309                      | _ => raise Fail "bogus type for E_LoadNrrd"                      | _ => raise Fail "bogus type for E_LoadNrrd"
310                    (* end case *))                    (* end case *))
311                  | AST.E_Coerce{dstTy, e=AST.E_Lit(Literal.Int n), ...} => (case cvtTy dstTy
312                       of SimpleTypes.T_Tensor[] => (stms, S.E_Lit(Literal.Real(RealLit.fromInt n)))
313                        | _ => raise Fail "impossible: bad coercion"
314                      (* end case *))
315                | AST.E_Coerce{srcTy, dstTy, e} => let                | AST.E_Coerce{srcTy, dstTy, e} => let
316                    val (stms, x) = simplifyExpToVar (errStrm, e, stms)                    val (stms, x) = simplifyExpToVar (errStrm, e, stms)
317                    val dstTy = cvtTy dstTy                    val dstTy = cvtTy dstTy
# Line 395  Line 347 
347              f (exps, [], stms)              f (exps, [], stms)
348            end            end
349    
350      fun simplifyStrand (errStrm, AST.Strand{name, params, state, initM, updateM, stabilizeM}) = let    (* simplify a block and then prune unreachable and dead code *)
351        fun simplifyAndPruneBlock errStrm blk =
352              DeadCode.eliminate (simplifyBlock (errStrm, blk))
353    
354        fun simplifyStrand (errStrm, strand) = let
355              val AST.Strand{name, params, state, stateInit, initM, updateM, stabilizeM} = strand
356            val params' = cvtVars params            val params' = cvtVars params
357            fun simplifyState ([], xs, stms) = (List.rev xs, mkBlock stms)            fun simplifyState ([], xs, stms) = (List.rev xs, mkBlock stms)
358              | simplifyState ((x, optE) :: r, xs, stms) = let              | simplifyState ((x, optE) :: r, xs, stms) = let
# Line 406  Line 363 
363                      | SOME e => let                      | SOME e => let
364                          val (stms, e') = simplifyExp (errStrm, e, stms)                          val (stms, e') = simplifyExp (errStrm, e, stms)
365                          in                          in
366                            simplifyState (r, x'::xs, S.S_Var(x', SOME e') :: stms)                            simplifyState (r, x'::xs, S.S_Assign(x', e') :: stms)
367                          end                          end
368                    (* end case *)                    (* end case *)
369                  end                  end
# Line 417  Line 374 
374                  params = params',                  params = params',
375                  state = xs,                  state = xs,
376                  stateInit = stm,                  stateInit = stm,
377                  initM = Option.map (simplifyBlock errStrm) initM,                  initM = Option.map (simplifyAndPruneBlock errStrm) initM,
378                  updateM = simplifyBlock errStrm updateM,                  updateM = simplifyAndPruneBlock errStrm updateM,
379                  stabilizeM = Option.map (simplifyBlock errStrm) stabilizeM                  stabilizeM = Option.map (simplifyAndPruneBlock errStrm) stabilizeM
380                }                }
381            end            end
382    
     fun simplifyCreate (errStrm, AST.C_Grid(dim, stm)) = S.C_Grid(dim, simplifyBlock errStrm stm)  
       | simplifyCreate (errStrm, AST.C_Collection stm) = S.C_Collection(simplifyBlock errStrm stm)  
   
383      fun transform (errStrm, prog) = let      fun transform (errStrm, prog) = let
384            val AST.Program{            val AST.Program{
385                    props, const_dcls, input_dcls, globals, init, strand, create, update                    props, const_dcls, input_dcls, globals, globInit, strand, create, init, update
386                  } = prog                  } = prog
387            val consts' = ref[]            val consts' = ref[]
388            val constInit = ref[]            val constInit = ref[]
# Line 446  Line 400 
400            fun simplifyInputDcl ((x, NONE), desc) = let            fun simplifyInputDcl ((x, NONE), desc) = let
401                  val x' = cvtVar x                  val x' = cvtVar x
402                  val init = (case SimpleVar.typeOf x'                  val init = (case SimpleVar.typeOf x'
403                         of SimpleTypes.T_Image{dim, shape} => let                         of STy.T_Image{dim, shape} => let
404                              val info = ImageInfo.mkInfo(dim, shape)                              val info = ImageInfo.mkInfo(dim, shape)
405                              in                              in
406                                S.Image info                                S.Image info
# Line 455  Line 409 
409                        (* end case *))                        (* end case *))
410                  val inp = S.INP{                  val inp = S.INP{
411                          var = x',                          var = x',
412                            name = Var.nameOf x,
413                            ty =  apiTypeOf x',
414                          desc = desc,                          desc = desc,
415                          init = init                          init = init
416                        }                        }
# Line 463  Line 419 
419                  end                  end
420              | simplifyInputDcl ((x, SOME(AST.E_LoadNrrd(tvs, nrrd, ty))), desc) = let              | simplifyInputDcl ((x, SOME(AST.E_LoadNrrd(tvs, nrrd, ty))), desc) = let
421                  val x' = cvtVar x                  val x' = cvtVar x
               (* load the nrrd proxy here *)  
                 val info = NrrdInfo.getInfo (errStrm, nrrd)  
422                  val init = (case SimpleVar.typeOf x'                  val init = (case SimpleVar.typeOf x'
423                         of SimpleTypes.T_Sequence(_, NONE) => S.LoadSeq nrrd                         of SimpleTypes.T_Sequence(_, NONE) => S.LoadSeq nrrd
424                          | SimpleTypes.T_Image{dim, shape} => inputImage(errStrm, nrrd, dim, shape)                          | SimpleTypes.T_Image{dim, shape} => (
425                                case NrrdInfo.getInfo (errStrm, nrrd)
426                                 of SOME info => (case ImageInfo.fromNrrd(info, dim, shape)
427                                       of NONE => (
428                                            Error.error (errStrm, [
429                                                "proxy nrrd file \"", nrrd,
430                                                "\" does not have expected type"
431                                              ]);
432                                            S.Image(ImageInfo.mkInfo(dim, shape)))
433                                        | SOME info => S.Proxy(nrrd, info)
434                                      (* end case *))
435                                  | NONE => (
436                                      Error.warning (errStrm, [
437                                          "proxy nrrd file \"", nrrd, "\" does not exist"
438                                        ]);
439                                      S.Image(ImageInfo.mkInfo(dim, shape)))
440                                (* end case *))
441                          | _ => raise Fail "impossible"                          | _ => raise Fail "impossible"
442                        (* end case *))                        (* end case *))
443                  val inp = S.INP{                  val inp = S.INP{
444                          var = x',                          var = x',
445                            name = Var.nameOf x,
446                            ty = apiTypeOf x',
447                          desc = desc,                          desc = desc,
448                          init = init                          init = init
449                        }                        }
# Line 483  Line 455 
455                  val (stms, e') = simplifyExp (errStrm, e, [])                  val (stms, e') = simplifyExp (errStrm, e, [])
456                  val inp = S.INP{                  val inp = S.INP{
457                          var = x',                          var = x',
458                            name = Var.nameOf x,
459                            ty = apiTypeOf x',
460                          desc = desc,                          desc = desc,
461                          init = S.ConstExpr                          init = S.ConstExpr
462                        }                        }
# Line 506  Line 480 
480              | simplifyGlobalDcl (AST.D_Func(f, params, body)) = let              | simplifyGlobalDcl (AST.D_Func(f, params, body)) = let
481                  val f' = cvtVar f                  val f' = cvtVar f
482                  val params' = cvtVars params                  val params' = cvtVars params
483                  val body' = pruneUnreachableCode (simplifyBlock errStrm body)                  val body' = simplifyAndPruneBlock errStrm body
484                  in                  in
485                    funcs := S.Func{f=f', params=params', body=body'} :: !funcs                    funcs := S.Func{f=f', params=params', body=body'} :: !funcs
486                  end                  end
487            in            val () = (
488              List.app simplifyConstDcl const_dcls;              List.app simplifyConstDcl const_dcls;
489              List.app simplifyInputDcl input_dcls;              List.app simplifyInputDcl input_dcls;
490              List.app simplifyGlobalDcl globals;                  List.app simplifyGlobalDcl globals)
491            (* make the global-initialization block *)
492              val globInit = (case globInit
493                     of SOME stm => mkBlock (simplifyStmt (errStrm, stm, !globalInit))
494                      | NONE => mkBlock (!globalInit)
495                    (* end case *))
496              in
497              S.Program{              S.Program{
498                  props = props,                  props = props,
499                  consts = List.rev(!consts'),                  consts = List.rev(!consts'),
500                  inputs = List.rev(!inputs'),                  inputs = List.rev(!inputs'),
501                  constInit = mkBlock (!constInit),                  constInit = mkBlock (!constInit),
502                  globals = List.rev(!globals'),                  globals = List.rev(!globals'),
503                  init = mkBlock (!globalInit),                  globInit = globInit,
504                  funcs = List.rev(!funcs),                  funcs = List.rev(!funcs),
505                  strand = simplifyStrand (errStrm, strand),                  strand = simplifyStrand (errStrm, strand),
506                  create = simplifyCreate (errStrm, create),                  create = Create.map (simplifyAndPruneBlock errStrm) create,
507                  update = Option.map (simplifyBlock errStrm) update                  init = Option.map (simplifyAndPruneBlock errStrm) init,
508                    update = Option.map (simplifyAndPruneBlock errStrm) update
509                }                }
510            end            end
511    

Legend:
Removed from v.3465  
changed lines
  Added in v.4113

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