{"id":"8e6d946a-af32-45c3-a2f4-f53eef4e8427","arxiv_id":"2411.13617","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Derives computable maximum-norm a posteriori error bounds for the extrapolated Euler method with FEM for linear parabolic equations, for any extrapolation order L.","lead":"This paper proves maximum-norm a posteriori error bounds for parabolic equations solved by finite elements in space and Richardson-extrapolated backward Euler in time. The result extends an earlier second-order analysis to extrapolation of arbitrary order, a step toward reliable adaptive time stepping.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 1's second constant is off by a factor of 4 (κ1′τ_j^2/24 should be κ1′τ_j^2/6), so Theorem 1's η_t is not a guaranteed bound as written; the theorem also mishandles the L=2 case.","rationale":"The reader identified Assumption 1 as the weakest assumption, and I agree that the unproved generality of the elliptic estimator limits applicability. However, the more immediate blocker is internal: Lemma 1's stated second bound is algebraically wrong by a factor of four, and the min structure of μ_{j,i} makes this an actual invalidation rather than a harmless typo. The theorem also uses an undefined η_ψ/η_Ψ and the L=2 limit of 1/(2L−5) is inconsistent with the derivation in §4.4.2. These are not disagreements with external consensus; they are localized correctness failures. The approach may be salvageable by correcting the constant, adding a special L=2 case, and cleaning up the theorem statement, but the central claim as stated is not established. Hence the verdict should move from CONDITIONAL to REJECT for this version, with the expectation that a corrected version could be reconsidered.","tokens_in":10172,"tokens_out":18137,"duration_ms":179295,"concrete_test":"Independently re-derive the last step of Lemma 1 with κ1=0, κ1′=1, γ=0 (so φ1≡1): the proof's own displayed bound yields ‖v‖∞ τ_j^2/6, while the lemma states ‖v‖∞ τ_j^2/24. If the discrepancy reproduces, Theorem 1's η_t is not a guaranteed bound and the central claim is invalid as written.","verdict_should_be":"REJECT","load_bearing_attack":"In the proof of Lemma 1, after integration by parts the proof reaches the bound |∫_{I_j}⟨v,G(T−s)⟩P_i(ξ_j(s))ds| ≤ ‖v‖∞ σ_j [ (κ1/τ_j)∫_{I_j} (t_j−s)(s−t_{j−1})/(T−s) ds + (κ1′/τ_j)∫_{I_j} (t_j−s)(s−t_{j−1}) ds ]. Since ∫_{t_{j−1}}^{t_j} (t_j−s)(s−t_{j−1}) ds = τ_j^3/6, the second term equals κ1′ τ_j^2/6. The lemma states κ1′ τ_j^2/24 instead. Because μ_{j,i} is defined as a minimum of this estimate and the first estimate, the erroneous 24 can make μ_{j,i} four times too small, and Theorem 1's η_t, built from μ_{j,L−2} and μ_{j,L−1}, is then not a guaranteed upper bound. This is an internal algebraic inconsistency, independent of external Assumption 1. Two textual inconsistencies reinforce the problem: Theorem 1 announces η_ψ + η_Ψ although only η_t is defined; and for L=2 the coefficient 1/(2L−5)=−1 contradicts the P0 coefficient +Δ0ψ_j+2Δ1u_j^h/τ_j in (21), so the theorem needs a separate L=2 case.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a framework for guaranteed maximum-norm a posteriori error bounds for linear parabolic initial--boundary value problems discretised by the L-step Richardson extrapolation of the backward Euler method in time and conforming FEM in space. The final-time error is represented via the parabolic Green's function; the spatial error is controlled by elliptic reconstructions together with a generic maximum-norm elliptic estimator (Assumption 1), and the temporal error is controlled by a piecewise polynomial reconstruction of degree L-1 built from integrated Legendre polynomials, chosen so that the residual psi + partial_t u_h collapses to its two highest Legendre modes (Eq. (22)). Lemma 1 bounds the resulting mode contributions, and Theorem 1 assembles the claimed computable bound eta_M = eta_init + eta_f + eta_t + eta_ell. The overall strategy is coherent and the displayed computations in the proof of Lemma 1 are essentially correct, but the printed statements of Lemma 1 and Theorem 1 contain defects that, as printed, invalidate the guaranteed-bound claim and require correction.","tokens_in":10513,"tokens_out":26508,"duration_ms":233440,"significance":"If corrected, the paper gives a useful and reusable extension of the authors' second-order analysis [19] to extrapolation of arbitrary order L >= 2. The Legendre-bubble reconstruction with its recursive coefficient choice (Section 4.4.2) is explicit, elegant, and implementable; the estimator contains no fitted parameters and is derived from the PDE, the scheme, and two clearly stated external assumptions, so the approach is genuinely a posteriori and not circular. The main theorem is, however, conditional on an elliptic estimator whose verified instances in the literature are narrower than the paper's advertised scope, which limits the practical reach but not the value of the framework once the defects below are fixed. The manuscript contains no numerical experiments and would be materially strengthened by one simple test problem with reported effectivity indices.","major_comments":[{"comment":"The constant multiplying kappa_1' in the second bound of Lemma 1 is printed as kappa_1' tau_j^2/24, both in the lemma statement and in the displayed bound (24), but the proof's own computation gives kappa_1' tau_j^2/6: after the substitution ||P_i'||_\\infty = i(i+1)/2, the bound contains (kappa_1'/tau_j) \\int_{I_j} (t_j-s)(s-t_{j-1}) ds, and since this integral equals tau_j^3/6, the surviving constant is kappa_1' tau_j^2/6. Because mu_{j,i} is defined as the minimum of the two candidate bounds, the too-small printed constant can make mu_{j,i}, and hence the temporal estimator eta_t in Theorem 1, four times smaller than the true guaranteed bound. This error is load-bearing: with the printed constant, Theorem 1 does not assert a valid upper bound. The first term in the same displayed bound should also be made consistent between the lemma statement and the proof, which yields (kappa_1/tau_j) \\int ... (the typeset formula may be read as kappa_1 tau_j \\int ...).","section":"Section 4.4.2, Lemma 1 and Eq. (24)"},{"comment":"The displayed definition of eta_M reads eta_M := eta_init + eta_f + eta_psi + eta_Psi + eta_ell, but the symbols eta_psi and eta_Psi are never defined; the subsequent list defines only eta_init, eta_f, eta_t, and eta_ell. As printed, the equality defining eta_M is unverifiable. The intended temporal term is evidently eta_t, and the display should say so explicitly.","section":"Section 4.4.4, Theorem 1"},{"comment":"The formula for eta_t contains the coefficient 1/(2L-5) in front of Delta_{L-2} psi_j. For L = 2 this coefficient equals -1, whereas the P0 coefficient of psi + partial_t u_h for L = 2 in Eq. (21) is +Delta_0 psi_j + 2 Delta_1 u_j^h / tau_j; moreover, the reduction to the two-mode form (22) is derived only for L > 2, since the recursion defining Delta_2 psi_j does not exist when L = 2. Because the paper explicitly claims L >= 2, Theorem 1 either needs a separate statement for L = 2 based on (21), or should be restricted to L >= 3 with the L = 2 case referred to the authors' earlier treatment [19].","section":"Section 4.4.4, Theorem 1 (L = 2 case)"},{"comment":"The Introduction states that 'any finite-element method might be used' in space, but every eta_ell term of Theorem 1 rests on Assumption 1, which postulates a maximum-norm elliptic estimator for the discrete problem (4) with, in general, quadrature-perturbed forms a_h(.,.) and (.,.)_h. The works cited in Section 4.1 establish such estimators only in restricted settings (specific operators, quasiuniform or anisotropic meshes, and in several cases without quadrature), and the paper proves no such estimator for the general FEM with quadrature described in Section 3. The advertised computable, guaranteed bound is therefore contingent on a premise that is not verified for the method actually analyzed. The authors should either delimit the theorem's scope to FEMs for which Assumption 1 is known to hold or state concrete verified instances, including quadrature if possible.","section":"Section 4.1, Assumption 1; Section 3"}],"minor_comments":[{"comment":"In the definition of eta_f, the summation is printed as \\sum_{m}^{j=1}; it should read \\sum_{j=1}^{M} sigma_j \\int_{I_j} ||(f - \\hat f)(s)||_{\\infty,\\Omega} ds.","section":"Section 4.4.4"},{"comment":"The formula defining the extrapolation weights alpha_ell through the matrix B is garbled in the printed algorithm; please state the weights explicitly, for example as the solution of the Vandermonde system \\sum_{\\ell=1}^{L} alpha_ell / \\ell^k = \\delta_{k0} for k = 0,...,L-1, so that the extrapolation coefficients are unambiguous.","section":"Section 3"},{"comment":"The manuscript contains no numerical experiments. In view of the paper's claim that the bound is computable and of the approximations described in Remark 3, at least one academic test problem with reported effectivity indices would substantially support the practical claims; I recommend adding one during revision.","section":"Numerical examples"},{"comment":"In expansion (18), the coefficients Delta_i phi_j for i >= 2 are introduced only verbally; a sentence stating that they are the Legendre coefficients of the piecewise polynomial phi (or equivalently the coefficients determined by interpolation at the L-2 interior points) would make the reconstruction fully explicit, especially since these coefficients are later determined through (19).","section":"Section 4.4.2"},{"comment":"Typos such as 'differ to some extend' (Introduction) and 'finite dimentional' (Section 3) should be corrected in the revision.","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the defects identified in Lemma 1 and Theorem 1 appear to be genuine slips rather than structural flaws; the proof of Lemma 1 contains the computation that yields the corrected tau_j^2/6 constant, and the L = 2 case is already worked out in Eq. (21), so a revision can plausibly close the gaps. The manuscript builds substantially on the authors' own prior work ([12,13,17,18,19]), but the citations are to published, relevant results and are appropriate for the elliptic-reconstruction programme. The main editorial question is whether to require numerical experiments for a 'computable bound' paper; I would encourage at least one example but would not make it a condition for acceptance. The paper fits the journal's scope well."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the core idea is real: the paper extends the elliptic reconstruction framework to Richardson extrapolation of arbitrary order L, and the systematic construction of degree L-1 temporal reconstructions from nodal values via integrated Legendre polynomials is a genuine contribution, not in the authors' earlier L=2 paper. Second, as printed, Lemma 1 has a factor error that breaks the guarantee in Theorem 1's η_t, and the L=2 case is mishandled.\n\nWhat is actually new and good: the recursive definition of the auxiliary unknowns Δ_i ψ and Δ_i u is clever and nontrivial. The proof structure is coherent, the paper honestly lists its external assumptions, and no fitted parameters or manufactured predictions enter. If the algebraic issues are fixed, this gives a useful computable maximum-norm a posteriori bound for a wide class of parabolic FEM discretisations.\n\nSoft spots, in order of severity. (1) In Lemma 1, after integration by parts the proof reaches 1/τ_j ∫ φ_1(T-s)(t_j-s)(s-t_{j-1}) ds. Since that integral is τ_j^3/6, the second constant should be κ_1' τ_j^2/6, not κ_1' τ_j^2/24. Because μ_{j,i} is defined as a minimum, the too-small printed constant can make μ_{j,i} four times too small, so η_t is not a guaranteed upper bound as written. This is load-bearing, not cosmetic. (2) Theorem 1's formula for η_t uses 1/(2L-5), which for L=2 equals -1 and does not match the P0 coefficient in equation (21); a separate L=2 case is needed. Also, the theorem announces η_ψ + η_Ψ while only η_t is defined. (3) Assumption 1 is a real condition: the elliptic estimator is cited from the literature in special settings (quasiuniform meshes, certain operators), but the paper states it for general FEM with quadrature on arbitrary meshes. That makes the main theorem conditional on an assumption that may not hold in the stated generality. Minor point: there are no numerical experiments, so nothing tests sharpness or verifies the constants.\n\nWho this is for: numerical analysts working on adaptive methods for parabolic problems, especially those interested in maximum-norm control. The framework is worth engaging with, but not in its current form. I would send it to a serious referee—the construction is valuable and the errors seem fixable—but the referee should be asked to verify Lemma 1 and the L=2 case carefully. I would not cite it until the factor error and edge case are corrected.","headline":"The arbitrary-order reconstruction machinery is a genuine extension, but a factor error in Lemma 1 makes the printed η_t not a guaranteed bound until fixed.","tokens_in":11031,"tokens_out":3056,"would_cite":false,"duration_ms":31229,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M15","65M50","65M60"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes a fully computable maximum-norm a posteriori error bound for the L-step Richardson-extrapolated backward Euler method combined with FEM, valid at the final time for linear parabolic equations.","keywords":["parabolic problems","maximum-norm a posteriori error estimates","backward Euler","Richardson extrapolation","elliptic reconstructions","Green's function","finite element method"],"falsifier":"Take a concrete linear parabolic problem with known exact solution, run the L-step extrapolated Euler-FEM scheme, evaluate each term of $\\eta_M$, and compare the true error $\\|u(T)-u_M^h\\|_{\\infty,\\Omega}$ with the estimator; because Theorem 1 claims a guaranteed upper bound, any single computed example with true error larger than $\\eta_M$ would refute the bound.","tokens_in":9971,"feed_emoji":"📐","tokens_out":8278,"duration_ms":77343,"temperature":0.7,"pith_summary":"The paper derives guaranteed, computable upper bounds on the pointwise (maximum-norm) error at the final time for linear parabolic equations solved by Richardson extrapolation of the backward Euler method in time and finite elements in space. The central claim is that, although the extrapolation produces high-order approximations only at temporal mesh points, one can still build a piecewise-polynomial-in-time reconstruction whose residual controls the error, and all ingredients in the bound are evaluated from the computed discrete solution and given data. This matters because maximum-norm a posteriori bounds of this type have previously been available only for low-order or special time discretisations; the paper extends them to extrapolations of arbitrary order. If the bound is correct, an adaptive solver can certify the pointwise accuracy of its final output without knowing the exact solution.","feed_headline":"Extrapolated Euler-FEM solvers get pointwise error certificates","feed_subtitle":"The bound needs no exact solution and works for any order of the time extrapolation.","key_machinery":"The machinery is the elliptic reconstruction $R(t)$, defined on each time slice as the solution of $LR(t)=\\hat f(t)+\\psi(t)$, where $\\psi$ measures the deviation of the discrete equation from the continuous one; the finite element solution $u_h(t)$ is then an approximation of $R(t)$, and the elliptic estimator from Assumption 1 bounds $R-u_h$ pointwise. The second key ingredient is the temporal reconstruction: on each interval the numerical solution is represented in a basis built from integrated Legendre polynomials $N_i$, and the free coefficients are chosen so that $\\psi+\\partial_t u_h$ reduces to a linear combination of the two highest Legendre polynomials on each interval. Lemma 1 then bounds the Green's-function-weighted time integrals of such terms using the assumed decay bounds on $G$ and its derivative.","core_discovery":"The paper's central result is Theorem 1: for the L-step Richardson extrapolation method applied to the backward Euler discretisation, combined with any finite-element discretisation that satisfies the paper's assumptions, the error at the final time obeys $$\\|u(T)-u_M^h\\|_{\\infty,\\$\\Omega$} \\le \\eta_M := \\eta_{\\mathrm{init}} + \\eta_f + \\eta_t + \\eta_{\\mathrm{ell}},$$ where each term is assembled from the computed approximations, the interpolated source term, elliptic reconstructions, an available elliptic a posteriori estimator, and constants depending only on the Green's function of the parabolic operator. The bound is fully computable from the numerical solution and the data, with no reference to the unknown exact solution. The proof uses the Green's function representation of the error and an elliptic reconstruction of the discrete solution, with the temporal reconstruction chosen so that the residual term $\\psi+\\partial_t u_h$ collapses to a sum of two Legendre polynomials on each time interval.","pith_inferences":["The paper does not report numerical experiments; a natural check is to compute the effectivity ratio on problems with known solutions to see how far the estimator is from one.","Because the Legendre-bubble construction is the only method-dependent step, the same reconstruction may extend to Richardson extrapolation of other one-step methods, though the paper does not make that claim.","The estimator's per-interval terms suggest a time-adaptive strategy driven by the temporal contribution and local Green's-function weights, an algorithmic step not developed here.","If sharper Green's-function constants became available for specific operators, the bound would tighten automatically; the paper takes the constants as input data."],"forward_implications":["For any extrapolation order L ≥ 2, the final-time solution carries a pointwise error certificate computed from the discrete solution and data alone.","The four additive terms in the estimator separate initial-data, source-interpolation, temporal, and spatial contributions, so an adaptive procedure can target the dominant source of error.","Evaluating the estimate requires only the mesh-point values of the extrapolated solution, auxiliary extrapolation values, the interpolated source, and one elliptic estimator per time level; no additional parabolic solves are needed.","For L = 2 the estimator reduces to the previously known second-order framework, which the present construction generalises systematically.","Any improvement in the available maximum-norm elliptic estimators immediately improves the parabolic bound without changing the temporal analysis."],"supporting_citations":[{"why":"Introduces the elliptic reconstruction idea used to convert the discrete solution into a continuous function whose residual is controlled pointwise.","marker":"[20]"},{"why":"Supplies the Green's function representation of the error that is the starting point of the bound.","marker":"[8]"},{"why":"Provides the class of parabolic operators for which the assumed Green's function decay bounds hold.","marker":"[13]"},{"why":"Establishes the previous unified framework for second-order time discretisations, which the present arbitrary-order construction extends.","marker":"[19]"},{"why":"An available maximum-norm elliptic a posteriori estimator that satisfies Assumption 1 in a concrete setting.","marker":"[21]"},{"why":"Another maximum-norm elliptic estimator, for singularly perturbed reaction-diffusion problems, giving a further setting where Assumption 1 holds.","marker":"[3]"}],"fun_headline_variants":["Pointwise error bounds for parabolic FEM without exact solution","A posteriori max-norm error control for extrapolated Euler","Computable error certificates for parabolic time-FEM","Legendre time reconstructions enable pointwise error bounds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's bound depends on the unproved assumption that a computable maximum-norm error estimator exists for the chosen finite-element discretisation of the elliptic problem, including the quadrature and mesh actually used.","fun_headline_variants_meta":{"raw":{"variants":["Pointwise error bounds for parabolic FEM without exact solution","A posteriori max-norm error control for extrapolated Euler","Computable error certificates for parabolic time-FEM","Legendre time reconstructions enable pointwise error bounds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000547,"raw_usage":{"total_tokens":2549,"prompt_tokens":817,"completion_tokens":1732,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":433,"completion_tokens_details":{"reasoning_tokens":1668}},"tokens_in":433,"tokens_out":1732,"duration_ms":23901,"temperature":1.0,"reasoning_tokens":1668,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:48:31.389077+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a concrete linear parabolic problem with known exact solution, run the L-step extrapolated Euler-FEM scheme, evaluate each term of $\\eta_M$, and compare the true error $\\|u(T)-u_M^h\\|_{\\infty,\\Omega}$ with the estimator; because Theorem 1 claims a guaranteed upper bound, any single computed example with true error larger than $\\eta_M$ would refute the bound.","supporting_citations":[{"cited_title":"Makridakis and R","cited_arxiv_id":null,"evidence_quote":"Introduces the elliptic reconstruction idea used to convert the discrete solution into a continuous function whose residual is controlled pointwise."},{"cited_title":"Friedman","cited_arxiv_id":null,"evidence_quote":"Supplies the Green's function representation of the error that is the starting point of the bound."},{"cited_title":"Kopteva and T","cited_arxiv_id":null,"evidence_quote":"Provides the class of parabolic operators for which the assumed Green's function decay bounds hold."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the previous unified framework for second-order time discretisations, which the present arbitrary-order construction extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"An available maximum-norm elliptic a posteriori estimator that satisfies Assumption 1 in a concrete setting."},{"cited_title":"Demlow and N","cited_arxiv_id":null,"evidence_quote":"Another maximum-norm elliptic estimator, for singularly perturbed reaction-diffusion problems, giving a further setting where Assumption 1 holds."}],"review_version":1}