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