pith. sign in

arxiv: 2606.03141 · v2 · pith:JFONW7HQnew · submitted 2026-06-02 · 🧮 math.NA · cs.NA

Owner-selected bubble transforms and coefficient-robust Schwarz preconditioners for variable-degree hp finite elements

Situan Li , Weiying Zheng This is my paper

Pith reviewed 2026-06-28 09:25 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords hp finite elementsSchwarz preconditionersbubble transformsvariable-degreeadditive Schwarzreaction-diffusionNitsche discretizationfitted interfaces
0
0 comments X

The pith

Owner-selected Falk-Winther bubble transforms produce L2- and H1-stable decompositions for variable-degree hp finite elements with constants independent of mesh size and degrees.

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

This paper constructs space decompositions and Schwarz preconditioners for hp finite element discretizations where polynomial degrees can vary from element to element. The key step is an owner-selected bubble transform that splits the space into components stable in L2 and H1 norms with bounds that ignore both h and p variations. These decompositions lead to preconditioners for reaction-diffusion equations and fitted-interface problems that remain effective regardless of how degrees are distributed or how coefficients jump locally. A sympathetic reader would care because such robustness removes the need to adjust solver parameters when adapting the mesh or polynomial order in high-order simulations. The constructions use minimal-degree owners for general cases and coefficient-adapted ones when chains of elements satisfy certain conditions.

Core claim

On conforming simplicial meshes, an owner-selected Falk--Winther bubble transform gives L²- and H¹-stable components with constants independent of the mesh size, the local polynomial degrees, and the degree distribution. Minimal-degree owners preserve arbitrary variable-degree spaces with p_K≥1, while coefficient-adapted owners yield weighted estimates under local chain conditions. Combined with a weighted continuous piecewise affine extraction, this gives hp-uniform Schwarz preconditioners for conforming reaction-diffusion problems with locally comparable coefficients, and a coefficient-weighted conforming variant in the uniform-degree case. For three-dimensional fitted-interface problems,

What carries the argument

The owner-selected Falk--Winther bubble transform, which assigns each bubble function to a chosen owner element to produce L2- and H1-stable splitting of the finite element space independent of local degrees.

Load-bearing premise

The stability constants remain independent of degree distribution only if minimal-degree or coefficient-adapted owners can always be selected while preserving the spaces and satisfying local chain conditions, plus a common-degree condition at interfaces.

What would settle it

A sequence of refined meshes with increasing local degree variation where the computed condition number of the vertex-patch preconditioner grows unbounded would show the claimed independence does not hold.

Figures

Figures reproduced from arXiv: 2606.03141 by Situan Li, Weiying Zheng.

Figure 4
Figure 4. Figure 4 [PITH_FULL_IMAGE:figures/full_fig_p023_4.png] view at source ↗
Figure 6
Figure 6. Figure 6 [PITH_FULL_IMAGE:figures/full_fig_p034_6.png] view at source ↗
Figure 6
Figure 6. Figure 6 [PITH_FULL_IMAGE:figures/full_fig_p037_6.png] view at source ↗
read the original abstract

We construct $h$- and $p$-robust, degree-preserving space decompositions and additive Schwarz preconditioners for variable-degree $hp$ finite element discretizations of conforming reaction-diffusion and fitted-interface problems. On conforming simplicial meshes, an owner-selected Falk--Winther bubble transform gives $L^2$- and $H^1$-stable components with constants independent of the mesh size, the local polynomial degrees, and the degree distribution. Minimal-degree owners preserve arbitrary variable-degree spaces with $p_K\ge1$, while coefficient-adapted owners yield weighted estimates under local chain conditions. Combined with a weighted continuous piecewise affine extraction, this gives $hp$-uniform Schwarz preconditioners for conforming reaction-diffusion problems with locally comparable coefficients, and a coefficient-weighted conforming variant in the uniform-degree case. For three-dimensional fitted-interface problems, we use a symmetric Nitsche discretization on a tetrahedral mesh fitted to a piecewise planar interface. Surface jump components are lifted into the side selected by the penalty scaling using patch-level $p$-robust trace liftings. The conforming remainder is decomposed by the low-order extraction and a weighted one-sided bubble transform. Grouping the resulting components by vertices yields a practical vertex-patch Schwarz preconditioner whose condition number is independent of the mesh size, local polynomial degrees, diffusion contrast, and coefficient magnitudes under a common-degree condition on interface-touching tetrahedra. Numerical experiments for pure diffusion problems support the theory and suggest robustness beyond the common-degree assumption.

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.

Referee Report

2 major / 2 minor

Summary. The paper constructs owner-selected Falk--Winther bubble transforms to obtain L²- and H¹-stable space decompositions for variable-degree hp finite element discretizations of conforming reaction-diffusion problems and symmetric Nitsche discretizations of fitted-interface problems on tetrahedral meshes. Minimal-degree owners handle arbitrary variable-degree spaces (p_K ≥ 1), while coefficient-adapted owners are used under local chain conditions; for interfaces, surface jumps are lifted via p-robust trace liftings and the conforming remainder is decomposed with a weighted one-sided bubble transform. These yield vertex-patch additive Schwarz preconditioners whose condition numbers are independent of h, local degrees, degree distribution, and coefficient contrast under a common-degree condition on interface-touching tetrahedra. Numerical experiments support the theory and indicate possible robustness beyond the common-degree assumption.

Significance. If the stability results hold, the work provides a technically substantive advance in the construction of parameter-robust preconditioners for hp methods with variable degrees and interfaces. The owner-selection mechanism for bubble transforms, combined with weighted continuous piecewise affine extraction, offers a systematic way to achieve h- and p-uniform bounds while preserving the finite element space; the numerical validation for pure diffusion problems adds concrete support. Such tools are valuable for efficient solvers in adaptive high-order simulations of heterogeneous media.

major comments (2)
  1. [Abstract] Abstract: The opening claim that the owner-selected Falk--Winther bubble transform yields L²- and H¹-stable components with constants independent of mesh size, local polynomial degrees, and the degree distribution is stated without immediate qualification, yet the subsequent paragraph on fitted-interface problems explicitly requires the common-degree condition on all tetrahedra touching the interface for the vertex-patch preconditioner. This condition is load-bearing because the weighted one-sided bubble transform and patch grouping rely on matching degrees to control trace liftings and Nitsche jumps; without it the uniform bound may fail. The abstract should foreground the precise scope of the independence claim.
  2. [fitted-interface preconditioner section] Section describing the fitted-interface preconditioner (near the end of the abstract and the corresponding theory): The manuscript correctly notes that experiments suggest robustness beyond the common-degree assumption, but provides no additional analysis, counterexample, or relaxed theorem. Because the trace-lifting and coefficient-weighted extraction arguments explicitly invoke matching polynomial degrees on interface-touching elements, the gap between the proven statement and the suggested broader applicability is a load-bearing limitation for the robustness conclusion.
minor comments (2)
  1. [Preliminaries / notation] The distinction between minimal-degree owners (for arbitrary variable-degree spaces) and coefficient-adapted owners (under local chain conditions) is introduced in the abstract but would benefit from an early explicit definition and a small illustrative diagram in the preliminaries.
  2. [Numerical experiments] Numerical experiments section: Adding a short table or paragraph that reports condition numbers for at least one configuration violating the common-degree condition would make the claim of suggested robustness beyond the assumption more concrete and falsifiable.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback. We address the two major comments below and will make corresponding revisions to clarify the scope of the results.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The opening claim that the owner-selected Falk--Winther bubble transform yields L²- and H¹-stable components with constants independent of mesh size, local polynomial degrees, and the degree distribution is stated without immediate qualification, yet the subsequent paragraph on fitted-interface problems explicitly requires the common-degree condition on all tetrahedra touching the interface for the vertex-patch preconditioner. This condition is load-bearing because the weighted one-sided bubble transform and patch grouping rely on matching degrees to control trace liftings and Nitsche jumps; without it the uniform bound may fail. The abstract should foreground the precise scope of the independence claim.

    Authors: We agree that the abstract would benefit from foregrounding the scope of the claims. The first paragraph addresses the general reaction-diffusion case (no common-degree restriction), while the interface case carries the additional assumption. We will revise the abstract to state the independence claims with the appropriate qualification from the outset, distinguishing the two settings explicitly. revision: yes

  2. Referee: [fitted-interface preconditioner section] Section describing the fitted-interface preconditioner (near the end of the abstract and the corresponding theory): The manuscript correctly notes that experiments suggest robustness beyond the common-degree assumption, but provides no additional analysis, counterexample, or relaxed theorem. Because the trace-lifting and coefficient-weighted extraction arguments explicitly invoke matching polynomial degrees on interface-touching elements, the gap between the proven statement and the suggested broader applicability is a load-bearing limitation for the robustness conclusion.

    Authors: We acknowledge the gap: the analysis relies on the common-degree condition for the trace liftings and weighted extraction, and we provide neither a relaxed theorem nor a counterexample. The numerical suggestion of robustness beyond the assumption is observational only. We will revise the abstract and conclusion to remove any implication of proven robustness without the condition, stating only the proven bound under common degree and noting that experiments indicate the assumption may not be necessary in practice. revision: yes

Circularity Check

0 steps flagged

No significant circularity; constructions are independent of inputs

full rationale

The paper defines owner-selected Falk-Winther bubble transforms and proves L2/H1 stability plus Schwarz preconditioner bounds directly from standard finite-element trace and extension operators on simplicial meshes. The common-degree condition on interface tetrahedra is stated explicitly as a hypothesis required for the vertex-patch bound, not derived from the result itself. No parameters are fitted to data and then relabeled as predictions, no self-citations are invoked to justify uniqueness or ansatzes, and the derivations do not reduce by construction to the inputs they claim to establish. The theory is therefore self-contained against external benchmarks in finite-element analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The paper builds on standard finite element space properties and introduces new transform constructions; no free parameters, invented entities, or ad-hoc axioms beyond domain assumptions are apparent from the abstract.

axioms (3)
  • standard math Conforming simplicial meshes support Falk-Winther bubble transforms with the stated stability
    Invoked for the core decomposition on conforming meshes.
  • domain assumption Local chain conditions hold for coefficient-adapted owners
    Required for weighted estimates in reaction-diffusion case.
  • domain assumption Common-degree condition on interface-touching tetrahedra
    Needed for the interface preconditioner robustness.

pith-pipeline@v0.9.1-grok · 5796 in / 1345 out tokens · 30178 ms · 2026-06-28T09:25:21.346053+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Stable Triangle Projections for Variable-Degree Tetrahedral Spaces and Uniform IPDG Preconditioning

    math.NA 2026-06 unverdicted novelty 6.0

    Develops stable variable-degree projections on triangles and tetrahedra to obtain h-p-uniform preconditioners for conforming and symmetric IPDG Laplace operators whose constants depend only on local degree bounds, mes...

Reference graph

Works this paper leans on

36 extracted references · 17 canonical work pages · cited by 1 Pith paper

  1. [1]

    Ainsworth and S

    M. Ainsworth and S. Jiang , Preconditioning the mass matrix for high order finite element approximation on triangles , SIAM Journal on Numerical Analysis, 57 (2019), pp. 355--377

    2019

  2. [2]

    A384--A414

    height 2pt depth -1.6pt width 23pt, Preconditioning the mass matrix for high order finite element approximation on tetrahedra , SIAM Journal on Scientific Computing, 43 (2021), pp. A384--A414

    2021

  3. [3]

    1465--1491

    height 2pt depth -1.6pt width 23pt, Uniform substructuring preconditioners for high order FEM on triangles and the influence of nodal basis functions , SIAM Journal on Numerical Analysis, 62 (2024), pp. 1465--1491

    2024

  4. [4]

    Babu s ka, A

    I. Babu s ka, A. Craig, J. Mandel, and J. Pitk \"a ranta , Efficient preconditioning for the p-version finite element method in two dimensions , SIAM Journal on Numerical Analysis, 28 (1991), pp. 624--661

    1991

  5. [5]

    Burman, J

    E. Burman, J. Guzm \'a n, M. A. S \'a nchez, and M. Sarkis , Robust flux error estimation of an unfitted Nitsche method for high-contrast interface problems , IMA Journal of Numerical Analysis, 38 (2018), pp. 646--668

    2018

  6. [6]

    R. S. Falk and R. Winther , The bubble transform: A new tool for analysis of finite element methods , Foundations of Computational Mathematics, 16 (2016), pp. 297--328

    2016

  7. [7]

    Grisvard , Elliptic Problems in Nonsmooth Domains , vol

    P. Grisvard , Elliptic Problems in Nonsmooth Domains , vol. 24 of Monographs and Studies in Mathematics, Pitman, Boston, 1985

    1985

  8. [8]

    Guo and W

    B. Guo and W. Cao , An additive Schwarz method for the h - p version of the finite element method in three dimensions , SIAM Journal on Numerical Analysis, 35 (1998), pp. 632--654

    1998

  9. [9]

    Hansbo and P

    A. Hansbo and P. Hansbo , An unfitted finite element method, based on Nitsche 's method, for elliptic interface problems , Computer Methods in Applied Mechanics and Engineering, 191 (2002), pp. 5537--5552

    2002

  10. [10]

    Karkulik, J

    M. Karkulik, J. M. Melenk, and A. Rieder , On interpolation spaces of piecewise polynomials on mixed meshes , ESAIM: Mathematical Modelling and Numerical Analysis, 59 (2025), pp. 231--264

    2025

  11. [11]

    Li and W

    S. Li and W. Zheng , p -robust trace liftings for discrete harmonic extensions and boundary-preserving hp interpolation on tetrahedral meshes , 2026. arXiv preprint arXiv:2606.02086

    Pith/arXiv arXiv 2026

  12. [12]

    L. F. Pavarino , Additive Schwarz methods for the p -version finite element method , Numerische Mathematik, 66 (1994), pp. 493--515

    1994

  13. [13]

    Sch \"o berl, J

    J. Sch \"o berl, J. M. Melenk, C. Pechstein, and S. Zaglmayr , Additive Schwarz preconditioning for p -version triangular and tetrahedral finite elements , IMA Journal of Numerical Analysis, 28 (2008), pp. 1--24

    2008

  14. [14]

    and Winther, Ragnar , year = 2016, month = feb, journal =

    Falk, Richard S. and Winther, Ragnar , year = 2016, month = feb, journal =. The Bubble Transform: A New Tool for Analysis of Finite Element Methods , shorttitle =. doi:10.1007/s10208-015-9252-1 , urldate =

    work page doi:10.1007/s10208-015-9252-1 2016

  15. [15]

    SIAM Journal on Numerical Analysis , volume =

    Polynomial Liftings on a Tetrahedron and Applications to the H-p Version of the Finite Element Method in Three Dimensions , author =. SIAM Journal on Numerical Analysis , volume =. doi:10.1137/S0036142994267552 , urldate =

    work page doi:10.1137/s0036142994267552

  16. [16]

    Part I , author =

    Polynomial Extension Operators. Part I , author =. SIAM Journal on Numerical Analysis , volume =. doi:10.1137/070698786 , langid =

    work page doi:10.1137/070698786

  17. [17]

    SIAM Journal on Numerical Analysis , volume =

    hp-Interpolation of Nonsmooth Functions and an Application to hp-a Posteriori Error Estimation , author =. SIAM Journal on Numerical Analysis , volume =. doi:10.1137/S0036142903432930 , urldate =

    work page doi:10.1137/s0036142903432930

  18. [18]

    Mathematics of Computation , volume =

    On Commuting p -Version Projection-Based Interpolation on Tetrahedra , author =. Mathematics of Computation , volume =. doi:10.1090/mcom/3454 , langid =

    work page doi:10.1090/mcom/3454

  19. [19]

    Approximation by Finite Element Functions Using Local Regularization , author =. RAIRO. Analyse Num

  20. [20]

    ESAIM: Mathematical Modelling and Numerical Analysis , volume =

    On Interpolation Spaces of Piecewise Polynomials on Mixed Meshes , author =. ESAIM: Mathematical Modelling and Numerical Analysis , volume =. doi:10.1051/m2an/2024069 , langid =

    work page doi:10.1051/m2an/2024069

  21. [21]

    Elliptic Problems in Nonsmooth Domains , author =

  22. [22]

    p - and hp -Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics , author =

  23. [23]

    doi:10.1007/b137868 , langid =

    Domain Decomposition Methods: Algorithms and Theory , author =. doi:10.1007/b137868 , langid =

    work page doi:10.1007/b137868

  24. [24]

    SIAM Journal on Numerical Analysis , volume =

    Efficient Preconditioning for the P-Version Finite Element Method in Two Dimensions , author =. SIAM Journal on Numerical Analysis , volume =. doi:10.1137/0728034 , urldate =. 2157721 , eprinttype =

    work page doi:10.1137/0728034

  25. [25]

    , year = 1994, journal =

    Pavarino, Luca F. , year = 1994, journal =. Additive. doi:10.1007/BF01385709 , langid =

    work page doi:10.1007/bf01385709 1994

  26. [26]

    An Additive

    Guo, Benqi and Cao, Weiming , year = 1998, journal =. An Additive. doi:10.1137/S0036142996299435 , langid =

    work page doi:10.1137/s0036142996299435 1998

  27. [27]

    Additive

    Sch. Additive. IMA Journal of Numerical Analysis , volume =. doi:10.1093/imanum/drl046 , langid =

    work page doi:10.1093/imanum/drl046

  28. [28]

    SIAM Journal on Numerical Analysis , volume =

    Preconditioning the Mass Matrix for High Order Finite Element Approximation on Triangles , author =. SIAM Journal on Numerical Analysis , volume =. doi:10.1137/18M1182450 , langid =

    work page doi:10.1137/18m1182450

  29. [29]

    SIAM Journal on Scientific Computing , volume =

    Preconditioning the Mass Matrix for High Order Finite Element Approximation on Tetrahedra , author =. SIAM Journal on Scientific Computing , volume =. doi:10.1137/20M1333018 , langid =

    work page doi:10.1137/20m1333018

  30. [30]

    Uniform Substructuring Preconditioners for High Order

    Ainsworth, Mark and Jiang, Shuai , year = 2024, journal =. Uniform Substructuring Preconditioners for High Order. doi:10.1137/23M1561920 , langid =

    work page doi:10.1137/23m1561920 2024

  31. [31]

    An Unfitted Finite Element Method, Based on

    Hansbo, Anita and Hansbo, Peter , year = 2002, journal =. An Unfitted Finite Element Method, Based on. doi:10.1016/S0045-7825(02)00524-8 , langid =

    work page doi:10.1016/s0045-7825(02)00524-8 2002

  32. [32]

    Robust Flux Error Estimation of an Unfitted

    Burman, Erik and Guzm. Robust Flux Error Estimation of an Unfitted. IMA Journal of Numerical Analysis , volume =. doi:10.1093/imanum/drx017 , langid =

    work page doi:10.1093/imanum/drx017

  33. [33]

    Foundations of Computational Mathematics , volume =

    Stable Liftings of Polynomial Traces on Tetrahedra , author =. Foundations of Computational Mathematics , volume =. doi:10.1007/s10208-024-09670-x , langid =

    work page doi:10.1007/s10208-024-09670-x

  34. [34]

    Li, Situan and Zheng, Weiying , year = 2026, note =

    2026

  35. [35]

    Li, Situan and Zheng, Weiying , year = 2026, note =. Uniform

    2026

  36. [36]

    Stable Triangle Projections for Variable-Degree Tetrahedral Spaces and Uniform

    Li, Situan and Zheng, Weiying , year = 2026, note =. Stable Triangle Projections for Variable-Degree Tetrahedral Spaces and Uniform

    2026