“…, Zu,t+1 , Zχ,t+1 , Z∇χ,t+1 , γu,t+1 , γ χ,t+1 , γ∇χ,t+1 ← performNonlinearSolve else end 54 Kelley [1995] and Kelley [2003] but we utilize an initialization procedure called Homotopy or the continuation method detailed in section 3.6 of Kincaid and Cheney [2009] where we break each increment into pseudo-time steps such that…”