Locally-averaged McCormick relaxations for discretization-regularized inverse problems
Barbara Kaltenbacher, Paul Manns
arXiv:2601.01995·math.NA·Published 2026-01-05·Updated 2026-03-19
In this paper, by means of a standard model problem, we devise an approach to computing approximate dual bounds for use in global optimization of coefficient identification in partial differential equations (PDEs) by, e.g., (spatial) branch-and-bound methods. Linearization is achieved by a McCormick relaxation (that is, replacing the bilinear PDE term by a linear one and adding inequality constraints), combined with local averaging to reduce the number of inequalities. Optimization-based bound tightening allows us to tighten the relaxation and thus reduce the induced error. Combining this with a quantification of the discretization error and the propagated noise, we prove that the resulting discretization regularizes the inverse problem, thus leading to an overall convergent scheme. Numerical experiments illustrate the theoretical findings.
TopicsGenerative Models & Discovery
Tagspartial-differential-equations
arXiv categoriesmath.NA
arXiv abstract pagePDF