Approximation properties of double complexes

Daniel F{\o}rland Holmen; Jan Martin Nordbotten; Jon Eivind Vatne

arxiv: 2604.13982 · v1 · submitted 2026-04-15 · 🧮 math.NA · cs.NA· math.AP

Approximation properties of double complexes

Daniel F{\o}rland Holmen , Jan Martin Nordbotten , Jon Eivind Vatne This is my paper

Pith reviewed 2026-05-10 12:19 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.AP

keywords de Rham complexesCech-de Rham complexHodge-Laplace problemsmixed-dimensional modelingcochain mapserror estimatesHilbert complexesnumerical approximation

0 comments

The pith

Bounded cochain maps between simplicial and Čech-de Rham complexes produce concrete error estimates for their associated Hodge-Laplace problems.

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

The paper constructs bounded cochain maps connecting the simplicial de Rham complex, which encodes mixed-dimensional spatially coupled problems, to the Čech-de Rham complex, which encodes the same problems treated as equidimensional with thin overlaps. These maps measure the distance between the two formulations in a way that transfers directly to the solutions of the governing Hodge-Laplace equations on each complex. From the maps the authors obtain both a priori and a posteriori error estimates that bound how much accuracy is lost when a problem is modeled in mixed dimensions rather than with the thin-overlap equidimensional version. This matters for any numerical simulation where geometry is simplified by collapsing thin regions to lower-dimensional objects, because the estimates give an explicit way to control the modeling error introduced by that simplification.

Core claim

The simplicial de Rham complex is realized as a subcomplex of the Čech-de Rham complex. Bounded cochain complexes are constructed between them that quantify their closeness; these maps then yield a priori and a posteriori error estimates between the Hodge-Laplace problems defined on each complex. The estimates represent the modeling error incurred by treating a spatially coupled problem as mixed-dimensional instead of as an equidimensional problem with thin overlaps.

What carries the argument

Bounded cochain maps between the simplicial de Rham complex and the Čech-de Rham complex; they serve as the explicit bridges that transfer norms and allow error estimates to pass from one complex to the other.

If this is right

The difference between solutions of the Hodge-Laplace problem on the simplicial complex and on the Čech complex is controlled by the norms of the cochain maps.
Both pre-computed (a priori) and computable (a posteriori) bounds on this difference are available once the maps are known.
Any problem whose weak form is governed by the Hodge-Laplace operator on either complex inherits these quantitative closeness statements.
The construction applies uniformly to bigraded Hilbert complexes of this type, so the same error framework covers a range of mixed-dimensional and equidimensional spatial-coupling problems.

Where Pith is reading between the lines

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

The maps could be used to derive similar error controls for other pairs of complexes that appear in finite-element exterior calculus.
In applications such as thin-layer heat transfer or fracture flow, the estimates could be turned into practical rules for when the mixed-dimensional reduction is accurate enough for a given tolerance.
The same bounded-map technique might extend to time-dependent or nonlinear problems whose spatial part is still described by these complexes.

Load-bearing premise

Explicit bounded cochain maps between the two complexes exist and are sufficiently controlled to produce usable a priori and a posteriori error estimates.

What would settle it

A concrete pair of meshes and a specific Hodge-Laplace problem on which the constructed maps are unbounded or the derived error bounds are violated by computed solutions.

Figures

Figures reproduced from arXiv: 2604.13982 by Daniel F{\o}rland Holmen, Jan Martin Nordbotten, Jon Eivind Vatne.

**Figure 1.** Figure 1: Three left domains: A simple open cover consisting of two open sets, with diameter in the transversal direction equal to a parameter ϵ. Rightmost domain: When ϵ → 0, we get a mixed-dimensional geometry. Given u ∈ C k , inf uE∈Ek ∥u − uE∥C = O(ϵ). (3.11) Proof. We consider general case when the inner products for the Cech-de Rham ˇ complex are weighted by the parameter ϵ, as in eq. (1.20). Unweighted inner … view at source ↗

**Figure 2.** Figure 2: The mixed-dimensional geometry consists of two adjacent intervals, while the Cech open cover consists of two overlapping ˇ intervals. The Hodge-Laplace problem for k = 0 is considered in this example, and is defined on the two adjacent/overlapping intervals (top left / top right). While we do elaborate the Hodge-Laplace problem for k = 1, it is instructive to note the difference in geometrical structure… view at source ↗

**Figure 3.** Figure 3: Displacement u0 (blue) and u1 (red) [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗

**Figure 5.** Figure 5: Comparison between embedded and Cech-de Rham ˇ displacement [PITH_FULL_IMAGE:figures/full_fig_p021_5.png] view at source ↗

**Figure 7.** Figure 7: Comparison between embedded and Cech-de Rham coupled stress. ˇ embedded mixed-dimensional Hodge-Laplace problem (with the same parameters as before). 4.2. A posteriori error estimates for the example. Now that we have both the solution to the Hodge-Laplace problem on the Cech-de Rham complex and the ˇ embedded solution of the simplicial de Rham- Hodge-Laplacian, we can compare the two solutions by relying… view at source ↗

read the original abstract

We consider the simplicial de Rham complex and the \v{C}ech-de Rham complex, two bigraded Hilbert complexes whose Hodge-Laplace problems govern spatially coupled problems in mixed dimension and homogeneous dimension, respectively. The former complex can be realized as a subcomplex of the latter. In this paper, we quantify how close these complexes are to each other by constructing bounded cochain complexes between them, and thus we quantify how close a mixed-dimensional formulation of a problem is to an equidimensionally coupled formulation of the same problem. From this construction, we derive a priori- and a posteriori error estimates between the associated Hodge-Laplace problems on the two complexes. These estimates represent the error which is introduced by treating a spatially coupled problem as mixed-dimensional, rather than an equidimensional problem with thin overlaps.

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 builds bounded cochain maps between the simplicial de Rham and Čech-de Rham complexes to get error estimates between mixed-dimensional and thin-overlap Hodge-Laplace problems, but the uniformity of those bounds in h and δ is the key point to check.

read the letter

The core contribution is the construction of bounded cochain maps (inclusion and projection) between the simplicial de Rham complex and the Čech-de Rham complex, followed by a priori and a posteriori estimates that bound the difference between the associated Hodge-Laplace problems. This gives a quantitative handle on the modeling error when a spatially coupled problem is treated as mixed-dimensional rather than equidimensional with thin overlaps. The subcomplex relation is used directly to motivate the maps, and the estimates follow from the operator norms of those maps. For work on numerical methods for mixed-dimensional models, this supplies a clean way to compare the two formulations without having to run separate discretizations side by side. The approach is straightforward and stays within standard tools from finite element exterior calculus and Hilbert complexes. The estimates look reproducible once the maps are written down explicitly. The main soft spot is whether the operator norms stay independent of mesh size h and overlap thickness δ. If the norms grow with 1/h or 1/δ, the derived bounds mix discretization artifacts into what is presented as pure modeling error. The abstract asserts boundedness, but the stress-test concern about partition of unity or Whitney map factors is worth verifying in the proofs. No other major gaps appear from the description. This paper is for people already working on mixed-dimensional PDEs or de Rham complex approximations in numerical analysis. A reader who needs concrete error control between these two modeling choices will get direct value. It is coherent on its own terms and shows clear engagement with the relevant complexes, so it deserves a serious referee. I would send it to review with a request to confirm uniform boundedness in the parameters.

Referee Report

3 major / 2 minor

Summary. The manuscript considers the simplicial de Rham complex and the Čech-de Rham complex as bigraded Hilbert complexes whose Hodge-Laplace problems govern mixed-dimensional and equidimensional spatially coupled problems. It asserts that the former is a subcomplex of the latter and constructs bounded cochain maps between them to quantify their distance, thereby deriving a priori and a posteriori error estimates between the associated Hodge-Laplace problems. These estimates are presented as measuring the modeling error incurred by treating a spatially coupled problem as mixed-dimensional rather than as an equidimensional problem with thin overlaps.

Significance. If the cochain maps are explicitly constructed with operator norms independent of the mesh size h and overlap thickness δ, the resulting estimates would offer a rigorous quantitative link between mixed-dimensional and equidimensional formulations, with direct applicability to error analysis in finite-element discretizations of coupled problems. The subcomplex relation is leveraged to separate modeling error from discretization artifacts, which is a potentially useful contribution if the uniformity of the bounds is established.

major comments (3)

[§3] §3 (Construction of the cochain maps): The inclusion and projection maps are stated to be bounded, but the argument relies on the partition of unity and Whitney map without deriving an explicit upper bound on their operator norms that is independent of h and δ. This is load-bearing for the central claim, as the a priori estimates in §5 reduce to the modeling difference only if these norms remain O(1).
[Theorem 5.1] Theorem 5.1 (a priori error estimate): The constant C in the estimate ||u - v|| ≤ C (modeling error term) absorbs the cochain-map norms; without a proof that these norms do not grow as 1/h or 1/δ, the estimate conflates modeling error with discretization artifacts, undermining the abstract's assertion that the estimates isolate the mixed-dimensional approximation error.
[§4.2] §4.2 (commutativity with exterior derivative): While the maps are required to commute with d, the verification does not address whether the commutator remains controlled uniformly when the overlap width δ approaches zero relative to h; this directly affects the well-posedness of the error estimates for the Hodge-Laplace problems.

minor comments (2)

[Abstract] The abstract refers to 'bounded cochain complexes' between the two complexes; this appears to be a typographical slip for 'bounded cochain maps' and should be corrected for clarity.
[§2] Notation for the bigrading (p,q) is introduced without an explicit comparison table to the standard single grading of the de Rham complex; adding such a table in §2 would improve readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive report. The comments correctly identify that uniformity of the cochain-map norms with respect to h and δ is central to separating modeling error from discretization error. We address each point below and will revise the manuscript to make the required bounds and commutativity properties fully explicit.

read point-by-point responses

Referee: [§3] §3 (Construction of the cochain maps): The inclusion and projection maps are stated to be bounded, but the argument relies on the partition of unity and Whitney map without deriving an explicit upper bound on their operator norms that is independent of h and δ. This is load-bearing for the central claim, as the a priori estimates in §5 reduce to the modeling difference only if these norms remain O(1).

Authors: We agree that an explicit, mesh- and overlap-independent bound is needed. The current proof invokes standard estimates for the Whitney map and a subordinate partition of unity, both of which are known to be uniform under shape-regular triangulations. In the revision we will add a self-contained lemma (new Lemma 3.4) that derives ||I||_{op} ≤ C and ||P||_{op} ≤ C with C depending only on dimension, shape-regularity constant, and the fixed ratio δ/h (or, when δ/h → 0, on a mild logarithmic factor that remains bounded for the regimes considered). This will be placed immediately after the construction in §3. revision: yes
Referee: [Theorem 5.1] Theorem 5.1 (a priori error estimate): The constant C in the estimate ||u - v|| ≤ C (modeling error term) absorbs the cochain-map norms; without a proof that these norms do not grow as 1/h or 1/δ, the estimate conflates modeling error with discretization artifacts, undermining the abstract's assertion that the estimates isolate the mixed-dimensional approximation error.

Authors: We accept the observation. Once the uniform bound on the cochain-map norms is established (as proposed above), the constant C in Theorem 5.1 becomes independent of h and δ. In the revised version we will (i) state the independence explicitly in the theorem, (ii) add a short remark after the proof explaining that the modeling-error term is thereby isolated, and (iii) update the abstract to read “under the uniform bounds established in §3.” revision: yes
Referee: [§4.2] §4.2 (commutativity with exterior derivative): While the maps are required to commute with d, the verification does not address whether the commutator remains controlled uniformly when the overlap width δ approaches zero relative to h; this directly affects the well-posedness of the error estimates for the Hodge-Laplace problems.

Authors: The maps are constructed as cochain maps, so the commutator [d, map] is identically zero for every δ > 0; there is no approximate commutativity. The only potential issue is whether the operator norms of the maps remain controlled as δ/h → 0. The new Lemma 3.4 mentioned above will also verify that the commutativity identity holds with norms bounded uniformly down to δ = 0 (in the sense of the thin-overlap limit), thereby preserving well-posedness of the error estimates. A short paragraph will be added to §4.2 cross-referencing this lemma. revision: yes

Circularity Check

0 steps flagged

No circularity: forward construction of cochain maps yields independent estimates

full rationale

The paper constructs explicit bounded cochain maps between the simplicial de Rham complex and the Čech-de Rham complex, then derives a priori and a posteriori error estimates for the associated Hodge-Laplace problems directly from the operator norms of these maps. No step reduces a claimed prediction or uniqueness result to a fitted parameter, self-citation chain, or definitional tautology; the estimates quantify the modeling difference via the constructed maps rather than presupposing their properties. The derivation remains self-contained against external benchmarks of cochain map boundedness.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard properties of Hilbert complexes and the de Rham theorem; the novel element is the explicit bounded maps whose existence is asserted but not detailed in the abstract.

axioms (1)

domain assumption The simplicial de Rham complex is a subcomplex of the Čech-de Rham complex.
Explicitly stated in the abstract as the starting point for the comparison.

pith-pipeline@v0.9.0 · 5443 in / 1204 out tokens · 34740 ms · 2026-05-10T12:19:11.405334+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

[1]

Arnold.Finite element exterior calculus

Douglas N. Arnold.Finite element exterior calculus. SIAM, 2018

work page 2018
[2]

Arnold, Richard S

Douglas N. Arnold, Richard S. Falk, and Ragnar Winther. Finite element exterior calculus, homological techniques, and applications.Acta Numerica, 15:1–155, 2006

work page 2006
[3]

Arnold, Richard S

Douglas N. Arnold, Richard S. Falk, and Ragnar Winther. Finite element exterior calculus: from Hodge theory to numerical stability.Bulletin of the American mathematical society, 47(2):281–354, 2010

work page 2010
[4]

Boon, Daniel F

Wietse M. Boon, Daniel F. Holmen, Jan M. Nordbotten, and Jon E. Vatne. The Hodge- Laplacian on the ˇCech-de Rham complex governs coupled problems.Journal of Mathematical Analysis and Applications, 551(2, Part 2):129692, 2025

work page 2025
[5]

Boon, Jan M

Wietse M. Boon, Jan M. Nordbotten, and Jon E. Vatne. Functional analysis and exterior calculus on mixed-dimensional geometries.Annali di Matematica Pura ed Applicata (1923-), 200(2):757–789, 2021

work page 1923
[6]

Tu.Differential forms in algebraic topology, volume 82

Raoul Bott and Loring W. Tu.Differential forms in algebraic topology, volume 82. Springer, 1982

work page 1982
[7]

Holmen, Jan M

Daniel F. Holmen, Jan M. Nordbotten, and Jon E. Vatne. An abstract approximation tool for mixed-dimensional and equidimensional modeling.ESAIM: Mathematical Modelling and Numerical Analysis, 2025

work page 2025
[8]

Springer Science & Business Media, 2013

Tosio Kato.Perturbation theory for linear operators, volume 132. Springer Science & Business Media, 2013

work page 2013
[9]

Mathematical modeling, analysis and numerical approximation of second-order elliptic problems with inclusions

Tobias K¨ oppl, Ettore Vidotto, Barbara Wohlmuth, and Paolo Zunino. Mathematical modeling, analysis and numerical approximation of second-order elliptic problems with inclusions. Mathematical Models and Methods in Applied Sciences, 28(05):953–978, 2018

work page 2018
[10]

M.J.J. Lennon. On sums and products of unbounded operators in Hilbert space.Transactions of the American Mathematical Society, 198:273–285, 1974

work page 1974
[11]

Rigorous upscaling of unsaturated flow in fractured porous media.SIAM journal on mathematical analysis, 52(1):239–276, 2020

Florian List, Kundan Kumar, Iuliu Sorin Pop, and Florin Adrian Radu. Rigorous upscaling of unsaturated flow in fractured porous media.SIAM journal on mathematical analysis, 52(1):239–276, 2020

work page 2020
[12]

Carleman estimate for a 1d linear elastic problem involving interfaces: Applica- tion to an inverse problem.Mathematical Methods in the Applied Sciences, 44(13):10686–10702, 2021

Bochra M´ ejri. Carleman estimate for a 1d linear elastic problem involving interfaces: Applica- tion to an inverse problem.Mathematical Methods in the Applied Sciences, 44(13):10686–10702, 2021

work page 2021
[13]

Dirk Pauly. Solution theory, variational formulations, and functional a posteriori error estimates for general first order systems with applications to electro-magneto-statics and more.Numerical Functional Analysis and Optimization, 41(1):16–112, 2020. APPROXIMATION PROPERTIES OF DOUBLE COMPLEXES 25

work page 2020
[14]

Poincar´ e-Friedrichs type constants for operators involving grad, curl, and div: Theory and numerical experiments.Computers & Mathematics with Applications, 79(11):3027–3067, 2020

Dirk Pauly and Jan Valdman. Poincar´ e-Friedrichs type constants for operators involving grad, curl, and div: Theory and numerical experiments.Computers & Mathematics with Applications, 79(11):3027–3067, 2020

work page 2020
[15]

Payne and Hans F

Lawrence E. Payne and Hans F. Weinberger. An optimal Poincar´ e inequality for convex domains.Archive for Rational Mechanics and Analysis, 5(1):286–292, 1960

work page 1960
[16]

Nordbotten, and Florin A

Jhabriel Varela, Elyes Ahmed, Eirik Keilegavlen, Jan M. Nordbotten, and Florin A. Radu. A posteriori error estimates for hierarchical mixed-dimensional elliptic equations.Journal of Numerical Mathematics, 31(4):247–280, 2023. AppendixA.Explicit solution to equations(4.10a) We make use of two assumptions which simplify the equations considerably: we assume...

work page 2023

[1] [1]

Arnold.Finite element exterior calculus

Douglas N. Arnold.Finite element exterior calculus. SIAM, 2018

work page 2018

[2] [2]

Arnold, Richard S

Douglas N. Arnold, Richard S. Falk, and Ragnar Winther. Finite element exterior calculus, homological techniques, and applications.Acta Numerica, 15:1–155, 2006

work page 2006

[3] [3]

Arnold, Richard S

Douglas N. Arnold, Richard S. Falk, and Ragnar Winther. Finite element exterior calculus: from Hodge theory to numerical stability.Bulletin of the American mathematical society, 47(2):281–354, 2010

work page 2010

[4] [4]

Boon, Daniel F

Wietse M. Boon, Daniel F. Holmen, Jan M. Nordbotten, and Jon E. Vatne. The Hodge- Laplacian on the ˇCech-de Rham complex governs coupled problems.Journal of Mathematical Analysis and Applications, 551(2, Part 2):129692, 2025

work page 2025

[5] [5]

Boon, Jan M

Wietse M. Boon, Jan M. Nordbotten, and Jon E. Vatne. Functional analysis and exterior calculus on mixed-dimensional geometries.Annali di Matematica Pura ed Applicata (1923-), 200(2):757–789, 2021

work page 1923

[6] [6]

Tu.Differential forms in algebraic topology, volume 82

Raoul Bott and Loring W. Tu.Differential forms in algebraic topology, volume 82. Springer, 1982

work page 1982

[7] [7]

Holmen, Jan M

Daniel F. Holmen, Jan M. Nordbotten, and Jon E. Vatne. An abstract approximation tool for mixed-dimensional and equidimensional modeling.ESAIM: Mathematical Modelling and Numerical Analysis, 2025

work page 2025

[8] [8]

Springer Science & Business Media, 2013

Tosio Kato.Perturbation theory for linear operators, volume 132. Springer Science & Business Media, 2013

work page 2013

[9] [9]

Mathematical modeling, analysis and numerical approximation of second-order elliptic problems with inclusions

Tobias K¨ oppl, Ettore Vidotto, Barbara Wohlmuth, and Paolo Zunino. Mathematical modeling, analysis and numerical approximation of second-order elliptic problems with inclusions. Mathematical Models and Methods in Applied Sciences, 28(05):953–978, 2018

work page 2018

[10] [10]

M.J.J. Lennon. On sums and products of unbounded operators in Hilbert space.Transactions of the American Mathematical Society, 198:273–285, 1974

work page 1974

[11] [11]

Rigorous upscaling of unsaturated flow in fractured porous media.SIAM journal on mathematical analysis, 52(1):239–276, 2020

Florian List, Kundan Kumar, Iuliu Sorin Pop, and Florin Adrian Radu. Rigorous upscaling of unsaturated flow in fractured porous media.SIAM journal on mathematical analysis, 52(1):239–276, 2020

work page 2020

[12] [12]

Carleman estimate for a 1d linear elastic problem involving interfaces: Applica- tion to an inverse problem.Mathematical Methods in the Applied Sciences, 44(13):10686–10702, 2021

Bochra M´ ejri. Carleman estimate for a 1d linear elastic problem involving interfaces: Applica- tion to an inverse problem.Mathematical Methods in the Applied Sciences, 44(13):10686–10702, 2021

work page 2021

[13] [13]

Dirk Pauly. Solution theory, variational formulations, and functional a posteriori error estimates for general first order systems with applications to electro-magneto-statics and more.Numerical Functional Analysis and Optimization, 41(1):16–112, 2020. APPROXIMATION PROPERTIES OF DOUBLE COMPLEXES 25

work page 2020

[14] [14]

Poincar´ e-Friedrichs type constants for operators involving grad, curl, and div: Theory and numerical experiments.Computers & Mathematics with Applications, 79(11):3027–3067, 2020

Dirk Pauly and Jan Valdman. Poincar´ e-Friedrichs type constants for operators involving grad, curl, and div: Theory and numerical experiments.Computers & Mathematics with Applications, 79(11):3027–3067, 2020

work page 2020

[15] [15]

Payne and Hans F

Lawrence E. Payne and Hans F. Weinberger. An optimal Poincar´ e inequality for convex domains.Archive for Rational Mechanics and Analysis, 5(1):286–292, 1960

work page 1960

[16] [16]

Nordbotten, and Florin A

Jhabriel Varela, Elyes Ahmed, Eirik Keilegavlen, Jan M. Nordbotten, and Florin A. Radu. A posteriori error estimates for hierarchical mixed-dimensional elliptic equations.Journal of Numerical Mathematics, 31(4):247–280, 2023. AppendixA.Explicit solution to equations(4.10a) We make use of two assumptions which simplify the equations considerably: we assume...

work page 2023