From 0b50b65d81e3739491e9164ef3190d4f565398b1 Mon Sep 17 00:00:00 2001 From: David Minton Date: Wed, 26 May 2021 23:08:41 -0400 Subject: [PATCH] Added overflow 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 b81d6a15b..1addfeb25 100644 --- a/src/util/util_solve_linear_system.f90 +++ b/src/util/util_solve_linear_system.f90 @@ -92,7 +92,7 @@ function solve_wbs(u) result(x) ! solve with backward substitution n = size(u, 1) if (allocated(x)) deallocate(x) if (.not.allocated(x)) allocate(x(n)) - if (any(abs(u) < tiny(1._DP))) then + if (any(abs(u) < tiny(1._DP)) .or. any(abs(u) > huge(1._DP))) then x(:) = 0._DP return end if