2 /* @(#)w_j1.c 5.1 93/09/24 */
4 * ====================================================
5 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
7 * Developed at SunPro, a Sun Microsystems, Inc. business.
8 * Permission to use, copy, modify, and distribute this
9 * software is freely granted, provided that this notice
11 * ====================================================
21 #ifndef _DOUBLE_IS_32BITS
24 double j1(double x) /* wrapper j1 */
26 double j1(x) /* wrapper j1 */
31 return __ieee754_j1(x);
36 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;
42 exc.arg1 = exc.arg2 = x;
44 if (_LIB_VERSION == _POSIX_)
46 else if (!matherr(&exc)) {
58 double y1(double x) /* wrapper y1 */
60 double y1(x) /* wrapper y1 */
65 return __ieee754_y1(x);
70 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;
76 SET_HIGH_WORD(inf,0x7ff00000); /* set inf to infinite */
78 /* y1(0) = -inf or y1(x<0) = NaN */
79 exc.type = DOMAIN; /* should be SING for IEEE */
82 exc.arg1 = exc.arg2 = x;
83 if (_LIB_VERSION == _SVID_)
86 exc.retval = -HUGE_VAL;
87 if (_LIB_VERSION == _POSIX_)
89 else if (!matherr(&exc)) {
101 exc.arg1 = exc.arg2 = x;
103 if (_LIB_VERSION == _POSIX_)
105 else if (!matherr(&exc)) {
116 #endif /* defined(_DOUBLE_IS_32BITS) */