many parallel applications of a sequential transform in repa
Here’s a possible modified version of your code, using delayed arrays: {-# LANGUAGE TypeOperators #-} import Data.Array.Repa as R mulM :: (Num a, Source r a) => Array r DIM2 a -> Array D DIM2 a mulM arr = traverse arr id mulM’ where mulM’ _ idx@(i’ :. i) = sumAllS $ extract (Z:.0) (Z:.(i+1)) … Read more