From 72948ae0fd1374534c0fa356c39619cfa1947adc Mon Sep 17 00:00:00 2001
From: Bruce Mitchener <bruce.mitchener@gmail.com>
Date: Wed, 23 Dec 2020 14:23:26 +0700
Subject: [PATCH] Clean up medleyfp.h. (#127)

* Include <ieeefp.h> for Solaris.
* Remove sections for sparc and Sun i386 for non-Solaris.
* Have Linux use the same code as FreeBSD and macOS.
* Make the FreeBSD / Linux / macOS code path the default, if nothing
  special is defined since it is the modern way.
* Remove fpstatus_() usage from the ancient days.
---
 inc/medleyfp.h | 41 +++++------------------------------------
 1 file changed, 5 insertions(+), 36 deletions(-)

diff --git a/inc/medleyfp.h b/inc/medleyfp.h
index cfa4b20..10d782c 100755
--- a/inc/medleyfp.h
+++ b/inc/medleyfp.h
@@ -27,7 +27,6 @@
 #ifdef FLTINT
 volatile extern int  FP_error;
 
-
 /*  Note that a compiler may very likely move code around the arithmetic
     operation, causing this test (set by an interrupt handler) to be
     incorrect.  For example, the Sun SPARC compiler with -O2 makes
@@ -37,49 +36,19 @@ volatile extern int  FP_error;
 #define FPCLEAR         FP_error = 0;
 #define FPTEST(result)  FP_error
 
-#else
-
-/*  fpstatus_ is a FORTRAN library routine (in libc) which
-    can be called to determine floating point status results.
-    Documented in the Sun manual, "Floating Point Programmer's Guide",
-    (Rev. A 19-Sep-86), pg. 34, it does *not* exist in libc for the
-    SPARC.
-
-    For sparc, should also check for isnan?  Don't know what isnormal
-    & issubnormal do (these are sunos4.0 only)
- */
-#if defined(OS5)
+#elif defined(OS5)
+#include <ieeefp.h>
 #define FPCLEAR
 #define FPTEST(result) (!finite(result))
 
-#elif defined(sparc) || defined(I386)
-#define FPCLEAR
-#define FPTEST(result) (isinf(result) || isnan(result))
-
-#elif defined(LINUX)
-#include <math.h>
-#define FPCLEAR
-#define FPTEST(result) ((!finite(result)) || isnan(result))
-
-#elif defined(MACOSX) || defined(FREEBSD)
-#include <math.h>
-#define FPCLEAR
-#define FPTEST(result) (!isfinite(result))
-
 #elif defined(DOS)
 #include <i32.h>
 #define FPCLEAR
 #define FPTEST(result) (_getrealerror() & ( I87_ZERO_DIVIDE | I87_OVERFLOW | I87_UNDERFLOW))
 
 #else
-static int constant0 = 0;
-unsigned int fpstatus_();
+#include <math.h>
+#define FPCLEAR
+#define FPTEST(result) (!isfinite(result))
 
-#define FPCLEAR         fpstatus_(&constant0);
-#define FPTEST(result) (fpstatus_(&constant0) & 0xF0)
-
-#endif
 #endif /* FLTINT */
-
-
-