Return to mneg.c CVS log

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

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

1.1     ! root        1: #include "mprec.h"
        !             2: #include <assert.h>
        !             3: 
        !             4: 
        !             5: /*
        !             6:  *     Mneg sets the mint pointed to by "b" to negative 1 times the
        !             7:  *     mint pointed to by "a".  Note that "a" == "b" is permissable.
        !             8:  */
        !             9: 
        !            10: void
        !            11: mneg(a, b)
        !            12: mint *a, *b;
        !            13: {
        !            14:        register char *ap, *rp;
        !            15:        register unsigned count;
        !            16:        mint res;
        !            17:        int mifl;
        !            18: 
        !            19:        /* allocate result space */
        !            20:        mifl = ispos(a);
        !            21:        res.len = a->len;
        !            22:        if (mifl)
        !            23:                ++res.len;
        !            24:        res.val = (char *)mpalc(res.len);
        !            25: 
        !            26:        /* negate and copy */
        !            27:        ap = a->val;
        !            28:        rp = res.val;
        !            29:        count = a->len;
        !            30:        while (count-- > 0)
        !            31:                *rp++ = NEFL - *ap++;
        !            32:        if (mifl)
        !            33:                *rp = NEFL;
        !            34:        ++*res.val;
        !            35:        assert(ap == a->val + a->len);
        !            36:        assert(rp == res.val + res.len - (mifl ? 1 : 0));
        !            37:        norm(&res);
        !            38: 
        !            39:        /* replace old value of b with res */
        !            40:        mpfree(b->val);
        !            41:        *b = res;
        !            42: }

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