Global error estimators for parametric monotone nonlinearities and neural approximations
Abstract
We construct computable error estimators, which double as loss functions for neural networks, for a class of parametric nonlinear partial differential equations with a monotonicity property, and prove that they are globally reliable and efficient. The value of such a loss function is bounded above and below by the squared error in the natural trial norm, for every trial function, not merely for those near the exact solution; this global property rests on monotonicity. The construction rests on splitting the nonlinear operator into a linear part and a strongly monotone closure, and on measuring the linear residual in a discrete dual norm. Since the closure contributes a dual norm that admits no closed form when the trial norm is stronger than an L2 norm, the estimator is built around a computable surrogate for it, required only to satisfy a pairing bound and a Lipschitz bound. The main theorem then yields two-sided bounds with explicit constants and covers two instances. A first-order system least-squares estimator on conforming trial spaces is the first instance studied: its two-sided bound holds on the whole trial space and therefore applies to arbitrary approximations, including those not in a discrete finite element space. The second instance is a discontinuous Petrov-Galerkin estimator on trial spaces of finite element functions for each parameter value, built with broken test spaces, for which no conformity is required and the dual norm is computed by independent element-local problems. All assumptions are verified for a model class of nonlinear fluxes, with the constants tracked explicitly in terms of the parameter range. Being computable for an arbitrary input, both estimators serve as variationally correct loss functions for neural network approximations of parameter-to-solution maps.
Related
- Locally computable error estimators for conforming approximations of interface problems cannot be robust
- Error Bounds for Statistical Estimators in BTL Model with Parametric Multivariate Utility Functions
- Variationally correct loss functions for nonlinear parametric problems and their neural approximations
- Learning Parametric Monotone Games
- Computational study of Proper Orthogonal Decomposition methods for parametric approximations
BibTeX
@misc{castillo2026global,
title = {{Global error estimators for parametric monotone nonlinearities and neural approximations}},
author = {Pablo Cortés Castillo and Wolfgang Dahmen and Jay Gopalakrishnan},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.05767},
eprint = {2610.05767},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.05767},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Global error estimators for parametric monotone nonlinearities and neural approximations
AU - Pablo Cortés Castillo
AU - Wolfgang Dahmen
AU - Jay Gopalakrishnan
PY - 2026
DA - 2026/10/05
PB - arXiv
DO - 10.48550/arxiv.2610.05767
UR - https://arxiv.org/abs/2610.05767
AB - We construct computable error estimators, which double as loss functions for neural networks, for a class of parametric nonlinear partial differential equations with a monotonicity property, and prove that they are globally reliable and efficient. The value of such a loss function is bounded above and below by the squared error in the natural trial norm, for every trial function, not merely for those near the exact solution; this global property rests on monotonicity. The construction rests on splitting the nonlinear operator into a linear part and a strongly monotone closure, and on measuring the linear residual in a discrete dual norm. Since the closure contributes a dual norm that admits no closed form when the trial norm is stronger than an L_2 norm, the estimator is built around a computable surrogate for it, required only to satisfy a pairing bound and a Lipschitz bound. The main theorem then yields two-sided bounds with explicit constants and covers two instances. A first-order system least-squares estimator on conforming trial spaces is the first instance studied: its two-sided bound holds on the whole trial space and therefore applies to arbitrary approximations, including those not in a discrete finite element space. The second instance is a discontinuous Petrov-Galerkin estimator on trial spaces of finite element functions for each parameter value, built with broken test spaces, for which no conformity is required and the dual norm is computed by independent element-local problems. All assumptions are verified for a model class of nonlinear fluxes, with the constants tracked explicitly in terms of the parameter range. Being computable for an arbitrary input, both estimators serve as variationally correct loss functions for neural network approximations of parameter-to-solution maps.
ER -