summaryrefslogtreecommitdiff
path: root/lib/libm/fminl.c
blob: 087b87b51cac87f96b547ba24b133e80707d1ffd (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
#include <float.h> // for LDBL_MANT_DIG, LDBL_MAX_EXP
#include <math.h>  // for signbit, fminl, isnan

#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double fminl(long double x, long double y)
{
	return fmin(x, y);
}
#else
long double fminl(long double x, long 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) ? x : y;
	return x < y ? x : y;
}
#endif