diff --git a/src/util/util_solve_linear_system.f90 b/src/util/util_solve_linear_system.f90
index b44f9168f..b81d6a15b 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 .or. (abs(qx) > huge(x)) .or. (abs(qx) < tiny(x))) then
+   if (lerr .or. (any(abs(qx) > huge(x))) .or. (any(abs(qx) < tiny(x)))) then
       x = 0.0_DP
       write(*,*) 'fpe in util_solve_linear_system'
    else