From 3d7aceb99bb35bb6b547f83d31ccba29713713ff Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Wed, 11 Nov 2020 20:57:17 +0100
Subject: [PATCH] avoid underflow/division by zero

---
 src/math.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/math.f90 b/src/math.f90
index b01cf9e26..1fc9bd7eb 100644
--- a/src/math.f90
+++ b/src/math.f90
@@ -1200,8 +1200,8 @@ subroutine selfTest
 
   if(any(dNeq(math_exp33(math_I3,0),math_I3))) &
     error stop 'math_exp33(math_I3,1)'
-  if(any(dNeq(math_exp33(math_I3,256),exp(1.0_pReal)*math_I3))) &
-    error stop 'math_exp33(math_I3,256)'
+  if(any(dNeq(math_exp33(math_I3,128),exp(1.0_pReal)*math_I3))) &
+    error stop 'math_exp33(math_I3,128)'
 
   call random_number(v9)
   if(any(dNeq(math_33to9(math_9to33(v9)),v9))) &