Back to index

glibc  2.9
w_j1f.c
Go to the documentation of this file.
00001 /* w_j1f.c -- float version of w_j1.c.
00002  * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
00003  */
00004 
00005 /*
00006  * ====================================================
00007  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
00008  *
00009  * Developed at SunPro, a Sun Microsystems, Inc. business.
00010  * Permission to use, copy, modify, and distribute this
00011  * software is freely granted, provided that this notice
00012  * is preserved.
00013  * ====================================================
00014  */
00015 
00016 #if defined(LIBM_SCCS) && !defined(lint)
00017 static char rcsid[] = "$NetBSD: w_j1f.c,v 1.3 1995/05/10 20:49:17 jtc Exp $";
00018 #endif
00019 
00020 /*
00021  * wrapper of j1f,y1f
00022  */
00023 
00024 #include <math.h>
00025 #include "math_private.h"
00026 
00027 #ifdef __STDC__
00028        float j1f(float x)          /* wrapper j1f */
00029 #else
00030        float j1f(x)                /* wrapper j1f */
00031        float x;
00032 #endif
00033 {
00034 #ifdef _IEEE_LIBM
00035        return __ieee754_j1f(x);
00036 #else
00037        float z;
00038        z = __ieee754_j1f(x);
00039        if(_LIB_VERSION == _IEEE_ || __isnanf(x) ) return z;
00040        if(fabsf(x)>(float)X_TLOSS) {
00041               /* j1(|x|>X_TLOSS) */
00042                return (float)__kernel_standard((double)x,(double)x,136);
00043        } else
00044            return z;
00045 #endif
00046 }
00047 
00048 #ifdef __STDC__
00049        float y1f(float x)          /* wrapper y1f */
00050 #else
00051        float y1f(x)                /* wrapper y1f */
00052        float x;
00053 #endif
00054 {
00055 #ifdef _IEEE_LIBM
00056        return __ieee754_y1f(x);
00057 #else
00058        float z;
00059        z = __ieee754_y1f(x);
00060        if(_LIB_VERSION == _IEEE_ || __isnanf(x) ) return z;
00061         if(x <= (float)0.0){
00062                 if(x==(float)0.0)
00063                     /* d= -one/(x-x); */
00064                     return (float)__kernel_standard((double)x,(double)x,110);
00065                 else
00066                     /* d = zero/(x-x); */
00067                     return (float)__kernel_standard((double)x,(double)x,111);
00068         }
00069        if(x>(float)X_TLOSS) {
00070               /* y1(x>X_TLOSS) */
00071                return (float)__kernel_standard((double)x,(double)x,137);
00072        } else
00073            return z;
00074 #endif
00075 }