Mercurial > hg > octave-nkf

diff liboctave/lo-specfun.cc @ 3225:7aae2c3636a7

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1998-12-04 23:20:12 by jwe]
author jwe
date Fri, 04 Dec 1998 23:20:26 +0000
parents 3deb1105fbc1
children 7da18459c08b
line wrap: on
line diff
--- a/liboctave/lo-specfun.cc
+++ b/liboctave/lo-specfun.cc
@@ -612,6 +612,9 @@
 
   F77_FCN (zairy, ZAIRY) (zr, zi, id, kode, ar, ai, nz, ierr);
 
+  if (zi == 0.0 && (! scaled || zr > 0.0))
+    ai = 0.0;
+
   return bessel_return_value (Complex (ar, ai), ierr);
 }
 
@@ -630,6 +633,9 @@
 
   F77_FCN (zbiry, ZBIRY) (zr, zi, id, kode, ar, ai, ierr);
 
+  if (zi == 0.0 && (! scaled || zr > 0.0))
+    ai = 0.0;
+
   return bessel_return_value (Complex (ar, ai), ierr);
 }