HMatrix矩阵双双积

Question

f1 :: Mesh -> Matrix Double
f1 me = knx where
  hx :: Double
  (hx , _) = h me
  a, knx :: Matrix Double
  a = fromLists [[1,2], [3,4]] 
  knx = hx * a 
  -- knx = 2 * a

我不明白为什么在上面的函数中，变2工作，而乘以hx = 0.5则不工作。OTOH，将一个Matrix Double与一个函数外部的Double相乘，效果是应该的。

Couldn't match expected type ‘Matrix Double’
            with actual type ‘Double’
In the first argument of ‘(*)’, namely ‘hx’
In the expression: hx * a
Failed, modules loaded: none.

我真的很困惑。任何指点都欢迎！

Answer 1

在HMatrix中，scale :: Container c e => e -> c e -> c e做它在标签上说的事情(将c e中的e乘以第一个e)。下面是一些使用示例：https://hackage.haskell.org/package/hmatrix-0.16.1.4/docs/src/Data-Packed-Internal-Numeric.html

应该注意的是，scale x通过fromList考虑x为单例列表来构造容器类型。

如果至少将常见的算术运算重载，这样公式可能类似于它们的数学对应，这将是非常方便的。我不确定定义函数同义词(例如(.*) = scale )是个好主意，还是只会增加一层复杂性。有什么想法吗？