1 2 /* 3 * Copyright (c) 1998, 2001, Oracle and/or its affiliates. All rights reserved. 4 * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER. 5 * 6 * This code is free software; you can redistribute it and/or modify it 7 * under the terms of the GNU General Public License version 2 only, as 8 * published by the Free Software Foundation. Oracle designates this 9 * particular file as subject to the "Classpath" exception as provided 10 * by Oracle in the LICENSE file that accompanied this code. 11 * 12 * This code is distributed in the hope that it will be useful, but WITHOUT 13 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 14 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 15 * version 2 for more details (a copy is included in the LICENSE file that 16 * accompanied this code). 17 * 18 * You should have received a copy of the GNU General Public License version 19 * 2 along with this work; if not, write to the Free Software Foundation, 20 * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. 21 * 22 * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA 23 * or visit www.oracle.com if you need additional information or have any 24 * questions. 25 */ 26 27 /* 28 * wrapper of j1,y1 29 */ 30 31 #include "fdlibm.h" 32 33 #ifdef __STDC__ 34 double j1(double x) /* wrapper j1 */ 35 #else 36 double j1(x) /* wrapper j1 */ 37 double x; 38 #endif 39 { 40 #ifdef _IEEE_LIBM 41 return __ieee754_j1(x); 42 #else 43 double z; 44 z = __ieee754_j1(x); 45 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z; 46 if(fabs(x)>X_TLOSS) { 47 return __kernel_standard(x,x,36); /* j1(|x|>X_TLOSS) */ 48 } else 49 return z; 50 #endif 51 } 52 53 #ifdef __STDC__ 54 double y1(double x) /* wrapper y1 */ 55 #else 56 double y1(x) /* wrapper y1 */ 57 double x; 58 #endif 59 { 60 #ifdef _IEEE_LIBM 61 return __ieee754_y1(x); 62 #else 63 double z; 64 z = __ieee754_y1(x); 65 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z; 66 if(x <= 0.0){ 67 if(x==0.0) 68 /* d= -one/(x-x); */ 69 return __kernel_standard(x,x,10); 70 else 71 /* d = zero/(x-x); */ 72 return __kernel_standard(x,x,11); 73 } 74 if(x>X_TLOSS) { 75 return __kernel_standard(x,x,37); /* y1(x>X_TLOSS) */ 76 } else 77 return z; 78 #endif 79 }