Mercurial > hg > octave-nkf

diff scripts/linear-algebra/rref.m @ 7795:df9519e9990c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Handle single precision eps values
author David Bateman <dbateman@free.fr>
date Mon, 12 May 2008 22:57:11 +0200
parents a1dbe9d80eee
children eb63fbe60fab
line wrap: on
line diff
--- a/scripts/linear-algebra/rref.m
+++ b/scripts/linear-algebra/rref.m
@@ -44,7 +44,11 @@
   [rows, cols] = size (A);
 
   if (nargin < 2)
-    tolerance = eps * max (rows, cols) * norm (A, inf);
+    if (isa (A, "single"))
+      tolerance = eps ("single") * max (rows, cols) * norm (A, inf ("single"));
+    else
+      tolerance = eps * max (rows, cols) * norm (A, inf);
+    endif
   endif
 
   used = zeros (1, cols);