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
