Use a copy! instead of copy in DIOM

JuliaSmoothOptimizers · Oct 14, 2018 · 3ab0c65 · 3ab0c65
1 parent 1e5569b
commit 3ab0c65
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/diom.jl b/src/diom.jl
@@ -156,7 +156,7 @@ function diom{T <: Number}(A :: AbstractLinearOperator, b :: AbstractVector{T};
 
     # Update x_old and residual norm.
     if !p[next_pos]
-      x_old = copy(x)
+      copy!(x_old, x)
       # ‖ b - Axₘ ‖ = hₘ₊₁.ₘ * |ξₘ / uₘ.ₘ| without pivoting
       rNorm = H[1] * abs(ξ / H[2])
     else