[common-math-cvs] r4 - trunk/common-math/numerics/linear-algebra

[mantoniotti at common-lisp.net]

Author: mantoniotti
Date: Sat Sep 23 18:23:10 2006
New Revision: 4

Modified:
   trunk/common-math/numerics/linear-algebra/inversion.lisp
   trunk/common-math/numerics/linear-algebra/lu-decomposition.lisp
Log:
Minor changes.

Modified: trunk/common-math/numerics/linear-algebra/inversion.lisp
==============================================================================
--- trunk/common-math/numerics/linear-algebra/inversion.lisp	(original)
+++ trunk/common-math/numerics/linear-algebra/inversion.lisp	Sat Sep 23 18:23:10 2006
@@ -20,10 +20,15 @@
   (inverse (matrix-data a)))
 
 
+(defmethod inverse ((a array)) (error "Unimplemented ~S" a))
 
 
 (defmethod inverse/double-float ((a matrix))
   (inverse/double-float (matrix-data a)))
 
 
+(defmethod inverse/double-float ((a array)) (error "Unimplemented ~S" a))
+
+
+
 ;;; end of file -- inversion.lisp --

Modified: trunk/common-math/numerics/linear-algebra/lu-decomposition.lisp
==============================================================================
--- trunk/common-math/numerics/linear-algebra/lu-decomposition.lisp	(original)
+++ trunk/common-math/numerics/linear-algebra/lu-decomposition.lisp	Sat Sep 23 18:23:10 2006
@@ -243,9 +243,8 @@
           do (cond ((not (cl:zerop ii))
                     (loop for j of-type fixnum from ii below (1- i)
                           do (decf sum (* (aref lu i j) (aref b i)))))
-                   ((not (cl:zerop sum))
-                    (setf (aref b i) sum)))
-          )
+                   ((not (cl:zerop sum)) (setf ii i)))
+          do (setf (aref b i) sum))
     (loop for i of-type fixnum from (1- n) downto 0
           for sum of-type number = (aref b i)
           do (loop for j of-type fixnum from (1+ i) below n