1*f504f610SAugustin Cavalier #include <math.h> 2*f504f610SAugustin Cavalier fdim(double x,double y)3*f504f610SAugustin Cavalierdouble fdim(double x, double y) 4*f504f610SAugustin Cavalier { 5*f504f610SAugustin Cavalier if (isnan(x)) 6*f504f610SAugustin Cavalier return x; 7*f504f610SAugustin Cavalier if (isnan(y)) 8*f504f610SAugustin Cavalier return y; 9*f504f610SAugustin Cavalier return x > y ? x - y : 0; 10*f504f610SAugustin Cavalier } 11