pith. sign in

arxiv: 2606.31829 · v1 · pith:Q7L5PURTnew · submitted 2026-06-30 · 🧮 math.NA · cs.NA

Complementary families of approximating polynomials with applications to finite element methods applied to differential equations of arbitrary even spatial order

Peter K Moore This is my paper

Pith reviewed 2026-07-01 04:00 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords finite element methodC^m continuityHermite splinesLegendre polynomialsbubble functionsinterpolation errorsuperconvergence
0
0 comments X

The pith

Complementary families of polynomials generate C^m finite element basis functions of order p >= 2m+2 for arbitrary m >= 0.

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

The paper shows how two related families of polynomials can be combined to build finite element spaces that enforce C^m continuity at element interfaces while using polynomials of degree p at least 2m+2. The first family supplies the nodal basis functions via explicit Hermite splines defined on the unit interval. The second family supplies the interior bubble functions by taking (m+1)th derivatives of Legendre polynomials of degree p-m-1 and weighting them with binomial factors at the endpoints. Together the families produce an interpolant whose error formula extends known results for m=0 and m=1, and they deliver a superconvergence statement plus asymptotic equivalence between the interpolant and the finite-element solution for linear problems measured in the H^{m+1} norm.

Core claim

Complementary families of polynomials generate C^m finite element basis functions of order p >= 2m+2 for arbitrary m >= 0, with explicit formulas for the nodal and bubble functions, a relationship between the families, an interpolation error formula extending the m=0 and m=1 cases, and a superconvergence result along with asymptotic equivalence of the interpolant and finite element solution in the linear case in H^{m+1}.

What carries the argument

Complementary families consisting of Hermite splines (nodal basis ensuring C^m continuity) and polynomials obtained from (m+1)th derivatives of Legendre polynomials of degree p-m-1 multiplied by binomial powers (bubble basis on the interior).

If this is right

  • Explicit formulas for both families permit direct coding of the basis for any m without deriving new expressions case by case.
  • The interpolation error formula holds for every m, not only the classical cases m=0 and m=1.
  • In the linear setting the finite-element solution and the interpolant become asymptotically equivalent in the H^{m+1} norm.
  • Superconvergence occurs at the nodes in the H^{m+1} norm for the linear model problem.

Where Pith is reading between the lines

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

  • The explicit relationship between the two families may allow one family to be generated from the other, reducing storage or assembly cost in code.
  • The construction is stated for even-order differential equations; the same polynomial families could be tested on odd-order or mixed-order problems to see whether the continuity requirement changes.
  • Because the bubble functions are tied to ultraspherical polynomials, quadrature rules already known for those polynomials might accelerate the assembly of stiffness matrices.

Load-bearing premise

The second family produces functions linearly independent from the Hermite nodal basis that span the required interior space while preserving C^m continuity when elements are assembled.

What would settle it

Direct verification, for a chosen small m and p >= 2m+2, that the assembled basis functions are C^m continuous at nodes and that the computed interpolation error matches the stated formula.

Figures

Figures reproduced from arXiv: 2606.31829 by Peter K Moore.

Figure 1
Figure 1. Figure 1: Plot of (231) (solid curve) as a function of [PITH_FULL_IMAGE:figures/full_fig_p029_1.png] view at source ↗
read the original abstract

Complementary families of polynomials are introduced to generate $C^m$ finite element basis functions of order $p \geq 2m+2$ for arbitrary $m \ge 0$. One family consists of the Hermite splines that serve as the nodal basis functions by ensuring $C^m$ continuity across element boundaries. Explicit formulas for these splines for any $m \ge 0$ are presented on the canonical interval $[0,1]$. The second family is derived on the interval $[-1,1]$ from derivatives of order $m+1$ of the Legendre polynomials of degree $p-m-1$ multiplied by binomial powers of degree $m+1$ at -1 and 1, respectively. These polynomials, related to the ultraspherical polynomials, serve as the interior or bubble basis functions. A relationship between the two families of polynomials is demonstrated. For a particular $m$ and $p$, an interpolant is constructed using these basis pairs together with the roots of the related ultraspherical polynomial and the interval endpoints. A formula for the interpolation error that extends the results for $m=0$ and $m=1$ is given. To prove the formula extensions of the Lagrange interpolants are introduced. A superconvergence result along with the related asymptotic equivalence of the interpolant and finite element solution is proved in the linear case in $H^{m+1}$. Computational results demonstrate the theory for a model problem.

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

0 major / 3 minor

Summary. The paper introduces complementary families of polynomials to generate C^m finite element basis functions of order p ≥ 2m+2 for arbitrary m ≥ 0. One family consists of explicit Hermite spline nodal basis functions on [0,1] ensuring C^m continuity. The second family comprises bubble functions on [-1,1] obtained from (m+1)th derivatives of Legendre polynomials of degree p-m-1 multiplied by (x+1)^{m+1}(x-1)^{m+1}, related to ultraspherical polynomials. A relationship between the families is shown. An interpolant is constructed using these bases together with ultraspherical roots and endpoints; an interpolation error formula extending the m=0,1 cases is derived via extensions of Lagrange interpolants. A superconvergence result and asymptotic equivalence of the interpolant and finite element solution are proved in the linear case in H^{m+1}. Computational results for a model problem are presented.

Significance. If the explicit constructions, dimension counts, linear independence arguments, error formula, and superconvergence proof hold, the work supplies a systematic, explicit basis construction for high-order C^m elements applicable to even-order differential equations, extending known cases for m=0,1. The parameter-free derivation from Legendre polynomials, the matching of bubble space dimension p-2m-1, and the use of orthogonality for independence are strengths that support reproducibility and potential use in FEM codes. The superconvergence claim in H^{m+1} would be a notable addition if verified.

minor comments (3)
  1. [Section on relationship between families] §3 (or wherever the relationship between the two families is stated): the claimed algebraic relationship between the Hermite nodal functions and the bubble family should be stated as an explicit identity or recurrence rather than left as 'demonstrated' without the formula.
  2. [Section on interpolation error formula] The interpolation error formula extension: clarify whether the proof via extended Lagrange interpolants applies directly for all m or requires induction; the current description leaves the inductive step implicit.
  3. [Computational results] Computational results section: specify the exact values of m and p used in the model problem tests and confirm that the observed rates match the predicted superconvergence order in H^{m+1}.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary and significance assessment of the manuscript, as well as the recommendation for minor revision. No specific major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper derives the nodal (Hermite) family via explicit formulas on [0,1] and the bubble family via (m+1)th derivatives of Legendre polynomials of degree p-m-1 multiplied by (x+1)^{m+1}(x-1)^{m+1} on [-1,1], with the relationship between families and the interpolation error formula proved directly from these constructions and extensions of Lagrange interpolants. Dimension matching (p-2m-1 bubbles), linear independence, and C^m continuity follow from standard properties of orthogonal polynomials and differentiation kernels, which are external facts rather than self-referential. No parameter is fitted and then relabeled as a prediction, no uniqueness theorem is imported from self-citation, and no ansatz is smuggled via prior work. The superconvergence claim is proved for the linear case in H^{m+1}. The derivation chain is self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities; the constructions rely on standard properties of Legendre polynomials and binomial coefficients, which are treated as background.

pith-pipeline@v0.9.1-grok · 5796 in / 1306 out tokens · 54123 ms · 2026-07-01T04:00:53.644262+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

31 extracted references · 22 canonical work pages

  1. [1]

    Abramowitz and I.A

    M. Abramowitz and I.A. Stegun, Eds.,Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 10th Edition, Dept of Commerce, 1972

    1972

  2. [2]

    Adjerid,A posteriori error estimates for fourth-order elliptic problems, Comput

    S.A. Adjerid,A posteriori error estimates for fourth-order elliptic problems, Comput. Methods Appl. Mech. Engrg.191 (2002), nos. 23–24, 2539-2559.https://doi.org/10.1016/S0045-7825(02)00403-6

    work page doi:10.1016/s0045-7825(02)00403-6 2002

  3. [3]

    ´Alvarez-Caudevilla, J.D

    P. ´Alvarez-Caudevilla, J.D. Evans and V.A. Galaktionov,The Cauchy problem for tenth-order thin film equation I. Bi- furcation of oscillatory fundamental solutions, Mediterr. J. Math.27(2013), no. 4, 1761-1792.https://doi.org/10.1007/ s00009-013-0263-3

    2013

  4. [4]

    ´Alvarez-Caudevilla, J.D

    P. ´Alvarez-Caudevilla, J.D. Evans and V.A. Galaktionov,The Cauchy problem for tenth-order thin film equation II. Oscillatory source-type and fundamental similarity solutions, Discrete Contin. Dyn. Syst.35(2015), no. 3, 807-827. https://doi.org/10.3934/dcds.2015.35.807

    work page doi:10.3934/dcds.2015.35.807 2015

  5. [5]

    Arora, I

    S. Arora, I. Kaur,An efficient scheme for numerical solution of Burgers’ equation using quintic Hermite interpolating polynomials, Arab. J. Math.5(2016), no. 1, 23-34.https://doi.org/10.1007/s40065-015-0137-6

    work page doi:10.1007/s40065-015-0137-6 2016

  6. [6]

    Babuˇ ska and G

    I. Babuˇ ska and G. S¨ oderlind,On Roundoff error growth in elliptic problems, ACM Trans. Math. Soft.44(2018), no. 3, 1-22.https://doi.org/10.1145/3134444

    work page doi:10.1145/3134444 2018

  7. [7]

    Birkhoff, M.H

    G. Birkhoff, M.H. Schultz, and R.S. Varga,Piecewise Hermite interpolation in one and two variables with applications to partial differential equations, Numer. Math.11(1968), no. 3, 232-256.https://doi.org/10.1007/BF02161845 30

    work page doi:10.1007/bf02161845 1968

  8. [8]

    Alicante, A

    W.K.F. Brito, C.D.C.D. Maia, A.V. Mendonca,Bending analysis of elastically connected Euler–Bernoulli double-beam system using the direct boundary element method, Appl. Math. Model.74(2019), 387-408.https://doi.org/10.1016/j. apm.2019.04.049

    work page doi:10.1016/j 2019

  9. [9]

    Bukal and B

    M. Bukal and B. Muha,Justification of a nonlinear sixth-order thin-film equation as a reduced model for a fluid-structure interaction problem, Nonlinearity35(2022), no. 8, 4695-4726.https://doi.org/10.1088/1361-6544/ac7d89

    work page doi:10.1088/1361-6544/ac7d89 2022

  10. [10]

    Budd, V.A

    C.J. Budd, V.A. Galaktionov, J.F. Williams,Self-similar blow-up in higher-order semilinear parabolic equations, SIAM J. Appl. Math.64(2004), no. 5, 1775-1809.https://doi.org/10.1137/S003613990241552X

    work page doi:10.1137/s003613990241552x 2004

  11. [11]

    Galaktionov, E

    V.A. Galaktionov, E. Mitidieri, S.I. Pohozaev,Global sign-changing solutions of a higher order semilinear heat equation in the subcritical Fujita range, Adv. Nonlin. Stud.12(2012), no. 3, 569–596.https://doi.org/10.1515/ans-2012-0308

    work page doi:10.1515/ans-2012-0308 2012

  12. [12]

    Ghoul, V.T

    T.-E. Ghoul, V.T. Nguyen, H. Zaag,Construction of type 1 blowup solutions for a higher order semilinear parabolic equation, Adv. Nonlinear Anal.9(2020), no. 1, 388-412.https://doi.org/10.1515/anona-2020-0006

    work page doi:10.1515/anona-2020-0006 2020

  13. [13]

    Goodrich, T

    J. Goodrich, T. Hagstrom, and J. Lorenz,Hermite methods for hyperbolic initial-value problems, Math. Comp.75(2005), no. 254, 595-630.https://doi.org/10.1090/S0025-5718-05-01808-9

    work page doi:10.1090/s0025-5718-05-01808-9 2005

  14. [14]

    Hanson,Integrating Feynman-Kac equations using Hermite quintic finite elements, White Paper (2008) 1-20

    R. Hanson,Integrating Feynman-Kac equations using Hermite quintic finite elements, White Paper (2008) 1-20

    2008

  15. [15]

    Johansson et al.,mpmath: a Python library for arbitrary-precision floating-point arithmetic, Version 1.3.0 (2023)http: //mpmath.org/

    F. Johansson et al.,mpmath: a Python library for arbitrary-precision floating-point arithmetic, Version 1.3.0 (2023)http: //mpmath.org/

    2023

  16. [16]

    Juncheng, L

    L. Juncheng, L. Chengzhi, Y. Yeqing,The C3 piecewise seventh-degree Hermite interpolating splines with parameters, J. Computer-Aided Design and Computer Graphics29(2017), no. 12, 2204-2215.10.3724/SP.J.1089.2017.16427

    work page doi:10.3724/sp.j.1089.2017.16427 2017

  17. [17]

    Khakalo and J

    S. Khakalo and J. Niiranen,Form II of Mindlin’s second strain gradient theory of elasticity with a simplification: for materials and structures from nano- to macro-scales, Eur. J. Mech. A Solids71(2018), 292-319.https://doi.org/10. 1016/j.euromechsol.2018.02.013

    2018

  18. [18]

    Kumari and V.K

    A. Kumari and V.K. Kukreja,Survey of Hermite interpolating polynomials for the solution of differential equations, Mathematics11(2023), no. 14, Art. 3157, 28 pp.https://doi.org/10.3390/math11143157

    work page doi:10.3390/math11143157 2023

  19. [19]

    Liu,Qualitative properties for a sixth-order thin film equation, Math

    C. Liu,Qualitative properties for a sixth-order thin film equation, Math. Model. Anal.15(2010), no. 4, 457-471.https: //doi.org/10.3846/1392-6292.2010.15.457-471

    work page doi:10.3846/1392-6292.2010.15.457-471 2010

  20. [20]

    Guo, and J

    L.Lu, X. Guo, and J. Zhao,A unified nonlocal strain gradient model for nanobeams and the importance of higher order terms, Int. J. Eng. Sci.119(2017), 265-277.https://doi.org/10.1016/j.ijengsci.2017.06.024

    work page doi:10.1016/j.ijengsci.2017.06.024 2017

  21. [21]

    Mindlin,Second gradient of strain and surface-tension in linear elasticity, Int

    R.D. Mindlin,Second gradient of strain and surface-tension in linear elasticity, Int. J. Solids Struct.1(1965), no.4, 417–438.https://doi.org/10.1016/0020-7683(65)90006-5

    work page doi:10.1016/0020-7683(65)90006-5 1965

  22. [22]

    Moore,Interpolation error-based a posteriori error estimation for two-point boundary value problems and parabolic equations in one space dimension, Numer

    P.K. Moore,Interpolation error-based a posteriori error estimation for two-point boundary value problems and parabolic equations in one space dimension, Numer. Math.90(2001), no. 1, 149-177.https://doi.org/10.1007/s002110100297

    work page doi:10.1007/s002110100297 2001

  23. [23]

    Moore,Applications of Lobatto polynomials to an adaptive finite element method: a posteriori error estimates for hp-adaptivity and grid-to-grid interpolation, Numer

    P.K. Moore,Applications of Lobatto polynomials to an adaptive finite element method: a posteriori error estimates for hp-adaptivity and grid-to-grid interpolation, Numer. Math.94(2003), no. 2, 367-401.https://doi.org/10.1007/ s00211-002-0424-y

    2003

  24. [24]

    Moore, M

    P.K. Moore, M. Rangelova,A posteriori error estimation for hp-adaptivity for fourth-order equations, Math. Comp.79 (2010), no. 270, 677-705.https://doi.org/10.1090/S0025-5718-09-02290-X

    work page doi:10.1090/s0025-5718-09-02290-x 2010

  25. [25]

    Olver and A

    S. Olver and A. Townsend,A fast and well-conditioned spectral method,SIAM Review55(2013), no. 3, 462-489.https: //doi.org/10.1137/120865458

    work page doi:10.1137/120865458 2013

  26. [26]

    Rainville,Special Functions, Chelsea Publishing Co., New York, 1960

    E.D. Rainville,Special Functions, Chelsea Publishing Co., New York, 1960

    1960

  27. [27]

    Siddiqi and G

    S.S. Siddiqi and G. Akram,Solutions of tenth-order boundary value problems using eleventh degree spline, Appl. Math. Comput.185(2007), no. 1, 115-127.https://doi.org/10.1016/j.amc.2006.05.070

    work page doi:10.1016/j.amc.2006.05.070 2007

  28. [28]

    Spitzbart,A generalization of Hermite’s interpolation formula, Amer

    A. Spitzbart,A generalization of Hermite’s interpolation formula, Amer. Math. Mon.67(1960), no. 1, 42-46.https: //doi.org/10.1080/00029890.1960.11989446

    work page doi:10.1080/00029890.1960.11989446 1960

  29. [29]

    Stoer and R

    J. Stoer and R. Bulirsch,Introduction to Numerical Analysis 2nd Edition, Springer-Verlag, New York, 1993

    1993

  30. [30]

    Voit,Some lower bounds for distances of roots of classical orthogonal polynomials, preprint, 2022, arXiv:2109.14999

    M. Voit,Some lower bounds for distances of roots of classical orthogonal polynomials, preprint, 2022, arXiv:2109.14999

    work page arXiv 2022

  31. [31]

    Wolfram Research, Inc.,Mathematica, Version 14.1, Champaign, IL (2024).https://www.wolfram.com/mathematica. 31

    2024