Return to mult.c CVS log

Up to [Research Unix] / researchv10no / libmp

Annotation of researchv10no/libmp/mult.c, revision 1.1

1.1     ! root        1: #include "mp.h"
        !             2: #undef unfast
        !             3: mult(a,b,c) mint *a,*b,*c;
        !             4: {      mint x,y,z;
        !             5:        int sign;
        !             6:        sign = 1;
        !             7:        x.val=a->val;
        !             8:        y.val=b->val;
        !             9:        z.len=0;
        !            10:        if(a->len<0)
        !            11:        {       x.len= -a->len;
        !            12:                sign= -sign;
        !            13:        }
        !            14:        else    x.len=a->len;
        !            15:        if(b->len<0)
        !            16:        {       y.len= -b->len;
        !            17:                sign= -sign;
        !            18:        }
        !            19:        else    y.len=b->len;
        !            20:        if(x.len == 0 || y.len == 0)
        !            21:                z.len = 0;
        !            22:        else if(x.len<y.len) m_mult(&y,&x,&z);
        !            23:        else m_mult(&x,&y,&z);
        !            24:        xfree(c);
        !            25:        if(sign<0) c->len= -z.len;
        !            26:        else c->len=z.len;
        !            27:        c->val=z.val;
        !            28: }
        !            29: #define S2 x=a->val[j];
        !            30: #define S3 x=x*b->val[i-j];
        !            31: #ifdef unfast
        !            32: #define S4 tradd(&carry,&sum,x);
        !            33: #else
        !            34: #define S4 sum.xx += x; if(sum.yy.high & 0100000) {sum.yy.high &= 077777; carry++;}
        !            35: #endif
        !            36: #define S5 c->val[i]=sum.yy.low&077777;
        !            37: #define S6 sum.xx=sum.xx>>15;
        !            38: #define S7 sum.yy.high=carry;
        !            39: m_mult(a,b,c)
        !            40: register mint *a,*b,*c;
        !            41: {      register long x;
        !            42:        union {long xx; struct half yy;} sum;
        !            43:        int carry;
        !            44:        int i,j;
        !            45:        c->val=xalloc(a->len+b->len,"m_mult");
        !            46:        sum.xx=0;
        !            47:        for(i=0;i<b->len;i++)
        !            48:        {       carry=0;
        !            49:                for(j=0;j<i+1;j++)
        !            50:                {       S2
        !            51:                        S3
        !            52:                        S4
        !            53:                }
        !            54:                S5
        !            55:                S6
        !            56:                S7
        !            57:        }
        !            58:        for(;i<a->len;i++)
        !            59:        {       carry=0;
        !            60:                for(j=i-b->len+1;j<i+1;j++)
        !            61:                {       S2
        !            62:                        S3
        !            63:                        S4
        !            64:                }
        !            65:                S5
        !            66:                S6
        !            67:                S7
        !            68:        }
        !            69:        for(;i<a->len+b->len;i++)
        !            70:        {       carry=0;
        !            71:                for(j=i-b->len+1;j<a->len;j++)
        !            72:                {       S2
        !            73:                        S3
        !            74:                        S4
        !            75:                }
        !            76:                S5
        !            77:                S6
        !            78:                S7
        !            79:        }
        !            80:        if(i == 0 || c->val[i-1]!=0)
        !            81:                c->len=a->len+b->len;
        !            82:        else    c->len=a->len+b->len-1;
        !            83:        return;
        !            84: }
        !            85: tradd(a,b,c)
        !            86: long c;
        !            87: int *a;
        !            88: register union {long xx; struct half yy;} *b;
        !            89: {
        !            90:        b->xx += c;
        !            91:        if(b->yy.high&0100000)
        !            92:        {       b->yy.high= b->yy.high&077777;
        !            93:                *a += 1;
        !            94:        }
        !            95:        return;
        !            96: }

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