use remake instead

Author: jClugstor

ext/LinearSolveForwardDiffExt.jl

@@ -144,8 +144,8 @@ function SciMLBase.init(
     ∂_A = partial_vals(A)
     ∂_b = partial_vals(b)
 
-    primal_prob = LinearProblem(new_A, new_b, u0 = new_u0)
-    #remake(prob; A = new_A, b = new_b, u0 = new_u0)
+    #primal_prob = LinearProblem(new_A, new_b, u0 = new_u0)
+    primal_prob = remake(prob; A = new_A, b = new_b, u0 = new_u0)
 
     if get_dual_type(prob.A) !== nothing
         dual_type = get_dual_type(prob.A)