Home | History | Annotate | Download | only in libm
      1 /*-
      2  * Copyright (c) 2003 Mike Barcroft <mike (at) FreeBSD.org>
      3  * Copyright (c) 2002 David Schultz <das (at) FreeBSD.ORG>
      4  * All rights reserved.
      5  *
      6  * Redistribution and use in source and binary forms, with or without
      7  * modification, are permitted provided that the following conditions
      8  * are met:
      9  * 1. Redistributions of source code must retain the above copyright
     10  *    notice, this list of conditions and the following disclaimer.
     11  * 2. Redistributions in binary form must reproduce the above copyright
     12  *    notice, this list of conditions and the following disclaimer in the
     13  *    documentation and/or other materials provided with the distribution.
     14  *
     15  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
     16  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
     17  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
     18  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
     19  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
     20  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
     21  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
     22  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
     23  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
     24  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
     25  * SUCH DAMAGE.
     26  *
     27  * $FreeBSD: src/lib/libc/include/fpmath.h,v 1.3 2005/02/06 03:23:31 das Exp $
     28  */
     29 
     30 #include <endian.h>
     31 #include "_fpmath.h"
     32 
     33 union IEEEf2bits {
     34 	float	f;
     35 	struct {
     36 #if __BYTE_ORDER == __LITTLE_ENDIAN
     37 		unsigned int	man	:23;
     38 		unsigned int	exp	:8;
     39 		unsigned int	sign	:1;
     40 #else /* _BIG_ENDIAN */
     41 		unsigned int	sign	:1;
     42 		unsigned int	exp	:8;
     43 		unsigned int	man	:23;
     44 #endif
     45 	} bits;
     46 };
     47 
     48 #define	DBL_MANH_SIZE	20
     49 #define	DBL_MANL_SIZE	32
     50 
     51 union IEEEd2bits {
     52 	double	d;
     53 	struct {
     54 /* #ifdef __ARMEB__ */
     55 #if (__BYTE_ORDER == __BIG_ENDIAN) || (defined(__arm__) && !defined(__VFP_FP__))
     56 		unsigned int	manh	:20;
     57 		unsigned int	exp	:11;
     58 		unsigned int	sign	:1;
     59 		unsigned int	manl	:32;
     60 #elif  __BYTE_ORDER == __LITTLE_ENDIAN
     61 		unsigned int	manl	:32;
     62 		unsigned int	manh	:20;
     63 		unsigned int	exp	:11;
     64 		unsigned int	sign	:1;
     65 #elif __BYTE_ORDER == __BIG_ENDIAN
     66 		unsigned int	sign	:1;
     67 		unsigned int	exp	:11;
     68 		unsigned int	manh	:20;
     69 		unsigned int	manl	:32;
     70 #endif
     71 	} bits;
     72 };
     73 
     74 /*
     75  * The BSD "long double" functions are broken when sizeof(long double) == sizeof(double).
     76  * Android works around those cases by replacing the broken functions with our own trivial stubs
     77  * that call the regular "double" function.
     78  */
     79 #define __fpclassifyl __broken__fpclassify
     80 #define __isfinitel __broken__isfinitel
     81 #define __isinfl __broken__isinfl
     82 #define __isnanl __broken__isnanl
     83 #define __isnormall __broken__isnormall
     84 #define __signbitl __broken_signbitl
     85