Return to mult.c CVS log

Up to [Research Unix] / researchv10no / cmd / lfactor

Annotation of researchv10no/cmd/lfactor/mult.c, revision 1.1

1.1     ! root        1: #include "mp.h"
        !             2: 
        !             3: mint
        !             4: mult(a,b)
        !             5: mint a, b;
        !             6: {
        !             7:        mint z;
        !             8:        mint c;
        !             9:        mint mtemp;
        !            10:        int sign;
        !            11:        double a3, a2, a1, a0;
        !            12:        double b1, b0;
        !            13:        double temp1, temp2;
        !            14: 
        !            15:        if((a.high==0) && (b.high==0)){
        !            16:                c.high = 0.;
        !            17:                c.low = a.low * b.low;
        !            18:                if((c.low<e16) && (c.low > -e16)){
        !            19:                        return(c);
        !            20:                }
        !            21:        }
        !            22: 
        !            23:        sign = 1;
        !            24: 
        !            25:        if((a.high<0) || (a.low<0)){
        !            26:                sign = -sign;
        !            27:                a.high = -a.high;
        !            28:                a.low = -a.low;
        !            29:        }
        !            30:        if((b.high<0) || (b.low<0)){
        !            31:                sign = -sign;
        !            32:                b.high = -b.high;
        !            33:                b.low = -b.low;
        !            34:        }
        !            35:        if(b.high != 0){
        !            36:                mtemp = a;
        !            37:                a = b;
        !            38:                b = mtemp;
        !            39:        }
        !            40:        if(b.high!=0) ofl("mult");
        !            41: 
        !            42:        modf(a.low/1e8 , &a1);
        !            43:        a0 = a.low - a1*1e8;
        !            44:        modf(b.low/1e8 , &b1);
        !            45:        b0 = b.low - b1*1e8;
        !            46: 
        !            47:        modf(a.high/1e8 , &a3);
        !            48:        a2 = a.high - a3*1e8;
        !            49: 
        !            50:        temp1 = a0*b1 + a1*b0;
        !            51:        modf(temp1/1e8 , &temp2);
        !            52:        z.low = a0*b0 + (temp1 - temp2*1e8)*1e8;
        !            53:        z.high = a1*b1 + a2*b0 + (a2*b1 + a3*b0)*1e8 + temp2;
        !            54:        while(z.low >= e16){
        !            55:                z.low = z.low - e16;
        !            56:                z.high = z.high + 1.;
        !            57:        }
        !            58:        if(z.high >= e16)
        !            59:                ofl("mult");
        !            60:        if(a3*b1!= 0)
        !            61:                ofl("mult");
        !            62: 
        !            63:        if(sign<0){
        !            64:                z.high = -z.high;
        !            65:                z.low = -z.low;
        !            66:        }
        !            67:        return(z);
        !            68: }
        !            69: 
        !            70: mint
        !            71: smult(a,b)
        !            72: mint a;
        !            73: int b;
        !            74: {
        !            75:        return(mult(a, itom(b)));
        !            76: }
        !            77: ofl(s)
        !            78: char *s;
        !            79: {
        !            80:        printf("%s: overflow\n",s);
        !            81:        abort();
        !            82: }

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