vấn đề ghc với hmatrix

Question

Tôi đang cố gắng sử dụng hmatrix và hmatrix-gsl-stats. Nhưng tôi nhận được kết quả tuyệt vời kỳ lạ. Ai có thể tái sản xuất những gì tôi thấy không?

import Numeric.LinearAlgebra 
import Numeric.GSL.Fitting.Linear 

main :: IO() 
main = do 
    let y = fromList [1,2,3,4,5] 
    let x = trans $ fromLists [[1,2,3,4,5], [1,1,1,1,1]] 
    let x3 = fromLists [[1,1], [2,1], [3,1], [4,1], [5,1]] 
    print $ x == x3 
    print $ multifit x y 
    print $ multifit x3 y 

Vì vậy, tôi chỉ cần biên dịch này và chạy nó, và tôi nhận được:

True 
(fromList [0.6121951219512196,0.39756097560975645],(2><2) 
[ 0.5100317271465397, -0.46568114217727535 
, -0.46568114217727535,  0.82048582193139 ],27.27560975609756) 
(fromList [1.0000000000000002,-8.881784197001252e-16],(2><2) 
[ 2.1364982849735737e-32, -6.40949485492072e-32 
, -6.40949485492072e-32, 2.350148113470931e-31 ],6.409494854920721e-31) 

Am tôi thiếu một cái gì đó hoàn toàn rõ ràng?

Answer 1

Tôi nghĩ đây là vấn đề với cách chuyển vị ma trận được xử lý nội bộ. Khi bạn gọi trans, nó chỉ là một chút lộn xộn trong siêu dữ liệu ... rõ ràng là không đang được kiểm tra khi truyền ma trận qua FFI.

Nếu bạn thử điều này, multifit (reshape (cols x) $ flatten x) y, trên đường áp chót, bạn có thể buộc chuyển vị.