From 1ff59c415f8dc228547b6c7e43cef63d90687fad Mon Sep 17 00:00:00 2001
From: Liam Healy <liam@thinkpad.local>
Date: Mon, 18 Jan 2010 20:52:03 -0500
Subject: [PATCH] Correct default vector size in SV-solve for non-square

Correct default vector size in SV-solve for non-square matrices; patch
by Evgeniy Zhemchugov.
---
 linear-algebra/svd.lisp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/linear-algebra/svd.lisp b/linear-algebra/svd.lisp
index 9cbae06c..c88cd8b6 100644
--- a/linear-algebra/svd.lisp
+++ b/linear-algebra/svd.lisp
@@ -101,7 +101,7 @@
    references for details).")
 
 (defmfun SV-solve
-    (U S V b &optional (x (make-marray 'double-float :dimensions (dim0 U))))
+    (U S V b &optional (x (make-marray 'double-float :dimensions (dim0 V))))
   "gsl_linalg_SV_solve"
   (((mpointer U) :pointer) ((mpointer V) :pointer)
    ((mpointer S) :pointer)
-- 
GitLab