1*1e651e1eSRoland Levillain 2*1e651e1eSRoland Levillain /* @(#)w_j1.c 1.3 95/01/18 */ 3*1e651e1eSRoland Levillain /* 4*1e651e1eSRoland Levillain * ==================================================== 5*1e651e1eSRoland Levillain * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 6*1e651e1eSRoland Levillain * 7*1e651e1eSRoland Levillain * Developed at SunSoft, a Sun Microsystems, Inc. business. 8*1e651e1eSRoland Levillain * Permission to use, copy, modify, and distribute this 9*1e651e1eSRoland Levillain * software is freely granted, provided that this notice 10*1e651e1eSRoland Levillain * is preserved. 11*1e651e1eSRoland Levillain * ==================================================== 12*1e651e1eSRoland Levillain */ 13*1e651e1eSRoland Levillain 14*1e651e1eSRoland Levillain /* 15*1e651e1eSRoland Levillain * wrapper of j1,y1 16*1e651e1eSRoland Levillain */ 17*1e651e1eSRoland Levillain 18*1e651e1eSRoland Levillain #include "fdlibm.h" 19*1e651e1eSRoland Levillain 20*1e651e1eSRoland Levillain #ifdef __STDC__ ieee_j1(double x)21*1e651e1eSRoland Levillain double ieee_j1(double x) /* wrapper j1 */ 22*1e651e1eSRoland Levillain #else 23*1e651e1eSRoland Levillain double ieee_j1(x) /* wrapper j1 */ 24*1e651e1eSRoland Levillain double x; 25*1e651e1eSRoland Levillain #endif 26*1e651e1eSRoland Levillain { 27*1e651e1eSRoland Levillain #ifdef _IEEE_LIBM 28*1e651e1eSRoland Levillain return __ieee754_j1(x); 29*1e651e1eSRoland Levillain #else 30*1e651e1eSRoland Levillain double z; 31*1e651e1eSRoland Levillain z = __ieee754_j1(x); 32*1e651e1eSRoland Levillain if(_LIB_VERSION == _IEEE_ || ieee_isnan(x) ) return z; 33*1e651e1eSRoland Levillain if(ieee_fabs(x)>X_TLOSS) { 34*1e651e1eSRoland Levillain return __kernel_standard(x,x,36); /* ieee_j1(|x|>X_TLOSS) */ 35*1e651e1eSRoland Levillain } else 36*1e651e1eSRoland Levillain return z; 37*1e651e1eSRoland Levillain #endif 38*1e651e1eSRoland Levillain } 39*1e651e1eSRoland Levillain 40*1e651e1eSRoland Levillain #ifdef __STDC__ ieee_y1(double x)41*1e651e1eSRoland Levillain double ieee_y1(double x) /* wrapper y1 */ 42*1e651e1eSRoland Levillain #else 43*1e651e1eSRoland Levillain double ieee_y1(x) /* wrapper y1 */ 44*1e651e1eSRoland Levillain double x; 45*1e651e1eSRoland Levillain #endif 46*1e651e1eSRoland Levillain { 47*1e651e1eSRoland Levillain #ifdef _IEEE_LIBM 48*1e651e1eSRoland Levillain return __ieee754_y1(x); 49*1e651e1eSRoland Levillain #else 50*1e651e1eSRoland Levillain double z; 51*1e651e1eSRoland Levillain z = __ieee754_y1(x); 52*1e651e1eSRoland Levillain if(_LIB_VERSION == _IEEE_ || ieee_isnan(x) ) return z; 53*1e651e1eSRoland Levillain if(x <= 0.0){ 54*1e651e1eSRoland Levillain if(x==0.0) 55*1e651e1eSRoland Levillain /* d= -one/(x-x); */ 56*1e651e1eSRoland Levillain return __kernel_standard(x,x,10); 57*1e651e1eSRoland Levillain else 58*1e651e1eSRoland Levillain /* d = zero/(x-x); */ 59*1e651e1eSRoland Levillain return __kernel_standard(x,x,11); 60*1e651e1eSRoland Levillain } 61*1e651e1eSRoland Levillain if(x>X_TLOSS) { 62*1e651e1eSRoland Levillain return __kernel_standard(x,x,37); /* ieee_y1(x>X_TLOSS) */ 63*1e651e1eSRoland Levillain } else 64*1e651e1eSRoland Levillain return z; 65*1e651e1eSRoland Levillain #endif 66*1e651e1eSRoland Levillain } 67