fix non-deterministic behavior by forcing instantiations (#284)

MadNLP · Dec 11, 2023 · 3685435 · 3685435
1 parent 805018a
commit 3685435
Showing 1 changed file with 7 additions and 8 deletions.
diff --git a/src/nlpmodels.jl b/src/nlpmodels.jl
@@ -247,10 +247,10 @@ function create_callback(
 
     x0   = get_x0(nlp)
 
-    con_buffer = similar(x0, m)
-    grad_buffer = similar(x0, n)
-    jac_buffer = similar(x0, nnzj)
-    hess_buffer = similar(x0, nnzh)
+    con_buffer = similar(x0, m)     ; fill!(con_buffer, zero(T))
+    grad_buffer = similar(x0, n)    ; fill!(grad_buffer, zero(T))
+    jac_buffer = similar(x0, nnzj)  ; fill!(jac_buffer, zero(T))
+    hess_buffer = similar(x0, nnzh) ; fill!(hess_buffer, zero(T))
 
     jac_I = similar(x0, Int, nnzj)
     jac_J = similar(x0, Int, nnzj)
@@ -274,7 +274,6 @@ function create_callback(
 
     equality_handler = opt.equality_treatment()
 
-
     return SparseCallback(
         nlp,
         n,m,nnzj,nnzh,
@@ -304,9 +303,9 @@ function create_callback(
     m = get_ncon(nlp)
 
     x0   = similar(get_x0(nlp))
-    con_buffer = similar(x0, m)
-    jac_buffer = similar(x0, m, n)
-    grad_buffer = similar(x0, n)
+    con_buffer = similar(x0, m) ;    fill!(con_buffer, zero(T))
+    jac_buffer = similar(x0, m, n) ; fill!(jac_buffer, zero(T))
+    grad_buffer = similar(x0, n) ;   fill!(grad_buffer, zero(T))
     obj_scale = Ref(one(T))
     con_scale = similar(x0, m) ; fill!(con_scale, one(T))