| 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 | ||||