shithub: npe

ref: 1a54d1f0b13124bfbc60a4ecf3ed1f15869c8e3d
dir: /libnpe/fmax.c/

View raw version 
#include <math.h>
/* taken from musl */

#define _sign(d) (*((u64int*)d) & (1ULL<<63))

double
fmax(double x, double y)
{
	if(isNaN(x))
		return y;
	if(isNaN(y))
		return x;
	if(_sign(&x) != _sign(&y))
		return _sign(&x) ? y : x;
	return x < y ? y : x;
}