#include <math.h>
#include <private/ieee754.h>

int __fpclassify(double x) {
  if (ieee754exp64(x) == 0) {
    if (ieee754man64(x) == 0x0) {
      return FP_ZERO;
    } else {
      return FP_SUBNORMAL;
    }
  }
  if (ieee754exp64(x) == 0x7FF) {
    if (ieee754man64(x) == 0) {
      return FP_INFINITE;
    } else {
      return FP_NAN;
    }
  }
  return FP_NORMAL;
}