From 79de855bb1c162df046a7a6a803f7a582e71f2fb Mon Sep 17 00:00:00 2001 From: David Minton Date: Wed, 26 May 2021 23:00:33 -0400 Subject: [PATCH] Added floating point over/underflow check --- src/util/util_solve_linear_system.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/util/util_solve_linear_system.f90 b/src/util/util_solve_linear_system.f90 index 73dc49764..b44f9168f 100644 --- a/src/util/util_solve_linear_system.f90 +++ b/src/util/util_solve_linear_system.f90 @@ -30,7 +30,7 @@ function util_solve_linear_system_d(A,b,n,lerr) result(x) call ieee_get_flag(ieee_usual, fpe_flag) lerr = any(fpe_flag) - if (lerr) then + if (lerr .or. (abs(qx) > huge(x)) .or. (abs(qx) < tiny(x))) then x = 0.0_DP write(*,*) 'fpe in util_solve_linear_system' else