GHC 9.6.2、hmatrix与Vector 0.13.0.0编译类型错误排查及修复

问题原因确认 #

你遇到的类型不匹配问题确实是vector 0.13.0.0的变更导致的。hmatrix中的Vector R本质是对Data.Vector.Storable.Vector Double的newtype封装，在vector 0.13版本之前，GHC允许这种名义类型的隐式强制转换，但新版本严格执行类型规则，不再允许这类隐式转换，因此直接混用hmatrix的Vector和vector库的V.Vector会触发类型错误。

代码修改方案 #

方案一：使用hmatrix原生函数（推荐） #

hmatrix本身提供了完整的向量操作API，完全可以替代vector库的函数，避免类型冲突：

import Numeric.LinearAlgebra

nnz :: Vector R -> Int
nnz = length . nonZero  -- nonZero返回所有非零元素组成的列表，取长度即非零元素数量

avg :: Vector R -> R
avg v = sumElements v / fromIntegral (dim v)  -- sumElements求和，dim获取向量长度

方案二：显式转换类型（保留vector库操作） #

如果必须使用vector库的函数，需要显式将hmatrix的Vector R转换为V.Vector Double，可以用hmatrix提供的unpack函数直接提取底层vector：

import Numeric.LinearAlgebra
import qualified Data.Vector.Storable as V

nnz :: Vector R -> Int
nnz v = V.length $ V.filter (/= 0) (unpack v)

avg :: Vector R -> R
avg v = let vec = unpack v
        in V.sum vec / fromIntegral (V.length vec)

也可以用convert函数（需要依赖convertible库）完成类型转换：

import Numeric.LinearAlgebra
import qualified Data.Vector.Storable as V
import Data.Convertible (convert)

nnz :: Vector R -> Int
nnz v = V.length $ V.filter (/= 0) (convert v :: V.Vector Double)

avg :: Vector R -> R
avg v = let vec = convert v :: V.Vector Double
        in V.sum vec / fromIntegral (V.length vec)

内容的提问来源于stack exchange，提问作者George Co