blob: c244fbcb0ac60ddf6830ffd24e9d33de64d5dcc5 [file] [log] [blame]
#include <math.h>
double fmax(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) ? y : x;
return x < y ? y : x;
}