{-# LANGUAGE BangPatterns #-}
{-# LANGUAGE CPP #-}
{-# LANGUAGE FlexibleContexts #-}
{-# LANGUAGE FlexibleInstances #-}
{-# LANGUAGE MultiParamTypeClasses #-}
{-# LANGUAGE ScopedTypeVariables #-}
{-# LANGUAGE TypeFamilies #-}
{-# LANGUAGE UndecidableInstances #-}
-- |
-- Module      : Data.Massiv.Array.Delayed.Pull
-- Copyright   : (c) Alexey Kuleshevich 2018-2021
-- License     : BSD3
-- Maintainer  : Alexey Kuleshevich <lehins@yandex.ru>
-- Stability   : experimental
-- Portability : non-portable
--
module Data.Massiv.Array.Delayed.Pull
  ( D(..)
  , Array(..)
  , delay
  , eqArrays
  , compareArrays
  , imap
  , liftArray2'
  , liftArray2M
  , unsafeExtract
  , unsafeSlice
  , unsafeInnerSlice
  ) where

import           Control.Applicative
import qualified Data.Foldable as F
import           Data.Massiv.Array.Ops.Fold.Internal as A
import           Data.Massiv.Core.Common as A
import           Data.Massiv.Core.List (L, showArrayList, showsArrayPrec)
import           Data.Massiv.Core.Operations
import           Data.Massiv.Vector.Stream as S (steps)
import           GHC.Base (build)
import           Prelude hiding (zipWith)

#include "massiv.h"


-- | Delayed representation.
data D = D deriving Int -> D -> ShowS
[D] -> ShowS
D -> String
(Int -> D -> ShowS) -> (D -> String) -> ([D] -> ShowS) -> Show D
forall a.
(Int -> a -> ShowS) -> (a -> String) -> ([a] -> ShowS) -> Show a
showList :: [D] -> ShowS
$cshowList :: [D] -> ShowS
show :: D -> String
$cshow :: D -> String
showsPrec :: Int -> D -> ShowS
$cshowsPrec :: Int -> D -> ShowS
Show


data instance Array D ix e = DArray { Array D ix e -> Comp
dComp :: !Comp
                                    , Array D ix e -> Sz ix
dSize :: !(Sz ix)
                                    , Array D ix e -> ix -> e
dIndex :: ix -> e }

instance (Ragged L ix e, Show e) => Show (Array D ix e) where
  showsPrec :: Int -> Array D ix e -> ShowS
showsPrec = (Array D ix e -> Array D ix e) -> Int -> Array D ix e -> ShowS
forall r r' ix e.
(Ragged L ix e, Load r ix e, Load r' ix e, Source r' e, Show e) =>
(Array r ix e -> Array r' ix e) -> Int -> Array r ix e -> ShowS
showsArrayPrec Array D ix e -> Array D ix e
forall a. a -> a
id
  showList :: [Array D ix e] -> ShowS
showList = [Array D ix e] -> ShowS
forall arr. Show arr => [arr] -> ShowS
showArrayList

instance Index ix => Shape D ix where
  maxLinearSize :: Array D ix e -> Maybe Sz1
maxLinearSize = Sz1 -> Maybe Sz1
forall a. a -> Maybe a
Just (Sz1 -> Maybe Sz1)
-> (Array D ix e -> Sz1) -> Array D ix e -> Maybe Sz1
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Int -> Sz1
forall ix. ix -> Sz ix
SafeSz (Int -> Sz1) -> (Array D ix e -> Int) -> Array D ix e -> Sz1
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Array D ix e -> Int
forall ix r e. (Index ix, Size r) => Array r ix e -> Int
elemsCount
  {-# INLINE maxLinearSize #-}

instance Size D where
  size :: Array D ix e -> Sz ix
size = Array D ix e -> Sz ix
forall ix e. Array D ix e -> Sz ix
dSize
  {-# INLINE size #-}
  unsafeResize :: Sz ix' -> Array D ix e -> Array D ix' e
unsafeResize !Sz ix'
sz !Array D ix e
arr =
    Comp -> Sz ix' -> (ix' -> e) -> Array D ix' e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array D ix e -> Comp
forall ix e. Array D ix e -> Comp
dComp Array D ix e
arr) Sz ix'
sz ((ix' -> e) -> Array D ix' e) -> (ix' -> e) -> Array D ix' e
forall a b. (a -> b) -> a -> b
$ \ !ix'
ix ->
      Array D ix e -> ix -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array D ix e
arr (Sz ix -> Int -> ix
forall ix. Index ix => Sz ix -> Int -> ix
fromLinearIndex (Array D ix e -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array D ix e
arr) (Sz ix' -> ix' -> Int
forall ix. Index ix => Sz ix -> ix -> Int
toLinearIndex Sz ix'
sz ix'
ix))
  {-# INLINE unsafeResize #-}

instance Strategy D where
  setComp :: Comp -> Array D ix e -> Array D ix e
setComp Comp
c Array D ix e
arr = Array D ix e
R:ArrayDixe ix e
arr { dComp :: Comp
dComp = Comp
c }
  {-# INLINE setComp #-}
  getComp :: Array D ix e -> Comp
getComp = Array D ix e -> Comp
forall ix e. Array D ix e -> Comp
dComp
  {-# INLINE getComp #-}

instance Source D e where
  unsafeIndex :: Array D ix e -> ix -> e
unsafeIndex = INDEX_CHECK("(Source D ix e).unsafeIndex", size, dIndex)
  {-# INLINE unsafeIndex #-}

  unsafeOuterSlice :: Array D ix e -> Sz (Lower ix) -> Int -> Array D (Lower ix) e
unsafeOuterSlice !Array D ix e
arr !Sz (Lower ix)
szL !Int
i = Comp -> Sz (Lower ix) -> (Lower ix -> e) -> Array D (Lower ix) e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array D ix e -> Comp
forall ix e. Array D ix e -> Comp
dComp Array D ix e
arr) Sz (Lower ix)
szL (Array D ix e -> ix -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array D ix e
arr (ix -> e) -> (Lower ix -> ix) -> Lower ix -> e
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Int -> Lower ix -> ix
forall ix. Index ix => Int -> Lower ix -> ix
consDim Int
i)
  {-# INLINE unsafeOuterSlice #-}

  unsafeLinearSlice :: Int -> Sz1 -> Array D ix e -> Array D Int e
unsafeLinearSlice !Int
o !Sz1
sz Array D ix e
arr =
    Comp -> Sz1 -> (Int -> e) -> Array D Int e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array D ix e -> Comp
forall ix e. Array D ix e -> Comp
dComp Array D ix e
arr) Sz1
sz ((Int -> e) -> Array D Int e) -> (Int -> e) -> Array D Int e
forall a b. (a -> b) -> a -> b
$ \ !Int
i -> Array D ix e -> ix -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array D ix e
arr (Sz ix -> Int -> ix
forall ix. Index ix => Sz ix -> Int -> ix
fromLinearIndex (Array D ix e -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array D ix e
arr) (Int
i Int -> Int -> Int
forall a. Num a => a -> a -> a
+ Int
o))
  {-# INLINE unsafeLinearSlice #-}


-- | /O(1)/ - Extract a portion of an array. Staring index and new size are
-- not validated.
unsafeExtract :: (Source r e, Index ix) => ix -> Sz ix -> Array r ix e -> Array D ix e
unsafeExtract :: ix -> Sz ix -> Array r ix e -> Array D ix e
unsafeExtract !ix
sIx !Sz ix
newSz !Array r ix e
arr =
  Comp -> Sz ix -> (ix -> e) -> Array D ix e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array r ix e -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r ix e
arr) Sz ix
newSz (Array r ix e -> ix -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r ix e
arr (ix -> e) -> (ix -> ix) -> ix -> e
forall b c a. (b -> c) -> (a -> b) -> a -> c
. (Int -> Int -> Int) -> ix -> ix -> ix
forall ix. Index ix => (Int -> Int -> Int) -> ix -> ix -> ix
liftIndex2 Int -> Int -> Int
forall a. Num a => a -> a -> a
(+) ix
sIx)
{-# INLINE unsafeExtract #-}

-- | /O(1)/ - Take a slice out of an array from within
unsafeSlice :: (Source r e, Index ix, Index (Lower ix), MonadThrow m) =>
  Array r ix e -> ix -> Sz ix -> Dim -> m (Array D (Lower ix) e)
unsafeSlice :: Array r ix e -> ix -> Sz ix -> Dim -> m (Array D (Lower ix) e)
unsafeSlice Array r ix e
arr ix
start cut :: Sz ix
cut@(SafeSz ix
cutSz) Dim
dim = do
  Lower ix
newSz <- ix -> Dim -> m (Lower ix)
forall (m :: * -> *) ix.
(MonadThrow m, Index ix) =>
ix -> Dim -> m (Lower ix)
dropDimM ix
cutSz Dim
dim
  Array D (Lower ix) e -> m (Array D (Lower ix) e)
forall (m :: * -> *) a. Monad m => a -> m a
return (Array D (Lower ix) e -> m (Array D (Lower ix) e))
-> Array D (Lower ix) e -> m (Array D (Lower ix) e)
forall a b. (a -> b) -> a -> b
$ Sz (Lower ix) -> Array D ix e -> Array D (Lower ix) e
forall r ix ix' e.
(Size r, Index ix, Index ix') =>
Sz ix' -> Array r ix e -> Array r ix' e
unsafeResize (Lower ix -> Sz (Lower ix)
forall ix. ix -> Sz ix
SafeSz Lower ix
newSz) (ix -> Sz ix -> Array r ix e -> Array D ix e
forall r e ix.
(Source r e, Index ix) =>
ix -> Sz ix -> Array r ix e -> Array D ix e
unsafeExtract ix
start Sz ix
cut Array r ix e
arr)
{-# INLINE unsafeSlice #-}

-- | /O(1)/ - Take a slice out of an array from the inside
unsafeInnerSlice ::
     (Source r e, Index ix) => Array r ix e -> Sz (Lower ix) -> Int -> Array D (Lower ix) e
unsafeInnerSlice :: Array r ix e -> Sz (Lower ix) -> Int -> Array D (Lower ix) e
unsafeInnerSlice !Array r ix e
arr Sz (Lower ix)
szL !Int
i = Comp -> Sz (Lower ix) -> (Lower ix -> e) -> Array D (Lower ix) e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array r ix e -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r ix e
arr) Sz (Lower ix)
szL (Array r ix e -> ix -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r ix e
arr (ix -> e) -> (Lower ix -> ix) -> Lower ix -> e
forall b c a. (b -> c) -> (a -> b) -> a -> c
. (Lower ix -> Int -> ix
forall ix. Index ix => Lower ix -> Int -> ix
`snocDim` Int
i))
{-# INLINE unsafeInnerSlice #-}



instance (Eq e, Index ix) => Eq (Array D ix e) where
  == :: Array D ix e -> Array D ix e -> Bool
(==) = (e -> e -> Bool) -> Array D ix e -> Array D ix e -> Bool
forall ix r1 e1 r2 e2.
(Index ix, Source r1 e1, Source r2 e2) =>
(e1 -> e2 -> Bool) -> Array r1 ix e1 -> Array r2 ix e2 -> Bool
eqArrays e -> e -> Bool
forall a. Eq a => a -> a -> Bool
(==)
  {-# INLINE (==) #-}

instance (Ord e, Index ix) => Ord (Array D ix e) where
  compare :: Array D ix e -> Array D ix e -> Ordering
compare = (e -> e -> Ordering) -> Array D ix e -> Array D ix e -> Ordering
forall ix r1 e1 r2 e2.
(Index ix, Source r1 e1, Source r2 e2) =>
(e1 -> e2 -> Ordering)
-> Array r1 ix e1 -> Array r2 ix e2 -> Ordering
compareArrays e -> e -> Ordering
forall a. Ord a => a -> a -> Ordering
compare
  {-# INLINE compare #-}

instance Functor (Array D ix) where
  fmap :: (a -> b) -> Array D ix a -> Array D ix b
fmap a -> b
f (DArray c sz g) = Comp -> Sz ix -> (ix -> b) -> Array D ix b
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray Comp
c Sz ix
sz (a -> b
f (a -> b) -> (ix -> a) -> ix -> b
forall b c a. (b -> c) -> (a -> b) -> a -> c
. ix -> a
g)
  {-# INLINE fmap #-}
  <$ :: a -> Array D ix b -> Array D ix a
(<$) a
e (DArray c sz _) = Comp -> Sz ix -> (ix -> a) -> Array D ix a
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray Comp
c Sz ix
sz (a -> ix -> a
forall a b. a -> b -> a
const a
e)
  {-# INLINE (<$) #-}


instance Index ix => Applicative (Array D ix) where
  pure :: a -> Array D ix a
pure = a -> Array D ix a
forall r ix e. Load r ix e => e -> Array r ix e
singleton
  {-# INLINE pure #-}
  <*> :: Array D ix (a -> b) -> Array D ix a -> Array D ix b
(<*>) = ((a -> b) -> a -> b)
-> Array D ix (a -> b) -> Array D ix a -> Array D ix b
forall ix r1 a r2 b e.
(HasCallStack, Index ix, Source r1 a, Source r2 b) =>
(a -> b -> e) -> Array r1 ix a -> Array r2 ix b -> Array D ix e
liftArray2' (a -> b) -> a -> b
forall a. a -> a
id
  {-# INLINE (<*>) #-}
#if MIN_VERSION_base(4,10,0)
  liftA2 :: (a -> b -> c) -> Array D ix a -> Array D ix b -> Array D ix c
liftA2 = (a -> b -> c) -> Array D ix a -> Array D ix b -> Array D ix c
forall ix r1 a r2 b e.
(HasCallStack, Index ix, Source r1 a, Source r2 b) =>
(a -> b -> e) -> Array r1 ix a -> Array r2 ix b -> Array D ix e
liftArray2'
  {-# INLINE liftA2 #-}
#endif

-- | Row-major sequential folding over a Delayed array.
instance Index ix => Foldable (Array D ix) where
  fold :: Array D ix m -> m
fold = Array D ix m -> m
forall e ix r.
(Monoid e, Index ix, Source r e) =>
Array r ix e -> e
A.fold
  {-# INLINE fold #-}
  foldMap :: (a -> m) -> Array D ix a -> m
foldMap = (a -> m) -> Array D ix a -> m
forall ix r e m.
(Index ix, Source r e, Monoid m) =>
(e -> m) -> Array r ix e -> m
A.foldMono
  {-# INLINE foldMap #-}
  foldl :: (b -> a -> b) -> b -> Array D ix a -> b
foldl = (b -> a -> b) -> b -> Array D ix a -> b
forall ix r e a.
(Index ix, Source r e) =>
(a -> e -> a) -> a -> Array r ix e -> a
lazyFoldlS
  {-# INLINE foldl #-}
  foldl' :: (b -> a -> b) -> b -> Array D ix a -> b
foldl' = (b -> a -> b) -> b -> Array D ix a -> b
forall ix r e a.
(Index ix, Source r e) =>
(a -> e -> a) -> a -> Array r ix e -> a
foldlS
  {-# INLINE foldl' #-}
  foldr :: (a -> b -> b) -> b -> Array D ix a -> b
foldr = (a -> b -> b) -> b -> Array D ix a -> b
forall ix r e b.
(Index ix, Source r e) =>
(e -> b -> b) -> b -> Array r ix e -> b
foldrFB
  {-# INLINE foldr #-}
  foldr' :: (a -> b -> b) -> b -> Array D ix a -> b
foldr' = (a -> b -> b) -> b -> Array D ix a -> b
forall ix r e b.
(Index ix, Source r e) =>
(e -> b -> b) -> b -> Array r ix e -> b
foldrS
  {-# INLINE foldr' #-}
  null :: Array D ix a -> Bool
null (DArray _ sz _) = Sz ix -> Int
forall ix. Index ix => Sz ix -> Int
totalElem Sz ix
sz Int -> Int -> Bool
forall a. Eq a => a -> a -> Bool
== Int
0
  {-# INLINE null #-}
  length :: Array D ix a -> Int
length = Sz ix -> Int
forall ix. Index ix => Sz ix -> Int
totalElem (Sz ix -> Int) -> (Array D ix a -> Sz ix) -> Array D ix a -> Int
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Array D ix a -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size
  {-# INLINE length #-}
  elem :: a -> Array D ix a -> Bool
elem a
e = (a -> Bool) -> Array D ix a -> Bool
forall ix r e.
(Index ix, Source r e) =>
(e -> Bool) -> Array r ix e -> Bool
A.any (a
e a -> a -> Bool
forall a. Eq a => a -> a -> Bool
==)
  {-# INLINE elem #-}
  toList :: Array D ix a -> [a]
toList Array D ix a
arr = (forall b. (a -> b -> b) -> b -> b) -> [a]
forall a. (forall b. (a -> b -> b) -> b -> b) -> [a]
build (\ a -> b -> b
c b
n -> (a -> b -> b) -> b -> Array D ix a -> b
forall ix r e b.
(Index ix, Source r e) =>
(e -> b -> b) -> b -> Array r ix e -> b
foldrFB a -> b -> b
c b
n Array D ix a
arr)
  {-# INLINE toList #-}


instance Index ix => Load D ix e where
  makeArray :: Comp -> Sz ix -> (ix -> e) -> Array D ix e
makeArray = Comp -> Sz ix -> (ix -> e) -> Array D ix e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray
  {-# INLINE makeArray #-}
  iterArrayLinearST_ :: Scheduler s () -> Array D ix e -> (Int -> e -> ST s ()) -> ST s ()
iterArrayLinearST_ !Scheduler s ()
scheduler !Array D ix e
arr =
    Scheduler s ()
-> Int -> (Int -> e) -> (Int -> e -> ST s ()) -> ST s ()
forall s (m :: * -> *) b.
MonadPrimBase s m =>
Scheduler s () -> Int -> (Int -> b) -> (Int -> b -> m ()) -> m ()
splitLinearlyWith_ Scheduler s ()
scheduler (Array D ix e -> Int
forall ix r e. (Index ix, Size r) => Array r ix e -> Int
elemsCount Array D ix e
arr) (Array D ix e -> Int -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> Int -> e
unsafeLinearIndex Array D ix e
arr)
  {-# INLINE iterArrayLinearST_ #-}

instance Index ix => StrideLoad D ix e

instance Index ix => Stream D ix e where
  toStream :: Array D ix e -> Steps Id e
toStream = Array D ix e -> Steps Id e
forall r ix e (m :: * -> *).
(Monad m, Index ix, Source r e) =>
Array r ix e -> Steps m e
S.steps
  {-# INLINE toStream #-}
  toStreamIx :: Array D ix e -> Steps Id (ix, e)
toStreamIx = Array D ix (ix, e) -> Steps Id (ix, e)
forall r ix e (m :: * -> *).
(Monad m, Index ix, Source r e) =>
Array r ix e -> Steps m e
S.steps (Array D ix (ix, e) -> Steps Id (ix, e))
-> (Array D ix e -> Array D ix (ix, e))
-> Array D ix e
-> Steps Id (ix, e)
forall b c a. (b -> c) -> (a -> b) -> a -> c
. (ix -> e -> (ix, e)) -> Array D ix e -> Array D ix (ix, e)
forall r ix e a.
(Index ix, Source r e) =>
(ix -> e -> a) -> Array r ix e -> Array D ix a
imap (,)
  {-# INLINE toStreamIx #-}

-- | Map an index aware function over an array
--
-- @since 0.1.0
imap :: forall r ix e a. (Index ix, Source r e) => (ix -> e -> a) -> Array r ix e -> Array D ix a
imap :: (ix -> e -> a) -> Array r ix e -> Array D ix a
imap ix -> e -> a
f !Array r ix e
arr = Comp -> Sz ix -> (ix -> a) -> Array D ix a
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array r ix e -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r ix e
arr) (Array r ix e -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r ix e
arr) (\ !ix
ix -> ix -> e -> a
f ix
ix (Array r ix e -> ix -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r ix e
arr ix
ix))
{-# INLINE imap #-}


instance Num e => FoldNumeric D e where
  unsafeDotProduct :: Array D ix e -> Array D ix e -> e
unsafeDotProduct = Array D ix e -> Array D ix e -> e
forall e ix r.
(Num e, Index ix, Source r e) =>
Array r ix e -> Array r ix e -> e
defaultUnsafeDotProduct
  {-# INLINE unsafeDotProduct #-}
  powerSumArray :: Array D ix e -> Int -> e
powerSumArray = Array D ix e -> Int -> e
forall ix r e.
(Index ix, Source r e, Num e) =>
Array r ix e -> Int -> e
defaultPowerSumArray
  {-# INLINE powerSumArray #-}
  foldArray :: (e -> e -> e) -> e -> Array D ix e -> e
foldArray = (e -> e -> e) -> e -> Array D ix e -> e
forall ix r e.
(Index ix, Source r e) =>
(e -> e -> e) -> e -> Array r ix e -> e
defaultFoldArray
  {-# INLINE foldArray #-}

instance Num e => Numeric D e where
  unsafeLiftArray :: (e -> e) -> Array D ix e -> Array D ix e
unsafeLiftArray e -> e
f Array D ix e
arr = Array D ix e
R:ArrayDixe ix e
arr {dIndex :: ix -> e
dIndex = e -> e
f (e -> e) -> (ix -> e) -> ix -> e
forall b c a. (b -> c) -> (a -> b) -> a -> c
. Array D ix e -> ix -> e
forall ix e. Array D ix e -> ix -> e
dIndex Array D ix e
arr}
  {-# INLINE unsafeLiftArray #-}
  unsafeLiftArray2 :: (e -> e -> e) -> Array D ix e -> Array D ix e -> Array D ix e
unsafeLiftArray2 e -> e -> e
f Array D ix e
a1 Array D ix e
a2 = -- TODO: possibly use the first size, it is unsafe anyways.
    Comp -> Sz ix -> (ix -> e) -> Array D ix e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array D ix e -> Comp
forall ix e. Array D ix e -> Comp
dComp Array D ix e
a1 Comp -> Comp -> Comp
forall a. Semigroup a => a -> a -> a
<> Array D ix e -> Comp
forall ix e. Array D ix e -> Comp
dComp Array D ix e
a2) (ix -> Sz ix
forall ix. ix -> Sz ix
SafeSz ((Int -> Int -> Int) -> ix -> ix -> ix
forall ix. Index ix => (Int -> Int -> Int) -> ix -> ix -> ix
liftIndex2 Int -> Int -> Int
forall a. Ord a => a -> a -> a
min (Sz ix -> ix
forall ix. Sz ix -> ix
unSz (Array D ix e -> Sz ix
forall ix e. Array D ix e -> Sz ix
dSize Array D ix e
a1)) (Sz ix -> ix
forall ix. Sz ix -> ix
unSz (Array D ix e -> Sz ix
forall ix e. Array D ix e -> Sz ix
dSize Array D ix e
a2)))) ((ix -> e) -> Array D ix e) -> (ix -> e) -> Array D ix e
forall a b. (a -> b) -> a -> b
$ \ix
i ->
      e -> e -> e
f (Array D ix e -> ix -> e
forall ix e. Array D ix e -> ix -> e
dIndex Array D ix e
a1 ix
i) (Array D ix e -> ix -> e
forall ix e. Array D ix e -> ix -> e
dIndex Array D ix e
a2 ix
i)
  {-# INLINE unsafeLiftArray2 #-}


instance Floating e => NumericFloat D e



-- | /O(1)/ Conversion from a source array to `D` representation.
delay :: (Index ix, Source r e) => Array r ix e -> Array D ix e
delay :: Array r ix e -> Array D ix e
delay Array r ix e
arr = Comp -> Sz ix -> (ix -> e) -> Array D ix e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array r ix e -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r ix e
arr) (Array r ix e -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r ix e
arr) (Array r ix e -> ix -> e
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r ix e
arr)
{-# INLINE [1] delay #-}

{-# RULES
"delay" [~1] forall (arr :: Array D ix e) . delay arr = arr
 #-}

-- | Compute array equality by applying a comparing function to each
-- element. Empty arrays are always equal, regardless of their size.
--
-- @since 0.5.7
eqArrays :: (Index ix, Source r1 e1, Source r2 e2) =>
            (e1 -> e2 -> Bool) -> Array r1 ix e1 -> Array r2 ix e2 -> Bool
eqArrays :: (e1 -> e2 -> Bool) -> Array r1 ix e1 -> Array r2 ix e2 -> Bool
eqArrays e1 -> e2 -> Bool
f Array r1 ix e1
arr1 Array r2 ix e2
arr2 =
  let sz1 :: Sz ix
sz1 = Array r1 ix e1 -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r1 ix e1
arr1
      sz2 :: Sz ix
sz2 = Array r2 ix e2 -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r2 ix e2
arr2
   in (Sz ix
sz1 Sz ix -> Sz ix -> Bool
forall a. Eq a => a -> a -> Bool
== Sz ix
sz2 Bool -> Bool -> Bool
&&
       Bool -> Bool
not
         ((Bool -> Bool) -> Array D ix Bool -> Bool
forall ix r e.
(Index ix, Source r e) =>
(e -> Bool) -> Array r ix e -> Bool
A.any
            Bool -> Bool
not
            (Comp -> Sz ix -> (ix -> Bool) -> Array D ix Bool
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array r1 ix e1 -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r1 ix e1
arr1 Comp -> Comp -> Comp
forall a. Semigroup a => a -> a -> a
<> Array r2 ix e2 -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r2 ix e2
arr2) (Array r1 ix e1 -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r1 ix e1
arr1) ((ix -> Bool) -> Array D ix Bool)
-> (ix -> Bool) -> Array D ix Bool
forall a b. (a -> b) -> a -> b
$ \ix
ix ->
               e1 -> e2 -> Bool
f (Array r1 ix e1 -> ix -> e1
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r1 ix e1
arr1 ix
ix) (Array r2 ix e2 -> ix -> e2
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r2 ix e2
arr2 ix
ix)))) Bool -> Bool -> Bool
||
      (Sz ix -> Bool
forall ix. Index ix => Sz ix -> Bool
isZeroSz Sz ix
sz1 Bool -> Bool -> Bool
&& Sz ix -> Bool
forall ix. Index ix => Sz ix -> Bool
isZeroSz Sz ix
sz2)
{-# INLINE eqArrays #-}

-- | Compute array ordering by applying a comparing function to each element.
-- The exact ordering is unspecified so this is only intended for use in maps and the like where
-- you need an ordering but do not care about which one is used.
--
-- @since 0.5.7
compareArrays :: (Index ix, Source r1 e1, Source r2 e2) =>
       (e1 -> e2 -> Ordering) -> Array r1 ix e1 -> Array r2 ix e2 -> Ordering
compareArrays :: (e1 -> e2 -> Ordering)
-> Array r1 ix e1 -> Array r2 ix e2 -> Ordering
compareArrays e1 -> e2 -> Ordering
f Array r1 ix e1
arr1 Array r2 ix e2
arr2 =
  Sz ix -> Sz ix -> Ordering
forall a. Ord a => a -> a -> Ordering
compare (Array r1 ix e1 -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r1 ix e1
arr1) (Array r2 ix e2 -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r2 ix e2
arr2) Ordering -> Ordering -> Ordering
forall a. Semigroup a => a -> a -> a
<>
  Array D ix Ordering -> Ordering
forall e ix r.
(Monoid e, Index ix, Source r e) =>
Array r ix e -> e
A.fold
    (Comp -> Sz ix -> (ix -> Ordering) -> Array D ix Ordering
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray (Array r1 ix e1 -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r1 ix e1
arr1 Comp -> Comp -> Comp
forall a. Semigroup a => a -> a -> a
<> Array r2 ix e2 -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r2 ix e2
arr2) (Array r1 ix e1 -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r1 ix e1
arr1) ((ix -> Ordering) -> Array D ix Ordering)
-> (ix -> Ordering) -> Array D ix Ordering
forall a b. (a -> b) -> a -> b
$ \ix
ix ->
       e1 -> e2 -> Ordering
f (Array r1 ix e1 -> ix -> e1
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r1 ix e1
arr1 ix
ix) (Array r2 ix e2 -> ix -> e2
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r2 ix e2
arr2 ix
ix))
{-# INLINE compareArrays #-}

-- | Same as `liftArray2M`, but throws an imprecise exception on mismatched
-- sizes.
--
-- @since 1.0.0
liftArray2'
  :: (HasCallStack, Index ix, Source r1 a, Source r2 b)
  => (a -> b -> e) -> Array r1 ix a -> Array r2 ix b -> Array D ix e
liftArray2' :: (a -> b -> e) -> Array r1 ix a -> Array r2 ix b -> Array D ix e
liftArray2' a -> b -> e
f Array r1 ix a
arr1 Array r2 ix b
arr2 = Either SomeException (Array D ix e) -> Array D ix e
forall a. HasCallStack => Either SomeException a -> a
throwEither (Either SomeException (Array D ix e) -> Array D ix e)
-> Either SomeException (Array D ix e) -> Array D ix e
forall a b. (a -> b) -> a -> b
$ (a -> b -> e)
-> Array r1 ix a
-> Array r2 ix b
-> Either SomeException (Array D ix e)
forall ix r1 a r2 b (m :: * -> *) e.
(Index ix, Source r1 a, Source r2 b, MonadThrow m) =>
(a -> b -> e) -> Array r1 ix a -> Array r2 ix b -> m (Array D ix e)
liftArray2M a -> b -> e
f Array r1 ix a
arr1 Array r2 ix b
arr2
{-# INLINE liftArray2' #-}


-- | Similar to `Data.Massiv.Array.zipWith`, except dimensions of both arrays
-- have to be the same, otherwise it throws `SizeMismatchException`.
--
-- @since 1.0.0
liftArray2M
  :: (Index ix, Source r1 a, Source r2 b, MonadThrow m)
  => (a -> b -> e) -> Array r1 ix a -> Array r2 ix b -> m (Array D ix e)
liftArray2M :: (a -> b -> e) -> Array r1 ix a -> Array r2 ix b -> m (Array D ix e)
liftArray2M a -> b -> e
f !Array r1 ix a
arr1 !Array r2 ix b
arr2
  | Sz ix
sz1 Sz ix -> Sz ix -> Bool
forall a. Eq a => a -> a -> Bool
== Sz ix
sz2 =
    Array D ix e -> m (Array D ix e)
forall (f :: * -> *) a. Applicative f => a -> f a
pure (Array D ix e -> m (Array D ix e))
-> Array D ix e -> m (Array D ix e)
forall a b. (a -> b) -> a -> b
$
    Comp -> Sz ix -> (ix -> e) -> Array D ix e
forall ix e. Comp -> Sz ix -> (ix -> e) -> Array D ix e
DArray
      (Array r1 ix a -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r1 ix a
arr1 Comp -> Comp -> Comp
forall a. Semigroup a => a -> a -> a
<> Array r2 ix b -> Comp
forall r ix e. Strategy r => Array r ix e -> Comp
getComp Array r2 ix b
arr2)
      Sz ix
sz1
      (\ !ix
ix -> a -> b -> e
f (Array r1 ix a -> ix -> a
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r1 ix a
arr1 ix
ix) (Array r2 ix b -> ix -> b
forall r e ix. (Source r e, Index ix) => Array r ix e -> ix -> e
unsafeIndex Array r2 ix b
arr2 ix
ix))
  | Sz ix -> Bool
forall ix. Index ix => Sz ix -> Bool
isZeroSz Sz ix
sz1 Bool -> Bool -> Bool
&& Sz ix -> Bool
forall ix. Index ix => Sz ix -> Bool
isZeroSz Sz ix
sz2 = Array D ix e -> m (Array D ix e)
forall (f :: * -> *) a. Applicative f => a -> f a
pure Array D ix e
forall r ix e. Load r ix e => Array r ix e
A.empty
  | Bool
otherwise = SizeException -> m (Array D ix e)
forall (m :: * -> *) e a. (MonadThrow m, Exception e) => e -> m a
throwM (SizeException -> m (Array D ix e))
-> SizeException -> m (Array D ix e)
forall a b. (a -> b) -> a -> b
$ Sz ix -> Sz ix -> SizeException
forall ix. Index ix => Sz ix -> Sz ix -> SizeException
SizeMismatchException (Array r1 ix a -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r1 ix a
arr1) (Array r2 ix b -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r2 ix b
arr2)
  where
    sz1 :: Sz ix
sz1 = Array r1 ix a -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r1 ix a
arr1
    sz2 :: Sz ix
sz2 = Array r2 ix b -> Sz ix
forall r ix e. Size r => Array r ix e -> Sz ix
size Array r2 ix b
arr2
{-# INLINE liftArray2M #-}