1#include <math.h>
2#include <stdint.h>
3
4double fabs(double x)
5{
6	union {double f; uint64_t i;} u = {x};
7	u.i &= -1ULL/2;
8	return u.f;
9}
10