From ndpolyspace at gmail.com  Mon Jul 12 23:52:51 2021
From: ndpolyspace at gmail.com (Richard Holcombe)
Date: Mon, 12 Jul 2021 19:52:51 -0400
Subject: A question about cholesky-solve
Message-ID: <CADEd6UzS8cqCSUb=4OeOLMM-XoiakBi_bUzy3drt44zRJaPWJw@mail.gmail.com>

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.common-lisp.net/pipermail/gsll-devel/attachments/20210712/eb879b25/attachment.html>