  Solving nonlinear systems is an important problem. Numerical continuation
methods efficiently solve certain nonlinear systems. The Asymptotic Numerical
Method (ANM) is a powerful continuation method that usually converges faster
than Newtonian methods. ANM explores the landscape of the function by following
a parameterized solution curve approximated with a high-order power series.
Although ANM has successfully solved a few graphics and engineering problems,
prior to our work, applying ANM to new problems required significant effort
because the standard ANM assumes quadratic functions, while manually deriving
the power series expansion for nonquadratic systems is a tedious and
challenging task.
  This paper presents a novel solver, SANM, that applies ANM to solve
symbolically represented nonlinear systems. SANM solves such systems in a fully
automated manner. SANM also extends ANM to support many nonquadratic operators,
including intricate ones such as singular value decomposition. Furthermore,
SANM generalizes ANM to support the implicit homotopy form. Moreover, SANM
achieves high computing performance via optimized system design and
implementation.
  We deploy SANM to solve forward and inverse elastic force equilibrium
problems and controlled mesh deformation problems with a few constitutive
models. Our results show that SANM converges faster than Newtonian solvers,
requires little programming effort for new problems, and delivers comparable or
better performance than a hand-coded, specialized ANM solver. While we
demonstrate on mesh deformation problems, SANM is generic and potentially
applicable to many tasks.
