OpenCores
URL https://opencores.org/ocsvn/or1k/or1k/trunk

Subversion Repositories or1k

[/] [or1k/] [trunk/] [newlib/] [newlib/] [libm/] [math/] [wf_exp.c] - Blame information for rev 1765

Details | Compare with Previous | View Log

Line No. Rev Author Line
1 39 lampret
/* wf_exp.c -- float version of w_exp.c.
2
 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3
 */
4
 
5
/*
6
 * ====================================================
7
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8
 *
9
 * Developed at SunPro, a Sun Microsystems, Inc. business.
10
 * Permission to use, copy, modify, and distribute this
11
 * software is freely granted, provided that this notice
12
 * is preserved.
13
 * ====================================================
14
 */
15
 
16
/*
17
 * wrapper expf(x)
18
 */
19
 
20
#include "fdlibm.h"
21
#include <errno.h>
22
 
23
#ifdef __STDC__
24
static const float
25
#else
26
static float
27
#endif
28
o_threshold=  8.8721679688e+01,  /* 0x42b17180 */
29
u_threshold= -1.0397208405e+02;  /* 0xc2cff1b5 */
30
 
31
#ifdef __STDC__
32
        float expf(float x)             /* wrapper expf */
33
#else
34
        float expf(x)                   /* wrapper expf */
35
        float x;
36
#endif
37
{
38
#ifdef _IEEE_LIBM
39
        return __ieee754_expf(x);
40
#else
41
        float z;
42
        struct exception exc;
43
        z = __ieee754_expf(x);
44
        if(_LIB_VERSION == _IEEE_) return z;
45
        if(finitef(x)) {
46
            if(x>o_threshold) {
47
                /* expf(finite) overflow */
48
#ifndef HUGE_VAL
49
#define HUGE_VAL inf
50
                double inf = 0.0;
51
 
52
                SET_HIGH_WORD(inf,0x7ff00000);  /* set inf to infinite */
53
#endif
54
                exc.type = OVERFLOW;
55
                exc.name = "expf";
56 56 joel
                exc.err = 0;
57
                exc.arg1 = exc.arg2 = (double)x;
58 39 lampret
                if (_LIB_VERSION == _SVID_)
59
                  exc.retval = HUGE;
60
                else
61
                  exc.retval = HUGE_VAL;
62
                if (_LIB_VERSION == _POSIX_)
63
                  errno = ERANGE;
64
                else if (!matherr(&exc)) {
65
                        errno = ERANGE;
66
                }
67
                if (exc.err != 0)
68
                   errno = exc.err;
69
                return exc.retval;
70
            } else if(x<u_threshold) {
71
                /* expf(finite) underflow */
72
                exc.type = UNDERFLOW;
73
                exc.name = "expf";
74 56 joel
                exc.err = 0;
75
                exc.arg1 = exc.arg2 = (double)x;
76 39 lampret
                exc.retval = 0.0;
77
                if (_LIB_VERSION == _POSIX_)
78
                  errno = ERANGE;
79
                else if (!matherr(&exc)) {
80
                        errno = ERANGE;
81
                }
82
                if (exc.err != 0)
83
                   errno = exc.err;
84
                return exc.retval;
85
            }
86
        }
87
        return z;
88
#endif
89
}
90
 
91
#ifdef _DOUBLE_IS_32BITS
92
 
93
#ifdef __STDC__
94
        double exp(double x)
95
#else
96
        double exp(x)
97
        double x;
98
#endif
99
{
100
        return (double) expf((float) x);
101
}
102
 
103
#endif /* defined(_DOUBLE_IS_32BITS) */

powered by: WebSVN 2.1.0

© copyright 1999-2024 OpenCores.org, equivalent to Oliscience, all rights reserved. OpenCores®, registered trademark.