# GIPC: Fast and stable Gauss-Newton optimization of IPC barrier energy

arXiv: https://arxiv.org/abs/2308.09400 · HTML: https://mlsanity.com/paper/2308.09400/

## Key facts

- Cited by 7 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 18 Aug 2023.

## Built on by (most cited first, 7 of 7)

- [DICE: End-to-end Deformation Capture of Hand-Face Interactions from a Single Image](https://mlsanity.com/paper/2406.17988/) (2024)
- [WorldParticle: Unified World Simulation of Lagrangian Particle Dynamics via Transformer](https://mlsanity.com/paper/2605.15305/) (2026)
- [ClothTransformer: Unified Latent-Space Transformers for Scalable Cloth Simulation](https://mlsanity.com/paper/2605.27852/) (2026)
- [4DCodeBench: Benchmarking Agents on Inverse Graphics of Dynamic Scenes](https://mlsanity.com/paper/2610.03715/) (2026)
- [Tac2Real: Reliable and GPU Visuotactile Simulation for Online Reinforcement Learning and Zero-Shot Real-World Deployment](https://mlsanity.com/paper/2603.28475/) (2026)
- [Adaptive Volumetric Mechanical Property Fields Invariant to Resolution](https://mlsanity.com/paper/2606.18231/) (2026)
- [An Efficient Multilevel Preconditioned Nonlinear Conjugate Gradient Method for Incremental Potential Contact](https://mlsanity.com/paper/2604.19892/) (2026)

## Abstract (arXiv)

Barrier functions are crucial for maintaining an intersection and inversion free simulation trajectory but existing methods which directly use distance can restrict implementation design and performance. We present an approach to rewriting the barrier function for arriving at an efficient and robust approximation of its Hessian. The key idea is to formulate a simplicial geometric measure of contact using mesh boundary elements, from which analytic eigensystems are derived and enhanced with filtering and stiffening terms that ensure robustness with respect to the convergence of a Project-Newton solver. A further advantage of our rewriting of the barrier function is that it naturally caters to the notorious case of nearly-parallel edge-edge contacts for which we also present a novel analytic eigensystem. Our approach is thus well suited for standard second order unconstrained optimization strategies for resolving contacts, minimizing nonlinear nonconvex functions where the Hessian may be indefinite. The efficiency of our eigensystems alone yields a 3x speedup over the standard IPC barrier formulation. We further apply our analytic proxy eigensystems to produce an entirely GPU-based implementation of IPC with significant further acceleration.