Added overflow check

MintonGroup · May 27, 2021 · 0b50b65 · 0b50b65
1 parent c3f51cb
commit 0b50b65
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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