Plan 9 from Bell Labs’s /usr/web/sources/plan9/sys/src/ape/lib/ap/math/hypot.c

Copyright © 2021 Plan 9 Foundation.
Distributed under the MIT License.
Download the Plan 9 distribution.


#include <math.h>
/*
 * sqrt(a^2 + b^2)
 *	(but carefully)
 */

double
hypot(double a, double b)
{
	double t;

	if(a < 0)
		a = -a;
	if(b < 0)
		b = -b;
	if(a > b) {
		t = a;
		a = b;
		b = t;
	}
	if(b == 0) 
		return 0;
	a /= b;
	/*
	 * pathological overflow possible
	 * in the next line.
	 */
	return b * sqrt(1 + a*a);
}

Bell Labs OSI certified Powered by Plan 9

(Return to Plan 9 Home Page)

Copyright © 2021 Plan 9 Foundation. All Rights Reserved.
Comments to [email protected].