Mercurial > hg > octave-lyh

diff scripts/sparse/normest.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 cadc73247d65
line wrap: on
line diff
--- a/scripts/sparse/normest.m
+++ b/scripts/sparse/normest.m
@@ -32,8 +32,14 @@
   if (nargin < 2)
     tol = 1e-6;
   endif
-  if (tol < eps)
-    tol = eps
+  if (isa (A, "single"))
+    if (tol < eps ("single"))
+      tol = eps ("single");
+    endif
+  else
+    if (tol < eps)
+      tol = eps
+    endif
   endif
   if (ndims(A) != 2)
     error ("normest: A must be a matrix");