Mercurial > hg > octave-nkf

diff liboctave/array/dDiagMatrix.cc @ 18549:16b0cd465ecd

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Handle special case of 0 for pinv with Diagonal matrices. * CDiagMatrix.cc (pseudo_inverse), dDiagMatrix.cc (pseudo_inverse), fCDiagMatrix.cc (pseudo_inverse), fDiagMatrix.cc (pseudo_inverse): Check for special case where element is 0 to avoid a division by zero error.
author Rik <rik@octave.org>
date Mon, 17 Feb 2014 10:04:27 -0800
parents c5a101de2d88
children 1fa5bba16218
line wrap: on
line diff
--- a/liboctave/array/dDiagMatrix.cc
+++ b/liboctave/array/dDiagMatrix.cc
@@ -302,10 +302,11 @@
 
   for (octave_idx_type i = 0; i < len; i++)
     {
-      if (std::abs (elem (i, i)) < tol)
+      double val = std::abs (elem (i, i));
+      if (val < tol || val == 0.0)
         retval.elem (i, i) = 0.0;
       else
-        retval.elem (i, i) = 1.0 / elem (i, i);
+        retval.elem (i, i) = 1.0 / val;
     }
 
   return retval;