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