I am running sbcl 2.1.6 on linux. I installed gsll using quicklisp I picked an arbitrary linear system to solve:
(gsl:cholesky-solve (grid:make-foreign-array 'double-float :dimensions '(3 3) :initial-contents '((1 2 3)(4 5 6)(12 42 65))) (grid:make-foreign-array 'double-float :dimensions '(3) :initial-contents '(4 5 3)) t)
and it comes back with: #m(6.002328994082840d0 -0.534239053254438d0 0.011218934911243d0)
but the correct answer is [12.524;-28.048;15.857]
I can't figure out the number it is giving me. What could I be missing?
Richard