vx32

Local 9vx git repository for patches.
git clone git://r-36.net/vx32
Log | Files | Refs

w_jnf.c (1006B)

      1 /* w_jnf.c -- float version of w_jn.c.
      2  * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
      3  */
      4 
      5 /*
      6  * ====================================================
      7  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
      8  *
      9  * Developed at SunPro, a Sun Microsystems, Inc. business.
     10  * Permission to use, copy, modify, and distribute this
     11  * software is freely granted, provided that this notice
     12  * is preserved.
     13  * ====================================================
     14  */
     15 
     16 #ifndef lint
     17 static char rcsid[] = "$FreeBSD: src/lib/msun/src/w_jnf.c,v 1.8 2002/05/28 18:15:04 alfred Exp $";
     18 #endif
     19 
     20 #include "math.h"
     21 #include "math_private.h"
     22 
     23 float
     24 jnf(int n, float x)	/* wrapper jnf */
     25 {
     26 #ifdef _IEEE_LIBM
     27 	return __ieee754_jnf(n,x);
     28 #else
     29 	float z;
     30 	z = __ieee754_jnf(n,x);
     31 	if(_LIB_VERSION == _IEEE_ || isnanf(x) ) return z;
     32 	if(fabsf(x)>(float)X_TLOSS) {
     33 	    /* jn(|x|>X_TLOSS,n) */
     34 	    return (float)__kernel_standard((double)n,(double)x,138);
     35 	} else
     36 	    return z;
     37 #endif
     38 }