RTAI API: base/math/w_j1.c Source File

00001 /* @(#)w_j1.c 5.1 93/09/24 */ 00002 /* 00003 * ==================================================== 00004 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 00005 * 00006 * Developed at SunPro, a Sun Microsystems, Inc. business. 00007 * Permission to use, copy, modify, and distribute this 00008 * software is freely granted, provided that this notice 00009 * is preserved. 00010 * ==================================================== 00011 */ 00012 00013 #if defined(LIBM_SCCS) && !defined(lint) 00014 static char rcsid[] = "$NetBSD: w_j1.c,v 1.6 1995/05/10 20:49:15 jtc Exp $"; 00015 #endif 00016 00017 /* 00018 * wrapper of j1,y1 00019 */ 00020 00021 #include "math.h" 00022 #include "mathP.h" 00023 00024 #ifdef __STDC__ 00025 double j1(double x) /* wrapper j1 */ 00026 #else 00027 double j1(x) /* wrapper j1 */ 00028 double x; 00029 #endif 00030 { 00031 #ifdef _IEEE_LIBM 00032 return __ieee754_j1(x); 00033 #else 00034 double z; 00035 z = __ieee754_j1(x); 00036 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z; 00037 if(fabs(x)>X_TLOSS) { 00038 return __kernel_standard(x,x,36); /* j1(|x|>X_TLOSS) */ 00039 } else 00040 return z; 00041 #endif 00042 } 00043 00044 #ifdef __STDC__ 00045 double y1(double x) /* wrapper y1 */ 00046 #else 00047 double y1(x) /* wrapper y1 */ 00048 double x; 00049 #endif 00050 { 00051 #ifdef _IEEE_LIBM 00052 return __ieee754_y1(x); 00053 #else 00054 double z; 00055 z = __ieee754_y1(x); 00056 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z; 00057 if(x <= 0.0){ 00058 if(x==0.0) 00059 /* d= -one/(x-x); */ 00060 return __kernel_standard(x,x,10); 00061 else 00062 /* d = zero/(x-x); */ 00063 return __kernel_standard(x,x,11); 00064 } 00065 if(x>X_TLOSS) { 00066 return __kernel_standard(x,x,37); /* y1(x>X_TLOSS) */ 00067 } else 00068 return z; 00069 #endif 00070 }