Return to norm.c CVS log

Up to [Research Unix] / researchv10no / libj

Annotation of researchv10no/libj/norm.c, revision 1.1

1.1     ! root        1: norm(x,y,z)
        !             2: {
        !             3:        return (sqrt(x*x + y*y + z*z));
        !             4: }
        !             5: 
        !             6: sqrtryz(x,y,z)
        !             7: {
        !             8:        register long sumsq;
        !             9: 
        !            10:        sumsq = x*x - y*y - z*z;
        !            11:        if (sumsq <= 0)
        !            12:                return 0;
        !            13:        return ( sqrt(sumsq) );
        !            14: }
        !            15: 
        !            16: #define MAXROOT 0xb504
        !            17: sqrt(x)
        !            18:        register long x;
        !            19: {
        !            20:        register long high=MAXROOT;
        !            21:        register long low=0;
        !            22:        register long current=MAXROOT/2;
        !            23:        if(x<=0)
        !            24:                return 0;
        !            25:        if(x>=MAXROOT*MAXROOT)
        !            26:                return(MAXROOT);
        !            27:        while(high>low+1){
        !            28:                if(current*current==x)
        !            29:                        return (current);
        !            30:                if(current*current>x)
        !            31:                        high=current;
        !            32:                else
        !            33:                        low=current;
        !            34:                current=(high+low)>>1;
        !            35:        }
        !            36:        return(current);
        !            37: }

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