Return to spow.c CVS log

Up to [MW Coherent from dump] / coherent / d / usr / lib / libmp

Annotation of coherent/d/usr/lib/libmp/spow.c, revision 1.1

1.1     ! root        1: #include "mprec.h"
        !             2: 
        !             3: 
        !             4: /*
        !             5:  *     Spow sets "b" to "a" raised to the "n" power.  Note that
        !             6:  *     a = b is acceptable.
        !             7:  */
        !             8: 
        !             9: void
        !            10: spow(a, n, b)
        !            11: register mint  *a, *b;
        !            12: register unsigned n;
        !            13: {
        !            14:        mint    ac;
        !            15: 
        !            16:        minit(&ac);
        !            17:        mcopy(a, &ac);
        !            18:        mcopy(mone, b);
        !            19:        if (n != 0) {
        !            20:                for (; n != 1; n >>= 1, mult(&ac, &ac, &ac))
        !            21:                        if (n % 2 != 0)
        !            22:                                mult(&ac, b, b);
        !            23:                mult(&ac, b, b);
        !            24:        }
        !            25:        mpfree(ac.val);
        !            26: }

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