SANM: A Symbolic Asymptotic Numerical Solver with Applications in Mesh Deformation

05/18/2021
by   Kai Jia, et al.
0

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.

READ FULL TEXT
research
09/30/2022

A high-order fast direct solver for surface PDEs

We introduce a fast direct solver for variable-coefficient elliptic part...
research
04/21/2022

Singularity-Avoiding Multi-Dimensional Root-Finder

We proposed in this paper a new method, which we named the W4 method, to...
research
01/22/2021

Implicit shock tracking for unsteady flows by the method of lines

A recently developed high-order implicit shock tracking (HOIST) framewor...
research
01/22/2020

An asymptotic preserving semi-implicit multiderivative solver

In this work we construct a multiderivative implicit-explicit (IMEX) sch...
research
12/24/2019

Implicit shock tracking using an optimization-based, r-adaptive, high-order discontinuous Galerkin method

A novel framework for resolving discontinuous solutions of conservation ...
research
05/12/2018

Fast Symbolic 3D Registration Solution

3D registration has always been performed invoking singular value decomp...

Please sign up or login with your details

Forgot password? Click here to reset