4d77c1151c
* reworked amd64/_fpmath.h and arm64/_fpmath.h to support 128-bit long doubles. * improved tests to cover long double cases * made modfl an alias for LP32 Tests pass on x86, x86_64, arm, arm64 and mips. Bug: 12921273 Change-Id: Ibe39acde57972447a8950fa45b1e702acc68ebeb
93 lines
2.9 KiB
C
93 lines
2.9 KiB
C
/*-
|
|
* Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org>
|
|
* Copyright (c) 2002 David Schultz <das@FreeBSD.ORG>
|
|
* All rights reserved.
|
|
*
|
|
* Redistribution and use in source and binary forms, with or without
|
|
* modification, are permitted provided that the following conditions
|
|
* are met:
|
|
* 1. Redistributions of source code must retain the above copyright
|
|
* notice, this list of conditions and the following disclaimer.
|
|
* 2. Redistributions in binary form must reproduce the above copyright
|
|
* notice, this list of conditions and the following disclaimer in the
|
|
* documentation and/or other materials provided with the distribution.
|
|
*
|
|
* THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
|
|
* ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
|
|
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
|
|
* ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
|
|
* FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
|
|
* DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
|
|
* OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
|
|
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
|
|
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
|
|
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
|
|
* SUCH DAMAGE.
|
|
*
|
|
* $FreeBSD: src/lib/libc/include/fpmath.h,v 1.3 2005/02/06 03:23:31 das Exp $
|
|
*/
|
|
|
|
#ifndef _FPMATH_
|
|
#define _FPMATH_
|
|
|
|
#include <endian.h>
|
|
#include "_fpmath.h"
|
|
|
|
union IEEEf2bits {
|
|
float f;
|
|
struct {
|
|
#if __BYTE_ORDER == __LITTLE_ENDIAN
|
|
unsigned int man :23;
|
|
unsigned int exp :8;
|
|
unsigned int sign :1;
|
|
#else /* _BIG_ENDIAN */
|
|
unsigned int sign :1;
|
|
unsigned int exp :8;
|
|
unsigned int man :23;
|
|
#endif
|
|
} bits;
|
|
};
|
|
|
|
#define DBL_MANH_SIZE 20
|
|
#define DBL_MANL_SIZE 32
|
|
|
|
union IEEEd2bits {
|
|
double d;
|
|
struct {
|
|
/* #ifdef __ARMEB__ */
|
|
#if (__BYTE_ORDER == __BIG_ENDIAN) || (defined(__arm__) && !defined(__VFP_FP__))
|
|
unsigned int manh :20;
|
|
unsigned int exp :11;
|
|
unsigned int sign :1;
|
|
unsigned int manl :32;
|
|
#elif __BYTE_ORDER == __LITTLE_ENDIAN
|
|
unsigned int manl :32;
|
|
unsigned int manh :20;
|
|
unsigned int exp :11;
|
|
unsigned int sign :1;
|
|
#elif __BYTE_ORDER == __BIG_ENDIAN
|
|
unsigned int sign :1;
|
|
unsigned int exp :11;
|
|
unsigned int manh :20;
|
|
unsigned int manl :32;
|
|
#endif
|
|
} bits;
|
|
};
|
|
|
|
/*
|
|
* The BSD "long double" functions are broken when sizeof(long double) == sizeof(double).
|
|
* Android works around those cases by replacing the broken functions with our own trivial stubs
|
|
* that call the regular "double" function.
|
|
*/
|
|
#ifndef __LP64__
|
|
|
|
#define __fpclassifyl __broken__fpclassify
|
|
#define __isfinitel __broken__isfinitel
|
|
#define __isinfl __broken__isinfl
|
|
#define __isnanl __broken__isnanl
|
|
#define __isnormall __broken__isnormall
|
|
#define __signbitl __broken_signbitl
|
|
|
|
#endif // __LP64__
|
|
|
|
#endif // _FPMATH_
|