Skip to content
No description or website provided.
C HTML Makefile
Find file
Latest commit fcd1725 Andreas Gampe committed with Gerrit Code Review Merge "Fdlibm: Silence warnings"
Failed to load latest commit information. Fdlibm: Silence warnings Fix x86 build (problems spotted by StrictMath tests).
MD5 auto import from //depot/cupcake/@135843
MODULE_LICENSE_BSD_LIKE auto import from //depot/cupcake/@135843
NOTICE Automated import from //branches/donutburger/...@140818,140818
changes auto import from //depot/cupcake/@135843
configure auto import from //depot/cupcake/@135843 auto import from //depot/cupcake/@135843
e_acos.c auto import from //depot/cupcake/@135843
e_acosh.c auto import from //depot/cupcake/@135843
e_asin.c auto import from //depot/cupcake/@135843
e_atan2.c auto import from //depot/cupcake/@135843
e_atanh.c auto import from //depot/cupcake/@135843
e_cosh.c auto import from //depot/cupcake/@135843
e_exp.c auto import from //depot/cupcake/@135843
e_fmod.c auto import from //depot/cupcake/@135843
e_gamma.c auto import from //depot/cupcake/@135843
e_gamma_r.c auto import from //depot/cupcake/@135843
e_hypot.c auto import from //depot/cupcake/@135843
e_j0.c auto import from //depot/cupcake/@135843
e_j1.c auto import from //depot/cupcake/@135843
e_jn.c auto import from //depot/cupcake/@135843
e_lgamma.c auto import from //depot/cupcake/@135843
e_lgamma_r.c auto import from //depot/cupcake/@135843
e_log.c auto import from //depot/cupcake/@135843
e_log10.c auto import from //depot/cupcake/@135843
e_pow.c auto import from //depot/cupcake/@135843
e_rem_pio2.c auto import from //depot/cupcake/@135843
e_remainder.c auto import from //depot/cupcake/@135843
e_scalb.c auto import from //depot/cupcake/@135843
e_sinh.c auto import from //depot/cupcake/@135843
e_sqrt.c auto import from //depot/cupcake/@135843
fdlibm.h [ENDIAN] Select endianess correctly for MIPS
index auto import from //depot/cupcake/@135843
index.html auto import from //depot/cupcake/@135843
k_cos.c auto import from //depot/cupcake/@135843
k_rem_pio2.c auto import from //depot/cupcake/@135843
k_sin.c auto import from //depot/cupcake/@135843
k_standard.c auto import from //depot/cupcake/@135843
k_tan.c auto import from //depot/cupcake/@135843 auto import from //depot/cupcake/@135843
readme auto import from //depot/cupcake/@135843
s_asinh.c auto import from //depot/cupcake/@135843
s_cbrt.c auto import from //depot/cupcake/@135843
s_ceil.c auto import from //depot/cupcake/@135843
s_copysign.c auto import from //depot/cupcake/@135843
s_cos.c auto import from //depot/cupcake/@135843
s_erf.c auto import from //depot/cupcake/@135843
s_expm1.c auto import from //depot/cupcake/@135843
s_fabs.c auto import from //depot/cupcake/@135843
s_finite.c auto import from //depot/cupcake/@135843
s_floor.c auto import from //depot/cupcake/@135843
s_frexp.c auto import from //depot/cupcake/@135843
s_ilogb.c auto import from //depot/cupcake/@135843
s_isnan.c auto import from //depot/cupcake/@135843
s_ldexp.c auto import from //depot/cupcake/@135843
s_lib_version.c auto import from //depot/cupcake/@135843
s_log1p.c auto import from //depot/cupcake/@135843
s_logb.c auto import from //depot/cupcake/@135843
s_matherr.c auto import from //depot/cupcake/@135843
s_modf.c auto import from //depot/cupcake/@135843
s_nextafter.c auto import from //depot/cupcake/@135843
s_rint.c auto import from //depot/cupcake/@135843
s_scalbn.c auto import from //depot/cupcake/@135843
s_signgam.c auto import from //depot/cupcake/@135843
s_significand.c auto import from //depot/cupcake/@135843
s_sin.c auto import from //depot/cupcake/@135843
s_tan.c auto import from //depot/cupcake/@135843
s_tanh.c auto import from //depot/cupcake/@135843
w_acos.c auto import from //depot/cupcake/@135843
w_acosh.c auto import from //depot/cupcake/@135843
w_asin.c auto import from //depot/cupcake/@135843
w_atan2.c auto import from //depot/cupcake/@135843
w_atanh.c auto import from //depot/cupcake/@135843
w_cosh.c auto import from //depot/cupcake/@135843
w_exp.c auto import from //depot/cupcake/@135843
w_fmod.c auto import from //depot/cupcake/@135843
w_gamma.c auto import from //depot/cupcake/@135843
w_gamma_r.c auto import from //depot/cupcake/@135843
w_hypot.c auto import from //depot/cupcake/@135843
w_j0.c auto import from //depot/cupcake/@135843
w_j1.c auto import from //depot/cupcake/@135843
w_jn.c auto import from //depot/cupcake/@135843
w_lgamma.c auto import from //depot/cupcake/@135843
w_lgamma_r.c auto import from //depot/cupcake/@135843
w_log.c auto import from //depot/cupcake/@135843
w_log10.c auto import from //depot/cupcake/@135843
w_pow.c auto import from //depot/cupcake/@135843
w_remainder.c auto import from //depot/cupcake/@135843
w_scalb.c auto import from //depot/cupcake/@135843
w_sinh.c auto import from //depot/cupcake/@135843
w_sqrt.c auto import from //depot/cupcake/@135843


 	  * Announcing FDLIBM Version 5.3 *
	developed at Sun Microsystems, Inc. 

What's new in FDLIBM 5.3?

	To build FDLIBM, edit the supplied Makefile or create
	a local Makefile by running "sh configure" 
	using the supplied configure script contributed by Nelson Beebe


    1. e_pow.c incorrect results when 
	x is very close to -1.0 and y is very large, e.g.
  	pow(-1.0000000000000002e+00,4.5035996273704970e+15) = 0
  	pow(-9.9999999999999978e-01,4.5035996273704970e+15) = 0
	Correct results are close to -e and -1/e.

    2. k_tan.c error was > 1 ulp target for FDLIBM
	5.2: Worst error at least 1.45 ulp at
	tan(1.7765241907548024E+269) = 1.7733884462610958E+16
	5.3: Worst error 0.96 ulp


    3. Compiler failure on non-standard code
	Statements like
	            *(1+(int*)&t1) = 0;
	are not standard C and cause some optimizing compilers (e.g. GCC)
	to generate bad code under optimization.    These cases
	are to be addressed in the next release.
FDLIBM (Freely Distributable LIBM) is a C math library 
for machines that support IEEE 754 floating-point arithmetic. 
In this release, only double precision is supported.

FDLIBM is intended to provide a reasonably portable (see 
assumptions below), reference quality (below one ulp for
major functions like sin,cos,exp,log) math library 
(libm.a).  For a copy of FDLIBM, please see

FDLIBM (double precision version) assumes:
 a.  IEEE 754 style (if not precise compliance) arithmetic;
 b.  32 bit 2's complement integer arithmetic;
 c.  Each double precision floating-point number must be in IEEE 754 
     double format, and that each number can be retrieved as two 32-bit 
     integers through the using of pointer bashing as in the example 

     Example: let y = 2.0
	double fp number y: 	2.0
	IEEE double format:	0x4000000000000000

	Referencing y as two integers:
	*(int*)&y,*(1+(int*)&y) =	{0x40000000,0x0} (on sparc)
					{0x0,0x40000000} (on 386)

	Note: Four macros are defined in fdlibm.h to handle this kind of

	__HI(x)		the high part of a double x 
			(sign,exponent,the first 21 significant bits)
	__LO(x)		the least 32 significant bits of x
	__HIp(x)	same as __HI except that the argument is a pointer
			to a double
	__LOp(x)	same as __LO except that the argument is a pointer
			to a double
	To ensure obtaining correct ordering, one must define  __LITTLE_ENDIAN
	during compilation for little endian machine (like 386,486). The 
	default is big endian.

	If the behavior of pointer bashing is undefined, one may hack on the 
	macro in fdlibm.h.
  d. IEEE exceptions may trigger "signals" as is common in Unix

All exception cases in the FDLIBM functions will be mapped
to one of the following four exceptions:

   +-huge*huge, +-tiny*tiny,    +-1.0/0.0,	+-0.0/0.0
    (overflow)	(underflow)  (divided-by-zero) 	(invalid)

For example, ieee_log(0) is a singularity and is thus mapped to 
	-1.0/0.0 = -infinity.
That is, FDLIBM's log will compute -one/zero and return the
computed value.  On an IEEE machine, this will trigger the 
divided-by-zero exception and a negative infinity is returned by 

Similarly, ieee_exp(-huge) will be mapped to tiny*tiny to generate
an underflow signal. 

The default FDLIBM functions (compiled with -D_IEEE_LIBM flag)  
are in "IEEE spirit" (i.e., return the most reasonable result in 
floating-point arithmetic). If one wants FDLIBM to comply with
standards like SVID, X/OPEN, or POSIX/ANSI, then one can 
create a multi-standard compliant FDLIBM. In this case, each
function in FDLIBM is actually a standard compliant wrapper

File organization:
    1. For FDLIBM's kernel (internal) function,
		File name	Entry point
		k_sin.c		__kernel_sin
		k_tan.c		__kernel_tan
    2. For functions that have no standards conflict 
		File name	Entry point
		s_sin.c		sin
		s_erf.c		erf
    3. Ieee754 core functions
		File name	Entry point
		e_exp.c		__ieee754_exp
		e_sinh.c	__ieee754_sinh
    4. Wrapper functions
		File name	Entry point
		w_exp.c		exp
		w_sinh.c	sinh

Wrapper functions will twist the result of the ieee754 
function to comply to the standard specified by the value 
    if _LIB_VERSION = _IEEE_, return the ieee754 result;
    if _LIB_VERSION = _SVID_, return SVID result;
    if _LIB_VERSION = _XOPEN_, return XOPEN result;
    if _LIB_VERSION = _POSIX_, return POSIX/ANSI result.
(These are macros, see fdlibm.h for their definition.)

There are two types of libm.a. One is IEEE only, and the other is
multi-standard compliant (supports IEEE,XOPEN,POSIX/ANSI,SVID).

To create the IEEE only libm.a, use 
	    make "CFLAGS = -D_IEEE_LIBM"	 
This will create an IEEE libm.a, which is smaller in size, and 
somewhat faster.

To create a multi-standard compliant libm, use
    make "CFLAGS = -D_IEEE_MODE"   --- multi-standard fdlibm: default
					 to IEEE
    make "CFLAGS = -D_XOPEN_MODE"  --- multi-standard fdlibm: default
					 to X/OPEN
    make "CFLAGS = -D_POSIX_MODE"  --- multi-standard fdlibm: default
					 to POSIX/ANSI
    make "CFLAGS = -D_SVID3_MODE"  --- multi-standard fdlibm: default
					 to SVID

Here is how one makes a SVID compliant libm.
    Make the library by
		make "CFLAGS = -D_SVID3_MODE".
    The libm.a of FDLIBM will be multi-standard compliant and 
    _LIB_VERSION is initialized to the value _SVID_ . 

		double ieee_y0();
		printf("y0(1e300) = %1.20e\n",y0(1e300));

    % cc example1.c libm.a
    % a.out
    y0: TLOSS error
    ieee_y0(1e300) = 0.00000000000000000000e+00

It is possible to change the default standard in multi-standard 
fdlibm. Here is an example of how to do it:
	#include "fdlibm.h"	/* must include FDLIBM's fdlibm.h */
		double ieee_y0();
		printf("IEEE: ieee_y0(1e300) = %1.20e\n",y0(1e300));
		printf("XOPEN ieee_y0(1e300) = %1.20e\n",y0(1e300));
		printf("POSIX ieee_y0(1e300) = %1.20e\n",y0(1e300));
		printf("SVID  ieee_y0(1e300) = %1.20e\n",y0(1e300));

    % cc example2.c libm.a
    % a.out
      IEEE: ieee_y0(1e300) = -1.36813604503424810557e-151
      XOPEN ieee_y0(1e300) = 0.00000000000000000000e+00
      POSIX ieee_y0(1e300) = 0.00000000000000000000e+00
      y0: TLOSS error
      SVID  ieee_y0(1e300) = 0.00000000000000000000e+00

Note:	Here _LIB_VERSION is a global variable. If global variables 
	are forbidden, then one should modify fdlibm.h to change
	_LIB_VERSION to be a global constant. In this case, one
	may not change the value of _LIB_VERSION as in example2.

	Care must be taken when installing FDLIBM over existing
	All co-existing function prototypes must agree, otherwise
	users will encounter mysterious failures.

	So far, the only known likely conflict is the declaration 
	of the IEEE recommended function scalb:

		double ieee_scalb(double,double)	(1)	SVID3 defined
		double ieee_scalb(double,int)	(2)	IBM,DEC,...

	FDLIBM follows Sun definition and use (1) as default. 
	If one's existing libm.a uses (2), then one may raise
	the flags _SCALB_INT during the compilation of FDLIBM
	to get the correct function prototype.
	(E.g., make "CFLAGS = -D_IEEE_LIBM -D_SCALB_INT".)
	NOTE that if -D_SCALB_INT is raised, it won't be SVID3

Please send comments and bug reports to the electronic mail address
suggested by: 
		fdlibm-comments AT

Something went wrong with that request. Please try again.