module PrimitiveExtras.Folds
(
  indexCounts,
  unliftedArray,
  primMultiArray,
)
where

import PrimitiveExtras.Prelude hiding (fold, foldM)
import PrimitiveExtras.Types
import Control.Foldl
import qualified PrimitiveExtras.UnliftedArray as UA


unsafeIO :: (state -> input -> IO state) -> IO state -> (state -> IO output) -> Fold input output
unsafeIO :: (state -> input -> IO state)
-> IO state -> (state -> IO output) -> Fold input output
unsafeIO state -> input -> IO state
stepInIO IO state
initInIO state -> IO output
extractInIO =
  (state -> input -> state)
-> state -> (state -> output) -> Fold input output
forall a b x. (x -> a -> x) -> x -> (x -> b) -> Fold a b
Fold
    (\ !state
state input
input -> IO state -> state
forall a. IO a -> a
unsafeDupablePerformIO (state -> input -> IO state
stepInIO state
state input
input))
    (IO state -> state
forall a. IO a -> a
unsafeDupablePerformIO IO state
initInIO)
    (\ state
state -> let !output :: output
output = IO output -> output
forall a. IO a -> a
unsafePerformIO (state -> IO output
extractInIO state
state) in output
output)

foldMInUnsafeDupableIO :: FoldM IO input output -> Fold input output
foldMInUnsafeDupableIO :: FoldM IO input output -> Fold input output
foldMInUnsafeDupableIO (FoldM x -> input -> IO x
step IO x
init x -> IO output
extract) = (x -> input -> IO x)
-> IO x -> (x -> IO output) -> Fold input output
forall state input output.
(state -> input -> IO state)
-> IO state -> (state -> IO output) -> Fold input output
unsafeIO x -> input -> IO x
step IO x
init x -> IO output
extract

{-|
Given a size of the array,
construct a fold, which produces an array of index counts.
-}
indexCounts :: (Integral count, Prim count) => Int {-^ Array size -} -> Fold Int (PrimArray count)
indexCounts :: Int -> Fold Int (PrimArray count)
indexCounts Int
size = (MutablePrimArray RealWorld count
 -> Int -> IO (MutablePrimArray RealWorld count))
-> IO (MutablePrimArray RealWorld count)
-> (MutablePrimArray RealWorld count -> IO (PrimArray count))
-> Fold Int (PrimArray count)
forall state input output.
(state -> input -> IO state)
-> IO state -> (state -> IO output) -> Fold input output
unsafeIO MutablePrimArray RealWorld count
-> Int -> IO (MutablePrimArray RealWorld count)
forall (m :: * -> *) a.
(Prim a, PrimMonad m, Enum a) =>
MutablePrimArray (PrimState m) a
-> Int -> m (MutablePrimArray (PrimState m) a)
step IO (MutablePrimArray RealWorld count)
IO (MutablePrimArray (PrimState IO) count)
init MutablePrimArray RealWorld count -> IO (PrimArray count)
forall (m :: * -> *) a.
PrimMonad m =>
MutablePrimArray (PrimState m) a -> m (PrimArray a)
extract where
  init :: IO (MutablePrimArray (PrimState IO) count)
init = PrimArray count -> IO (MutablePrimArray (PrimState IO) count)
forall (m :: * -> *) a.
PrimMonad m =>
PrimArray a -> m (MutablePrimArray (PrimState m) a)
unsafeThawPrimArray (Int -> count -> PrimArray count
forall a. Prim a => Int -> a -> PrimArray a
replicatePrimArray Int
size count
0)
  step :: MutablePrimArray (PrimState m) a
-> Int -> m (MutablePrimArray (PrimState m) a)
step MutablePrimArray (PrimState m) a
mutable Int
i = do
    a
count <- MutablePrimArray (PrimState m) a -> Int -> m a
forall a (m :: * -> *).
(Prim a, PrimMonad m) =>
MutablePrimArray (PrimState m) a -> Int -> m a
readPrimArray MutablePrimArray (PrimState m) a
mutable Int
i
    MutablePrimArray (PrimState m) a -> Int -> a -> m ()
forall a (m :: * -> *).
(Prim a, PrimMonad m) =>
MutablePrimArray (PrimState m) a -> Int -> a -> m ()
writePrimArray MutablePrimArray (PrimState m) a
mutable Int
i (a -> a
forall a. Enum a => a -> a
succ a
count)
    MutablePrimArray (PrimState m) a
-> m (MutablePrimArray (PrimState m) a)
forall (m :: * -> *) a. Monad m => a -> m a
return MutablePrimArray (PrimState m) a
mutable
  extract :: MutablePrimArray (PrimState m) a -> m (PrimArray a)
extract = MutablePrimArray (PrimState m) a -> m (PrimArray a)
forall (m :: * -> *) a.
PrimMonad m =>
MutablePrimArray (PrimState m) a -> m (PrimArray a)
unsafeFreezePrimArray

{-|
This function is partial in the sense that it expects the
index vector of produced elements to be within the specified amount.
-}
unliftedArray :: PrimUnlifted element => Int {-^ Size of the array -} -> Fold (Int, element) (UnliftedArray element)
unliftedArray :: Int -> Fold (Int, element) (UnliftedArray element)
unliftedArray Int
size =
  (MutableUnliftedArray RealWorld element
 -> (Int, element) -> IO (MutableUnliftedArray RealWorld element))
-> IO (MutableUnliftedArray RealWorld element)
-> (MutableUnliftedArray RealWorld element
    -> IO (UnliftedArray element))
-> Fold (Int, element) (UnliftedArray element)
forall state input output.
(state -> input -> IO state)
-> IO state -> (state -> IO output) -> Fold input output
unsafeIO MutableUnliftedArray RealWorld element
-> (Int, element) -> IO (MutableUnliftedArray RealWorld element)
forall (f :: * -> *) a.
(PrimMonad f, PrimUnlifted a) =>
MutableUnliftedArray (PrimState f) a
-> (Int, a) -> f (MutableUnliftedArray (PrimState f) a)
step IO (MutableUnliftedArray RealWorld element)
IO (MutableUnliftedArray (PrimState IO) element)
init MutableUnliftedArray RealWorld element
-> IO (UnliftedArray element)
forall (m :: * -> *) a.
PrimMonad m =>
MutableUnliftedArray (PrimState m) a -> m (UnliftedArray a)
extract
  where
    step :: MutableUnliftedArray (PrimState f) a
-> (Int, a) -> f (MutableUnliftedArray (PrimState f) a)
step MutableUnliftedArray (PrimState f) a
mutable (Int
index, a
element) =
      MutableUnliftedArray (PrimState f) a -> Int -> a -> f ()
forall (m :: * -> *) a.
(PrimMonad m, PrimUnlifted a) =>
MutableUnliftedArray (PrimState m) a -> Int -> a -> m ()
writeUnliftedArray MutableUnliftedArray (PrimState f) a
mutable Int
index a
element f ()
-> MutableUnliftedArray (PrimState f) a
-> f (MutableUnliftedArray (PrimState f) a)
forall (f :: * -> *) a b. Functor f => f a -> b -> f b
$> MutableUnliftedArray (PrimState f) a
mutable
    init :: IO (MutableUnliftedArray (PrimState IO) element)
init =
      Int -> IO (MutableUnliftedArray (PrimState IO) element)
forall (m :: * -> *) a.
PrimMonad m =>
Int -> m (MutableUnliftedArray (PrimState m) a)
unsafeNewUnliftedArray Int
size
    extract :: MutableUnliftedArray (PrimState m) a -> m (UnliftedArray a)
extract =
      MutableUnliftedArray (PrimState m) a -> m (UnliftedArray a)
forall (m :: * -> *) a.
PrimMonad m =>
MutableUnliftedArray (PrimState m) a -> m (UnliftedArray a)
unsafeFreezeUnliftedArray

{-|
Having a priorly computed array of inner dimension sizes,
e.g., using the 'indexCounts' fold,
construct a fold over indexed elements into a multi-array of elements.

Thus it allows to construct it in two passes over the indexed elements.
-}
primMultiArray :: forall size element. (Integral size, Prim size, Prim element) => PrimArray size -> Fold (Int, element) (PrimMultiArray element)
primMultiArray :: PrimArray size -> Fold (Int, element) (PrimMultiArray element)
primMultiArray PrimArray size
sizeArray =
  (Product2
   (MutablePrimArray RealWorld size)
   (UnliftedArray (MutablePrimArray RealWorld element))
 -> (Int, element)
 -> IO
      (Product2
         (MutablePrimArray RealWorld size)
         (UnliftedArray (MutablePrimArray RealWorld element))))
-> IO
     (Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray (MutablePrimArray RealWorld element)))
-> (Product2
      (MutablePrimArray RealWorld size)
      (UnliftedArray (MutablePrimArray RealWorld element))
    -> IO (PrimMultiArray element))
-> Fold (Int, element) (PrimMultiArray element)
forall state input output.
(state -> input -> IO state)
-> IO state -> (state -> IO output) -> Fold input output
unsafeIO Product2
  (MutablePrimArray RealWorld size)
  (UnliftedArray (MutablePrimArray RealWorld element))
-> (Int, element)
-> IO
     (Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray (MutablePrimArray RealWorld element)))
forall (m :: * -> *) a a.
(PrimMonad m, Prim a, Prim a, Integral a) =>
Product2
  (MutablePrimArray (PrimState m) a)
  (UnliftedArray (MutablePrimArray (PrimState m) a))
-> (Int, a)
-> m (Product2
        (MutablePrimArray (PrimState m) a)
        (UnliftedArray (MutablePrimArray (PrimState m) a)))
step IO
  (Product2
     (MutablePrimArray RealWorld size)
     (UnliftedArray (MutablePrimArray RealWorld element)))
init Product2
  (MutablePrimArray RealWorld size)
  (UnliftedArray (MutablePrimArray RealWorld element))
-> IO (PrimMultiArray element)
extract
  where
    outerLength :: Int
outerLength = PrimArray size -> Int
forall a. Prim a => PrimArray a -> Int
sizeofPrimArray PrimArray size
sizeArray
    init :: IO
  (Product2
     (MutablePrimArray RealWorld size)
     (UnliftedArray (MutablePrimArray RealWorld element)))
init =
      MutablePrimArray RealWorld size
-> UnliftedArray (MutablePrimArray RealWorld element)
-> Product2
     (MutablePrimArray RealWorld size)
     (UnliftedArray (MutablePrimArray RealWorld element))
forall a b. a -> b -> Product2 a b
Product2 (MutablePrimArray RealWorld size
 -> UnliftedArray (MutablePrimArray RealWorld element)
 -> Product2
      (MutablePrimArray RealWorld size)
      (UnliftedArray (MutablePrimArray RealWorld element)))
-> IO (MutablePrimArray RealWorld size)
-> IO
     (UnliftedArray (MutablePrimArray RealWorld element)
      -> Product2
           (MutablePrimArray RealWorld size)
           (UnliftedArray (MutablePrimArray RealWorld element)))
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> IO (MutablePrimArray RealWorld size)
initIndexArray IO
  (UnliftedArray (MutablePrimArray RealWorld element)
   -> Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray (MutablePrimArray RealWorld element)))
-> IO (UnliftedArray (MutablePrimArray RealWorld element))
-> IO
     (Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray (MutablePrimArray RealWorld element)))
forall (f :: * -> *) a b. Applicative f => f (a -> b) -> f a -> f b
<*> IO (UnliftedArray (MutablePrimArray RealWorld element))
initMultiArray
      where
        initIndexArray :: IO (MutablePrimArray RealWorld size)
        initIndexArray :: IO (MutablePrimArray RealWorld size)
initIndexArray =
          PrimArray size -> IO (MutablePrimArray (PrimState IO) size)
forall (m :: * -> *) a.
PrimMonad m =>
PrimArray a -> m (MutablePrimArray (PrimState m) a)
unsafeThawPrimArray (Int -> size -> PrimArray size
forall a. Prim a => Int -> a -> PrimArray a
replicatePrimArray Int
outerLength size
0)
        initMultiArray :: IO (UnliftedArray (MutablePrimArray RealWorld element))
        initMultiArray :: IO (UnliftedArray (MutablePrimArray RealWorld element))
initMultiArray =
          Int
-> (Int -> IO (MutablePrimArray RealWorld element))
-> IO (UnliftedArray (MutablePrimArray RealWorld element))
forall a.
PrimUnlifted a =>
Int -> (Int -> IO a) -> IO (UnliftedArray a)
UA.generate Int
outerLength ((Int -> IO (MutablePrimArray RealWorld element))
 -> IO (UnliftedArray (MutablePrimArray RealWorld element)))
-> (Int -> IO (MutablePrimArray RealWorld element))
-> IO (UnliftedArray (MutablePrimArray RealWorld element))
forall a b. (a -> b) -> a -> b
$ \ Int
index -> do
            Int -> IO (MutablePrimArray (PrimState IO) element)
forall (m :: * -> *) a.
(PrimMonad m, Prim a) =>
Int -> m (MutablePrimArray (PrimState m) a)
newPrimArray (size -> Int
forall a b. (Integral a, Num b) => a -> b
fromIntegral (PrimArray size -> Int -> size
forall a. Prim a => PrimArray a -> Int -> a
indexPrimArray PrimArray size
sizeArray Int
index))
    step :: Product2
  (MutablePrimArray (PrimState m) a)
  (UnliftedArray (MutablePrimArray (PrimState m) a))
-> (Int, a)
-> m (Product2
        (MutablePrimArray (PrimState m) a)
        (UnliftedArray (MutablePrimArray (PrimState m) a)))
step (Product2 MutablePrimArray (PrimState m) a
indexArray UnliftedArray (MutablePrimArray (PrimState m) a)
multiArray) (Int
outerIndex, a
element) = do
      let innerArray :: MutablePrimArray (PrimState m) a
innerArray = UnliftedArray (MutablePrimArray (PrimState m) a)
-> Int -> MutablePrimArray (PrimState m) a
forall a. PrimUnlifted a => UnliftedArray a -> Int -> a
indexUnliftedArray UnliftedArray (MutablePrimArray (PrimState m) a)
multiArray Int
outerIndex 
      a
innerIndex <- MutablePrimArray (PrimState m) a -> Int -> m a
forall a (m :: * -> *).
(Prim a, PrimMonad m) =>
MutablePrimArray (PrimState m) a -> Int -> m a
readPrimArray MutablePrimArray (PrimState m) a
indexArray Int
outerIndex
      MutablePrimArray (PrimState m) a -> Int -> a -> m ()
forall a (m :: * -> *).
(Prim a, PrimMonad m) =>
MutablePrimArray (PrimState m) a -> Int -> a -> m ()
writePrimArray MutablePrimArray (PrimState m) a
indexArray Int
outerIndex (a -> a
forall a. Enum a => a -> a
succ a
innerIndex)
      MutablePrimArray (PrimState m) a -> Int -> a -> m ()
forall a (m :: * -> *).
(Prim a, PrimMonad m) =>
MutablePrimArray (PrimState m) a -> Int -> a -> m ()
writePrimArray MutablePrimArray (PrimState m) a
innerArray (a -> Int
forall a b. (Integral a, Num b) => a -> b
fromIntegral a
innerIndex) a
element
      Product2
  (MutablePrimArray (PrimState m) a)
  (UnliftedArray (MutablePrimArray (PrimState m) a))
-> m (Product2
        (MutablePrimArray (PrimState m) a)
        (UnliftedArray (MutablePrimArray (PrimState m) a)))
forall (m :: * -> *) a. Monad m => a -> m a
return (MutablePrimArray (PrimState m) a
-> UnliftedArray (MutablePrimArray (PrimState m) a)
-> Product2
     (MutablePrimArray (PrimState m) a)
     (UnliftedArray (MutablePrimArray (PrimState m) a))
forall a b. a -> b -> Product2 a b
Product2 MutablePrimArray (PrimState m) a
indexArray UnliftedArray (MutablePrimArray (PrimState m) a)
multiArray)
    extract :: Product2
  (MutablePrimArray RealWorld size)
  (UnliftedArray (MutablePrimArray RealWorld element))
-> IO (PrimMultiArray element)
extract (Product2 MutablePrimArray RealWorld size
_ UnliftedArray (MutablePrimArray RealWorld element)
multiArray) = do
      MutableUnliftedArray RealWorld (PrimArray element)
copied <- Int -> IO (MutableUnliftedArray (PrimState IO) (PrimArray element))
forall (m :: * -> *) a.
PrimMonad m =>
Int -> m (MutableUnliftedArray (PrimState m) a)
unsafeNewUnliftedArray Int
outerLength
      Int -> (Int -> IO ()) -> IO ()
forall (m :: * -> *) a.
Applicative m =>
Int -> (Int -> m a) -> m ()
forMFromZero_ Int
outerLength ((Int -> IO ()) -> IO ()) -> (Int -> IO ()) -> IO ()
forall a b. (a -> b) -> a -> b
$ \ Int
outerIndex -> do
        let mutableInnerArray :: MutablePrimArray RealWorld element
mutableInnerArray = UnliftedArray (MutablePrimArray RealWorld element)
-> Int -> MutablePrimArray RealWorld element
forall a. PrimUnlifted a => UnliftedArray a -> Int -> a
indexUnliftedArray UnliftedArray (MutablePrimArray RealWorld element)
multiArray Int
outerIndex
        PrimArray element
frozenInnerArray <- MutablePrimArray (PrimState IO) element -> IO (PrimArray element)
forall (m :: * -> *) a.
PrimMonad m =>
MutablePrimArray (PrimState m) a -> m (PrimArray a)
unsafeFreezePrimArray MutablePrimArray RealWorld element
MutablePrimArray (PrimState IO) element
mutableInnerArray
        MutableUnliftedArray (PrimState IO) (PrimArray element)
-> Int -> PrimArray element -> IO ()
forall (m :: * -> *) a.
(PrimMonad m, PrimUnlifted a) =>
MutableUnliftedArray (PrimState m) a -> Int -> a -> m ()
writeUnliftedArray MutableUnliftedArray RealWorld (PrimArray element)
MutableUnliftedArray (PrimState IO) (PrimArray element)
copied Int
outerIndex PrimArray element
frozenInnerArray
      UnliftedArray (PrimArray element)
result <- MutableUnliftedArray (PrimState IO) (PrimArray element)
-> IO (UnliftedArray (PrimArray element))
forall (m :: * -> *) a.
PrimMonad m =>
MutableUnliftedArray (PrimState m) a -> m (UnliftedArray a)
unsafeFreezeUnliftedArray MutableUnliftedArray RealWorld (PrimArray element)
MutableUnliftedArray (PrimState IO) (PrimArray element)
copied
      PrimMultiArray element -> IO (PrimMultiArray element)
forall (m :: * -> *) a. Monad m => a -> m a
return (PrimMultiArray element -> IO (PrimMultiArray element))
-> PrimMultiArray element -> IO (PrimMultiArray element)
forall a b. (a -> b) -> a -> b
$ UnliftedArray (PrimArray element) -> PrimMultiArray element
forall a. UnliftedArray (PrimArray a) -> PrimMultiArray a
PrimMultiArray (UnliftedArray (PrimArray element) -> PrimMultiArray element)
-> UnliftedArray (PrimArray element) -> PrimMultiArray element
forall a b. (a -> b) -> a -> b
$ UnliftedArray (PrimArray element)
result