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;