This is an implementation of SGP2011 paper "Deformable 3D Shape Registration Based on Local Similarity Transforms" and SGP2013 paper "Locally Injective Mappings". Generally speaking, ASAP is an extension of ARAP with scaling factor added, which makes local deformation scalable. We added a barier function to penalty the element inversion, whose energy will go up to infinity when the volume of the tetrahedra is close to zero. Levenberg-Marquardt algorithm is ultilized to optimize the non-linear system. In each linearized step, cholesky decomposition are used to solve the normal equation. LIBIGL is adopted to visualize the deformation procedure.
