Newton + LineSearch in Hot Loop: Reduce allocations

Okay, one more thing: It seems to be not as trivial to set the Jacobian by hand (I cannot resort to AD). I found this post but unfortunately NonlinearSolve.jl still tries to compute the Jacobian using AD.