# HG changeset patch # User John W. Eaton # Date 1326391246 18000 # Node ID ce731bb6a1911dfab99edfcbd448c1ebcdc11d40 # Parent 35ce1eab7400af0b1faf404e3b863b0036eaa93c besseli: don't change sign of result for negative integer orders * lo-specfun.cc (zbesi): Omit sign change for negative integer orders. diff -r 35ce1eab7400 -r ce731bb6a191 liboctave/lo-specfun.cc --- a/liboctave/lo-specfun.cc Thu Jan 12 12:10:33 2012 -0500 +++ b/liboctave/lo-specfun.cc Thu Jan 12 13:00:46 2012 -0500 @@ -857,8 +857,6 @@ // zbesi can overflow as z->0, and cause troubles for generic case below alpha = -alpha; Complex tmp = zbesi (z, alpha, kode, ierr); - if ((static_cast (alpha)) & 1) - tmp = - tmp; retval = bessel_return_value (tmp, ierr); } else