如何减少数组和返回索引?
在GNU Octave中这个代码 -
[e, ix] = min(X);
将返回最小元素和它的位置。 你如何在修复任意二进制函数?
这就是我想到的:
min x = z $ foldl' f (e,0,0) es
where
(e:es) = toList x
f (a,ix,r) b = let ix' = ix+1 in if a < b then (a,ix',r) else (b,ix',ix')
z (a,ix,r) = (a,r)
在上面的例子中,我们将repa 1D矩阵转换为列表并使用带有两个累加器的foldl'(来自Data.List) - 一个用于计算迭代(ix),另一个用于保存最小元素(r)的位置。 但使用repa的全部重点是使用数组,而不是列表!
在修复中,数组类型(foldS和foldP)有两个折叠 - 但它们只能使用类型(a - > a - > a)的函数 - 这意味着我无法将元组传递给它。 还有遍历,原则上可以将一维数组减少为一个标量数组:
min x = traverse x to0D min
where
to0D (Z:.i) = Z
min f (Z) = ??? -- how to get elements for comparison?
首先想到的是
[f (Z:.i) | i <- [1..n]], where n = ((Z:.i) -> i) $ extent x
但是这也会将数组转换为列表,而不是对数组进行计算。
我不是Repa的专家,但这似乎适用于一维阵列。 它可能可以适用于其他维度。
import Data.Array.Repa
indexed arr = traverse arr id (src idx@(Z :. i) -> (src idx, i))
minimize arr = foldP f h t
where
(Z :. n) = extent arr
arr' = indexed arr
h = arr' ! (Z :. 0)
t = extract (Z :. 1) (Z :. (n-1)) arr'
f min@(valMin, iMin) x@(val, i) | val < valMin = x
| otherwise = min
冒着复活僵尸帖子的风险,这里是我制定的任意维度解决方案(其中,评论这些评论,正是@哈马尔所建议的)。 由于怪异性质的foldAllP
,合并元组的时候,我作为一个单位元使用,也需要为最小你正在寻找找到阵列提供一个上限。
import Data.Array.Repa
import Data.Array.Repa.Eval
import Data.Array.Repa.Repr.Unboxed
minimize :: (Shape sh, Source r a, Elt a, Unbox a, Monad m, Ord a) => Array r sh a -> a -> m (a,sh)
minimize arr upperBound = do
-- Zip the array with its (Int representation of) indices
let zippedWithIndex = Data.Array.Repa.traverse arr id (f idx -> (f idx, toIndex (extent arr) idx))
-- Define a function that compares two tuple of (a,Int) on the value of the first entry
let minimumIndex t@(v,_) t'@(v',_) = if v < v' then t else t'
-- Do a parallel fold of all elements in the array
(min,idx) <- foldAllP minimumIndex (upperBound, error "Empty array") zippedWithIndex
-- convert the indice back from its Int representation
return (min, fromIndex (extent arr) idx)
当然,如果你的数组包含一个具有Bounded
实例的类型,你可以做出更方便的函数
minimize' arr = minimize arr maxBound
您可以在提示符下执行一些测试:
λ> let x = fromListUnboxed (ix2 2 2) [20, 30, 10, 40] :: Array U DIM2 Int
λ> minimize' x
(10,(Z :. 1) :. 0)
链接地址: http://www.djcxy.com/p/11481.html