pith. machine review for the scientific record. sign in

arxiv: 2605.11842 · v1 · submitted 2026-05-12 · 🧮 math.NA · cs.NA

Recognition: 2 theorem links

· Lean Theorem

Dynamics of the Longest-Edge Altitude Bisection Algorithm

Authors on Pith no claims yet

Pith reviewed 2026-05-13 05:17 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords longest-edge bisectionaltitude bisectiontriangle shape spacemesh refinementright-triangle geodesicdiscretizationtriangulation
0
0 comments X

The pith

Longest-edge altitude bisection collapses the entire triangle shape space onto the right-triangle geodesic in one step.

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

The paper studies the longest-edge altitude bisection scheme for refining triangulations, in which each triangle is split by the altitude dropped to its longest edge. Using the normalized shape space of triangles, the authors give an explicit geometric construction for the child triangles as intersections of rays from the longest-edge endpoints with the right-triangle geodesic. This construction immediately shows that one round of bisection maps every possible triangle shape onto that geodesic, while every point on the geodesic remains unchanged by further bisections. The work also supplies two-sided bounds on the rate at which the mesh diameter shrinks under repeated application of the rule.

Core claim

The longest-edge altitude bisection admits a geometric description in the normalized shape space where the left and right children are obtained by intersecting the right-triangle geodesic with rays from the longest edge endpoints. This characterization shows that the refinement dynamics collapse the entire shape space onto the right-triangle geodesic in a single step, with every point on this geodesic being a fixed point. Explicit mappings and two-sided bounds for the contraction of the mesh size are also derived.

What carries the argument

The normalized shape space of triangles together with its geodesic of right triangles, which serves as the immediate attractor under the altitude-bisection mapping.

Load-bearing premise

The normalized shape space accurately records the geometric change produced by altitude bisection on any triangle.

What would settle it

Pick any acute or obtuse triangle, apply the longest-edge altitude bisection once, compute the normalized shape parameters of both children, and verify whether both parameters lie exactly on the right-triangle geodesic.

Figures

Figures reproduced from arXiv: 2605.11842 by J\'er\^ome Michaud, Sergey Korotov.

Figure 1
Figure 1. Figure 1: Illustration of the LEAB process. Triangle [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Piecewise definition of the left and right LEAB maps [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Illustration of the LEAB process on an equilateral and a right triangle. [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
read the original abstract

We study a longest-edge based refinement scheme for triangulations, termed the longest-edge altitude bisection (LEAB), in which each triangle is subdivided by dropping the altitude from the vertex opposite to its longest edge. Using the normalized shape space of triangles introduced by Perdomo and Plaza in: Properties of triangulations obtained by the longest-edge bisection. \emph{Cent. Eur. J. Math.}, 12(12) (2014), 1796-1810, we show that LEAB admits a simple geometric description: the normalized left and right children of a triangle in focus are obtained by intersecting the geodesic of right triangles with rays issued from the endpoints of the longest edge and explicit formulas for the mappings are derived. This characterization implies an interesting observation that the associated refinement dynamics collapse the entire shape space onto the right-triangle geodesic in a single step and that every point on this geodesic is fixed. Two-sided bounds for the contraction of the mesh size (discretization parameter) are derived. Also, applications and limitations of the method are briefly discussed.

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 / 1 minor

Summary. The manuscript studies the longest-edge altitude bisection (LEAB) refinement scheme for triangulations. In the normalized shape space of triangles introduced by Perdomo and Plaza (2014), it derives explicit mappings for the normalized left and right children by intersecting rays from the endpoints of the longest edge with the geodesic of right triangles. This leads to the claim that the refinement dynamics collapse the entire shape space onto the right-triangle geodesic in a single step, with every point on the geodesic being a fixed point, and two-sided bounds are derived for the contraction of the mesh size (discretization parameter).

Significance. If the geometric characterization and collapse property hold over the claimed domain, the work supplies a clean, explicit description of the refinement action together with contraction bounds that could be useful for analyzing convergence rates in adaptive finite-element computations. The reduction to a one-dimensional geodesic and the fixed-point property are attractive features when valid.

major comments (2)
  1. [Abstract] Abstract and the geometric-description paragraph: the claim that the ray-intersection construction yields the normalized children for every triangle in the shape space is not justified for obtuse triangles. When the angle opposite the longest edge exceeds 90°, the altitude foot lies outside the edge segment, so the actual subdivision does not produce the two right triangles predicted by the intersection; the universal collapse and the fixed-point property therefore hold at most on the acute/right subset.
  2. [bounds derivation] The two-sided mesh-size contraction bounds inherit the same restriction. Because the bounds are stated for the full space and the derivation relies on the same ray-intersection maps, the bounds require an explicit domain restriction or a separate argument for obtuse cases.
minor comments (1)
  1. The 2014 reference is cited only by title and journal; a full bibliographic entry with page numbers should be supplied.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for identifying the important restriction to acute and right triangles that was not made explicit in the manuscript. We will revise the abstract, geometric description, and bounds sections to incorporate the domain limitation and clarify the applicability of the collapse and fixed-point properties.

read point-by-point responses
  1. Referee: [Abstract] Abstract and the geometric-description paragraph: the claim that the ray-intersection construction yields the normalized children for every triangle in the shape space is not justified for obtuse triangles. When the angle opposite the longest edge exceeds 90°, the altitude foot lies outside the edge segment, so the actual subdivision does not produce the two right triangles predicted by the intersection; the universal collapse and the fixed-point property therefore hold at most on the acute/right subset.

    Authors: We agree that the ray-intersection construction assumes the altitude foot lies inside the longest-edge segment, which holds precisely when the triangle is acute or right-angled at the base vertices. For obtuse triangles with the obtuse angle opposite the longest edge the foot lies outside the segment, so the children are not the two right triangles given by the intersection. Consequently the one-step collapse onto the right-triangle geodesic and the fixed-point property are valid only on the acute/right subset of the shape space. We will revise the abstract and the geometric-description paragraph to state the domain restriction explicitly and, if space allows, add a brief note on the distinct behavior for obtuse triangles. revision: yes

  2. Referee: [bounds derivation] The two-sided mesh-size contraction bounds inherit the same restriction. Because the bounds are stated for the full space and the derivation relies on the same ray-intersection maps, the bounds require an explicit domain restriction or a separate argument for obtuse cases.

    Authors: We concur that the two-sided contraction bounds for the discretization parameter are obtained from the same ray-intersection maps and therefore inherit the same domain restriction. We will revise the bounds section to restrict the statements of the theorems to acute and right triangles, update the accompanying text, and note the limitation for obtuse triangles. A separate argument for the obtuse case will be added if it can be derived concisely; otherwise the restriction will be stated as a limitation of the present analysis. revision: yes

Circularity Check

0 steps flagged

No circularity; derivation uses independent external shape space and explicit geometric mappings

full rationale

The paper adopts the normalized shape space of triangles from the independent 2014 reference by Perdomo and Plaza (distinct authors) and derives the LEAB child mappings via explicit ray intersections with the right-triangle geodesic. The collapse property, fixed-point behavior, and two-sided mesh-size contraction bounds follow directly from these mappings and the geometric construction without any reduction to fitted parameters, self-citations, or redefinition of inputs. All load-bearing steps are self-contained mathematical derivations within the adopted coordinate system.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the geometric properties of the normalized shape space defined in the 2014 Perdomo-Plaza paper; no free parameters, new entities, or ad-hoc axioms are introduced in the abstract.

axioms (1)
  • domain assumption The normalized shape space of triangles introduced by Perdomo and Plaza (2014) correctly represents the geometry of altitude bisection.
    Invoked to obtain the geodesic description and collapse result.

pith-pipeline@v0.9.0 · 5485 in / 1186 out tokens · 25394 ms · 2026-05-13T05:17:54.958440+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

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

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.

Reference graph

Works this paper leans on

14 extracted references · 14 canonical work pages · 1 internal anchor

  1. [1]

    E. L. Allgower and K. Georg.Numerical Continuation Methods: An Introduction, volume 13 ofSpringer Series in Computational Mathematics. Springer, Berlin, 1990. 4

  2. [2]

    Anderson.Hyperbolic geometry

    J. Anderson.Hyperbolic geometry. Springer Science & Business Media, 2005

  3. [3]

    Apel.Anisotropic Finite Elements: Local Estimates and Applications

    T. Apel.Anisotropic Finite Elements: Local Estimates and Applications. B. G. Teubner, Stuttgart, 1999

  4. [4]

    Babuˇ ska and A

    I. Babuˇ ska and A. K. Aziz. On the angle condition in the finite element method.SIAM Journal on Numerical Analysis, 13(2):214–226, 1976

  5. [5]

    R. E. Bank and R. K. Smith. A posteriori error estimates based on hierarchical bases.SIAM Journal on Numerical Analysis, 30(4):921–935, 1993

  6. [6]

    On the stability, complexity, and distribution of similarity classes of the longest edge bisection process for triangles

    D. Kalmanovich and Y. Solomon. On the stability, complexity, and distribution of similarity classes of the longest edge bisection process for triangles.arXiv preprint arXiv:2601.13663, 2026

  7. [7]

    W. F. Mitchell. Adaptive refinement for arbitrary finite-element spaces with hierarchical bases.Journal of Computational and Applied Mathematics, 36(1):65–78, 1991

  8. [8]

    Perdomo and A

    F. Perdomo and A. Plaza. A new proof of the degeneracy property of the longest-edgen-section refinement scheme for triangular meshes.Applied Mathematics and Computation, 219:2342–2344, 2012

  9. [9]

    Perdomo and A

    F. Perdomo and A. Plaza. Proving the non-degeneracy of the longest-edge trisection by a space of triangular shapes with hyperbolic metric.Applied Mathematics and Computation, 221:424–432, 2013

  10. [10]

    Perdomo and A

    F. Perdomo and A. Plaza. Properties of triangulations obtained by the longest-edge bisection.Central European Journal of Mathematics, 12(12):1796–1810, 2014

  11. [11]

    M. C. Rivara. Mesh refinement processes based on the generalized bisection of simplices.SIAM Journal on Numerical Analysis, 21:604–613, 1984

  12. [12]

    I. G. Rosenberg and F. Stenger. A lower bound on the angles of triangles constructed by bisecting the longest side.Mathematics of Computation, 29:390–395, 1975

  13. [13]

    M. Stynes. On faster convergence of the bisection method for all triangles.Mathematics of Computation, 35:1195–1201, 1980

  14. [14]

    C. T. Traxler. An algorithm for adaptive mesh refinement in two dimensions.Computing, 59(2):115–137, 1997. 5