{-# LANGUAGE FlexibleContexts #-}
{-# LANGUAGE TupleSections #-}
{-# LANGUAGE TypeFamilies #-}

-- | It is well known that fully parallel loops can always be
-- interchanged inwards with a sequential loop.  This module
-- implements that transformation.
--
-- This is also where we implement loop-switching (for branches),
-- which is semantically similar to interchange.
module Futhark.Pass.ExtractKernels.Interchange
  ( SeqLoop (..),
    interchangeLoops,
    Branch (..),
    interchangeBranch,
    WithAccStm (..),
    interchangeWithAcc,
  )
where

import Control.Monad.Identity
import Data.List (find)
import Data.Maybe
import Futhark.IR.SOACS
import Futhark.MonadFreshNames
import Futhark.Pass.ExtractKernels.Distribution
  ( KernelNest,
    LoopNesting (..),
    kernelNestLoops,
  )
import Futhark.Tools
import Futhark.Transform.Rename

-- | An encoding of a sequential do-loop with no existential context,
-- alongside its result pattern.
data SeqLoop = SeqLoop [Int] Pat [(FParam, SubExp)] (LoopForm SOACS) Body

seqLoopStm :: SeqLoop -> Stm
seqLoopStm :: SeqLoop -> Stm
seqLoopStm (SeqLoop [Int]
_ Pat
pat [(FParam, SubExp)]
merge LoopForm SOACS
form Body
body) =
  Pat -> StmAux (ExpDec SOACS) -> Exp SOACS -> Stm
forall rep. Pat rep -> StmAux (ExpDec rep) -> Exp rep -> Stm rep
Let Pat
pat (() -> StmAux ()
forall dec. dec -> StmAux dec
defAux ()) (Exp SOACS -> Stm) -> Exp SOACS -> Stm
forall a b. (a -> b) -> a -> b
$ [(FParam, SubExp)] -> LoopForm SOACS -> Body -> Exp SOACS
forall rep.
[(FParam rep, SubExp)] -> LoopForm rep -> BodyT rep -> ExpT rep
DoLoop [(FParam, SubExp)]
merge LoopForm SOACS
form Body
body

interchangeLoop ::
  (MonadBuilder m, LocalScope SOACS m) =>
  (VName -> Maybe VName) ->
  SeqLoop ->
  LoopNesting ->
  m SeqLoop
interchangeLoop :: (VName -> Maybe VName) -> SeqLoop -> LoopNesting -> m SeqLoop
interchangeLoop
  VName -> Maybe VName
isMapParameter
  (SeqLoop [Int]
perm Pat
loop_pat [(FParam, SubExp)]
merge LoopForm SOACS
form Body
body)
  (MapNesting PatT Type
pat StmAux ()
aux SubExp
w [(Param Type, VName)]
params_and_arrs) = do
    [(Param (TypeBase Shape Uniqueness), SubExp)]
merge_expanded <-
      Scope SOACS
-> m [(Param (TypeBase Shape Uniqueness), SubExp)]
-> m [(Param (TypeBase Shape Uniqueness), SubExp)]
forall rep (m :: * -> *) a.
LocalScope rep m =>
Scope rep -> m a -> m a
localScope ([Param Type] -> Scope SOACS
forall rep dec. (LParamInfo rep ~ dec) => [Param dec] -> Scope rep
scopeOfLParams ([Param Type] -> Scope SOACS) -> [Param Type] -> Scope SOACS
forall a b. (a -> b) -> a -> b
$ ((Param Type, VName) -> Param Type)
-> [(Param Type, VName)] -> [Param Type]
forall a b. (a -> b) -> [a] -> [b]
map (Param Type, VName) -> Param Type
forall a b. (a, b) -> a
fst [(Param Type, VName)]
params_and_arrs) (m [(Param (TypeBase Shape Uniqueness), SubExp)]
 -> m [(Param (TypeBase Shape Uniqueness), SubExp)])
-> m [(Param (TypeBase Shape Uniqueness), SubExp)]
-> m [(Param (TypeBase Shape Uniqueness), SubExp)]
forall a b. (a -> b) -> a -> b
$
        ((Param (TypeBase Shape Uniqueness), SubExp)
 -> m (Param (TypeBase Shape Uniqueness), SubExp))
-> [(Param (TypeBase Shape Uniqueness), SubExp)]
-> m [(Param (TypeBase Shape Uniqueness), SubExp)]
forall (t :: * -> *) (m :: * -> *) a b.
(Traversable t, Monad m) =>
(a -> m b) -> t a -> m (t b)
mapM (Param (TypeBase Shape Uniqueness), SubExp)
-> m (Param (TypeBase Shape Uniqueness), SubExp)
expand [(Param (TypeBase Shape Uniqueness), SubExp)]
[(FParam, SubExp)]
merge

    let loop_pat_expanded :: PatT Type
loop_pat_expanded =
          [PatElemT Type] -> PatT Type
forall dec. [PatElemT dec] -> PatT dec
Pat ([PatElemT Type] -> PatT Type) -> [PatElemT Type] -> PatT Type
forall a b. (a -> b) -> a -> b
$ (PatElemT Type -> PatElemT Type)
-> [PatElemT Type] -> [PatElemT Type]
forall a b. (a -> b) -> [a] -> [b]
map PatElemT Type -> PatElemT Type
expandPatElem ([PatElemT Type] -> [PatElemT Type])
-> [PatElemT Type] -> [PatElemT Type]
forall a b. (a -> b) -> a -> b
$ PatT Type -> [PatElemT Type]
forall dec. PatT dec -> [PatElemT dec]
patElems PatT Type
Pat
loop_pat
        new_params :: [Param Type]
new_params =
          [Attrs -> VName -> Type -> Param Type
forall dec. Attrs -> VName -> dec -> Param dec
Param Attrs
attrs VName
pname (Type -> Param Type) -> Type -> Param Type
forall a b. (a -> b) -> a -> b
$ TypeBase Shape Uniqueness -> Type
forall shape.
TypeBase shape Uniqueness -> TypeBase shape NoUniqueness
fromDecl TypeBase Shape Uniqueness
ptype | (Param Attrs
attrs VName
pname TypeBase Shape Uniqueness
ptype, SubExp
_) <- [(Param (TypeBase Shape Uniqueness), SubExp)]
[(FParam, SubExp)]
merge]
        new_arrs :: [VName]
new_arrs = ((Param (TypeBase Shape Uniqueness), SubExp) -> VName)
-> [(Param (TypeBase Shape Uniqueness), SubExp)] -> [VName]
forall a b. (a -> b) -> [a] -> [b]
map (Param (TypeBase Shape Uniqueness) -> VName
forall dec. Param dec -> VName
paramName (Param (TypeBase Shape Uniqueness) -> VName)
-> ((Param (TypeBase Shape Uniqueness), SubExp)
    -> Param (TypeBase Shape Uniqueness))
-> (Param (TypeBase Shape Uniqueness), SubExp)
-> VName
forall b c a. (b -> c) -> (a -> b) -> a -> c
. (Param (TypeBase Shape Uniqueness), SubExp)
-> Param (TypeBase Shape Uniqueness)
forall a b. (a, b) -> a
fst) [(Param (TypeBase Shape Uniqueness), SubExp)]
merge_expanded
        rettype :: [Type]
rettype = (Type -> Type) -> [Type] -> [Type]
forall a b. (a -> b) -> [a] -> [b]
map Type -> Type
forall shape u.
ArrayShape shape =>
TypeBase shape u -> TypeBase shape u
rowType ([Type] -> [Type]) -> [Type] -> [Type]
forall a b. (a -> b) -> a -> b
$ PatT Type -> [Type]
forall dec. Typed dec => PatT dec -> [Type]
patTypes PatT Type
loop_pat_expanded

    -- If the map consumes something that is bound outside the loop
    -- (i.e. is not a merge parameter), we have to copy() it.  As a
    -- small simplification, we just remove the parameter outright if
    -- it is not used anymore.  This might happen if the parameter was
    -- used just as the inital value of a merge parameter.
    (([Param Type]
params', [VName]
arrs'), Stms SOACS
pre_copy_stms) <-
      Builder SOACS ([Param Type], [VName])
-> m (([Param Type], [VName]), Stms SOACS)
forall (m :: * -> *) somerep rep a.
(MonadFreshNames m, HasScope somerep m, SameScope somerep rep) =>
Builder rep a -> m (a, Stms rep)
runBuilder (Builder SOACS ([Param Type], [VName])
 -> m (([Param Type], [VName]), Stms SOACS))
-> Builder SOACS ([Param Type], [VName])
-> m (([Param Type], [VName]), Stms SOACS)
forall a b. (a -> b) -> a -> b
$
        Scope SOACS
-> Builder SOACS ([Param Type], [VName])
-> Builder SOACS ([Param Type], [VName])
forall rep (m :: * -> *) a.
LocalScope rep m =>
Scope rep -> m a -> m a
localScope ([Param Type] -> Scope SOACS
forall rep dec. (LParamInfo rep ~ dec) => [Param dec] -> Scope rep
scopeOfLParams [Param Type]
new_params) (Builder SOACS ([Param Type], [VName])
 -> Builder SOACS ([Param Type], [VName]))
-> Builder SOACS ([Param Type], [VName])
-> Builder SOACS ([Param Type], [VName])
forall a b. (a -> b) -> a -> b
$
          [(Param Type, VName)] -> ([Param Type], [VName])
forall a b. [(a, b)] -> ([a], [b])
unzip ([(Param Type, VName)] -> ([Param Type], [VName]))
-> ([Maybe (Param Type, VName)] -> [(Param Type, VName)])
-> [Maybe (Param Type, VName)]
-> ([Param Type], [VName])
forall b c a. (b -> c) -> (a -> b) -> a -> c
. [Maybe (Param Type, VName)] -> [(Param Type, VName)]
forall a. [Maybe a] -> [a]
catMaybes ([Maybe (Param Type, VName)] -> ([Param Type], [VName]))
-> BuilderT SOACS (State VNameSource) [Maybe (Param Type, VName)]
-> Builder SOACS ([Param Type], [VName])
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> ((Param Type, VName)
 -> BuilderT SOACS (State VNameSource) (Maybe (Param Type, VName)))
-> [(Param Type, VName)]
-> BuilderT SOACS (State VNameSource) [Maybe (Param Type, VName)]
forall (t :: * -> *) (m :: * -> *) a b.
(Traversable t, Monad m) =>
(a -> m b) -> t a -> m (t b)
mapM (Param Type, VName)
-> BuilderT SOACS (State VNameSource) (Maybe (Param Type, VName))
copyOrRemoveParam [(Param Type, VName)]
params_and_arrs

    let lam :: LambdaT SOACS
lam = [LParam SOACS] -> Body -> [Type] -> LambdaT SOACS
forall rep. [LParam rep] -> BodyT rep -> [Type] -> LambdaT rep
Lambda ([Param Type]
params' [Param Type] -> [Param Type] -> [Param Type]
forall a. Semigroup a => a -> a -> a
<> [Param Type]
new_params) Body
body [Type]
rettype
        map_stm :: Stm
map_stm =
          Pat -> StmAux (ExpDec SOACS) -> Exp SOACS -> Stm
forall rep. Pat rep -> StmAux (ExpDec rep) -> Exp rep -> Stm rep
Let PatT Type
Pat
loop_pat_expanded StmAux ()
StmAux (ExpDec SOACS)
aux (Exp SOACS -> Stm) -> Exp SOACS -> Stm
forall a b. (a -> b) -> a -> b
$
            Op SOACS -> Exp SOACS
forall rep. Op rep -> ExpT rep
Op (Op SOACS -> Exp SOACS) -> Op SOACS -> Exp SOACS
forall a b. (a -> b) -> a -> b
$ SubExp -> [VName] -> ScremaForm SOACS -> SOAC SOACS
forall rep. SubExp -> [VName] -> ScremaForm rep -> SOAC rep
Screma SubExp
w ([VName]
arrs' [VName] -> [VName] -> [VName]
forall a. Semigroup a => a -> a -> a
<> [VName]
new_arrs) (LambdaT SOACS -> ScremaForm SOACS
forall rep. Lambda rep -> ScremaForm rep
mapSOAC LambdaT SOACS
lam)
        res :: Result
res = [VName] -> Result
varsRes ([VName] -> Result) -> [VName] -> Result
forall a b. (a -> b) -> a -> b
$ PatT Type -> [VName]
forall dec. PatT dec -> [VName]
patNames PatT Type
loop_pat_expanded
        pat' :: PatT Type
pat' = [PatElemT Type] -> PatT Type
forall dec. [PatElemT dec] -> PatT dec
Pat ([PatElemT Type] -> PatT Type) -> [PatElemT Type] -> PatT Type
forall a b. (a -> b) -> a -> b
$ [Int] -> [PatElemT Type] -> [PatElemT Type]
forall a. [Int] -> [a] -> [a]
rearrangeShape [Int]
perm ([PatElemT Type] -> [PatElemT Type])
-> [PatElemT Type] -> [PatElemT Type]
forall a b. (a -> b) -> a -> b
$ PatT Type -> [PatElemT Type]
forall dec. PatT dec -> [PatElemT dec]
patElems PatT Type
pat

    SeqLoop -> m SeqLoop
forall (m :: * -> *) a. Monad m => a -> m a
return (SeqLoop -> m SeqLoop) -> SeqLoop -> m SeqLoop
forall a b. (a -> b) -> a -> b
$
      [Int]
-> Pat -> [(FParam, SubExp)] -> LoopForm SOACS -> Body -> SeqLoop
SeqLoop [Int]
perm PatT Type
Pat
pat' [(Param (TypeBase Shape Uniqueness), SubExp)]
[(FParam, SubExp)]
merge_expanded LoopForm SOACS
form (Body -> SeqLoop) -> Body -> SeqLoop
forall a b. (a -> b) -> a -> b
$
        Stms SOACS -> Result -> Body
forall rep. Buildable rep => Stms rep -> Result -> Body rep
mkBody (Stms SOACS
pre_copy_stms Stms SOACS -> Stms SOACS -> Stms SOACS
forall a. Semigroup a => a -> a -> a
<> Stm -> Stms SOACS
forall rep. Stm rep -> Stms rep
oneStm Stm
map_stm) Result
res
    where
      free_in_body :: Names
free_in_body = Body -> Names
forall a. FreeIn a => a -> Names
freeIn Body
body

      copyOrRemoveParam :: (Param Type, VName)
-> BuilderT SOACS (State VNameSource) (Maybe (Param Type, VName))
copyOrRemoveParam (Param Type
param, VName
arr)
        | Bool -> Bool
not (Param Type -> VName
forall dec. Param dec -> VName
paramName Param Type
param VName -> Names -> Bool
`nameIn` Names
free_in_body) =
          Maybe (Param Type, VName)
-> BuilderT SOACS (State VNameSource) (Maybe (Param Type, VName))
forall (m :: * -> *) a. Monad m => a -> m a
return Maybe (Param Type, VName)
forall a. Maybe a
Nothing
        | Bool
otherwise =
          Maybe (Param Type, VName)
-> BuilderT SOACS (State VNameSource) (Maybe (Param Type, VName))
forall (m :: * -> *) a. Monad m => a -> m a
return (Maybe (Param Type, VName)
 -> BuilderT SOACS (State VNameSource) (Maybe (Param Type, VName)))
-> Maybe (Param Type, VName)
-> BuilderT SOACS (State VNameSource) (Maybe (Param Type, VName))
forall a b. (a -> b) -> a -> b
$ (Param Type, VName) -> Maybe (Param Type, VName)
forall a. a -> Maybe a
Just (Param Type
param, VName
arr)

      expandedInit :: [Char] -> SubExp -> m SubExp
expandedInit [Char]
_ (Var VName
v)
        | Just VName
arr <- VName -> Maybe VName
isMapParameter VName
v =
          SubExp -> m SubExp
forall (m :: * -> *) a. Monad m => a -> m a
return (SubExp -> m SubExp) -> SubExp -> m SubExp
forall a b. (a -> b) -> a -> b
$ VName -> SubExp
Var VName
arr
      expandedInit [Char]
param_name SubExp
se =
        [Char] -> Exp (Rep m) -> m SubExp
forall (m :: * -> *).
MonadBuilder m =>
[Char] -> Exp (Rep m) -> m SubExp
letSubExp ([Char]
param_name [Char] -> [Char] -> [Char]
forall a. Semigroup a => a -> a -> a
<> [Char]
"_expanded_init") (Exp (Rep m) -> m SubExp) -> Exp (Rep m) -> m SubExp
forall a b. (a -> b) -> a -> b
$
          BasicOp -> Exp (Rep m)
forall rep. BasicOp -> ExpT rep
BasicOp (BasicOp -> Exp (Rep m)) -> BasicOp -> Exp (Rep m)
forall a b. (a -> b) -> a -> b
$ Shape -> SubExp -> BasicOp
Replicate ([SubExp] -> Shape
forall d. [d] -> ShapeBase d
Shape [SubExp
w]) SubExp
se

      expand :: (Param (TypeBase Shape Uniqueness), SubExp)
-> m (Param (TypeBase Shape Uniqueness), SubExp)
expand (Param (TypeBase Shape Uniqueness)
merge_param, SubExp
merge_init) = do
        Param (TypeBase Shape Uniqueness)
expanded_param <-
          [Char]
-> TypeBase Shape Uniqueness
-> m (Param (TypeBase Shape Uniqueness))
forall (m :: * -> *) dec.
MonadFreshNames m =>
[Char] -> dec -> m (Param dec)
newParam ([Char]
param_name [Char] -> [Char] -> [Char]
forall a. Semigroup a => a -> a -> a
<> [Char]
"_expanded") (TypeBase Shape Uniqueness
 -> m (Param (TypeBase Shape Uniqueness)))
-> TypeBase Shape Uniqueness
-> m (Param (TypeBase Shape Uniqueness))
forall a b. (a -> b) -> a -> b
$
            TypeBase Shape Uniqueness
-> Shape -> Uniqueness -> TypeBase Shape Uniqueness
forall shape u_unused u.
ArrayShape shape =>
TypeBase shape u_unused -> shape -> u -> TypeBase shape u
arrayOf (Param (TypeBase Shape Uniqueness) -> TypeBase Shape Uniqueness
forall dec. DeclTyped dec => Param dec -> TypeBase Shape Uniqueness
paramDeclType Param (TypeBase Shape Uniqueness)
merge_param) ([SubExp] -> Shape
forall d. [d] -> ShapeBase d
Shape [SubExp
w]) (Uniqueness -> TypeBase Shape Uniqueness)
-> Uniqueness -> TypeBase Shape Uniqueness
forall a b. (a -> b) -> a -> b
$
              TypeBase Shape Uniqueness -> Uniqueness
forall shape. TypeBase shape Uniqueness -> Uniqueness
uniqueness (TypeBase Shape Uniqueness -> Uniqueness)
-> TypeBase Shape Uniqueness -> Uniqueness
forall a b. (a -> b) -> a -> b
$ Param (TypeBase Shape Uniqueness) -> TypeBase Shape Uniqueness
forall t. DeclTyped t => t -> TypeBase Shape Uniqueness
declTypeOf Param (TypeBase Shape Uniqueness)
merge_param
        SubExp
expanded_init <- [Char] -> SubExp -> m SubExp
expandedInit [Char]
param_name SubExp
merge_init
        (Param (TypeBase Shape Uniqueness), SubExp)
-> m (Param (TypeBase Shape Uniqueness), SubExp)
forall (m :: * -> *) a. Monad m => a -> m a
return (Param (TypeBase Shape Uniqueness)
expanded_param, SubExp
expanded_init)
        where
          param_name :: [Char]
param_name = VName -> [Char]
baseString (VName -> [Char]) -> VName -> [Char]
forall a b. (a -> b) -> a -> b
$ Param (TypeBase Shape Uniqueness) -> VName
forall dec. Param dec -> VName
paramName Param (TypeBase Shape Uniqueness)
merge_param

      expandPatElem :: PatElemT Type -> PatElemT Type
expandPatElem (PatElem VName
name Type
t) =
        VName -> Type -> PatElemT Type
forall dec. VName -> dec -> PatElemT dec
PatElem VName
name (Type -> PatElemT Type) -> Type -> PatElemT Type
forall a b. (a -> b) -> a -> b
$ Type -> SubExp -> Type
forall d.
ArrayShape (ShapeBase d) =>
TypeBase (ShapeBase d) NoUniqueness
-> d -> TypeBase (ShapeBase d) NoUniqueness
arrayOfRow Type
t SubExp
w

-- | Given a (parallel) map nesting and an inner sequential loop, move
-- the maps inside the sequential loop.  The result is several
-- statements - one of these will be the loop, which will then contain
-- statements with @map@ expressions.
interchangeLoops ::
  (MonadFreshNames m, HasScope SOACS m) =>
  KernelNest ->
  SeqLoop ->
  m (Stms SOACS)
interchangeLoops :: KernelNest -> SeqLoop -> m (Stms SOACS)
interchangeLoops KernelNest
nest SeqLoop
loop = do
  (SeqLoop
loop', Stms SOACS
stms) <-
    Builder SOACS SeqLoop -> m (SeqLoop, Stms SOACS)
forall (m :: * -> *) somerep rep a.
(MonadFreshNames m, HasScope somerep m, SameScope somerep rep) =>
Builder rep a -> m (a, Stms rep)
runBuilder (Builder SOACS SeqLoop -> m (SeqLoop, Stms SOACS))
-> Builder SOACS SeqLoop -> m (SeqLoop, Stms SOACS)
forall a b. (a -> b) -> a -> b
$
      (SeqLoop -> LoopNesting -> Builder SOACS SeqLoop)
-> SeqLoop -> [LoopNesting] -> Builder SOACS SeqLoop
forall (t :: * -> *) (m :: * -> *) b a.
(Foldable t, Monad m) =>
(b -> a -> m b) -> b -> t a -> m b
foldM ((VName -> Maybe VName)
-> SeqLoop -> LoopNesting -> Builder SOACS SeqLoop
forall (m :: * -> *).
(MonadBuilder m, LocalScope SOACS m) =>
(VName -> Maybe VName) -> SeqLoop -> LoopNesting -> m SeqLoop
interchangeLoop VName -> Maybe VName
isMapParameter) SeqLoop
loop ([LoopNesting] -> Builder SOACS SeqLoop)
-> [LoopNesting] -> Builder SOACS SeqLoop
forall a b. (a -> b) -> a -> b
$
        [LoopNesting] -> [LoopNesting]
forall a. [a] -> [a]
reverse ([LoopNesting] -> [LoopNesting]) -> [LoopNesting] -> [LoopNesting]
forall a b. (a -> b) -> a -> b
$ KernelNest -> [LoopNesting]
kernelNestLoops KernelNest
nest
  Stms SOACS -> m (Stms SOACS)
forall (m :: * -> *) a. Monad m => a -> m a
return (Stms SOACS -> m (Stms SOACS)) -> Stms SOACS -> m (Stms SOACS)
forall a b. (a -> b) -> a -> b
$ Stms SOACS
stms Stms SOACS -> Stms SOACS -> Stms SOACS
forall a. Semigroup a => a -> a -> a
<> Stm -> Stms SOACS
forall rep. Stm rep -> Stms rep
oneStm (SeqLoop -> Stm
seqLoopStm SeqLoop
loop')
  where
    isMapParameter :: VName -> Maybe VName
isMapParameter VName
v =
      ((Param Type, VName) -> VName)
-> Maybe (Param Type, VName) -> Maybe VName
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
fmap (Param Type, VName) -> VName
forall a b. (a, b) -> b
snd (Maybe (Param Type, VName) -> Maybe VName)
-> Maybe (Param Type, VName) -> Maybe VName
forall a b. (a -> b) -> a -> b
$
        ((Param Type, VName) -> Bool)
-> [(Param Type, VName)] -> Maybe (Param Type, VName)
forall (t :: * -> *) a. Foldable t => (a -> Bool) -> t a -> Maybe a
find ((VName -> VName -> Bool
forall a. Eq a => a -> a -> Bool
== VName
v) (VName -> Bool)
-> ((Param Type, VName) -> VName) -> (Param Type, VName) -> Bool
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Param Type -> VName
forall dec. Param dec -> VName
paramName (Param Type -> VName)
-> ((Param Type, VName) -> Param Type)
-> (Param Type, VName)
-> VName
forall b c a. (b -> c) -> (a -> b) -> a -> c
. (Param Type, VName) -> Param Type
forall a b. (a, b) -> a
fst) ([(Param Type, VName)] -> Maybe (Param Type, VName))
-> [(Param Type, VName)] -> Maybe (Param Type, VName)
forall a b. (a -> b) -> a -> b
$
          (LoopNesting -> [(Param Type, VName)])
-> [LoopNesting] -> [(Param Type, VName)]
forall (t :: * -> *) a b. Foldable t => (a -> [b]) -> t a -> [b]
concatMap LoopNesting -> [(Param Type, VName)]
loopNestingParamsAndArrs ([LoopNesting] -> [(Param Type, VName)])
-> [LoopNesting] -> [(Param Type, VName)]
forall a b. (a -> b) -> a -> b
$ KernelNest -> [LoopNesting]
kernelNestLoops KernelNest
nest

data Branch = Branch [Int] Pat SubExp Body Body (IfDec (BranchType SOACS))

branchStm :: Branch -> Stm
branchStm :: Branch -> Stm
branchStm (Branch [Int]
_ Pat
pat SubExp
cond Body
tbranch Body
fbranch IfDec (BranchType SOACS)
ret) =
  Pat -> StmAux (ExpDec SOACS) -> Exp SOACS -> Stm
forall rep. Pat rep -> StmAux (ExpDec rep) -> Exp rep -> Stm rep
Let Pat
pat (() -> StmAux ()
forall dec. dec -> StmAux dec
defAux ()) (Exp SOACS -> Stm) -> Exp SOACS -> Stm
forall a b. (a -> b) -> a -> b
$ SubExp -> Body -> Body -> IfDec (BranchType SOACS) -> Exp SOACS
forall rep.
SubExp
-> BodyT rep -> BodyT rep -> IfDec (BranchType rep) -> ExpT rep
If SubExp
cond Body
tbranch Body
fbranch IfDec (BranchType SOACS)
ret

interchangeBranch1 ::
  (MonadBuilder m) =>
  Branch ->
  LoopNesting ->
  m Branch
interchangeBranch1 :: Branch -> LoopNesting -> m Branch
interchangeBranch1
  (Branch [Int]
perm Pat
branch_pat SubExp
cond Body
tbranch Body
fbranch (IfDec [BranchType SOACS]
ret IfSort
if_sort))
  (MapNesting PatT Type
pat StmAux ()
aux SubExp
w [(Param Type, VName)]
params_and_arrs) = do
    let ret' :: [TypeBase (ShapeBase (Ext SubExp)) NoUniqueness]
ret' = (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
 -> TypeBase (ShapeBase (Ext SubExp)) NoUniqueness)
-> [TypeBase (ShapeBase (Ext SubExp)) NoUniqueness]
-> [TypeBase (ShapeBase (Ext SubExp)) NoUniqueness]
forall a b. (a -> b) -> [a] -> [b]
map (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
-> Ext SubExp -> TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
forall d.
ArrayShape (ShapeBase d) =>
TypeBase (ShapeBase d) NoUniqueness
-> d -> TypeBase (ShapeBase d) NoUniqueness
`arrayOfRow` SubExp -> Ext SubExp
forall a. a -> Ext a
Free SubExp
w) [TypeBase (ShapeBase (Ext SubExp)) NoUniqueness]
[BranchType SOACS]
ret
        pat' :: PatT Type
pat' = [PatElemT Type] -> PatT Type
forall dec. [PatElemT dec] -> PatT dec
Pat ([PatElemT Type] -> PatT Type) -> [PatElemT Type] -> PatT Type
forall a b. (a -> b) -> a -> b
$ [Int] -> [PatElemT Type] -> [PatElemT Type]
forall a. [Int] -> [a] -> [a]
rearrangeShape [Int]
perm ([PatElemT Type] -> [PatElemT Type])
-> [PatElemT Type] -> [PatElemT Type]
forall a b. (a -> b) -> a -> b
$ PatT Type -> [PatElemT Type]
forall dec. PatT dec -> [PatElemT dec]
patElems PatT Type
pat

        ([Param Type]
params, [VName]
arrs) = [(Param Type, VName)] -> ([Param Type], [VName])
forall a b. [(a, b)] -> ([a], [b])
unzip [(Param Type, VName)]
params_and_arrs
        lam_ret :: [Type]
lam_ret = [Int] -> [Type] -> [Type]
forall a. [Int] -> [a] -> [a]
rearrangeShape [Int]
perm ([Type] -> [Type]) -> [Type] -> [Type]
forall a b. (a -> b) -> a -> b
$ (Type -> Type) -> [Type] -> [Type]
forall a b. (a -> b) -> [a] -> [b]
map Type -> Type
forall shape u.
ArrayShape shape =>
TypeBase shape u -> TypeBase shape u
rowType ([Type] -> [Type]) -> [Type] -> [Type]
forall a b. (a -> b) -> a -> b
$ PatT Type -> [Type]
forall dec. Typed dec => PatT dec -> [Type]
patTypes PatT Type
pat

        branch_pat' :: PatT Type
branch_pat' =
          [PatElemT Type] -> PatT Type
forall dec. [PatElemT dec] -> PatT dec
Pat ([PatElemT Type] -> PatT Type) -> [PatElemT Type] -> PatT Type
forall a b. (a -> b) -> a -> b
$ (PatElemT Type -> PatElemT Type)
-> [PatElemT Type] -> [PatElemT Type]
forall a b. (a -> b) -> [a] -> [b]
map ((Type -> Type) -> PatElemT Type -> PatElemT Type
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
fmap (Type -> SubExp -> Type
forall d.
ArrayShape (ShapeBase d) =>
TypeBase (ShapeBase d) NoUniqueness
-> d -> TypeBase (ShapeBase d) NoUniqueness
`arrayOfRow` SubExp
w)) ([PatElemT Type] -> [PatElemT Type])
-> [PatElemT Type] -> [PatElemT Type]
forall a b. (a -> b) -> a -> b
$ PatT Type -> [PatElemT Type]
forall dec. PatT dec -> [PatElemT dec]
patElems PatT Type
Pat
branch_pat

        mkBranch :: Body -> m Body
mkBranch Body
branch = (Body -> m Body
forall rep (m :: * -> *).
(Renameable rep, MonadFreshNames m) =>
Body rep -> m (Body rep)
renameBody (Body -> m Body) -> m Body -> m Body
forall (m :: * -> *) a b. Monad m => (a -> m b) -> m a -> m b
=<<) (m Body -> m Body) -> m Body -> m Body
forall a b. (a -> b) -> a -> b
$ do
          let lam :: LambdaT SOACS
lam = [LParam SOACS] -> Body -> [Type] -> LambdaT SOACS
forall rep. [LParam rep] -> BodyT rep -> [Type] -> LambdaT rep
Lambda [Param Type]
[LParam SOACS]
params Body
branch [Type]
lam_ret
              res :: Result
res = [VName] -> Result
varsRes ([VName] -> Result) -> [VName] -> Result
forall a b. (a -> b) -> a -> b
$ PatT Type -> [VName]
forall dec. PatT dec -> [VName]
patNames PatT Type
branch_pat'
              map_stm :: Stm
map_stm = Pat -> StmAux (ExpDec SOACS) -> Exp SOACS -> Stm
forall rep. Pat rep -> StmAux (ExpDec rep) -> Exp rep -> Stm rep
Let PatT Type
Pat
branch_pat' StmAux ()
StmAux (ExpDec SOACS)
aux (Exp SOACS -> Stm) -> Exp SOACS -> Stm
forall a b. (a -> b) -> a -> b
$ Op SOACS -> Exp SOACS
forall rep. Op rep -> ExpT rep
Op (Op SOACS -> Exp SOACS) -> Op SOACS -> Exp SOACS
forall a b. (a -> b) -> a -> b
$ SubExp -> [VName] -> ScremaForm SOACS -> SOAC SOACS
forall rep. SubExp -> [VName] -> ScremaForm rep -> SOAC rep
Screma SubExp
w [VName]
arrs (ScremaForm SOACS -> SOAC SOACS) -> ScremaForm SOACS -> SOAC SOACS
forall a b. (a -> b) -> a -> b
$ LambdaT SOACS -> ScremaForm SOACS
forall rep. Lambda rep -> ScremaForm rep
mapSOAC LambdaT SOACS
lam
          Body -> m Body
forall (m :: * -> *) a. Monad m => a -> m a
return (Body -> m Body) -> Body -> m Body
forall a b. (a -> b) -> a -> b
$ Stms SOACS -> Result -> Body
forall rep. Buildable rep => Stms rep -> Result -> Body rep
mkBody (Stm -> Stms SOACS
forall rep. Stm rep -> Stms rep
oneStm Stm
map_stm) Result
res

    Body
tbranch' <- Body -> m Body
mkBranch Body
tbranch
    Body
fbranch' <- Body -> m Body
mkBranch Body
fbranch
    Branch -> m Branch
forall (m :: * -> *) a. Monad m => a -> m a
return (Branch -> m Branch) -> Branch -> m Branch
forall a b. (a -> b) -> a -> b
$
      [Int]
-> Pat
-> SubExp
-> Body
-> Body
-> IfDec (BranchType SOACS)
-> Branch
Branch [Int
0 .. PatT Type -> Int
forall dec. PatT dec -> Int
patSize PatT Type
pat Int -> Int -> Int
forall a. Num a => a -> a -> a
-Int
1] PatT Type
Pat
pat' SubExp
cond Body
tbranch' Body
fbranch' (IfDec (BranchType SOACS) -> Branch)
-> IfDec (BranchType SOACS) -> Branch
forall a b. (a -> b) -> a -> b
$
        [TypeBase (ShapeBase (Ext SubExp)) NoUniqueness]
-> IfSort -> IfDec (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness)
forall rt. [rt] -> IfSort -> IfDec rt
IfDec [TypeBase (ShapeBase (Ext SubExp)) NoUniqueness]
ret' IfSort
if_sort

interchangeBranch ::
  (MonadFreshNames m, HasScope SOACS m) =>
  KernelNest ->
  Branch ->
  m (Stms SOACS)
interchangeBranch :: KernelNest -> Branch -> m (Stms SOACS)
interchangeBranch KernelNest
nest Branch
loop = do
  (Branch
loop', Stms SOACS
stms) <-
    Builder SOACS Branch -> m (Branch, Stms SOACS)
forall (m :: * -> *) somerep rep a.
(MonadFreshNames m, HasScope somerep m, SameScope somerep rep) =>
Builder rep a -> m (a, Stms rep)
runBuilder (Builder SOACS Branch -> m (Branch, Stms SOACS))
-> Builder SOACS Branch -> m (Branch, Stms SOACS)
forall a b. (a -> b) -> a -> b
$ (Branch -> LoopNesting -> Builder SOACS Branch)
-> Branch -> [LoopNesting] -> Builder SOACS Branch
forall (t :: * -> *) (m :: * -> *) b a.
(Foldable t, Monad m) =>
(b -> a -> m b) -> b -> t a -> m b
foldM Branch -> LoopNesting -> Builder SOACS Branch
forall (m :: * -> *).
MonadBuilder m =>
Branch -> LoopNesting -> m Branch
interchangeBranch1 Branch
loop ([LoopNesting] -> Builder SOACS Branch)
-> [LoopNesting] -> Builder SOACS Branch
forall a b. (a -> b) -> a -> b
$ [LoopNesting] -> [LoopNesting]
forall a. [a] -> [a]
reverse ([LoopNesting] -> [LoopNesting]) -> [LoopNesting] -> [LoopNesting]
forall a b. (a -> b) -> a -> b
$ KernelNest -> [LoopNesting]
kernelNestLoops KernelNest
nest
  Stms SOACS -> m (Stms SOACS)
forall (m :: * -> *) a. Monad m => a -> m a
return (Stms SOACS -> m (Stms SOACS)) -> Stms SOACS -> m (Stms SOACS)
forall a b. (a -> b) -> a -> b
$ Stms SOACS
stms Stms SOACS -> Stms SOACS -> Stms SOACS
forall a. Semigroup a => a -> a -> a
<> Stm -> Stms SOACS
forall rep. Stm rep -> Stms rep
oneStm (Branch -> Stm
branchStm Branch
loop')

data WithAccStm
  = WithAccStm [Int] Pat [(Shape, [VName], Maybe (Lambda, [SubExp]))] Lambda

withAccStm :: WithAccStm -> Stm
withAccStm :: WithAccStm -> Stm
withAccStm (WithAccStm [Int]
_ Pat
pat [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs LambdaT SOACS
lam) =
  Pat -> StmAux (ExpDec SOACS) -> Exp SOACS -> Stm
forall rep. Pat rep -> StmAux (ExpDec rep) -> Exp rep -> Stm rep
Let Pat
pat (() -> StmAux ()
forall dec. dec -> StmAux dec
defAux ()) (Exp SOACS -> Stm) -> Exp SOACS -> Stm
forall a b. (a -> b) -> a -> b
$ [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
-> LambdaT SOACS -> Exp SOACS
forall rep. [WithAccInput rep] -> Lambda rep -> ExpT rep
WithAcc [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs LambdaT SOACS
lam

interchangeWithAcc1 ::
  (MonadBuilder m, Rep m ~ SOACS) =>
  WithAccStm ->
  LoopNesting ->
  m WithAccStm
interchangeWithAcc1 :: WithAccStm -> LoopNesting -> m WithAccStm
interchangeWithAcc1
  (WithAccStm [Int]
perm Pat
_withacc_pat [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs LambdaT SOACS
acc_lam)
  (MapNesting PatT Type
map_pat StmAux ()
map_aux SubExp
w [(Param Type, VName)]
params_and_arrs) = do
    [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs' <- ((Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))
 -> m (Shape, [VName], Maybe (LambdaT SOACS, [SubExp])))
-> [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
-> m [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
forall (t :: * -> *) (m :: * -> *) a b.
(Traversable t, Monad m) =>
(a -> m b) -> t a -> m (t b)
mapM (Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))
-> m (Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))
onInput [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs
    [Param Type]
lam_params' <- [Param Type] -> m [Param Type]
forall (m :: * -> *) a.
MonadFreshNames m =>
[Param a] -> m [Param a]
newAccLamParams ([Param Type] -> m [Param Type]) -> [Param Type] -> m [Param Type]
forall a b. (a -> b) -> a -> b
$ LambdaT SOACS -> [LParam SOACS]
forall rep. LambdaT rep -> [LParam rep]
lambdaParams LambdaT SOACS
acc_lam
    Param Type
iota_p <- [Char] -> Type -> m (Param Type)
forall (m :: * -> *) dec.
MonadFreshNames m =>
[Char] -> dec -> m (Param dec)
newParam [Char]
"iota_p" (Type -> m (Param Type)) -> Type -> m (Param Type)
forall a b. (a -> b) -> a -> b
$ PrimType -> Type
forall shape u. PrimType -> TypeBase shape u
Prim PrimType
int64
    LambdaT SOACS
acc_lam' <- SubExp -> LambdaT SOACS -> m (LambdaT SOACS)
trLam (VName -> SubExp
Var (Param Type -> VName
forall dec. Param dec -> VName
paramName Param Type
iota_p)) (LambdaT SOACS -> m (LambdaT SOACS))
-> (m Result -> m (LambdaT SOACS)) -> m Result -> m (LambdaT SOACS)
forall (m :: * -> *) b c a.
Monad m =>
(b -> m c) -> (a -> m b) -> a -> m c
<=< [LParam (Rep m)] -> m Result -> m (Lambda (Rep m))
forall (m :: * -> *).
MonadBuilder m =>
[LParam (Rep m)] -> m Result -> m (Lambda (Rep m))
mkLambda [Param Type]
[LParam (Rep m)]
lam_params' (m Result -> m (LambdaT SOACS)) -> m Result -> m (LambdaT SOACS)
forall a b. (a -> b) -> a -> b
$ do
      let acc_params :: [Param Type]
acc_params = Int -> [Param Type] -> [Param Type]
forall a. Int -> [a] -> [a]
drop ([(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))] -> Int
forall (t :: * -> *) a. Foldable t => t a -> Int
length [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs) [Param Type]
lam_params'
          orig_acc_params :: [Param Type]
orig_acc_params = Int -> [Param Type] -> [Param Type]
forall a. Int -> [a] -> [a]
drop ([(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))] -> Int
forall (t :: * -> *) a. Foldable t => t a -> Int
length [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs) ([Param Type] -> [Param Type]) -> [Param Type] -> [Param Type]
forall a b. (a -> b) -> a -> b
$ LambdaT SOACS -> [LParam SOACS]
forall rep. LambdaT rep -> [LParam rep]
lambdaParams LambdaT SOACS
acc_lam
      VName
iota_w <-
        [Char] -> Exp (Rep m) -> m VName
forall (m :: * -> *).
MonadBuilder m =>
[Char] -> Exp (Rep m) -> m VName
letExp [Char]
"acc_inter_iota" (Exp SOACS -> m VName)
-> (BasicOp -> Exp SOACS) -> BasicOp -> m VName
forall b c a. (b -> c) -> (a -> b) -> a -> c
. BasicOp -> Exp SOACS
forall rep. BasicOp -> ExpT rep
BasicOp (BasicOp -> m VName) -> BasicOp -> m VName
forall a b. (a -> b) -> a -> b
$
          SubExp -> SubExp -> SubExp -> IntType -> BasicOp
Iota SubExp
w (IntType -> Integer -> SubExp
intConst IntType
Int64 Integer
0) (IntType -> Integer -> SubExp
intConst IntType
Int64 Integer
1) IntType
Int64
      let ([Param Type]
params, [VName]
arrs) = [(Param Type, VName)] -> ([Param Type], [VName])
forall a b. [(a, b)] -> ([a], [b])
unzip [(Param Type, VName)]
params_and_arrs
          maplam_ret :: [Type]
maplam_ret = LambdaT SOACS -> [Type]
forall rep. LambdaT rep -> [Type]
lambdaReturnType LambdaT SOACS
acc_lam
          maplam :: LambdaT SOACS
maplam = [LParam SOACS] -> Body -> [Type] -> LambdaT SOACS
forall rep. [LParam rep] -> BodyT rep -> [Type] -> LambdaT rep
Lambda (Param Type
iota_p Param Type -> [Param Type] -> [Param Type]
forall a. a -> [a] -> [a]
: [Param Type]
orig_acc_params [Param Type] -> [Param Type] -> [Param Type]
forall a. [a] -> [a] -> [a]
++ [Param Type]
params) (LambdaT SOACS -> Body
forall rep. LambdaT rep -> BodyT rep
lambdaBody LambdaT SOACS
acc_lam) [Type]
maplam_ret
      StmAux () -> m Result -> m Result
forall (m :: * -> *) anyrep a.
MonadBuilder m =>
StmAux anyrep -> m a -> m a
auxing StmAux ()
map_aux (m Result -> m Result)
-> (Exp SOACS -> m Result) -> Exp SOACS -> m Result
forall b c a. (b -> c) -> (a -> b) -> a -> c
. ([SubExp] -> Result) -> m [SubExp] -> m Result
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
fmap [SubExp] -> Result
subExpsRes (m [SubExp] -> m Result)
-> (Exp SOACS -> m [SubExp]) -> Exp SOACS -> m Result
forall b c a. (b -> c) -> (a -> b) -> a -> c
. [Char] -> Exp (Rep m) -> m [SubExp]
forall (m :: * -> *).
MonadBuilder m =>
[Char] -> Exp (Rep m) -> m [SubExp]
letTupExp' [Char]
"withacc_inter" (Exp SOACS -> m Result) -> Exp SOACS -> m Result
forall a b. (a -> b) -> a -> b
$
        Op SOACS -> Exp SOACS
forall rep. Op rep -> ExpT rep
Op (Op SOACS -> Exp SOACS) -> Op SOACS -> Exp SOACS
forall a b. (a -> b) -> a -> b
$ SubExp -> [VName] -> ScremaForm SOACS -> SOAC SOACS
forall rep. SubExp -> [VName] -> ScremaForm rep -> SOAC rep
Screma SubExp
w (VName
iota_w VName -> [VName] -> [VName]
forall a. a -> [a] -> [a]
: (Param Type -> VName) -> [Param Type] -> [VName]
forall a b. (a -> b) -> [a] -> [b]
map Param Type -> VName
forall dec. Param dec -> VName
paramName [Param Type]
acc_params [VName] -> [VName] -> [VName]
forall a. [a] -> [a] -> [a]
++ [VName]
arrs) (LambdaT SOACS -> ScremaForm SOACS
forall rep. Lambda rep -> ScremaForm rep
mapSOAC LambdaT SOACS
maplam)
    let pat :: PatT Type
pat = [PatElemT Type] -> PatT Type
forall dec. [PatElemT dec] -> PatT dec
Pat ([PatElemT Type] -> PatT Type) -> [PatElemT Type] -> PatT Type
forall a b. (a -> b) -> a -> b
$ [Int] -> [PatElemT Type] -> [PatElemT Type]
forall a. [Int] -> [a] -> [a]
rearrangeShape [Int]
perm ([PatElemT Type] -> [PatElemT Type])
-> [PatElemT Type] -> [PatElemT Type]
forall a b. (a -> b) -> a -> b
$ PatT Type -> [PatElemT Type]
forall dec. PatT dec -> [PatElemT dec]
patElems PatT Type
map_pat
    WithAccStm -> m WithAccStm
forall (f :: * -> *) a. Applicative f => a -> f a
pure (WithAccStm -> m WithAccStm) -> WithAccStm -> m WithAccStm
forall a b. (a -> b) -> a -> b
$ [Int]
-> Pat
-> [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
-> LambdaT SOACS
-> WithAccStm
WithAccStm [Int]
perm PatT Type
Pat
pat [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs' LambdaT SOACS
acc_lam'
    where
      newAccLamParams :: [Param a] -> m [Param a]
newAccLamParams [Param a]
ps = do
        let ([Param a]
cert_ps, [Param a]
acc_ps) = Int -> [Param a] -> ([Param a], [Param a])
forall a. Int -> [a] -> ([a], [a])
splitAt ([Param a] -> Int
forall (t :: * -> *) a. Foldable t => t a -> Int
length [Param a]
ps Int -> Int -> Int
forall a. Integral a => a -> a -> a
`div` Int
2) [Param a]
ps
        -- Should not rename the certificates.
        [Param a]
acc_ps' <- [Param a] -> (Param a -> m (Param a)) -> m [Param a]
forall (t :: * -> *) (m :: * -> *) a b.
(Traversable t, Monad m) =>
t a -> (a -> m b) -> m (t b)
forM [Param a]
acc_ps ((Param a -> m (Param a)) -> m [Param a])
-> (Param a -> m (Param a)) -> m [Param a]
forall a b. (a -> b) -> a -> b
$ \(Param Attrs
attrs VName
v a
t) ->
          Attrs -> VName -> a -> Param a
forall dec. Attrs -> VName -> dec -> Param dec
Param Attrs
attrs (VName -> a -> Param a) -> m VName -> m (a -> Param a)
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> [Char] -> m VName
forall (m :: * -> *). MonadFreshNames m => [Char] -> m VName
newVName (VName -> [Char]
baseString VName
v) m (a -> Param a) -> m a -> m (Param a)
forall (f :: * -> *) a b. Applicative f => f (a -> b) -> f a -> f b
<*> a -> m a
forall (f :: * -> *) a. Applicative f => a -> f a
pure a
t
        [Param a] -> m [Param a]
forall (f :: * -> *) a. Applicative f => a -> f a
pure ([Param a] -> m [Param a]) -> [Param a] -> m [Param a]
forall a b. (a -> b) -> a -> b
$ [Param a]
cert_ps [Param a] -> [Param a] -> [Param a]
forall a. Semigroup a => a -> a -> a
<> [Param a]
acc_ps'

      num_accs :: Int
num_accs = [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))] -> Int
forall (t :: * -> *) a. Foldable t => t a -> Int
length [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
inputs
      acc_certs :: [VName]
acc_certs = (Param Type -> VName) -> [Param Type] -> [VName]
forall a b. (a -> b) -> [a] -> [b]
map Param Type -> VName
forall dec. Param dec -> VName
paramName ([Param Type] -> [VName]) -> [Param Type] -> [VName]
forall a b. (a -> b) -> a -> b
$ Int -> [Param Type] -> [Param Type]
forall a. Int -> [a] -> [a]
take Int
num_accs ([Param Type] -> [Param Type]) -> [Param Type] -> [Param Type]
forall a b. (a -> b) -> a -> b
$ LambdaT SOACS -> [LParam SOACS]
forall rep. LambdaT rep -> [LParam rep]
lambdaParams LambdaT SOACS
acc_lam
      onArr :: VName -> m VName
onArr VName
v =
        VName -> m VName
forall (f :: * -> *) a. Applicative f => a -> f a
pure (VName -> m VName)
-> (Maybe (Param Type, VName) -> VName)
-> Maybe (Param Type, VName)
-> m VName
forall b c a. (b -> c) -> (a -> b) -> a -> c
. VName
-> ((Param Type, VName) -> VName)
-> Maybe (Param Type, VName)
-> VName
forall b a. b -> (a -> b) -> Maybe a -> b
maybe VName
v (Param Type, VName) -> VName
forall a b. (a, b) -> b
snd (Maybe (Param Type, VName) -> m VName)
-> Maybe (Param Type, VName) -> m VName
forall a b. (a -> b) -> a -> b
$
          ((Param Type, VName) -> Bool)
-> [(Param Type, VName)] -> Maybe (Param Type, VName)
forall (t :: * -> *) a. Foldable t => (a -> Bool) -> t a -> Maybe a
find ((VName -> VName -> Bool
forall a. Eq a => a -> a -> Bool
== VName
v) (VName -> Bool)
-> ((Param Type, VName) -> VName) -> (Param Type, VName) -> Bool
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Param Type -> VName
forall dec. Param dec -> VName
paramName (Param Type -> VName)
-> ((Param Type, VName) -> Param Type)
-> (Param Type, VName)
-> VName
forall b c a. (b -> c) -> (a -> b) -> a -> c
. (Param Type, VName) -> Param Type
forall a b. (a, b) -> a
fst) [(Param Type, VName)]
params_and_arrs
      onInput :: (Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))
-> m (Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))
onInput (Shape
shape, [VName]
arrs, Maybe (LambdaT SOACS, [SubExp])
op) =
        ([SubExp] -> Shape
forall d. [d] -> ShapeBase d
Shape [SubExp
w] Shape -> Shape -> Shape
forall a. Semigroup a => a -> a -> a
<> Shape
shape,,) ([VName]
 -> Maybe (LambdaT SOACS, [SubExp])
 -> (Shape, [VName], Maybe (LambdaT SOACS, [SubExp])))
-> m [VName]
-> m (Maybe (LambdaT SOACS, [SubExp])
      -> (Shape, [VName], Maybe (LambdaT SOACS, [SubExp])))
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> (VName -> m VName) -> [VName] -> m [VName]
forall (t :: * -> *) (m :: * -> *) a b.
(Traversable t, Monad m) =>
(a -> m b) -> t a -> m (t b)
mapM VName -> m VName
onArr [VName]
arrs m (Maybe (LambdaT SOACS, [SubExp])
   -> (Shape, [VName], Maybe (LambdaT SOACS, [SubExp])))
-> m (Maybe (LambdaT SOACS, [SubExp]))
-> m (Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))
forall (f :: * -> *) a b. Applicative f => f (a -> b) -> f a -> f b
<*> ((LambdaT SOACS, [SubExp]) -> m (LambdaT SOACS, [SubExp]))
-> Maybe (LambdaT SOACS, [SubExp])
-> m (Maybe (LambdaT SOACS, [SubExp]))
forall (t :: * -> *) (f :: * -> *) a b.
(Traversable t, Applicative f) =>
(a -> f b) -> t a -> f (t b)
traverse (LambdaT SOACS, [SubExp]) -> m (LambdaT SOACS, [SubExp])
forall (m :: * -> *) rep shape u b.
(MonadFreshNames m, LParamInfo rep ~ TypeBase shape u) =>
(LambdaT rep, b) -> m (LambdaT rep, b)
onOp Maybe (LambdaT SOACS, [SubExp])
op

      onOp :: (LambdaT rep, b) -> m (LambdaT rep, b)
onOp (LambdaT rep
op_lam, b
nes) = do
        -- We need to add an additional index parameter because we are
        -- extending the index space of the accumulator.
        Param (TypeBase shape u)
idx_p <- [Char] -> TypeBase shape u -> m (Param (TypeBase shape u))
forall (m :: * -> *) dec.
MonadFreshNames m =>
[Char] -> dec -> m (Param dec)
newParam [Char]
"idx" (TypeBase shape u -> m (Param (TypeBase shape u)))
-> TypeBase shape u -> m (Param (TypeBase shape u))
forall a b. (a -> b) -> a -> b
$ PrimType -> TypeBase shape u
forall shape u. PrimType -> TypeBase shape u
Prim PrimType
int64
        (LambdaT rep, b) -> m (LambdaT rep, b)
forall (f :: * -> *) a. Applicative f => a -> f a
pure (LambdaT rep
op_lam {lambdaParams :: [LParam rep]
lambdaParams = Param (TypeBase shape u)
idx_p Param (TypeBase shape u)
-> [Param (TypeBase shape u)] -> [Param (TypeBase shape u)]
forall a. a -> [a] -> [a]
: LambdaT rep -> [LParam rep]
forall rep. LambdaT rep -> [LParam rep]
lambdaParams LambdaT rep
op_lam}, b
nes)

      trType :: TypeBase shape u -> TypeBase shape u
      trType :: TypeBase shape u -> TypeBase shape u
trType (Acc VName
acc Shape
ispace [Type]
ts u
u)
        | VName
acc VName -> [VName] -> Bool
forall (t :: * -> *) a. (Foldable t, Eq a) => a -> t a -> Bool
`elem` [VName]
acc_certs =
          VName -> Shape -> [Type] -> u -> TypeBase shape u
forall shape u. VName -> Shape -> [Type] -> u -> TypeBase shape u
Acc VName
acc ([SubExp] -> Shape
forall d. [d] -> ShapeBase d
Shape [SubExp
w] Shape -> Shape -> Shape
forall a. Semigroup a => a -> a -> a
<> Shape
ispace) [Type]
ts u
u
      trType TypeBase shape u
t = TypeBase shape u
t

      trParam :: Param (TypeBase shape u) -> Param (TypeBase shape u)
      trParam :: Param (TypeBase shape u) -> Param (TypeBase shape u)
trParam = (TypeBase shape u -> TypeBase shape u)
-> Param (TypeBase shape u) -> Param (TypeBase shape u)
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
fmap TypeBase shape u -> TypeBase shape u
forall shape u. TypeBase shape u -> TypeBase shape u
trType

      trLam :: SubExp -> LambdaT SOACS -> m (LambdaT SOACS)
trLam SubExp
i (Lambda [LParam SOACS]
params Body
body [Type]
ret) =
        Scope SOACS -> m (LambdaT SOACS) -> m (LambdaT SOACS)
forall rep (m :: * -> *) a.
LocalScope rep m =>
Scope rep -> m a -> m a
localScope ([Param Type] -> Scope SOACS
forall rep dec. (LParamInfo rep ~ dec) => [Param dec] -> Scope rep
scopeOfLParams [Param Type]
[LParam SOACS]
params) (m (LambdaT SOACS) -> m (LambdaT SOACS))
-> m (LambdaT SOACS) -> m (LambdaT SOACS)
forall a b. (a -> b) -> a -> b
$
          [LParam SOACS] -> Body -> [Type] -> LambdaT SOACS
forall rep. [LParam rep] -> BodyT rep -> [Type] -> LambdaT rep
Lambda ((Param Type -> Param Type) -> [Param Type] -> [Param Type]
forall a b. (a -> b) -> [a] -> [b]
map Param Type -> Param Type
forall shape u.
Param (TypeBase shape u) -> Param (TypeBase shape u)
trParam [Param Type]
[LParam SOACS]
params) (Body -> [Type] -> LambdaT SOACS)
-> m Body -> m ([Type] -> LambdaT SOACS)
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> SubExp -> Body -> m Body
trBody SubExp
i Body
body m ([Type] -> LambdaT SOACS) -> m [Type] -> m (LambdaT SOACS)
forall (f :: * -> *) a b. Applicative f => f (a -> b) -> f a -> f b
<*> [Type] -> m [Type]
forall (f :: * -> *) a. Applicative f => a -> f a
pure ((Type -> Type) -> [Type] -> [Type]
forall a b. (a -> b) -> [a] -> [b]
map Type -> Type
forall shape u. TypeBase shape u -> TypeBase shape u
trType [Type]
ret)

      trBody :: SubExp -> Body -> m Body
trBody SubExp
i (Body BodyDec SOACS
dec Stms SOACS
stms Result
res) =
        Stms SOACS -> m Body -> m Body
forall rep a (m :: * -> *) b.
(Scoped rep a, LocalScope rep m) =>
a -> m b -> m b
inScopeOf Stms SOACS
stms (m Body -> m Body) -> m Body -> m Body
forall a b. (a -> b) -> a -> b
$ BodyDec SOACS -> Stms SOACS -> Result -> Body
forall rep. BodyDec rep -> Stms rep -> Result -> BodyT rep
Body BodyDec SOACS
dec (Stms SOACS -> Result -> Body)
-> m (Stms SOACS) -> m (Result -> Body)
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> (Stm -> m Stm) -> Stms SOACS -> m (Stms SOACS)
forall (t :: * -> *) (f :: * -> *) a b.
(Traversable t, Applicative f) =>
(a -> f b) -> t a -> f (t b)
traverse (SubExp -> Stm -> m Stm
trStm SubExp
i) Stms SOACS
stms m (Result -> Body) -> m Result -> m Body
forall (f :: * -> *) a b. Applicative f => f (a -> b) -> f a -> f b
<*> Result -> m Result
forall (f :: * -> *) a. Applicative f => a -> f a
pure Result
res

      trStm :: SubExp -> Stm -> m Stm
trStm SubExp
i (Let Pat
pat StmAux (ExpDec SOACS)
aux Exp SOACS
e) =
        Pat -> StmAux (ExpDec SOACS) -> Exp SOACS -> Stm
forall rep. Pat rep -> StmAux (ExpDec rep) -> Exp rep -> Stm rep
Let ((Type -> Type) -> PatT Type -> PatT Type
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
fmap Type -> Type
forall shape u. TypeBase shape u -> TypeBase shape u
trType PatT Type
Pat
pat) StmAux (ExpDec SOACS)
aux (Exp SOACS -> Stm) -> m (Exp SOACS) -> m Stm
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> SubExp -> Exp SOACS -> m (Exp SOACS)
trExp SubExp
i Exp SOACS
e

      trSOAC :: SubExp -> SOAC SOACS -> m (SOAC SOACS)
trSOAC SubExp
i = SOACMapper SOACS SOACS m -> SOAC SOACS -> m (SOAC SOACS)
forall (m :: * -> *) frep trep.
(Applicative m, Monad m) =>
SOACMapper frep trep m -> SOAC frep -> m (SOAC trep)
mapSOACM SOACMapper SOACS SOACS m
mapper
        where
          mapper :: SOACMapper SOACS SOACS m
mapper =
            SOACMapper Any Any m
forall (m :: * -> *) rep. Monad m => SOACMapper rep rep m
identitySOACMapper {mapOnSOACLambda :: LambdaT SOACS -> m (LambdaT SOACS)
mapOnSOACLambda = SubExp -> LambdaT SOACS -> m (LambdaT SOACS)
trLam SubExp
i}

      trExp :: SubExp -> Exp SOACS -> m (Exp SOACS)
trExp SubExp
i (WithAcc [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
acc_inputs LambdaT SOACS
lam) =
        [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
-> LambdaT SOACS -> Exp SOACS
forall rep. [WithAccInput rep] -> Lambda rep -> ExpT rep
WithAcc [(Shape, [VName], Maybe (LambdaT SOACS, [SubExp]))]
acc_inputs (LambdaT SOACS -> Exp SOACS) -> m (LambdaT SOACS) -> m (Exp SOACS)
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> SubExp -> LambdaT SOACS -> m (LambdaT SOACS)
trLam SubExp
i LambdaT SOACS
lam
      trExp SubExp
i (BasicOp (UpdateAcc VName
acc [SubExp]
is [SubExp]
ses)) = do
        Type
acc_t <- VName -> m Type
forall rep (m :: * -> *). HasScope rep m => VName -> m Type
lookupType VName
acc
        Exp SOACS -> m (Exp SOACS)
forall (f :: * -> *) a. Applicative f => a -> f a
pure (Exp SOACS -> m (Exp SOACS)) -> Exp SOACS -> m (Exp SOACS)
forall a b. (a -> b) -> a -> b
$ case Type
acc_t of
          Acc VName
cert Shape
_ [Type]
_ NoUniqueness
_
            | VName
cert VName -> [VName] -> Bool
forall (t :: * -> *) a. (Foldable t, Eq a) => a -> t a -> Bool
`elem` [VName]
acc_certs ->
              BasicOp -> Exp SOACS
forall rep. BasicOp -> ExpT rep
BasicOp (BasicOp -> Exp SOACS) -> BasicOp -> Exp SOACS
forall a b. (a -> b) -> a -> b
$ VName -> [SubExp] -> [SubExp] -> BasicOp
UpdateAcc VName
acc (SubExp
i SubExp -> [SubExp] -> [SubExp]
forall a. a -> [a] -> [a]
: [SubExp]
is) [SubExp]
ses
          Type
_ ->
            BasicOp -> Exp SOACS
forall rep. BasicOp -> ExpT rep
BasicOp (BasicOp -> Exp SOACS) -> BasicOp -> Exp SOACS
forall a b. (a -> b) -> a -> b
$ VName -> [SubExp] -> [SubExp] -> BasicOp
UpdateAcc VName
acc [SubExp]
is [SubExp]
ses
      trExp SubExp
i Exp SOACS
e = Mapper SOACS SOACS m -> Exp SOACS -> m (Exp SOACS)
forall (m :: * -> *) frep trep.
(Applicative m, Monad m) =>
Mapper frep trep m -> Exp frep -> m (Exp trep)
mapExpM Mapper SOACS SOACS m
mapper Exp SOACS
e
        where
          mapper :: Mapper SOACS SOACS m
mapper =
            Mapper Any Any m
forall (m :: * -> *) rep. Monad m => Mapper rep rep m
identityMapper
              { mapOnBody :: Scope SOACS -> Body -> m Body
mapOnBody = \Scope SOACS
scope -> Scope SOACS -> m Body -> m Body
forall rep (m :: * -> *) a.
LocalScope rep m =>
Scope rep -> m a -> m a
localScope Scope SOACS
scope (m Body -> m Body) -> (Body -> m Body) -> Body -> m Body
forall b c a. (b -> c) -> (a -> b) -> a -> c
. SubExp -> Body -> m Body
trBody SubExp
i,
                mapOnRetType :: RetType SOACS -> m (RetType SOACS)
mapOnRetType = TypeBase (ShapeBase (Ext SubExp)) Uniqueness
-> m (TypeBase (ShapeBase (Ext SubExp)) Uniqueness)
forall (f :: * -> *) a. Applicative f => a -> f a
pure (TypeBase (ShapeBase (Ext SubExp)) Uniqueness
 -> m (TypeBase (ShapeBase (Ext SubExp)) Uniqueness))
-> (TypeBase (ShapeBase (Ext SubExp)) Uniqueness
    -> TypeBase (ShapeBase (Ext SubExp)) Uniqueness)
-> TypeBase (ShapeBase (Ext SubExp)) Uniqueness
-> m (TypeBase (ShapeBase (Ext SubExp)) Uniqueness)
forall b c a. (b -> c) -> (a -> b) -> a -> c
. TypeBase (ShapeBase (Ext SubExp)) Uniqueness
-> TypeBase (ShapeBase (Ext SubExp)) Uniqueness
forall shape u. TypeBase shape u -> TypeBase shape u
trType,
                mapOnBranchType :: BranchType SOACS -> m (BranchType SOACS)
mapOnBranchType = TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
-> m (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness)
forall (f :: * -> *) a. Applicative f => a -> f a
pure (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
 -> m (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness))
-> (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
    -> TypeBase (ShapeBase (Ext SubExp)) NoUniqueness)
-> TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
-> m (TypeBase (ShapeBase (Ext SubExp)) NoUniqueness)
forall b c a. (b -> c) -> (a -> b) -> a -> c
. TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
-> TypeBase (ShapeBase (Ext SubExp)) NoUniqueness
forall shape u. TypeBase shape u -> TypeBase shape u
trType,
                mapOnFParam :: FParam -> m FParam
mapOnFParam = Param (TypeBase Shape Uniqueness)
-> m (Param (TypeBase Shape Uniqueness))
forall (f :: * -> *) a. Applicative f => a -> f a
pure (Param (TypeBase Shape Uniqueness)
 -> m (Param (TypeBase Shape Uniqueness)))
-> (Param (TypeBase Shape Uniqueness)
    -> Param (TypeBase Shape Uniqueness))
-> Param (TypeBase Shape Uniqueness)
-> m (Param (TypeBase Shape Uniqueness))
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Param (TypeBase Shape Uniqueness)
-> Param (TypeBase Shape Uniqueness)
forall shape u.
Param (TypeBase shape u) -> Param (TypeBase shape u)
trParam,
                mapOnLParam :: LParam SOACS -> m (LParam SOACS)
mapOnLParam = Param Type -> m (Param Type)
forall (f :: * -> *) a. Applicative f => a -> f a
pure (Param Type -> m (Param Type))
-> (Param Type -> Param Type) -> Param Type -> m (Param Type)
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Param Type -> Param Type
forall shape u.
Param (TypeBase shape u) -> Param (TypeBase shape u)
trParam,
                mapOnOp :: Op SOACS -> m (Op SOACS)
mapOnOp = SubExp -> SOAC SOACS -> m (SOAC SOACS)
trSOAC SubExp
i
              }

interchangeWithAcc ::
  (MonadFreshNames m, HasScope SOACS m) =>
  KernelNest ->
  WithAccStm ->
  m (Stms SOACS)
interchangeWithAcc :: KernelNest -> WithAccStm -> m (Stms SOACS)
interchangeWithAcc KernelNest
nest WithAccStm
withacc = do
  (WithAccStm
withacc', Stms SOACS
stms) <-
    Builder SOACS WithAccStm -> m (WithAccStm, Stms SOACS)
forall (m :: * -> *) somerep rep a.
(MonadFreshNames m, HasScope somerep m, SameScope somerep rep) =>
Builder rep a -> m (a, Stms rep)
runBuilder (Builder SOACS WithAccStm -> m (WithAccStm, Stms SOACS))
-> Builder SOACS WithAccStm -> m (WithAccStm, Stms SOACS)
forall a b. (a -> b) -> a -> b
$ (WithAccStm -> LoopNesting -> Builder SOACS WithAccStm)
-> WithAccStm -> [LoopNesting] -> Builder SOACS WithAccStm
forall (t :: * -> *) (m :: * -> *) b a.
(Foldable t, Monad m) =>
(b -> a -> m b) -> b -> t a -> m b
foldM WithAccStm -> LoopNesting -> Builder SOACS WithAccStm
forall (m :: * -> *).
(MonadBuilder m, Rep m ~ SOACS) =>
WithAccStm -> LoopNesting -> m WithAccStm
interchangeWithAcc1 WithAccStm
withacc ([LoopNesting] -> Builder SOACS WithAccStm)
-> [LoopNesting] -> Builder SOACS WithAccStm
forall a b. (a -> b) -> a -> b
$ [LoopNesting] -> [LoopNesting]
forall a. [a] -> [a]
reverse ([LoopNesting] -> [LoopNesting]) -> [LoopNesting] -> [LoopNesting]
forall a b. (a -> b) -> a -> b
$ KernelNest -> [LoopNesting]
kernelNestLoops KernelNest
nest
  Stms SOACS -> m (Stms SOACS)
forall (m :: * -> *) a. Monad m => a -> m a
return (Stms SOACS -> m (Stms SOACS)) -> Stms SOACS -> m (Stms SOACS)
forall a b. (a -> b) -> a -> b
$ Stms SOACS
stms Stms SOACS -> Stms SOACS -> Stms SOACS
forall a. Semigroup a => a -> a -> a
<> Stm -> Stms SOACS
forall rep. Stm rep -> Stms rep
oneStm (WithAccStm -> Stm
withAccStm WithAccStm
withacc')