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 (Unlifted element)
 -> (Int, element)
 -> IO (MutableUnliftedArray_ RealWorld element (Unlifted element)))
-> IO (MutableUnliftedArray_ RealWorld element (Unlifted element))
-> (MutableUnliftedArray_ RealWorld element (Unlifted 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 (Unlifted element)
-> (Int, element)
-> IO (MutableUnliftedArray_ RealWorld element (Unlifted element))
forall (f :: * -> *) a.
(PrimMonad f, PrimUnlifted a) =>
MutableUnliftedArray_ (PrimState f) a (Unlifted a)
-> (Int, a)
-> f (MutableUnliftedArray_ (PrimState f) a (Unlifted a))
step IO (MutableUnliftedArray_ RealWorld element (Unlifted element))
IO (MutableUnliftedArray (PrimState IO) element)
init MutableUnliftedArray_ RealWorld element (Unlifted element)
-> IO (UnliftedArray element)
forall (m :: * -> *) a.
PrimMonad m =>
MutableUnliftedArray_ (PrimState m) a (Unlifted a)
-> m (UnliftedArray_ a (Unlifted a))
extract
  where
    step :: MutableUnliftedArray_ (PrimState f) a (Unlifted a)
-> (Int, a)
-> f (MutableUnliftedArray_ (PrimState f) a (Unlifted a))
step MutableUnliftedArray_ (PrimState f) a (Unlifted a)
mutable (Int
index, a
element) =
      MutableUnliftedArray_ (PrimState f) a (Unlifted a)
-> Int -> a -> f ()
forall (m :: * -> *) a.
(PrimMonad m, PrimUnlifted a) =>
MutableUnliftedArray (PrimState m) a -> Int -> a -> m ()
writeUnliftedArray MutableUnliftedArray_ (PrimState f) a (Unlifted a)
mutable Int
index a
element f ()
-> MutableUnliftedArray_ (PrimState f) a (Unlifted a)
-> f (MutableUnliftedArray_ (PrimState f) a (Unlifted a))
forall (f :: * -> *) a b. Functor f => f a -> b -> f b
$> MutableUnliftedArray_ (PrimState f) a (Unlifted 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 (Unlifted a)
-> m (UnliftedArray_ a (Unlifted a))
extract =
      MutableUnliftedArray_ (PrimState m) a (Unlifted a)
-> m (UnliftedArray_ a (Unlifted a))
forall (m :: * -> *) a.
PrimMonad m =>
MutableUnliftedArray_ (PrimState m) a (Unlifted a)
-> m (UnliftedArray_ a (Unlifted 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) (MutableByteArray# RealWorld))
 -> (Int, element)
 -> IO
      (Product2
         (MutablePrimArray RealWorld size)
         (UnliftedArray_
            (MutablePrimArray RealWorld element)
            (MutableByteArray# RealWorld))))
-> IO
     (Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray_
           (MutablePrimArray RealWorld element)
           (MutableByteArray# RealWorld)))
-> (Product2
      (MutablePrimArray RealWorld size)
      (UnliftedArray_
         (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld))
    -> 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) (MutableByteArray# RealWorld))
-> (Int, element)
-> IO
     (Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray_
           (MutablePrimArray RealWorld element)
           (MutableByteArray# RealWorld)))
forall (m :: * -> *) a a.
(PrimMonad m, Prim a, Prim a, Integral a) =>
Product2
  (MutablePrimArray (PrimState m) a)
  (UnliftedArray_
     (MutablePrimArray (PrimState m) a)
     (MutableByteArray# (PrimState m)))
-> (Int, a)
-> m (Product2
        (MutablePrimArray (PrimState m) a)
        (UnliftedArray_
           (MutablePrimArray (PrimState m) a)
           (MutableByteArray# (PrimState m))))
step IO
  (Product2
     (MutablePrimArray RealWorld size)
     (UnliftedArray_
        (MutablePrimArray RealWorld element)
        (MutableByteArray# RealWorld)))
init Product2
  (MutablePrimArray RealWorld size)
  (UnliftedArray_
     (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld))
-> 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)
        (MutableByteArray# RealWorld)))
init =
      MutablePrimArray RealWorld size
-> UnliftedArray_
     (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld)
-> Product2
     (MutablePrimArray RealWorld size)
     (UnliftedArray_
        (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld))
forall a b. a -> b -> Product2 a b
Product2 (MutablePrimArray RealWorld size
 -> UnliftedArray_
      (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld)
 -> Product2
      (MutablePrimArray RealWorld size)
      (UnliftedArray_
         (MutablePrimArray RealWorld element)
         (MutableByteArray# RealWorld)))
-> IO (MutablePrimArray RealWorld size)
-> IO
     (UnliftedArray_
        (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld)
      -> Product2
           (MutablePrimArray RealWorld size)
           (UnliftedArray_
              (MutablePrimArray RealWorld element)
              (MutableByteArray# RealWorld)))
forall (f :: * -> *) a b. Functor f => (a -> b) -> f a -> f b
<$> IO (MutablePrimArray RealWorld size)
initIndexArray IO
  (UnliftedArray_
     (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld)
   -> Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray_
           (MutablePrimArray RealWorld element)
           (MutableByteArray# RealWorld)))
-> IO
     (UnliftedArray_
        (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld))
-> IO
     (Product2
        (MutablePrimArray RealWorld size)
        (UnliftedArray_
           (MutablePrimArray RealWorld element)
           (MutableByteArray# RealWorld)))
forall (f :: * -> *) a b. Applicative f => f (a -> b) -> f a -> f b
<*> IO
  (UnliftedArray_
     (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld))
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)
     (MutableByteArray# (PrimState m)))
-> (Int, a)
-> m (Product2
        (MutablePrimArray (PrimState m) a)
        (UnliftedArray_
           (MutablePrimArray (PrimState m) a)
           (MutableByteArray# (PrimState m))))
step (Product2 MutablePrimArray (PrimState m) a
indexArray UnliftedArray_
  (MutablePrimArray (PrimState m) a)
  (MutableByteArray# (PrimState m))
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)
  (MutableByteArray# (PrimState m))
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)
     (MutableByteArray# (PrimState m)))
-> m (Product2
        (MutablePrimArray (PrimState m) a)
        (UnliftedArray_
           (MutablePrimArray (PrimState m) a)
           (MutableByteArray# (PrimState m))))
forall (m :: * -> *) a. Monad m => a -> m a
return (MutablePrimArray (PrimState m) a
-> UnliftedArray_
     (MutablePrimArray (PrimState m) a)
     (MutableByteArray# (PrimState m))
-> Product2
     (MutablePrimArray (PrimState m) a)
     (UnliftedArray_
        (MutablePrimArray (PrimState m) a)
        (MutableByteArray# (PrimState m)))
forall a b. a -> b -> Product2 a b
Product2 MutablePrimArray (PrimState m) a
indexArray UnliftedArray_
  (MutablePrimArray (PrimState m) a)
  (MutableByteArray# (PrimState m))
multiArray)
    extract :: Product2
  (MutablePrimArray RealWorld size)
  (UnliftedArray_
     (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld))
-> IO (PrimMultiArray element)
extract (Product2 MutablePrimArray RealWorld size
_ UnliftedArray_
  (MutablePrimArray RealWorld element) (MutableByteArray# RealWorld)
multiArray) = do
      MutableUnliftedArray_ RealWorld (PrimArray element) ByteArray#
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) (MutableByteArray# RealWorld)
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) ByteArray#
MutableUnliftedArray (PrimState IO) (PrimArray element)
copied Int
outerIndex PrimArray element
frozenInnerArray
      UnliftedArray_ (PrimArray element) ByteArray#
result <- MutableUnliftedArray (PrimState IO) (PrimArray element)
-> IO (UnliftedArray (PrimArray element))
forall (m :: * -> *) a.
PrimMonad m =>
MutableUnliftedArray_ (PrimState m) a (Unlifted a)
-> m (UnliftedArray_ a (Unlifted a))
unsafeFreezeUnliftedArray MutableUnliftedArray_ RealWorld (PrimArray element) ByteArray#
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) ByteArray#
UnliftedArray (PrimArray element)
result