From e0244b0fa7d884b52b7d7155e577550f71140044 Mon Sep 17 00:00:00 2001
From: David Minton <daminton@purdue.edu>
Date: Wed, 26 May 2021 23:12:25 -0400
Subject: [PATCH] Reverted loop back to regular do instead of forall

---
 src/util/util_solve_linear_system.f90 | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/src/util/util_solve_linear_system.f90 b/src/util/util_solve_linear_system.f90
index 1addfeb25..502d9b35e 100644
--- a/src/util/util_solve_linear_system.f90
+++ b/src/util/util_solve_linear_system.f90
@@ -96,7 +96,10 @@ function solve_wbs(u) result(x) ! solve with backward substitution
       x(:) = 0._DP
       return
    end if
-   forall (i=n:1:-1) x(i) = (u(i, n + 1) - sum(u(i, i + 1:n) * x(i + 1:n))) / u(i, i)
+   call ieee_set_halting_mode(ieee_divide_by_zero, .false.)
+   do i = n, 1, -1 
+      x(i) = (u(i, n + 1) - sum(u(i, i + 1:n) * x(i + 1:n))) / u(i, i)
+   end do   
    return
 end function solve_wbs