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