如何在Haskell的Hmatrix中设置切片？

Question

Haskell的hmatrix允许您方便地 get slices：

m ?? (All, Take 3)

但是如何 set slices ，特别是非矩形的？例如，在Python的Numpy中，我会做：

>>> x = np.arange(12).reshape(3, 4)
>>> x
array([[ 0,  1,  2,  3],
       [ 4,  5,  6,  7],
       [ 8,  9, 10, 11]])
>>> x[[0,0], [1,3]] += x[[2,2], [1,3]]     # to ix (0, 1) add (2, 1)  and  to ix (0, 3) add (2, 3)
>>> x
array([[ 0, 10,  2, 14],
       [ 4,  5,  6,  7],
       [ 8,  9, 10, 11]])

我可以自己编写此辅助功能，但是，这是丑陋的，部分，并且需要我添加zipwith3m _，该_ 应该存在于向量库，但没有：

setSlice :: (Element t, Num t) => Matrix t -> (Extractor, Extractor) -> Vector t -> Matrix t
setSlice parent (Pos (V.map fromIntegral -> irows), Pos (V.map fromIntegral -> icols)) sub = runST $ do
  tparent <- thawMatrix parent
  zipWith3M_ (writeMatrix tparent) irows icols sub
  freezeMatrix tparent