A robust hybrid algorithm for computing multiple equilibrium solutions

A hybrid method is described that seeks to combine the efficiency of a quasi‐Newton method capable of locating stable and unstable equilibrium configurations with a robust homotopy method that is capable of tracking equilibrium paths with turning points while exploiting symmetry and sparsity of the Jacobian matrices. Numerical results are presented for a shallow arch problem.