Print this page
*** 69,79 ****
one = 1.0;
GENERIC
jn(int n, GENERIC x) {
int i, sgn;
! GENERIC a, b, temp;
GENERIC z, w, ox, on;
/* J(-n,x) = (-1)^n * J(n, x), J(n, -x) = (-1)^n * J(n, x)
* Thus, J(-n,x) = J(n,-x)
*/
--- 69,79 ----
one = 1.0;
GENERIC
jn(int n, GENERIC x) {
int i, sgn;
! GENERIC a, b, temp = 0;
GENERIC z, w, ox, on;
/* J(-n,x) = (-1)^n * J(n, x), J(n, -x) = (-1)^n * J(n, x)
* Thus, J(-n,x) = J(n,-x)
*/
*** 215,235 ****
GENERIC
yn(int n, GENERIC x) {
int i;
int sign;
! GENERIC a, b, temp, ox, on;
ox = x; on = (GENERIC)n;
if(isnan(x)) return x*x; /* + -> * for Cheetah */
! if (x <= zero)
! if(x==zero)
/* return -one/zero; */
return _SVID_libm_err((GENERIC)n,x,12);
! else
/* return zero/zero; */
return _SVID_libm_err((GENERIC)n,x,13);
if (!((int) _lib_version == libm_ieee ||
(__xpg6 & _C99SUSv3_math_errexcept) != 0)) {
if(x > X_TLOSS) return _SVID_libm_err(on,ox,39);
}
sign = 1;
--- 215,237 ----
GENERIC
yn(int n, GENERIC x) {
int i;
int sign;
! GENERIC a, b, temp = 0, ox, on;
ox = x; on = (GENERIC)n;
if(isnan(x)) return x*x; /* + -> * for Cheetah */
! if (x <= zero) {
! if(x==zero) {
/* return -one/zero; */
return _SVID_libm_err((GENERIC)n,x,12);
! } else {
/* return zero/zero; */
return _SVID_libm_err((GENERIC)n,x,13);
+ }
+ }
if (!((int) _lib_version == libm_ieee ||
(__xpg6 & _C99SUSv3_math_errexcept) != 0)) {
if(x > X_TLOSS) return _SVID_libm_err(on,ox,39);
}
sign = 1;