[unix-history] / usr / src / usr.bin / f77 / libF77 / cabs.c

/*
 *	"@(#)cabs.c	1.1"
 */

double cabs(real, imag)
double real, imag;
{
double temp, sqrt();

if(real < 0)
	real = -real;
if(imag < 0)
	imag = -imag;
if(imag > real){
	temp = real;
	real = imag;
	imag = temp;
}
if((real+imag) == real)
	return(real);

temp = imag/real;
temp = real*sqrt(1.0 + temp*temp);  /*overflow!!*/
return(temp);
}
Commit	Line	Data
af189a9d DW	1	/*
	2	* "@(#)cabs.c 1.1"
	3	*/
	4
	5	double cabs(real, imag)
	6	double real, imag;
	7	{
	8	double temp, sqrt();
	9
	10	if(real < 0)
	11	real = -real;
	12	if(imag < 0)
	13	imag = -imag;
	14	if(imag > real){
	15	temp = real;
	16	real = imag;
	17	imag = temp;
	18	}
	19	if((real+imag) == real)
	20	return(real);
	21
	22	temp = imag/real;
	23	temp = realsqrt(1.0 + temptemp); /overflow!!/
	24	return(temp);
	25	}