From eecf116894e1bbc5fdee0fec7984008f9f3b04b4 Mon Sep 17 00:00:00 2001
From: Jarrett Revels <jarrettrevels@gmail.com>
Date: Fri, 31 Aug 2018 10:52:34 -0400
Subject: [PATCH] Revert "Use broadcast expressions instead of some loops."

This reverts commit 19003af62af38a110fa6ec12d3b52cd3451e4a03.

See https://github.com/JuliaDiff/ForwardDiff.jl/pull/353#issuecomment-417688351
---
 src/apiutils.jl | 20 ++++++++++++++------
 src/gradient.jl |  5 +++--
 2 files changed, 17 insertions(+), 8 deletions(-)

diff --git a/src/apiutils.jl b/src/apiutils.jl
index 44b931df..eec39653 100644
--- a/src/apiutils.jl
+++ b/src/apiutils.jl
@@ -53,28 +53,36 @@ end
 
 function seed!(duals::AbstractArray{Dual{T,V,N}}, x,
                seed::Partials{N,V} = zero(Partials{N,V})) where {T,V,N}
-    duals .= Dual{T,V,N}.(x, Base.RefValue(seed))
+    for i in eachindex(duals)
+        duals[i] = Dual{T,V,N}(x[i], seed)
+    end
     return duals
 end
 
 function seed!(duals::AbstractArray{Dual{T,V,N}}, x,
                seeds::NTuple{N,Partials{N,V}}) where {T,V,N}
-    duals[1:N] .= Dual{T,V,N}.(x[1:N], seeds[1:N])
+    for i in 1:N
+        duals[i] = Dual{T,V,N}(x[i], seeds[i])
+    end
     return duals
 end
 
 function seed!(duals::AbstractArray{Dual{T,V,N}}, x, index,
                seed::Partials{N,V} = zero(Partials{N,V})) where {T,V,N}
     offset = index - 1
-    chunk = (1:N) .+ offset
-    duals[chunk] .= Dual{T,V,N}.(x[chunk], Base.RefValue(seed))
+    for i in 1:N
+        j = i + offset
+        duals[j] = Dual{T,V,N}(x[j], seed)
+    end
     return duals
 end
 
 function seed!(duals::AbstractArray{Dual{T,V,N}}, x, index,
                seeds::NTuple{N,Partials{N,V}}, chunksize = N) where {T,V,N}
     offset = index - 1
-    chunk = (1:chunksize) .+ offset
-    duals[chunk] .= Dual{T,V,N}.(x[chunk], seeds[1:chunksize])
+    for i in 1:chunksize
+        j = i + offset
+        duals[j] = Dual{T,V,N}(x[j], seeds[i])
+    end
     return duals
 end
diff --git a/src/gradient.jl b/src/gradient.jl
index 771bffa7..8721af6d 100644
--- a/src/gradient.jl
+++ b/src/gradient.jl
@@ -77,8 +77,9 @@ extract_gradient!(::Type{T}, result::AbstractArray, dual::Dual) where {T}= copyt
 
 function extract_gradient_chunk!(::Type{T}, result, dual, index, chunksize) where {T}
     offset = index - 1
-    chunk = (1:chunksize) .+ offset
-    result[chunk] .= ForwardDiff.partials(T, dual, 1:chunksize)
+    for i in 1:chunksize
+        result[i + offset] = partials(T, dual, i)
+    end
     return result
 end