Some error estimates for the DEC method in the plane

Miguel Angel Moreles; Rafael Herrera; Ruben Carrillo

arxiv: 1907.06762 · v1 · pith:RLP6IGRKnew · submitted 2019-07-15 · 🧮 math.NA · cs.NA

Some error estimates for the DEC method in the plane

Ruben Carrillo , Miguel Angel Moreles , Rafael Herrera This is my paper

Pith reviewed 2026-05-24 21:04 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords Discrete Exterior Calculusbox methodPoisson equationerror estimatesfinite element methodconvergence rates

0 comments

The pith

The Discrete Exterior Calculus method for the Poisson problem equals the box method on the same mesh.

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

This paper shows that the Discrete Exterior Calculus method applied to the Poisson equation in two dimensions is identical to the box method when using the same mesh. Because of this, the error analysis already developed for the box method transfers directly to DEC. The result is that DEC has the same error estimates as the finite element method with linear elements. This matters to readers who want to use DEC but need proven bounds on its accuracy rather than relying on numerical experiments alone.

Core claim

We show that the Discrete Exterior Calculus (DEC) method can be cast as the earlier box method for the Poisson problem in the plane. Consequently, error estimates are established, proving that the DEC method is comparable to the Finite Element Method with linear elements. We also discuss some virtues, others than convergence, of the DEC method.

What carries the argument

Casting the DEC discretization as the box method on identical meshes, which transfers all prior error analysis without modification.

If this is right

The same error estimates hold for DEC as for the box method.
DEC is comparable in accuracy to linear finite elements.
No additional error terms arise from the DEC formulation.

Where Pith is reading between the lines

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

The same approach might identify DEC with other methods in different settings.
Users of DEC can draw on the broader literature on box methods for theoretical results.
The unification suggests focusing on implementation differences rather than theoretical ones for this problem.

Load-bearing premise

The DEC discretization produces exactly the same linear system as the box method on the same mesh.

What would settle it

Explicitly constructing and comparing the linear systems from both methods on a single triangular mesh to check for equality.

read the original abstract

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

The paper recasts DEC for the 2D Poisson problem as the box method on the same mesh so that existing box-method error bounds transfer directly and match linear FEM rates.

read the letter

The central observation is that DEC for -Δu = f in the plane produces the same linear system as the box method, so prior error analysis carries over and shows DEC is comparable to linear elements. That equivalence is the new piece; the abstract frames the box method as earlier work and uses the link to supply estimates that were missing for DEC. The authors also note other practical features of DEC beyond convergence rates, which is a reasonable addition for readers already using the method. The argument stays non-circular by leaning on independent box-method theory. The load-bearing step is whether the two schemes really agree on the matrix and the load vector for general triangulations, including how dual cells and the source term are integrated. If the paper proves an exact match rather than an approximate one, the transfer works; any hidden scaling or boundary correction would require extra justification. The abstract alone leaves that detail open, but the claim is stated cleanly enough that a referee can check the derivation. This is narrow-scope numerical analysis aimed at people who already care about DEC or finite-volume schemes and want convergence rates without starting from scratch. It is not opening new domains, but the equivalence is the kind of concrete bridge that saves others time. The thinking is straightforward and the literature framing looks honest. Send it to peer review so the equivalence can be verified in detail.

Referee Report

2 major / 2 minor

Summary. The paper claims that the Discrete Exterior Calculus (DEC) discretization of the Poisson problem −Δu = f in the plane can be exactly recast as the box (finite-volume) method on the same primal triangulation. This equivalence permits the direct transfer of existing error estimates for the box method, establishing that DEC achieves the same O(h) convergence in the H1 norm and O(h²) in L2 as the linear finite-element method. The manuscript also discusses additional practical virtues of DEC beyond convergence rates.

Significance. If the claimed identity between the two linear systems holds without residual terms, the result supplies a short route to a priori error bounds for DEC without developing new analysis from scratch. It also positions DEC as a geometrically natural method whose convergence theory is already covered by classical box-method results, which may be useful for applications on unstructured meshes where DEC's built-in preservation of Stokes' theorem is advantageous.

major comments (2)

[§3] §3 (equivalence of the discrete operators): the central claim that the DEC stiffness matrix and load vector are identical to those of the box method is asserted but the explicit algebraic verification that the dual-cell integration of f and the discrete Hodge star applied to the gradient produce exactly the same entries (including on boundary dual cells) is not supplied. Without this identity, the error estimates do not transfer verbatim.
[Theorem 4.1] Theorem 4.1 (error estimate transfer): the statement that the DEC solution satisfies the same a priori bounds as the box method assumes the two schemes produce the identical linear system on every admissible triangulation; any difference in the treatment of the source term on primal edges incident to the boundary would introduce an additional consistency error that must be controlled separately.

minor comments (2)

[Abstract] The abstract states that error estimates are 'established' but the manuscript only transfers existing bounds; a short remark clarifying that no new analysis is performed would avoid overstatement.
[§2] Notation for the dual mesh and the averaging operator used to define the box-method right-hand side should be introduced once and used consistently when the equivalence is proved.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. Below we provide point-by-point responses to the major comments and indicate the revisions we will make.

read point-by-point responses

Referee: [§3] §3 (equivalence of the discrete operators): the central claim that the DEC stiffness matrix and load vector are identical to those of the box method is asserted but the explicit algebraic verification that the dual-cell integration of f and the discrete Hodge star applied to the gradient produce exactly the same entries (including on boundary dual cells) is not supplied. Without this identity, the error estimates do not transfer verbatim.

Authors: We acknowledge that an explicit verification of the algebraic identity would improve the clarity of the equivalence. In the revised manuscript, we will include a detailed derivation showing that the DEC discretization yields the same stiffness matrix and load vector as the box method, with particular attention to the boundary dual cells. This will confirm that the linear systems are identical. revision: yes
Referee: [Theorem 4.1] Theorem 4.1 (error estimate transfer): the statement that the DEC solution satisfies the same a priori bounds as the box method assumes the two schemes produce the identical linear system on every admissible triangulation; any difference in the treatment of the source term on primal edges incident to the boundary would introduce an additional consistency error that must be controlled separately.

Authors: Since the equivalence in §3 establishes that the two methods produce identical linear systems on any admissible triangulation, including the discretization of the source term, the error estimates transfer directly. We will revise the statement of Theorem 4.1 and its proof to explicitly reference this identity and note that no additional consistency error arises from boundary treatment. revision: yes

Circularity Check

0 steps flagged

No circularity; equivalence proof transfers independent error analysis

full rationale

The paper's central step is an explicit algebraic identification of the DEC discretization (primal-dual operators, integration over dual cells, and load vector) with the box/finite-volume scheme on the same triangulation. Once this identity is established by direct comparison of the resulting linear systems, the error estimates for the box method (treated as prior independent work) apply verbatim. No parameter is fitted and then relabeled as a prediction, no ansatz is smuggled via self-citation, and the uniqueness or correctness of the equivalence is not justified by citing the authors' own earlier results. The derivation therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based on abstract only; no free parameters, invented entities, or non-standard axioms are mentioned. The result relies on standard properties of the Poisson problem and the box method.

axioms (1)

domain assumption The box method possesses established a priori error estimates for the Poisson problem on suitable meshes in the plane.
Invoked when the paper states that error estimates follow 'consequently' from the casting.

pith-pipeline@v0.9.0 · 5567 in / 1268 out tokens · 38823 ms · 2026-05-24T21:04:26.802384+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We show that the Discrete Exterior Calculus (DEC) method can be cast as the earlier box method for the Poisson problem in the plane. Consequently, error estimates are established...
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The box method consists on finding uB∈P10(T) such that ā(uB,v̄)=(f,v̄)... stiffness matrix is identical to that arising from the Finite Element Method...

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.