Efficient parallel optimization for approximating CAD curves featuring super-convergence
We present an efficient, parallel, constrained optimization technique for approximating CAD curves with super-convergent rates. The optimization function is a disparity measure in terms of a piece-wise polynomial approximation and a curve re-parametrization. The constrained problem solves the disparity functional fixing the mesh element interfaces. We have numerical evidence that the constrained disparity preserves the original super-convergence: 2p order for planar curves and ⌊3/2(p-1)⌋ + 2 for 3D curves, p being the mesh polynomial degree. Our optimization scheme consists of a globalized Newton method with a nonmonotone line search, and a log barrier function preventing element inversion in the curve re-parameterization. Moreover, we introduce a Julia interface to the EGADS geometry kernel and a parallel optimization algorithm. We test the potential of our curve mesh generation tool on a computer cluster using several aircraft CAD models. We conclude that the solver is well-suited for parallel computing, producing super-convergent approximations to CAD curves.
READ FULL TEXT