1#include <math.h>
2#include <float.h>
3
4#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
5long double fdiml(long double x, long double y)
6{
7	return fdim(x, y);
8}
9#else
10long double fdiml(long double x, long double y)
11{
12	if (isnan(x))
13		return x;
14	if (isnan(y))
15		return y;
16	return x > y ? x - y : 0;
17}
18#endif
19