TOGEARI: Interaction-Space Preconditioning for Condensed Finite-Element Systems with IPC Contact
Abstract
TOGEARI constructs a compressed contact correction for condensed finite-element equations with incremental potential contact (IPC). A core-conditioned interaction matrix ranks contact combinations by their mechanical responses; a raw-factor proxy acquires only the retained responses. For a symmetric positive-definite core, the reference analysis gives the exact spectrum, scaled inverse error and optimal worst omitted interaction. The selected responses are maintained through signed contact changes using supported operator images and a small coarse matrix. A complementary Schur spectrum and perturbation bounds connect this construction to the full current linear equation, damped local Newton steps and compatible episode reuse. On a 325,260-coordinate finger/cup problem, eight selected directions reduce Krylov work from eleven columns to five, matching full contact. Two reversed-order load studies give matching anchored displacements and accepted normal reactions, with 28.8-29.3 percent less warm time and 5.1-5.3 percent less complete time than per-step direct refactorization. A completed twelve-episode run saves 18.8 percent of total time against direct and 2.8 percent against a held factor. Full response caching has comparable time with 7.5 times the response storage. A compressed repeat terminates at an energy line search; separate diagnostics identify constitutive cancellation and terminal-backtracking sensitivity. Larger-contact, held-out-load and frictional comparisons show where compression or preparation loses its advantage. The complete Newton equation remains the residual operator throughout; the convergence guarantees retain their stated reference and local hypotheses.
Related
- Quantum Domain Decomposition for Preconditioning the Finite Element Method
- Ideal interaction dynamics of defects and background fields in condensed matter
- Low-Order Refined Preconditioning for Spectral/hp Element Method for Complex, 3D Geometries
- Preconditioning Magnetic Systems in Kohn-Sham Density Functional Theory
- A unified variational framework for augmented Lagrangian preconditioning of saddle point systems: applications to mixed-dimensional and fictitious domain methods
BibTeX
@misc{liu2026togeari,
title = {{TOGEARI: Interaction-Space Preconditioning for Condensed Finite-Element Systems with IPC Contact}},
author = {Yanlin Liu and Chao Huang and Kaixiang Yao and Yao Shen},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2608.14162},
eprint = {2608.14162},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.14162},
note = {CC BY-SA 4.0},
}
RIS
TY - UNPB
TI - TOGEARI: Interaction-Space Preconditioning for Condensed Finite-Element Systems with IPC Contact
AU - Yanlin Liu
AU - Chao Huang
AU - Kaixiang Yao
AU - Yao Shen
PY - 2026
DA - 2026/09/11
PB - arXiv
DO - 10.48550/arxiv.2608.14162
UR - https://arxiv.org/abs/2608.14162
AB - TOGEARI constructs a compressed contact correction for condensed finite-element equations with incremental potential contact (IPC). A core-conditioned interaction matrix ranks contact combinations by their mechanical responses; a raw-factor proxy acquires only the retained responses. For a symmetric positive-definite core, the reference analysis gives the exact spectrum, scaled inverse error and optimal worst omitted interaction. The selected responses are maintained through signed contact changes using supported operator images and a small coarse matrix. A complementary Schur spectrum and perturbation bounds connect this construction to the full current linear equation, damped local Newton steps and compatible episode reuse. On a 325,260-coordinate finger/cup problem, eight selected directions reduce Krylov work from eleven columns to five, matching full contact. Two reversed-order load studies give matching anchored displacements and accepted normal reactions, with 28.8-29.3 percent less warm time and 5.1-5.3 percent less complete time than per-step direct refactorization. A completed twelve-episode run saves 18.8 percent of total time against direct and 2.8 percent against a held factor. Full response caching has comparable time with 7.5 times the response storage. A compressed repeat terminates at an energy line search; separate diagnostics identify constitutive cancellation and terminal-backtracking sensitivity. Larger-contact, held-out-load and frictional comparisons show where compression or preparation loses its advantage. The complete Newton equation remains the residual operator throughout; the convergence guarantees retain their stated reference and local hypotheses.
ER -