0 
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，将 aMatrix 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.

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

4

1 回答 1

1

在HMatrix, scale :: Container c e => e -> c e -> c e执行它在标签上所说的（将ea乘以c e第一个e）。这里有一些使用示例：https ://hackage.haskell.org/package/hmatrix-0.16.1.4/docs/src/Data-Packed-Internal-Numeric.html

需要注意的是，scale x通过考虑 xa 单例列表来构造 Container 类型，通过fromList.

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

于 2015-03-02T15:52:02.307 回答