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

Subversion Repositories or1k

[/] [or1k/] [trunk/] [newlib/] [newlib/] [libm/] [math/] [wf_scalb.c] - Blame information for rev 1775

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

Line No. Rev Author Line
1 39 lampret 
/* wf_scalb.c -- float version of w_scalb.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 scalbf(float x, float fn) is provide for
18 
 * passing various standard test suite. One
19 
 * should use scalbn() instead.
20 
 */
21 
 
22 
#include "fdlibm.h"
23 
#include <errno.h>
24 
 
25 
#ifdef __STDC__
26 
#ifdef _SCALB_INT
27 
        float scalbf(float x, int fn)           /* wrapper scalbf */
28 
#else
29 
        float scalbf(float x, float fn)         /* wrapper scalbf */
30 
#endif
31 
#else
32 
        float scalbf(x,fn)                      /* wrapper scalbf */
33 
#ifdef _SCALB_INT
34 
        float x; int fn;
35 
#else
36 
        float x,fn;
37 
#endif
38 
#endif
39 
{
40 
#ifdef _IEEE_LIBM
41 
        return __ieee754_scalbf(x,fn);
42 
#else
43 
        float z;
44 
#ifndef HUGE_VAL
45 
#define HUGE_VAL inf
46 
        double inf = 0.0;
47 
 
48 
        SET_HIGH_WORD(inf,0x7ff00000);  /* set inf to infinite */
49 
#endif
50 
        struct exception exc;
51 
        z = __ieee754_scalbf(x,fn);
52 
        if(_LIB_VERSION == _IEEE_) return z;
53 
        if(!(finitef(z)||isnanf(z))&&finitef(x)) {
54 
            /* scalbf overflow; SVID also returns +-HUGE_VAL */
55 
            exc.type = OVERFLOW;
56 
            exc.name = "scalbf";
57 56 joel 
            exc.err = 0;
58 
            exc.arg1 = (double)x;
59 
            exc.arg2 = (double)fn;
60 39 lampret 
            exc.retval = x > 0.0 ? HUGE_VAL : -HUGE_VAL;
61 
            if (_LIB_VERSION == _POSIX_)
62 
               errno = ERANGE;
63 
            else if (!matherr(&exc)) {
64 
               errno = ERANGE;
65 
            }
66 
            if (exc.err != 0)
67 
               errno = exc.err;
68 
            return exc.retval;
69 
        }
70 
        if(z==(float)0.0&&z!=x) {
71 
            /* scalbf underflow */
72 
            exc.type = UNDERFLOW;
73 
            exc.name = "scalbf";
74 56 joel 
            exc.err = 0;
75 
            exc.arg1 = (double)x;
76 
            exc.arg2 = (double)fn;
77 39 lampret 
            exc.retval = copysign(0.0,x);
78 
            if (_LIB_VERSION == _POSIX_)
79 
               errno = ERANGE;
80 
            else if (!matherr(&exc)) {
81 
               errno = ERANGE;
82 
            }
83 
            if (exc.err != 0)
84 
               errno = exc.err;
85 
            return exc.retval;
86 
        }
87 
#ifndef _SCALB_INT
88 
        if(!finitef(fn)) errno = ERANGE;
89 
#endif
90 
        return z;
91 
#endif
92 
}
93 
 
94 
#ifdef _DOUBLE_IS_32BITS
95 
 
96 
#ifdef __STDC__
97 
#ifdef _SCALB_INT
98 
        double scalb(double x, int fn)
99 
#else
100 
        double scalb(double x, double fn)
101 
#endif
102 
#else
103 
        double scalb(x, fn)
104 
#ifdef _SCALB_INT
105 
        double x; int fn;
106 
#else
107 
        double x,fn;
108 
#endif
109 
#endif
110 
{
111 
#ifdef _SCALB_INT
112 
        return (double) scalbf((float) x, fn);
113 
#else
114 
        return (double) scalbf((float) x, (float) fn);
115 
#endif
116 
}
117 
 
118 
#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.