{"id":"arxiv/2610.13591","url":"https://arxiv.org/abs/2610.13591","tags":["math.NA"],"score":0,"title":"Preconditioning proximal Galerkin: mesh and degree robust solvers for variational inequalities","domain":"math.NA","source":"arxiv","authors":["Ioannis P. A. Papadopoulos","Pablo D. Brubeck"],"reviews":0,"subject":"math.NA","summary":"The proximal Galerkin algorithm is a structure-preserving solver for variational inequalities that results in discretization-independent nonlinear iteration counts. However, in large three-dimensional simulations, the arising Newton saddle point systems cannot be feasibly solved by direct methods. To leverage an iterative linear solver, we propose a robust and scalable block preconditioner. Instead of using the usual dual Schur complement that acts on the latent variable, we design a preconditioner for the primal Schur complement. This leads to a block preconditioner that requires two…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-10-09","reactions":{}}