Return to abs.c CVS log

Up to [Research Unix] / researchv9 / cmd / cfront / libC / complex

Annotation of researchv9/cmd/cfront/libC/complex/abs.c, revision 1.1.1.1

1.1       root        1: 
                      2: #include "complex.h"
                      3: #include "const.h"
                      4: 
                      5: double abs(complex z)
                      6: {
                      7:        register double temp;
                      8:        register double r = z.re;
                      9:        register double i = z.im;
                     10:        
                     11:        if (r < 0) r = -r;
                     12:        if (i == 0) return r;
                     13: 
                     14:        if (i < 0) i = -i;
                     15:        if (r == 0) return i;
                     16: 
                     17:        if (i > r) { temp = r; r = i; i = temp; }
                     18:        temp = i/r;
                     19:        temp = r*sqrt(1.0 + temp*temp);  /*overflow!!*/
                     20:        return temp;
                     21: }
                     22: 
                     23: double norm(complex z)
                     24: {
                     25:        register double r = z.re;
                     26:        register double i = z.im;
                     27: 
                     28: #define        SAFE 6.521908912666391000e+18  /* 0.5*sqrt(GREATEST)*/
                     29: 
                     30:        if (r<SAFE && i<SAFE) return r*r+i*i;
                     31:        return GREATEST;
                     32: }

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