switched to broadcasting operator instead of loop

QuantEcon · Nov 26, 2014 · 8d502d1 · 8d502d1
1 parent 74e812c
commit 8d502d1
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/src/mc_tools.jl b/src/mc_tools.jl
@@ -86,9 +86,8 @@ function gth_solve{T<:Real}(A::AbstractMatrix{T})
             n = k
             break
         end
-        for i in k+1:n
-            A1[i, k] /= scale
-        end
+        A1[k+1:n, k] ./= scale
+
         for j in k+1:n, i in k+1:n
             A1[i, j] += A1[i, k] * A1[k, j]
         end