comment out initialize method

Author: vyudu

lib/SimpleImplicitDiscreteSolve/src/solve.jl

@@ -53,29 +53,29 @@ function _initialize_dae!(integrator, prob::ImplicitDiscreteProblem,
 end
 
 #### TODO: Implement real algorithm
-function _initialize_dae!(integrator, prob::ImplicitDiscreteProblem, alg::BrownFullBasicInit, isinplace::Val{true}) 
-    @unpack p, t, f = integrator
-    u0 = integrator.u
-    
-    nlequation! = (out, u, p) -> begin 
-        f(out, u, u0, p, t)
-    end
-
-    nlfunc = NonlinearFunction(nlequation!; jac_prototype = f.jac_prototype)
-    nlprob = NonlinearProblem(nlfunc, ifelse.(differential_vars, du, u), p)
-    nlsol = solve(nlprob, nlsolve; abstol = alg.abstol, reltol = integrator.opts.reltol)
-
-    @. du = ifelse(differential_vars, nlsol.u, du)
-    @. u = ifelse(differential_vars, u, nlsol.u)
-
-    recursivecopy!(integrator.uprev, integrator.u)
-    if alg_extrapolates(integrator.alg)
-        recursivecopy!(integrator.uprev2, integrator.uprev)
-    end
-
-    if nlsol.retcode != ReturnCode.Success
-        integrator.sol = SciMLBase.solution_new_retcode(integrator.sol,
-            ReturnCode.InitialFailure)
-    end
-    return
-end
+# function _initialize_dae!(integrator, prob::ImplicitDiscreteProblem, alg::BrownFullBasicInit, isinplace::Val{true}) 
+#     @unpack p, t, f = integrator
+#     u0 = integrator.u
+#     
+#     nlequation! = (out, u, p) -> begin 
+#         f(out, u, u0, p, t)
+#     end
+# 
+#     nlfunc = NonlinearFunction(nlequation!; jac_prototype = f.jac_prototype)
+#     nlprob = NonlinearProblem(nlfunc, ifelse.(differential_vars, du, u), p)
+#     nlsol = solve(nlprob, nlsolve; abstol = alg.abstol, reltol = integrator.opts.reltol)
+# 
+#     @. du = ifelse(differential_vars, nlsol.u, du)
+#     @. u = ifelse(differential_vars, u, nlsol.u)
+# 
+#     recursivecopy!(integrator.uprev, integrator.u)
+#     if alg_extrapolates(integrator.alg)
+#         recursivecopy!(integrator.uprev2, integrator.uprev)
+#     end
+# 
+#     if nlsol.retcode != ReturnCode.Success
+#         integrator.sol = SciMLBase.solution_new_retcode(integrator.sol,
+#             ReturnCode.InitialFailure)
+#     end
+#     return
+# end