1 /* $NetBSD: mathimpl.h,v 1.10 2011/11/02 02:34:56 christos Exp $ */ 2 /* 3 * Copyright (c) 1988, 1993 4 * The Regents of the University of California. 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 * 3. Neither the name of the University nor the names of its contributors 15 * may be used to endorse or promote products derived from this software 16 * without specific prior written permission. 17 * 18 * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND 19 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 20 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 21 * ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE 22 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 23 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 24 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 25 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 26 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 27 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 28 * SUCH DAMAGE. 29 * 30 * @(#)mathimpl.h 8.1 (Berkeley) 6/4/93 31 */ 32 #ifndef _NOIEEE_SRC_MATHIMPL_H_ 33 #define _NOIEEE_SRC_MATHIMPL_H_ 34 35 #include <sys/cdefs.h> 36 #include <math.h> 37 #include <stdint.h> 38 39 #if defined(__vax__) || defined(tahoe) 40 41 /* Deal with different ways to concatenate in cpp */ 42 #define cat3(a,b,c) a ## b ## c 43 44 /* Deal with vax/tahoe byte order issues */ 45 # ifdef __vax__ 46 # define cat3t(a,b,c) cat3(a,b,c) 47 # else 48 # define cat3t(a,b,c) cat3(a,c,b) 49 # endif 50 51 # define vccast(name) (cat3(__,name,x).d) 52 53 /* 54 * Define a constant to high precision on a Vax or Tahoe. 55 * 56 * Args are the name to define, the decimal floating point value, 57 * four 16-bit chunks of the float value in hex 58 * (because the vax and tahoe differ in float format!), the power 59 * of 2 of the hex-float exponent, and the hex-float mantissa. 60 * Most of these arguments are not used at compile time; they are 61 * used in a post-check to make sure the constants were compiled 62 * correctly. 63 * 64 * People who want to use the constant will have to do their own 65 * #define foo vccast(foo) 66 * since CPP cannot do this for them from inside another macro (sigh). 67 * We define "vccast" if this needs doing. 68 */ 69 #ifdef _LIBM_DECLARE 70 # define vc(name, value, x1,x2,x3,x4, bexp, xval) \ 71 const union { uint32_t l[2]; double d; } cat3(__,name,x) = { \ 72 .l = { [0] = cat3t(0x,x1,x2), [1] = cat3t(0x,x3,x4) } }; 73 #elif defined(_LIBM_STATIC) 74 # define vc(name, value, x1,x2,x3,x4, bexp, xval) \ 75 static const union { uint32_t l[2]; double d; } cat3(__,name,x) = { \ 76 .l = { [0] = cat3t(0x,x1,x2), [1] = cat3t(0x,x3,x4) } }; 77 #else 78 # define vc(name, value, x1,x2,x3,x4, bexp, xval) \ 79 extern const union { uint32_t l[2]; double d; } cat3(__,name,x); 80 #endif 81 # define ic(name, value, bexp, xval) 82 83 #else /* __vax__ or tahoe */ 84 85 /* Hooray, we have an IEEE machine */ 86 # undef vccast 87 # define vc(name, value, x1,x2,x3,x4, bexp, xval) 88 89 #ifdef _LIBM_DECLARE 90 # define ic(name, value, bexp, xval) \ 91 const double __CONCAT(__,name) = value; 92 #elif _LIBM_STATIC 93 # define ic(name, value, bexp, xval) \ 94 static const double __CONCAT(__,name) = value; 95 #else 96 # define ic(name, value, bexp, xval) \ 97 extern const double __CONCAT(__,name); 98 #endif 99 100 #endif /* defined(__vax__)||defined(tahoe) */ 101 102 #ifdef __vax__ 103 #include <machine/float.h> 104 #define _TINY DBL_EPSILON 105 #define _TINYER DBL_EPSILON 106 #define _HUGE DBL_MAX 107 #else 108 #define _TINY 1e-300 109 #define _TINYER 1e-308 110 #define _HUGE 1e+300 111 #endif 112 113 114 /* 115 * Functions internal to the math package, yet not static. 116 */ 117 extern double __exp__E(double, double); 118 extern double __log__L(double); 119 extern int infnan(int); 120 121 struct Double {double a, b;}; 122 double __exp__D(double, double); 123 struct Double __log__D(double); 124 125 #endif /* _NOIEEE_SRC_MATHIMPL_H_ */ 126