Return to madd.c CVS log

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

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

1.1     ! root        1: #include "mp.h"
        !             2: 
        !             3: mint
        !             4: madd(a,b)
        !             5: mint a,b;
        !             6: {
        !             7:        mint x;
        !             8:        int i;
        !             9: 
        !            10:        x.low = a.low + b.low;
        !            11:        x.high = a.high + b.high;
        !            12:        if(x.low > e16){
        !            13:                x.low -= e16;
        !            14:                x.high += 1.;
        !            15:        }else
        !            16:        if(x.low < -e16){
        !            17:                x.low += e16;
        !            18:                x.high -= 1.;
        !            19:        }
        !            20:        if(x.high >= e16) ofl("madd");
        !            21:        if((x.high>0) && (x.low<0)){
        !            22:                x.high -= 1.;
        !            23:                x.low += e16;
        !            24:        }
        !            25:        if((x.high<0)  && (x.low>0)){
        !            26:                x.high += 1.;
        !            27:                x.low -= e16;
        !            28:        }
        !            29:        return(x);
        !            30: }
        !            31: 
        !            32: mint
        !            33: msub(a,b)
        !            34: mint a,b;
        !            35: {
        !            36:        b.high = -b.high;
        !            37:        b.low = -b.low;
        !            38:        return(madd(a,b));
        !            39: }
        !            40: 
        !            41: int
        !            42: mcmp(a,b)
        !            43: mint a, b;
        !            44: {
        !            45:        if(a.high > b.high)
        !            46:                return(1);
        !            47:        if(a.high < b.high)
        !            48:                return(-1);
        !            49:        if(a.low > b.low)
        !            50:                return(1);
        !            51:        if(a.low < b.low)
        !            52:                return(-1);
        !            53:        return(0);
        !            54: }
        !            55: 
        !            56: 
        !            57: mint
        !            58: itom(i)
        !            59: int i;
        !            60: {
        !            61:        mint x;
        !            62: 
        !            63:        x.high = 0.;
        !            64:        x.low = i;
        !            65:        return(x);
        !            66: }
        !            67: mint
        !            68: mchs(a)
        !            69: mint a;
        !            70: {
        !            71:        a.high = -a.high;
        !            72:        a.low = -a.low;
        !            73:        return(a);
        !            74: }

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