|
|
1.1 ! root 1: double cabs(real, imag) ! 2: double real, imag; ! 3: { ! 4: double temp, sqrt(); ! 5: ! 6: if(real < 0) ! 7: real = -real; ! 8: if(imag < 0) ! 9: imag = -imag; ! 10: if(imag > real){ ! 11: temp = real; ! 12: real = imag; ! 13: imag = temp; ! 14: } ! 15: if((real+imag) == real) ! 16: return(real); ! 17: ! 18: temp = imag/real; ! 19: temp = real*sqrt(1.0 + temp*temp); /*overflow!!*/ ! 20: return(temp); ! 21: }
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.