Return to zabs.c CVS log

Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.lib / libF77

Annotation of cci/usr/src/usr.lib/libF77/zabs.c, revision 1.1

1.1     ! root        1: /*
        !             2:  *     "@(#)zabs.c     1.1"
        !             3:  */
        !             4: 
        !             5: /* THIS IS BASED ON TAHOE FP REPRESENTATION */
        !             6: #include <FP.h>
        !             7: 
        !             8: double zabs(real, imag)
        !             9: double real, imag;
        !            10: {
        !            11: double temp, sqrt();
        !            12: 
        !            13: if(real < 0)
        !            14:        *(long int *)&real ^= SIGN_BIT;
        !            15: if(imag < 0)
        !            16:        *(long int *)&imag ^= SIGN_BIT;
        !            17: if(imag > real){
        !            18:        temp = real;
        !            19:        real = imag;
        !            20:        imag = temp;
        !            21: }
        !            22: if(imag == 0.)         /* if((real+imag) == real) */
        !            23:        return(real);
        !            24: 
        !            25: temp = imag/real;
        !            26: temp = real*sqrt(1.0 + temp*temp);  /*overflow!!*/
        !            27: return(temp);
        !            28: }

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