@@ -52,30 +52,3 @@ function _initialize_dae!(integrator, prob::ImplicitDiscreteProblem,
5252    end 
5353end 
5454
55- # ### TODO : Implement real algorithm
56- #  function _initialize_dae!(integrator, prob::ImplicitDiscreteProblem, alg::BrownFullBasicInit, isinplace::Val{true}) 
57- #      @unpack p, t, f = integrator
58- #      u0 = integrator.u
59- #      
60- #      nlequation! = (out, u, p) -> begin 
61- #          f(out, u, u0, p, t)
62- #      end
63- #  
64- #      nlfunc = NonlinearFunction(nlequation!; jac_prototype = f.jac_prototype)
65- #      nlprob = NonlinearProblem(nlfunc, ifelse.(differential_vars, du, u), p)
66- #      nlsol = solve(nlprob, nlsolve; abstol = alg.abstol, reltol = integrator.opts.reltol)
67- #  
68- #      @. du = ifelse(differential_vars, nlsol.u, du)
69- #      @. u = ifelse(differential_vars, u, nlsol.u)
70- #  
71- #      recursivecopy!(integrator.uprev, integrator.u)
72- #      if alg_extrapolates(integrator.alg)
73- #          recursivecopy!(integrator.uprev2, integrator.uprev)
74- #      end
75- #  
76- #      if nlsol.retcode != ReturnCode.Success
77- #          integrator.sol = SciMLBase.solution_new_retcode(integrator.sol,
78- #              ReturnCode.InitialFailure)
79- #      end
80- #      return
81- #  end
0 commit comments