Return to cabs.c CVS log

Up to [Research Unix] / researchv10no / libF77

Annotation of researchv10no/libF77/cabs.c, revision 1.1.1.1

1.1       root        1: #ifdef KR_headers
                      2: extern double sqrt();
                      3: double f__cabs(real, imag) double real, imag;
                      4: #else
                      5: #undef abs
                      6: #include "math.h"
                      7: double f__cabs(double real, double imag)
                      8: #endif
                      9: {
                     10: double temp;
                     11: 
                     12: if(real < 0)
                     13:        real = -real;
                     14: if(imag < 0)
                     15:        imag = -imag;
                     16: if(imag > real){
                     17:        temp = real;
                     18:        real = imag;
                     19:        imag = temp;
                     20: }
                     21: if((real+imag) == real)
                     22:        return(real);
                     23: 
                     24: temp = imag/real;
                     25: temp = real*sqrt(1.0 + temp*temp);  /*overflow!!*/
                     26: return(temp);
                     27: }

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