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