Fully Discrete Pointwise Smoothing Error Estimates for Measure Valued Initial Data

Boris Vexler; Dmitriy Leykekhman; Jakob Wagner

arxiv: 2304.13694 · v2 · pith:5W7VGDN6new · submitted 2023-04-26 · 🧮 math.NA · cs.NA

Fully Discrete Pointwise Smoothing Error Estimates for Measure Valued Initial Data

Dmitriy Leykekhman , Boris Vexler , Jakob Wagner This is my paper

Pith reviewed 2026-05-24 09:17 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords parabolic equationsmeasure valued initial datafinite element methoddiscontinuous Galerkininterior error estimatespointwise estimatessmoothing estimates

0 comments

The pith

Interior L^∞ error estimates hold at final time for fully discrete parabolic schemes when measure initial data support is separated from the evaluation point.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper proves interior maximum-norm error bounds at the final time for approximations of parabolic equations whose initial data is a regular Borel measure. Time discretization uses discontinuous Galerkin of arbitrary degree and space discretization uses continuous linear or quadratic finite elements. The estimates rely on the fact that the heat kernel smooths the solution enough to control pointwise errors inside the domain when the measure is supported strictly away from the point of interest. A reader would care because the result justifies pointwise accuracy for singular sources such as Dirac measures without requiring global regularity. The work also extends earlier interior estimates to quadratic elements for square-integrable initial data.

Core claim

We show parabolic smoothing results for the continuous, semidiscrete and fully discrete problems. Our main results are interior L^∞ error estimates for the evaluation at the endtime, in cases where the initial data is supported in a subdomain. We additionally show interior L^∞ error estimates for L² initial data and quadratic finite elements.

What carries the argument

Separation of the support of the initial measure from the interior evaluation point, which activates parabolic smoothing to produce the regularity needed for the interior L^∞ bounds.

If this is right

Pointwise values at interior locations can be computed reliably at the final time even when the initial data is a singular measure.
The same interior estimates now cover quadratic finite elements for square-integrable initial data.
Smoothing properties carry over from the continuous problem to both the semidiscrete and fully discrete schemes.
The estimates remain valid for arbitrary polynomial degree in the discontinuous Galerkin time discretization.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The separation argument could be tested numerically by letting the support approach the evaluation point and observing when the interior bound breaks.
Similar interior estimates might hold for other time discretizations that preserve parabolic smoothing.
The technique could be applied to error analysis in optimal control problems whose controls are measures.

Load-bearing premise

The initial measure is supported in a subdomain strictly separated from the interior point where the error is evaluated.

What would settle it

Numerical computation of the fully discrete solution at an interior point with fixed separation distance between support and evaluation point, showing that the L^∞ error fails to converge at the rate predicted by the estimates.

read the original abstract

In this paper we analyze a homogeneous parabolic problem with initial data in the space of regular Borel measures. The problem is discretized in time with a discontinuous Galerkin scheme of arbitrary degree and in space with continuous finite elements of orders one or two. We show parabolic smoothing results for the continuous, semidiscrete and fully discrete problems. Our main results are interior $L^\infty$ error estimates for the evaluation at the endtime, in cases where the initial data is supported in a subdomain. In order to obtain these, we additionally show interior $L^\infty$ error estimates for $L^2$ initial data and quadratic finite elements, which extends the corresponding result previously established by the authors for linear finite elements.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Solid incremental extension of interior L^∞ estimates to quadratic elements and fully discrete DG+FE schemes for measure initial data.

read the letter

This paper extends the authors' earlier interior estimates to quadratic finite elements for L² data and to a fully discrete discontinuous Galerkin time plus continuous FE space scheme when the initial data is a measure. They first establish parabolic smoothing for the continuous, semidiscrete, and fully discrete problems, then derive the interior L^∞ bounds at final time under the assumption that the measure support is strictly separated from the evaluation subdomain. The quadratic extension and the measure-data fully discrete case are the actual new pieces. The approach is the standard one: use the smoothing property of the heat operator to recover enough regularity away from the initial support, then apply approximation theory. That part holds up. The central limitation is the separation requirement on the support. Without it the interior estimates do not apply, which is mathematically expected but restricts the result to localized singularities. The paper contains no numerical tests, which is normal for this style of analysis but leaves sharpness of constants unchecked. This is narrow technical work aimed at specialists already doing error analysis for parabolic problems with low-regularity data. Those readers will find the precise statements and the quadratic extension useful as a reference. The claims are grounded in established techniques and the extension is genuine, so the paper deserves peer review.

Referee Report

0 major / 2 minor

Summary. The paper analyzes the homogeneous parabolic problem with initial data in the space of regular Borel measures. It discretizes in time via discontinuous Galerkin of arbitrary degree and in space via continuous finite elements of order one or two. Parabolic smoothing estimates are established for the continuous, semidiscrete, and fully discrete problems. The central results are interior L^∞ error estimates at the final time when the initial measure is supported in a subdomain strictly separated from the evaluation region; an auxiliary result extends interior L^∞ estimates to L² initial data with quadratic elements.

Significance. If the derivations hold, the work supplies rigorous pointwise interior error bounds that exploit parabolic smoothing for singular (measure) data, extending prior linear-element results to quadratic elements and the fully discrete setting. Such estimates are useful for applications involving localized sources or point evaluations in parabolic models.

minor comments (2)

The abstract states that the initial-measure support is 'strictly separated' from the interior evaluation subdomain; the precise geometric condition (distance lower bound) should be stated explicitly in the main theorem statements.
Notation for the DG time-stepping polynomial degree and the finite-element mesh size should be introduced uniformly in the introduction and used consistently in all error statements.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary of our work on fully discrete interior L^∞ error estimates for parabolic problems with measure-valued initial data, the recognition of its significance, and the recommendation of minor revision. No specific major comments appear in the report.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper derives interior L^∞ error estimates for a parabolic problem with measure-valued initial data supported away from the evaluation subdomain, using parabolic smoothing properties of the continuous operator, standard DG time discretization, and continuous FE approximation theory for linear and quadratic elements. The extension to quadratic elements for L² data builds on the authors' prior result for linear elements, but this is a normal incremental mathematical extension rather than a load-bearing self-citation that reduces the new claims to unverified inputs. No derivation step reduces by construction to a fitted parameter, self-definition, or renaming of known results; the analysis remains self-contained against external benchmarks of parabolic regularity and FE error theory.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on standard existence and regularity theory for parabolic equations with measure data plus approximation properties of DG and FE spaces; no new free parameters or invented entities are introduced.

axioms (1)

domain assumption The homogeneous parabolic problem with regular Borel measure initial data admits solutions possessing the parabolic smoothing properties needed for interior L^∞ estimates.
Invoked to justify the continuous, semidiscrete, and fully discrete smoothing results stated in the abstract.

pith-pipeline@v0.9.0 · 5650 in / 1256 out tokens · 42471 ms · 2026-05-24T09:17:50.573254+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

[1]

R. A. A DAMS AND J. F OURNIER , Cone conditions and properties of Sobolev spaces , J. Math. Anal. Appl., 61 (1977), pp. 713–734

work page 1977
[2]

R. A. A DAMS AND J. J. F. F OURNIER , Sobolev spaces , vol. 140 of Pure and Applied Mathematics (Amsterdam), Else vier/Academic Press, Amsterdam, second ed., 2003

work page 2003
[3]

A SHYRALYEV AND P

A. A SHYRALYEV AND P. E. S OBOLEVSKII , W ell-Posedness of Parabolic Difference Equations, Birkh¨ auser Basel, 1994

work page 1994
[4]

S. C. B RENNER AND L. R. S COTT , The mathematical theory of ﬁnite element methods , vol. 15 of Texts in Applied Mathematics, Springer, New Y ork, third ed., 2008

work page 2008
[5]

C ASAS , B

E. C ASAS , B. V EXLER , AND E. Z UAZUA , Sparse initial data identiﬁcation for parabolic PDE and its ﬁnite element approximations , Math. Control Relat. Fields, 5 (2015), pp. 377–399

work page 2015
[6]

D RELICHMAN , R

I. D RELICHMAN , R. G. D UR ´AN , AND I. O JEA , A weighted setting for the numerical approximation of the po isson problem with singular sources, SIAM J. Numer. Anal., 58 (2020), pp. 590–606

work page 2020
[7]

E RIKSSON , C

K. E RIKSSON , C. J OHNSON , AND S. L ARSSON , Adaptive ﬁnite element methods for parabolic problems. VI. Analytic semigroups, SIAM J. Numer. Anal., 35 (1998), pp. 1315–1325

work page 1998
[8]

L. C. E VANS, Partial differential equations, no. 19 in Graduate studies in mathematics, American Mathem atical Society, second ed., 2010

work page 2010
[9]

H ANSBO , Strong stability and non-smooth data error estimates for di scretizations of linear parabolic problems, BIT, 42 (2002), pp

A. H ANSBO , Strong stability and non-smooth data error estimates for di scretizations of linear parabolic problems, BIT, 42 (2002), pp. 351–379

work page 2002
[10]

H ELL , A

T. H ELL , A. O STERMANN , AND M. S ANDBICHLER , Modiﬁcation of dimension-splitting methods—overcoming the order reduction due to corner singularities, IMA J. Numer. Anal., 35 (2015), pp. 1078–1091

work page 2015
[11]

L ESAINT AND P

P. L ESAINT AND P. R AVIART , On a ﬁnite element method for solving the neutron transport e quation, in Mathematical Aspects of Finite Elements in Partial Differential Equations, Academic Press, pp. 89– 123

work page
[12]

L EYKEKHMAN AND B

D. L EYKEKHMAN AND B. V EXLER , Pointwise best approximation results for Galerkin ﬁnite el ement solutions of parabolic problems , SIAM J. Numer. Anal., 54 (2016), pp. 1365–1384

work page 2016
[13]

Math., 135 (2017), pp

, Discrete maximal parabolic regularity for Galerkin ﬁnite e lement methods, Numer. Math., 135 (2017), pp. 923–952

work page 2017
[14]

L EYKEKHMAN , B

D. L EYKEKHMAN , B. V EXLER , AND J. W AGNER , Numerical analysis of sparse initial data identiﬁcation fo r parabolic problems with pointwise ﬁnal time observations, in preparation, (2023)

work page 2023
[15]

L EYKEKHMAN , B

D. L EYKEKHMAN , B. V EXLER , AND D. W ALTER , Numerical analysis of sparse initial data identiﬁcation fo r parabolic problems, ESAIM Math. Model. Numer. Anal., 54 (2020), pp. 1139–1180

work page 2020
[16]

M EIDNER , R

D. M EIDNER , R. R ANNACHER , AND B. V EXLER , A priori error estimates for ﬁnite element discretizations of parabolic optimization problems with pointwise state constraints in time , SIAM J. Control Optim., 49 (2011), pp. 1961–1997

work page 2011
[17]

E. B. S AFF AND R. S. V ARGA , On the zeros and poles of pad´ e approximants to eˆz, Numerische Mathematik, 25 (1975), pp. 1–14

work page 1975
[18]

A. H. S CHATZ , V. C. T HOM ´EE , AND L. B. W AHLBIN , Maximum norm stability and error estimates in parabolic ﬁnite element equations, Comm. Pure Appl. Math., 33 (1980), pp. 265–304

work page 1980
[19]

A. H. S CHATZ AND L. B. W AHLBIN , Interior maximum norm estimates for ﬁnite element methods , Math. Comp., 31 (1977), pp. 414–442

work page 1977
[20]

T HOM ´EE, Galerkin ﬁnite element methods for parabolic problems , vol

V. T HOM ´EE, Galerkin ﬁnite element methods for parabolic problems , vol. 25 of Springer Series in Computational Mathematics, S pringer-V erlag, Berlin, second ed., 2006

work page 2006
[21]

W ANNER , E

G. W ANNER , E. H AIRER , AND S. P. N ØRSETT , Order stars and stability theorems, BIT, 18 (1978), pp. 475–489

work page 1978

[1] [1]

R. A. A DAMS AND J. F OURNIER , Cone conditions and properties of Sobolev spaces , J. Math. Anal. Appl., 61 (1977), pp. 713–734

work page 1977

[2] [2]

R. A. A DAMS AND J. J. F. F OURNIER , Sobolev spaces , vol. 140 of Pure and Applied Mathematics (Amsterdam), Else vier/Academic Press, Amsterdam, second ed., 2003

work page 2003

[3] [3]

A SHYRALYEV AND P

A. A SHYRALYEV AND P. E. S OBOLEVSKII , W ell-Posedness of Parabolic Difference Equations, Birkh¨ auser Basel, 1994

work page 1994

[4] [4]

S. C. B RENNER AND L. R. S COTT , The mathematical theory of ﬁnite element methods , vol. 15 of Texts in Applied Mathematics, Springer, New Y ork, third ed., 2008

work page 2008

[5] [5]

C ASAS , B

E. C ASAS , B. V EXLER , AND E. Z UAZUA , Sparse initial data identiﬁcation for parabolic PDE and its ﬁnite element approximations , Math. Control Relat. Fields, 5 (2015), pp. 377–399

work page 2015

[6] [6]

D RELICHMAN , R

I. D RELICHMAN , R. G. D UR ´AN , AND I. O JEA , A weighted setting for the numerical approximation of the po isson problem with singular sources, SIAM J. Numer. Anal., 58 (2020), pp. 590–606

work page 2020

[7] [7]

E RIKSSON , C

K. E RIKSSON , C. J OHNSON , AND S. L ARSSON , Adaptive ﬁnite element methods for parabolic problems. VI. Analytic semigroups, SIAM J. Numer. Anal., 35 (1998), pp. 1315–1325

work page 1998

[8] [8]

L. C. E VANS, Partial differential equations, no. 19 in Graduate studies in mathematics, American Mathem atical Society, second ed., 2010

work page 2010

[9] [9]

H ANSBO , Strong stability and non-smooth data error estimates for di scretizations of linear parabolic problems, BIT, 42 (2002), pp

A. H ANSBO , Strong stability and non-smooth data error estimates for di scretizations of linear parabolic problems, BIT, 42 (2002), pp. 351–379

work page 2002

[10] [10]

H ELL , A

T. H ELL , A. O STERMANN , AND M. S ANDBICHLER , Modiﬁcation of dimension-splitting methods—overcoming the order reduction due to corner singularities, IMA J. Numer. Anal., 35 (2015), pp. 1078–1091

work page 2015

[11] [11]

L ESAINT AND P

P. L ESAINT AND P. R AVIART , On a ﬁnite element method for solving the neutron transport e quation, in Mathematical Aspects of Finite Elements in Partial Differential Equations, Academic Press, pp. 89– 123

work page

[12] [12]

L EYKEKHMAN AND B

D. L EYKEKHMAN AND B. V EXLER , Pointwise best approximation results for Galerkin ﬁnite el ement solutions of parabolic problems , SIAM J. Numer. Anal., 54 (2016), pp. 1365–1384

work page 2016

[13] [13]

Math., 135 (2017), pp

, Discrete maximal parabolic regularity for Galerkin ﬁnite e lement methods, Numer. Math., 135 (2017), pp. 923–952

work page 2017

[14] [14]

L EYKEKHMAN , B

D. L EYKEKHMAN , B. V EXLER , AND J. W AGNER , Numerical analysis of sparse initial data identiﬁcation fo r parabolic problems with pointwise ﬁnal time observations, in preparation, (2023)

work page 2023

[15] [15]

L EYKEKHMAN , B

D. L EYKEKHMAN , B. V EXLER , AND D. W ALTER , Numerical analysis of sparse initial data identiﬁcation fo r parabolic problems, ESAIM Math. Model. Numer. Anal., 54 (2020), pp. 1139–1180

work page 2020

[16] [16]

M EIDNER , R

D. M EIDNER , R. R ANNACHER , AND B. V EXLER , A priori error estimates for ﬁnite element discretizations of parabolic optimization problems with pointwise state constraints in time , SIAM J. Control Optim., 49 (2011), pp. 1961–1997

work page 2011

[17] [17]

E. B. S AFF AND R. S. V ARGA , On the zeros and poles of pad´ e approximants to eˆz, Numerische Mathematik, 25 (1975), pp. 1–14

work page 1975

[18] [18]

A. H. S CHATZ , V. C. T HOM ´EE , AND L. B. W AHLBIN , Maximum norm stability and error estimates in parabolic ﬁnite element equations, Comm. Pure Appl. Math., 33 (1980), pp. 265–304

work page 1980

[19] [19]

A. H. S CHATZ AND L. B. W AHLBIN , Interior maximum norm estimates for ﬁnite element methods , Math. Comp., 31 (1977), pp. 414–442

work page 1977

[20] [20]

T HOM ´EE, Galerkin ﬁnite element methods for parabolic problems , vol

V. T HOM ´EE, Galerkin ﬁnite element methods for parabolic problems , vol. 25 of Springer Series in Computational Mathematics, S pringer-V erlag, Berlin, second ed., 2006

work page 2006

[21] [21]

W ANNER , E

G. W ANNER , E. H AIRER , AND S. P. N ØRSETT , Order stars and stability theorems, BIT, 18 (1978), pp. 475–489

work page 1978