我试图在 Haskell 中获得一个可用版本的多维数组,与numpy
Python 和其他语言中的数组相当。我发现了其他关于如何为特定维度的数组编写自定义函数的问题,但我的目标更重要的是,得到Data.Vector
'sslice
函数的类似行为,这很直观并且可以完成其他语言的括号索引数组的工作。
Vector
的slice
函数有类型
V.slice :: Int -> Int -> V.Vector a -> V.Vector a
所以切片向量v
很简单
import Data.Vector as V
let v = V.fromList [1..10]
i = 1
j = 5
V.slice i j v
Repa
的slice
另一方面有类型
R.slice :: (R.Slice sl, R.Shape (R.FullShape sl), R.Source r e) => R.Array r (R.FullShape sl) e -> sl -> R.Array R.D (R.SliceShape sl) e
所以它需要一个 Repa 数组和一个形状并返回一个延迟数组。我知道 Repa 不将整数作为索引,但我正在寻找slice
任意维度的函数的一般用途,无论是使用Repa
's(Z :. i)
还是ixn
维度规范。
我不是在使用 寻找与维度相关的函数traverse
,并且不希望进入任何模板 haskell 来概括它,尽管如果slice
通常无法使用该函数,它就是它。
那么问题是:是否可以使用Repa
'slice
函数来获取多维数组的任意切片,例如在numpy
'sv[x1:x2,y1:y2] 或 C++Eigen
中matrix.block<p,q>(i,j)
?