#include double fdim(double x, double y) { if (isnan(x)) return x; if (isnan(y)) return y; return x > y ? x - y : 0; } #ifdef TEST int main(void) { return 0; } #endif