00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013 #if defined(LIBM_SCCS) && !defined(lint)
00014 static char rcsid[] = "$NetBSD: w_jn.c,v 1.6 1995/05/10 20:49:19 jtc Exp $";
00015 #endif
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043 #include "math.h"
00044 #include "mathP.h"
00045
00046 #ifdef __STDC__
00047 double jn(int n, double x)
00048 #else
00049 double jn(n,x)
00050 double x; int n;
00051 #endif
00052 {
00053 #ifdef _IEEE_LIBM
00054 return __ieee754_jn(n,x);
00055 #else
00056 double z;
00057 z = __ieee754_jn(n,x);
00058 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;
00059 if(fabs(x)>X_TLOSS) {
00060 return __kernel_standard((double)n,x,38);
00061 } else
00062 return z;
00063 #endif
00064 }
00065
00066 #ifdef __STDC__
00067 double yn(int n, double x)
00068 #else
00069 double yn(n,x)
00070 double x; int n;
00071 #endif
00072 {
00073 #ifdef _IEEE_LIBM
00074 return __ieee754_yn(n,x);
00075 #else
00076 double z;
00077 z = __ieee754_yn(n,x);
00078 if(_LIB_VERSION == _IEEE_ || isnan(x) ) return z;
00079 if(x <= 0.0){
00080 if(x==0.0)
00081
00082 return __kernel_standard((double)n,x,12);
00083 else
00084
00085 return __kernel_standard((double)n,x,13);
00086 }
00087 if(x>X_TLOSS) {
00088 return __kernel_standard((double)n,x,39);
00089 } else
00090 return z;
00091 #endif
00092 }