nsz | 9560b6b | 2012-03-13 19:51:14 +0100 | [diff] [blame] | 1 | #include <math.h> |
2 | #include <float.h> | ||||
Rich Felker | b69f695 | 2012-03-13 01:17:53 -0400 | [diff] [blame] | 3 | |
4 | #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 | ||||
5 | long double fdiml(long double x, long double y) | ||||
6 | { | ||||
7 | return fdim(x, y); | ||||
8 | } | ||||
9 | #else | ||||
10 | long 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 |