Repaで[i] = f（a [i-1]）をどのように計算しますか？

Question

Repa？の過去の値（すなわち、より小さなインデックス）に依存する配列を計算することは可能ですか？アレイの最初の部分（例えば、a[0]）が与えられる。 （配列の要素を示すにはCのような表記法を使用していますが、混同しないでください）

私はtutorialを読んですぐにハックをチェックしましたが、それを行う機能が見つかりませんでした。

（私はあなたがそれを並列化できないため、1次元配列で計算のこの種をやってするのrepAでローミングサービスをしないと思います。しかし、私はあなたが2次元以上の場合には、それを並列化することができると思います。）

EDIT： おそらく私はどのような種類のfを使用したいのかについて具体的に述べるべきです。 a[i]がスカラーの場合に並列化する方法がないので、a[i]はN dimベクトルです。 a[i]はベクトルに「展開」することができるため、より高い次元（行列など）である必要はありません。したがって、fは、R^NをR^Nにマッピングする関数です。

例ほとんどが、それはこのようなものだ：bはN DIMベクトルで

b = M a[i-1] 
a[i][j] = g(b)[j] 

、MはN Nの行列（スパースなし仮定）であり、そしてgは、いくつかの非線形関数です。そして、i=1,..N-1の値をa[0],g、Mと計算します。私の希望は、（1）このタイプの計算を並列化し、（2）bなどの中間変数を効率的に割り当てる（C言語のような言語で、再利用することができる、Repaまたは同様のライブラリは純粋さを壊さずに魔法のように行うことができます）。

Answer 1

編集：実際には、私は質問を誤解したと思う。マルチためにそれを拡張

R.computeUnboxedS $       -- force the vector to be "real" 
    R.traverse x         -- traverse the vector 
    (\ (Z :. i) -> (Z :. (i - 1)))     -- function to get the shape of the result 
    (\f (Z :. i) -> f (Z :. (i + 1)) - f (Z :. i)) -- actual "stencil" 

：

Prelude Data.Array.Repa R> let x = fromListUnboxed (Z :. 10 :: DIM1) [1..10] 
Prelude Data.Array.Repa R> R.computeUnboxedS $ R.traverse x (\ (Z :. i) -> (Z :. (i - 1))) (\f (Z :. i) -> f (Z :. (i + 1)) - f (Z :. i)) 
AUnboxed (Z :. 9) (fromList [1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0,1.0]) 

それを解剖：私は

が

あなたはtraversehttp://hackage.haskell.org/packages/archive/repa/3.2.1.1/doc/html/Data-Array-Repa.html#v:traverseを使用することができます...それは他の誰かのために便利だ場合には、ここに私の答えを残しておきます2次元アレイは些細でなければならない。

Answer 2

私はこれを行うためのRepaの方法を見ることができません。しかし、ベクターがある：Data.Vector.iterateNは、あなたが望むベクトルを構築する。次にData.Array.Repa.fromUnboxedをVectorからRepaに変換します。

iterateN :: Int -> (a -> a) -> a -> Vector aSource 

O（n）の値に関数をn回適用します。ゼロ要素は元の値です。