From ac201ff1ac0a893840672d06cac345b9a83392cc Mon Sep 17 00:00:00 2001
From: Jameson Nash <jameson@mit.edu>
Date: Tue, 22 May 2012 19:02:48 -0400
Subject: [PATCH] fpclassify definition

---
 Make.inc           |  2 +-
 src/s_fpclassify.c | 81 +++++++++++++++++++++++++++++++++++++++++++---
 2 files changed, 78 insertions(+), 5 deletions(-)

diff --git a/Make.inc b/Make.inc
index 4c996de..aeec0ab 100644
--- a/Make.inc
+++ b/Make.inc
@@ -1,5 +1,5 @@
 CC=gcc
-CFLAGS=-O2 -I. -I../include -I../ld128 -I../src -D__BSD_VISIBLE -Wno-implicit-function-declaration
+CFLAGS=-Wall -O2 -I. -I../include -I../ld128 -I../src -D__BSD_VISIBLE -Wno-implicit-function-declaration
 
 default: all
 
diff --git a/src/s_fpclassify.c b/src/s_fpclassify.c
index 6d9e329..03a1e98 100644
--- a/src/s_fpclassify.c
+++ b/src/s_fpclassify.c
@@ -29,7 +29,80 @@
 
 #include "fpmath.h"
 
-#warning "todo: fpclassify needs to be defined"
-int	__fpclassifyd(double);
-int	__fpclassifyf(float);
-int	__fpclassifyl(long double);
+//#define FP_INFINITE     0x01
+//#define FP_NAN          0x02
+//#define FP_NORMAL       0x04
+//#define FP_SUBNORMAL    0x08
+//#define FP_ZERO         0x10
+//#define fpclassify(x) \
+//    ((sizeof (x) == sizeof (float)) ? __fpclassifyf(x) \
+//    : (sizeof (x) == sizeof (double)) ? __fpclassifyd(x) \
+//    : __fpclassifyl(x))
+//
+
+int
+__fpclassifyd(double d)
+{
+	union IEEEd2bits u;
+
+	u.d = d;
+	if (u.bits.exp == 2047) {
+		if (u.bits.manl == 0 && u.bits.manh == 0) {
+			return FP_INFINITE;
+		} else {
+			return FP_NAN;
+		}
+	} else if (u.bits.exp != 0) {
+		return FP_NORMAL;
+	} else if (u.bits.manl == 0 && u.bits.manh == 0) {
+		return FP_ZERO;
+	} else {
+		return FP_SUBNORMAL;
+	}
+}
+		
+
+int
+__fpclassifyf(float f)
+{
+	union IEEEf2bits u;
+
+	u.f = f;
+	if (u.bits.exp == 255) {
+	   if (u.bits.man == 0) {
+		   return FP_INFINITE;
+	   } else {
+		   return FP_NAN;
+	   }
+	} else if (u.bits.exp != 0) {
+		return FP_NORMAL;
+	} else if (u.bits.man == 0) {
+		return FP_ZERO;
+	} else {
+		return FP_SUBNORMAL;
+	}
+}
+
+int
+__fpclassifyl(long double e)
+{
+	union IEEEl2bits u;
+
+	u.e = e;
+	mask_nbit_l(u);
+	if (u.bits.exp == 32767) {
+		if (u.bits.manl == 0 && u.bits.manh == 0) {
+			return FP_INFINITE;
+		} else {
+			return FP_NAN;
+		}
+	} else if (u.bits.exp != 0) {
+		return FP_NORMAL;
+	} else if (u.bits.manl == 0 && u.bits.manh == 0) {
+		return FP_ZERO;
+	} else {
+		return FP_SUBNORMAL;
+	}
+}
+
+