fpudebug

1 files changed, 18 insertions, 1 deletions

diff --git a/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c b/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c
index f4ae0f6..77a8847 100644
--- a/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c
+++ b/lib/Numeric/LinearAlgebra/LAPACK/lapack-aux.c
@@ -47,7 +47,24 @@
 //---------------------------------------
 void asm_finit() {
 #ifdef i386
-    asm("finit");
+
+    #if FPUDEBUG
+        uint val=0;
+        asm("fstsw"
+        : "=a" (val)
+        : "a"  (val)
+        );
+
+        val = (val%16384)/2048; // bits 13-11
+
+        if (val != 0) {
+            printf("Warning: FPU Stack: %d\n",val);
+            asm("finit");
+        }
+    #else
+        asm("finit");
+    #endif
+
 #endif
 }
 //---------------------------------------