Monolithic homotopy continuation with predictor based on higher derivatives

The predictor component of a monolithic homotopy continuation algorithm is augmented with higher derivative information for use as an efficient, robust, and scalable continuation algorithm suitable for application to large sparse systems of nonlinear algebraic equations. Convergence of the algorithm is established analytically, and efficiency studies are performed by applying the method to a practical computational aerodynamics problem.