On Sat, 15 Mar 2008, Liam Healy wrote:
If you update GSLL, you will now find the file ;; $Id: lu.lisp 36 2008-03-15 19:55:54Z lhealy $
with a defun at the end using the new definition of #'lu-decomp (defun invert-matrix (mat) "Invert the matrix." (letm ((mmat (matrix-double-float mat)) (dim (array-dimension mat 0)) (per (permutation dim)) (inv (matrix-double-float dim dim))) (lu-decomp mmat per) (lu-invert mmat per inv) (data inv)))
(also note that GSL matrices are now declared matrix-double-float). I hope this is what you're looking for.
Dear Liam,
With the latest svn checkout (36) sbcl 1.0.15 raises the following error:
There is no class named GSLL::MATRIX-DOUBLE-FLOAT. [Condition of type SIMPLE-ERROR]
Restarts: 0: [RETRY] Retry performing #<ASDF:LOAD-OP NIL {1003992C31}> on #<ASDF:CL-SOURCE-FILE "lu" {1003CB1B01}>. 1: [ACCEPT] Continue, treating #<ASDF:LOAD-OP NIL {1003992C31}> on #<ASDF:CL-SOURCE-FILE "lu" {1003CB1B01}> as having been successful. 2: [RETRY] Retry installation 3: [ABORT] Return to SLIME's top level. 4: [TERMINATE-THREAD] Terminate this thread (#<THREAD "repl-thread" {1003BC5301}>)
when compiling the lu code.
Jason