#include <math.h>

double fmin(double x, double y) {
  if (isnan(x))
    return y;
  if (isnan(y))
    return x;
  /* handle signed zeros, see C99 Annex F.9.9.2 */
  if (signbit(x) != signbit(y))
    return signbit(x) ? x : y;
  return x < y ? x : y;
}