Mercurial > hg > octave-nkf

diff scripts/general/quadl.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/general/quadl.m
+++ b/scripts/general/quadl.m
@@ -63,14 +63,19 @@
   if (nargin < 5)
     trace = []; 
   endif
+  if (isa (a, "single") || isa (b, "single"))
+    myeps = eps ("single");
+  else
+    myeps = eps;
+  endif
   if (isempty (tol))
-    tol = eps; 
+    tol = myeps; 
   endif
   if (isempty (trace))
     trace = 0; 
   endif
-  if (tol < eps)
-    tol = eps;
+  if (tol < myeps)
+    tol = myeps;
   endif
 
   m = (a+b)/2; 
@@ -119,7 +124,7 @@
   if (R > 0 && R < 1)
     tol = tol/R; 
   endif
-  is = s*abs(is)*tol/eps;
+  is = s*abs(is)*tol/myeps;
   if (is == 0)
     is = b-a;
   endif