Return to hypot.c-1 CVS log

Up to [Research Unix] / researchv10dc / libc / math

Annotation of researchv10dc/libc/math/hypot.c-1, revision 1.1

1.1     ! root        1: /*
        !             2:  * sqrt(a^2 + b^2)
        !             3:  *     (but carefully)
        !             4:  */
        !             5: 
        !             6: double sqrt();
        !             7: double
        !             8: hypot(a,b)
        !             9: double a,b;
        !            10: {
        !            11:        double t;
        !            12:        if(a<0) a = -a;
        !            13:        if(b<0) b = -b;
        !            14:        if(a > b) {
        !            15:                t = a;
        !            16:                a = b;
        !            17:                b = t;
        !            18:        }
        !            19:        if(b==0) return(0.);
        !            20:        a /= b;
        !            21:        /*
        !            22:         * pathological overflow possible
        !            23:         * in the next line.
        !            24:         */
        !            25:        return(b*sqrt(1. + a*a));
        !            26: }
        !            27: 
        !            28: struct complex
        !            29: {
        !            30:        double  r;
        !            31:        double  i;
        !            32: };
        !            33: 
        !            34: double
        !            35: cabs(arg)
        !            36: struct complex arg;
        !            37: {
        !            38:        double hypot();
        !            39: 
        !            40:        return(hypot(arg.r, arg.i));
        !            41: }

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