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

Subversion Repositories or1k

[/] [or1k/] [trunk/] [newlib/] [newlib/] [libm/] [math/] [wf_j0.c] - Blame information for rev 1774

Go to most recent revision | Details | Compare with Previous | View Log


/* wf_j0.c -- float version of w_j0.c.
 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
 */
 
/*
 * ====================================================
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
 *
 * Developed at SunPro, a Sun Microsystems, Inc. business.
 * Permission to use, copy, modify, and distribute this
 * software is freely granted, provided that this notice
 * is preserved.
 * ====================================================
 */
 
/*
 * wrapper j0f(float x), y0f(float x)
 */
 
#include "fdlibm.h"
#include <errno.h>
 
#ifdef __STDC__
        float j0f(float x)              /* wrapper j0f */
#else
        float j0f(x)                    /* wrapper j0f */
        float x;
#endif
{
#ifdef _IEEE_LIBM
        return __ieee754_j0f(x);
#else
        struct exception exc;
        float z = __ieee754_j0f(x);
        if(_LIB_VERSION == _IEEE_ || isnanf(x)) return z;
        if(fabsf(x)>(float)X_TLOSS) {
            /* j0f(|x|>X_TLOSS) */
            exc.type = TLOSS;
            exc.name = "j0f";
            exc.err = 0;
            exc.arg1 = exc.arg2 = (double)x;
            exc.retval = 0.0;
            if (_LIB_VERSION == _POSIX_)
               errno = ERANGE;
            else if (!matherr(&exc)) {
               errno = ERANGE;
            }
            if (exc.err != 0)
               errno = exc.err;
            return (float)exc.retval;
        } else
            return z;
#endif
}
 
#ifdef __STDC__
        float y0f(float x)              /* wrapper y0f */
#else
        float y0f(x)                    /* wrapper y0f */
        float x;
#endif
{
#ifdef _IEEE_LIBM
        return __ieee754_y0f(x);
#else
        float z;
        struct exception exc;
        z = __ieee754_y0f(x);
        if(_LIB_VERSION == _IEEE_ || isnanf(x) ) return z;
        if(x <= (float)0.0){
#ifndef HUGE_VAL 
#define HUGE_VAL inf
            double inf = 0.0;
 
            SET_HIGH_WORD(inf,0x7ff00000);      /* set inf to infinite */
#endif
            /* y0f(0) = -inf  or y0f(x<0) = NaN */
            exc.type = DOMAIN;  /* should be SING for IEEE y0f(0) */
            exc.name = "y0f";
            exc.err = 0;
            exc.arg1 = exc.arg2 = (double)x;
            if (_LIB_VERSION == _SVID_)
               exc.retval = -HUGE;
            else
               exc.retval = -HUGE_VAL;
            if (_LIB_VERSION == _POSIX_)
               errno = EDOM;
            else if (!matherr(&exc)) {
               errno = EDOM;
            }
            if (exc.err != 0)
               errno = exc.err;
            return (float)exc.retval;
        }
        if(x>(float)X_TLOSS) {
            /* y0f(x>X_TLOSS) */
            exc.type = TLOSS;
            exc.name = "y0f";
            exc.err = 0;
            exc.arg1 = exc.arg2 = (double)x;
            exc.retval = 0.0;
            if (_LIB_VERSION == _POSIX_)
                errno = ERANGE;
            else if (!matherr(&exc)) {
                errno = ERANGE;
            }
            if (exc.err != 0)
               errno = exc.err;
            return (float)exc.retval;
        } else
            return z;
#endif
}
 
#ifdef _DOUBLE_IS_32BITS
 
#ifdef __STDC__
        double j0(double x)
#else
        double j0(x)
        double x;
#endif
{
        return (double) j0f((float) x);
}
 
#ifdef __STDC__
        double y0(double x)
#else
        double y0(x)
        double x;
#endif
{
        return (double) y0f((float) x);
}
 
#endif /* defined(_DOUBLE_IS_32BITS) */

Line No.	Rev	Author	Line
1	39	lampret	`/* wf_j0.c -- float version of w_j0.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 j0f(float x), y0f(float x)`
18			`*/`
19
20			`#include "fdlibm.h"`
21			`#include <errno.h>`
22
23			`#ifdef __STDC__`
24			`float j0f(float x) /* wrapper j0f */`
25			`#else`
26			`float j0f(x) /* wrapper j0f */`
27			`float x;`
28			`#endif`
29			`{`
30			`#ifdef _IEEE_LIBM`
31			`return __ieee754_j0f(x);`
32			`#else`
33			`struct exception exc;`
34			`float z = __ieee754_j0f(x);`
35			`if(_LIB_VERSION == _IEEE_ \|\| isnanf(x)) return z;`
36			`if(fabsf(x)>(float)X_TLOSS) {`
37			`/* j0f(\|x\|>X_TLOSS) */`
38			`exc.type = TLOSS;`
39			`exc.name = "j0f";`
40	56	joel	`exc.err = 0;`
41			`exc.arg1 = exc.arg2 = (double)x;`
42	39	lampret	`exc.retval = 0.0;`
43			`if (_LIB_VERSION == _POSIX_)`
44			`errno = ERANGE;`
45			`else if (!matherr(&exc)) {`
46			`errno = ERANGE;`
47			`}`
48			`if (exc.err != 0)`
49			`errno = exc.err;`
50			`return (float)exc.retval;`
51			`} else`
52			`return z;`
53			`#endif`
54			`}`
55
56			`#ifdef __STDC__`
57			`float y0f(float x) /* wrapper y0f */`
58			`#else`
59			`float y0f(x) /* wrapper y0f */`
60			`float x;`
61			`#endif`
62			`{`
63			`#ifdef _IEEE_LIBM`
64			`return __ieee754_y0f(x);`
65			`#else`
66			`float z;`
67			`struct exception exc;`
68			`z = __ieee754_y0f(x);`
69			`if(_LIB_VERSION == _IEEE_ \|\| isnanf(x) ) return z;`
70			`if(x <= (float)0.0){`
71			`#ifndef HUGE_VAL`
72			`#define HUGE_VAL inf`
73			`double inf = 0.0;`
74
75			`SET_HIGH_WORD(inf,0x7ff00000); /* set inf to infinite */`
76			`#endif`
77			`/* y0f(0) = -inf or y0f(x<0) = NaN */`
78			`exc.type = DOMAIN; /* should be SING for IEEE y0f(0) */`
79			`exc.name = "y0f";`
80	56	joel	`exc.err = 0;`
81			`exc.arg1 = exc.arg2 = (double)x;`
82	39	lampret	`if (_LIB_VERSION == _SVID_)`
83			`exc.retval = -HUGE;`
84			`else`
85			`exc.retval = -HUGE_VAL;`
86			`if (_LIB_VERSION == _POSIX_)`
87			`errno = EDOM;`
88			`else if (!matherr(&exc)) {`
89			`errno = EDOM;`
90			`}`
91			`if (exc.err != 0)`
92			`errno = exc.err;`
93			`return (float)exc.retval;`
94			`}`
95			`if(x>(float)X_TLOSS) {`
96			`/* y0f(x>X_TLOSS) */`
97			`exc.type = TLOSS;`
98			`exc.name = "y0f";`
99	56	joel	`exc.err = 0;`
100			`exc.arg1 = exc.arg2 = (double)x;`
101	39	lampret	`exc.retval = 0.0;`
102			`if (_LIB_VERSION == _POSIX_)`
103			`errno = ERANGE;`
104			`else if (!matherr(&exc)) {`
105			`errno = ERANGE;`
106			`}`
107			`if (exc.err != 0)`
108			`errno = exc.err;`
109			`return (float)exc.retval;`
110			`} else`
111			`return z;`
112			`#endif`
113			`}`
114
115			`#ifdef _DOUBLE_IS_32BITS`
116
117			`#ifdef __STDC__`
118			`double j0(double x)`
119			`#else`
120			`double j0(x)`
121			`double x;`
122			`#endif`
123			`{`
124			`return (double) j0f((float) x);`
125			`}`
126
127			`#ifdef __STDC__`
128			`double y0(double x)`
129			`#else`
130			`double y0(x)`
131			`double x;`
132			`#endif`
133			`{`
134			`return (double) y0f((float) x);`
135			`}`
136
137			`#endif /* defined(_DOUBLE_IS_32BITS) */`

Browse

Tools

Subversion Repositories or1k

[/] [or1k/] [trunk/] [newlib/] [newlib/] [libm/] [math/] [wf_j0.c] - Blame information for rev 1774