如何使用hmatrix构建零矩阵？

尝试使用hmatrix来创建零芯片。出于某种原因，当我尝试这个命令行，它的工作原理：如何使用hmatrix构建零矩阵？

buildMatrix 2 3 (\(r,c) -> fromIntegral 0)

然而，当我试图做同样的事情在我的代码：

type Dim = (Int, Int) 

buildFull :: Matrix Double -> Vector Int -> Vector Int -> Dim -> Int 
buildFull matrix basic nonbasic (m, n) = do 
    -- Build mxn matrix of zeroes 
    let f = buildMatrix m n (\(r,c) -> fromIntegral 0) 
    m

失败：

Pivot.hs:23:17: 
    Ambiguous type variable `a0' in the constraints: 
     (Element a0) arising from a use of `buildMatrix' 
        at Pivot.hs:23:17-27 
     (Num a0) arising from a use of `fromIntegral' at Pivot.hs:23:44-55 
    Probable fix: add a type signature that fixes these type variable(s) 
    In the expression: buildMatrix m n (\ (r, c) -> fromIntegral 0) 
    In an equation for `f': 
     f = buildMatrix m n (\ (r, c) -> fromIntegral 0) 
    In the expression: 
     do { let f = buildMatrix m n (\ (r, c) -> ...); 
      m } 
Failed, modules loaded: none.

来源

2011-12-08 drozzy

type Dim = (Int, Int) 

buildFull :: Matrix Double -> Vector Int -> Vector Int -> Dim -> Int 
buildFull matrix basic nonbasic (m, n) = do 
    -- Build mxn matrix of zeroes 
    let f = buildMatrix m n (\(r,c) -> fromIntegral 0) 
    m

首先，使用do -notation，你需要一个单子返回类型，这样就不会编译甚至固定的暧昧元素类型（如我被@Carl提醒后，它会好起来的，而在这里只有一个表达式，因此不需要(>>=)或(>>)）。

关于元素类型，在let-binding中，没有办法找出使用哪种类型，fromIntegral是否应返回Double,Integer或其他什么。通常使用的类型可以通过上下文来推断，也可以通过它所使用的表达式来推断。这里没有使用f，所以没有上下文。因此，在这种情况下，你必须指定一个签名的类型，可以

let f :: Matrix Double 
    f = buildMatrix m n (const 0)

或

let f = buildMatrix m n (\_ -> (0 :: Double))

，如果你想要的元素O型是 Double

。

来源

2011-12-08 17:20:50

这并不完全正确，你需要一个单子结果类型使用做记号。事实上，除了在desugaring过程之后它们是很好的类型之外，符号不会对它所产生的表达式的形式引入任何约束。由于这只会解除“让...在......”，并且不以任何方式涉及>> =或>>运算符，所以这种情况不需要单点返回类型。 – Carl

嗯..我得到'非法签名模式：矩阵 - >双f 使用-XScopedTypeVariables来允许它错误，如果我做第一种方法。 – drozzy

另外'let f = buildMatrix mn 0 :: Double'似乎不工作：“无法匹配预期的类型'（Int，Int） - > a0' 与实际类型'Double' 在第三个参数'buildMatrix'，即'（0 :: Double）' 在表达式中：buildMatrix mn（0 :: Double）在'f'的等式中：f = buildMatrix mn（0 :: Double）失败，加载模块：没有。” – drozzy

将fromIntegral 0替换为0::Double。否则，你想建立的矩阵是不受约束的。在提示时，扩展的违约规则可能会为你解决这个问题。

来源

2011-12-08 17:17:35 sclv

您还可以使用konst从Numeric.Container：

import Numeric.LinearAlgebra 

m = konst 0 (2,3) :: Matrix Double 

v = konst 7 10 :: Vector (Complex Float)

来源

2011-12-08 18:12:45

如何使用hmatrix构建零矩阵？

回答

相关问题