Return to pow.c CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libm

Annotation of coherent/b/lib/libm/pow.c, revision 1.1

1.1     ! root        1: /*
        !             2:  * libm/pow.c
        !             3:  * Raise x to the power y.
        !             4:  */
        !             5: 
        !             6: #include <math.h>
        !             7: 
        !             8: #if    EMU87
        !             9: #include "emumath.h"
        !            10: #endif
        !            11: 
        !            12: double
        !            13: pow(x, y)
        !            14: double x, y;
        !            15: {
        !            16:        double r;
        !            17:        register unsigned s, i, e;
        !            18: 
        !            19:        s = 0;
        !            20:        i = 0;
        !            21:        if (x == 0.0) {
        !            22:                if (y <= 0.0)
        !            23:                        errno = EDOM;
        !            24:                return 0.0;
        !            25:        }
        !            26:        if ((r = modf(y, &r)) < 0.0)
        !            27:                r += 1.0;
        !            28:        if (x < 0.0) {
        !            29:                if (r != 0.0) {
        !            30:                        errno = EDOM;
        !            31:                        return 0.0;
        !            32:                }
        !            33:                x = -x;
        !            34:                if (((int) y) & 1)
        !            35:                        s = 1;
        !            36:        }
        !            37:        if (y < 0.0) {
        !            38:                y = -y;
        !            39:                i = 1;
        !            40:        }
        !            41:        if (r!=0.0 || y>16384.0)
        !            42:                r = _two(y*log10(x)*LOG10B2);
        !            43:        else {
        !            44:                r = 1.0;
        !            45:                for (e=y; e; e>>=1) {
        !            46:                        if (e&01)
        !            47:                                r *= x;
        !            48:                        x *= x;
        !            49:                }
        !            50:        }
        !            51:        if (i)
        !            52:                r = 1/r;
        !            53:        return s ? -r : r;
        !            54: }
        !            55: 
        !            56: /* end of libm/pow.c */

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.