X-Git-Url: http://git.shadowcat.co.uk/gitweb/gitweb.cgi?a=blobdiff_plain;f=numeric.c;h=f12e9ba2d3c7353c0516063fef8b8f83ad62c034;hb=0af80b6034aad516a126a9414dadccac4de7f9dc;hp=e7a5e8d7e69a0aa7d21088a0d76f10fe0033b36e;hpb=aa8b85def9003d4bba549b6080f6f4864f347f13;p=p5sagit%2Fp5-mst-13.2.git

diff --git a/numeric.c b/numeric.c
index e7a5e8d..f12e9ba 100644
--- a/numeric.c
+++ b/numeric.c
@@ -572,7 +572,19 @@ S_mulexp10(NV value, I32 exponent)
 	negative = 1;
 	exponent = -exponent;
     }
-#ifdef __VAX /* avoid %SYSTEM-F-FLTOVF_F sans VAXC$ESTABLISH */
+
+    /* Avoid %SYSTEM-F-FLTOVF_F sans VAXC$ESTABLISH.
+     * In VAX VMS we by default use the D_FLOAT double format,
+     * and that format does not have *easy* capabilities [1] for
+     * overflowing doubles 'silently' as IEEE fp does.  Therefore we
+     * need to detect early whether we would overflow (this is
+     * the behaviour of the native string-to-float conversion routines,
+     * and therefore the behaviour of native applications, too.)
+     *
+     * [1] VAXC$EXTABLISH is the capability but it is basically a signal
+     * handler setup routine, and one cannot return from a fp exception
+     * handler and except much anything useful. */
+#if defined(VMS) && !defined(__IEEE_FP)
 #  if defined(__DECC_VER) && __DECC_VER <= 50390006
     /* __F_FLT_MAX_10_EXP - 5 == 33 */
     if (!negative &&
@@ -580,11 +592,26 @@ S_mulexp10(NV value, I32 exponent)
         return NV_MAX;
 #  endif
 #endif
+
+    /* In UNICOS and in certain Cray models (such as T90) there is no
+     * IEEE fp, and no way at all from C to catch fp overflows gracefully.
+     * There is something you can do if you are willing to use some
+     * inline assembler: the instruction is called DFI-- but that will
+     * disable *all* floating point interrupts, a little bit too large
+     * a hammer.  Therefore we need to catch potential overflows before
+     * it's too late. */
+#if defined(_UNICOS) && defined(NV_MAX_10_EXP)
+    if (!negative &&
+	(log10(value) + exponent) >= NV_MAX_10_EXP)
+        return NV_MAX;
+#endif
+
     for (bit = 1; exponent; bit <<= 1) {
 	if (exponent & bit) {
 	    exponent ^= bit;
 	    result *= power;
 	}
+	/* Floating point exceptions are supposed to be turned off. */
 	power *= power;
     }
     return negative ? value / result : value * result;