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
