Return to pow_ci.c CVS log

Up to [CSRG BSD Unix] / 43BSDTahoe / usr.lib / libF77

Annotation of 43BSDTahoe/usr.lib/libF77/pow_ci.c, revision 1.1

1.1     ! root        1: /*
        !             2:  * Copyright (c) 1980 Regents of the University of California.
        !             3:  * All rights reserved.  The Berkeley software License Agreement
        !             4:  * specifies the terms and conditions for redistribution.
        !             5:  *
        !             6:  *     @(#)pow_ci.c    5.3     1/24/88
        !             7:  */
        !             8: 
        !             9: #include "complex"
        !            10: 
        !            11: #ifdef tahoe
        !            12: 
        !            13: #define        C_MULEQ(A,B)    \
        !            14:        t = (A).real * (B).real - (A).imag * (B).imag,\
        !            15:        (A).imag = (A).real * (B).imag + (A).imag * (B).real,\
        !            16:        (A).real = t    /* A *= B */
        !            17: 
        !            18: void
        !            19: pow_ci(p, a, b)        /* p = a**b  */
        !            20:        complex *p, *a;
        !            21:        long *b;
        !            22: {
        !            23:        register long n = *b;
        !            24:        register float t;
        !            25:        complex x;
        !            26: 
        !            27:        x = *a;
        !            28:        p->real = (float)1, p->imag = (float)0;
        !            29:        if (!n)
        !            30:                return;
        !            31:        if (n < 0) {
        !            32:                c_div(&x, p, a);
        !            33:                n = -n;
        !            34:        }
        !            35:        while (!(n&1)) {
        !            36:                C_MULEQ(x, x);
        !            37:                n >>= 1;
        !            38:        }
        !            39:        for (*p = x; --n > 0; C_MULEQ(*p, x))
        !            40:                while (!(n&1)) {
        !            41:                        C_MULEQ(x, x);
        !            42:                        n >>= 1;
        !            43:                }
        !            44: }
        !            45: 
        !            46: #else /* !tahoe */
        !            47: 
        !            48: extern void pow_zi();
        !            49: 
        !            50: void
        !            51: pow_ci(p, a, b)        /* p = a**b  */
        !            52:        complex *p, *a;
        !            53:        long *b;
        !            54: {
        !            55:        dcomplex p1, a1;
        !            56: 
        !            57:        a1.dreal = a->real;
        !            58:        a1.dimag = a->imag;
        !            59: 
        !            60:        pow_zi(&p1, &a1, b);
        !            61: 
        !            62:        p->real = p1.dreal;
        !            63:        p->imag = p1.dimag;
        !            64: }
        !            65: 
        !            66: #endif /* tahoe */

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