github.com/afumu/libc@v0.0.6/musl/src/math/fdiml.c (about)

     1  #include <math.h>
     2  #include <float.h>
     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