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

Subversion Repositories or1k

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

Details | Compare with Previous | View Log

Line No. Rev Author Line
1 39 lampret
/* wf_gamma.c -- float version of w_gamma.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
#include "fdlibm.h"
18 56 joel
#include <reent.h>
19 39 lampret
#include <errno.h>
20
 
21
#ifdef __STDC__
22
        float gammaf(float x)
23
#else
24
        float gammaf(x)
25
        float x;
26
#endif
27
{
28
#ifdef _IEEE_LIBM
29 56 joel
        return __ieee754_gammaf_r(x,&(_REENT->_new._reent._gamma_signgam));
30 39 lampret
#else
31
        float y;
32
        struct exception exc;
33 56 joel
        y = __ieee754_gammaf_r(x,&(_REENT->_new._reent._gamma_signgam));
34 39 lampret
        if(_LIB_VERSION == _IEEE_) return y;
35
        if(!finitef(y)&&finitef(x)) {
36
#ifndef HUGE_VAL 
37
#define HUGE_VAL inf
38
            double inf = 0.0;
39
 
40
            SET_HIGH_WORD(inf,0x7ff00000);      /* set inf to infinite */
41
#endif
42
            if(floorf(x)==x&&x<=(float)0.0) {
43
                /* gammaf(-integer) or gammaf(0) */
44
                exc.type = SING;
45
                exc.name = "gammaf";
46 56 joel
                exc.err = 0;
47
                exc.arg1 = exc.arg2 = (double)x;
48 39 lampret
                if (_LIB_VERSION == _SVID_)
49
                  exc.retval = HUGE;
50
                else
51
                  exc.retval = HUGE_VAL;
52
                if (_LIB_VERSION == _POSIX_)
53
                  errno = EDOM;
54
                else if (!matherr(&exc)) {
55
                  errno = EDOM;
56
                }
57
            } else {
58
                /* gammaf(finite) overflow */
59
                exc.type = OVERFLOW;
60
                exc.name = "gammaf";
61 56 joel
                exc.err = 0;
62
                exc.arg1 = exc.arg2 = (double)x;
63 39 lampret
                if (_LIB_VERSION == _SVID_)
64
                  exc.retval = HUGE;
65
                else
66
                  exc.retval = HUGE_VAL;
67
                if (_LIB_VERSION == _POSIX_)
68
                  errno = ERANGE;
69
                else if (!matherr(&exc)) {
70
                  errno = ERANGE;
71
                }
72
            }
73
            if (exc.err != 0)
74
               errno = exc.err;
75
            return (float)exc.retval;
76
        } else
77
            return y;
78
#endif
79
}
80
 
81
#ifdef _DOUBLE_IS_32BITS
82
 
83
#ifdef __STDC__
84
        double gamma(double x)
85
#else
86
        double gamma(x)
87
        double x;
88
#endif
89
{
90
        return (double) gammaf((float) x);
91
}
92
 
93
#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.