如何在Edward Kmett" Linear"中使用可变大小的向量?图书馆?

时间:2017-02-15 00:53:04

标签: haskell matrix vector linear-algebra ghc

我试图使用ekmett的线性库,而且我在Linear.V中遇到了可变长度向量的问题。如何使用dim函数来获取向量的大小?如何在由嵌套trace s组成的大方阵上使用V?我在这两种情况下都会遇到错误。

最小代码:

import qualified Data.Vector as Vector
import Linear.V (V(V), dim)
import Linear.Vector (outer)
import Linear.Matrix (trace)

v, w :: V n Double -- What do I do here?
v = V $ Vector.fromList [1..5]
w = V $ Vector.fromList [2, 3, 5, 7, 11]

d = dim v
m = outer v w
t = trace m

它给出了我不理解的这些错误:

• Ambiguous type variable ‘n0’ arising from a use of ‘dim’
  prevents the constraint ‘(Linear.V.Dim n0)’ from being solved.
  Probable fix: use a type annotation to specify what ‘n0’ should be.
  These potential instances exist:
    two instances involving out-of-scope types
    (use -fprint-potential-instances to see them all)
• In the expression: dim v
  In an equation for ‘d’: d = dim v

• Ambiguous type variable ‘n1’ arising from a use of ‘trace’
  prevents the constraint ‘(Linear.V.Dim n1)’ from being solved.
  Probable fix: use a type annotation to specify what ‘n1’ should be.
  These potential instances exist:
    two instances involving out-of-scope types
    (use -fprint-potential-instances to see them all)
• In the expression: trace m
  In an equation for ‘t’: t = trace m

1 个答案:

答案 0 :(得分:5)

因为Haskell不依赖于类型,所以它不能将类型级别提升到它可能仅在运行时获得的列表的长度。话虽如此,n的目的是你可以创建超出向量大小的多态代码(例如,你可以确保你没有采用具有不同长度的向量的点积) 。 但是如果要使用该信息,还需要在编译时明确指定实际向量的长度。

linear 为您提供的是fromVector,它在运行时检查您提供的向量是否与您指定的类型匹配。例如,

ghci> :set +t -XDataKinds -XOverloadedLists
ghci> import Linear
ghci> import Linear.V
ghci> fromVector [1,2,3] :: Maybe (V 3 Int)
Just (V {toVector = [1,2,3]})
it :: Maybe (V 3 Int)
ghci> fromVector [1,2,3] :: Maybe (V 2 Int)
Nothing
it :: Maybe (V 3 Int)

因此,在您的情况下,您可能应该执行以下操作:

ghci> Just v = fromVector [1..5]           :: Maybe (V 5 Double)
v :: V 5 Double
ghci> Just w = fromVector [2, 3, 5, 7, 11] :: Maybe (V 5 Double)
w :: V 5 Double
ghci> dim v
5
it :: Int
ghci> m = outer v w
m :: V 5 (V 5 Double)
ghci> trace m
<interactive>:44:1: error:
   • No instance for (Trace (V 5)) arising from a use of ‘trace’
   • In the expression: trace m
     In an equation for ‘it’: it = trace m

... annnnd是的 - 我认为最后的互动是一个错误(除非有人能看到我错过的东西)。 Trace (V 5)变体应该可以通过Dim n => Trace (V n)实例来满足,但由于某种原因它不是。

修改

正如@ user2407038所指出的那样,问题在于我上面提到的Dim n => Trace (V n)不是多边形的 - 它只适用于n :: *,而我们希望它适用于任何类型(特别是{{} 1}}在这种情况下)。这种限制没有任何理由,因此我们可以继续定义我们自己的实例版本

n :: Nat

我打开了issue

编辑2

现在问题已解决。我认为它应该可以进入下一版ghci> :set -XPolyKinds ghci> instance Dim n => Trace (V n) ghci> trace m 106.0

作为附注,我使用linear以便我可以编写类型级文字(类型-XDataKinds - 它们是特殊的并且硬连接到GHC中)和GHC.TypeLits.Nat所以我可以写-XOverloadedLists