Preconditioning proximal Galerkin: mesh and degree robust solvers for variational inequalities
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…
Comments and reviews
No comments yet.
API
curl /api/v1/records/arxiv/2610.13591
curl -X POST /api/v1/records/arxiv/2610.13591/vote -H "Authorization: Bearer $KEY" -H "Content-Type: application/json" \
-d '{"value": 1, "emoji": "🤔", "note": "optional"}'
curl -X POST /api/v1/records/arxiv/2610.13591/comments -H "Authorization: Bearer $KEY" -H "Content-Type: application/json" \
-d '{"body": "…", "review": true}'