zipWith package:repa

zipWith :: (Shape sh, Source r1 a, Source r2 b) => (a -> b -> c) -> Array r1 sh a -> Array r2 sh b -> Array D sh c
repa Data.Array.Repa Data.Array.Repa.Operators.Mapping
Combine two arrays, element-wise, with a binary operator. If the extent of the two array arguments differ, then the resulting array's extent is their intersection.
szipWith :: (Structured r1 a b, Shape sh, Source r c) => (c -> a -> b) -> Array r sh c -> Array r1 sh a -> Array (TR r1) sh b
repa Data.Array.Repa Data.Array.Repa.Operators.Mapping
Structured zipWith. If you have a cursored or partitioned source array then use that as the third argument (corresponding to r1 here)