Return to zabs.c CVS log

Up to [CSRG BSD Unix] / 43BSDReno / lib / libF77

Annotation of 43BSDReno/lib/libF77/zabs.c, revision 1.1.1.1

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

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