aboutsummaryrefslogtreecommitdiff
path: root/src/math/fmax.c
blob: ee52f3fc4f4de2f7777ba8a265a7c93fc531ec81 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
#include <math.h>

double fmax(double x, double y)
{
	if (isnan(x))
		return y;
	if (isnan(y))
		return x;
	/* handle signed zeros, see C99 Annex F.9.9.2 */
	if (signbit(x) != signbit(y))
		return signbit(x) ? y : x;
	return x < y ? y : x;
}


#ifdef TEST
int
main(void) {
	return 0;
}
#endif