Home | History | Annotate | Line # | Download | only in noieee_src
mathimpl.h revision 1.6
      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