Pith. sign in

REVIEW 2 major objections 2 minor 47 references

An Unconditionally Linearly Convergent ADMM Approach for the Allen-Cahn Equation with Flory-Huggins Potential

T0 review · 2 major / 2 minor · reviewed 2026-06-29 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read An ADMM solver for the Allen-Cahn equation with Flory-Huggins potential converges unconditionally and at a linear rate when the multiplier update step size lies in (0, (√5 + 1)/2).

desk verdict The paper extends an ADMM convergence result from Cahn-Hilliard to Allen-Cahn but the structural differences between the discretizations make the direct transfer of the alpha range and linear rate a point that needs explicit verification in the proofs. read the letter →

arxiv 2605.26623 v1 pith:R45KGALF submitted 2026-05-26 math.NA cs.NA

Peng Jiang , Shengtong Liang , Tiao Lu This is my paper
classification math.NAcs.NA
keywords Allen-CahnequationFlory-HugginspotentialADMMunconditionalconvergencelinearconvexsplittingphasefieldmodelsingularnonlinearsystem
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper proves that an alternating direction method of multipliers solver, applied after convex splitting discretization of the Allen-Cahn equation with logarithmic Flory-Huggins potential, converges for any time step when the multiplier step size α stays below the golden ratio and does so at a linear rate. This removes the need for artificial time-step restrictions or strict positivity conditions that usually constrain simulations of phase separation. The argument adapts an existing unconditional ADMM framework from the Cahn-Hilliard case to the singular nonlinear algebraic systems that arise here while preserving bound preservation and discrete energy stability. A reader would care because the result lets practitioners run stable, accurate long-time simulations of phase-field models without tuning parameters to satisfy hidden separation assumptions.

What carries the argument

The ADMM iteration with multiplier update step size α, applied to the singular nonlinear system obtained from convex splitting of the Allen-Cahn equation.

What would settle it

A concrete numerical run of the scheme on a one-dimensional problem with α slightly larger than (√5 + 1)/2 that either diverges or fails to exhibit a linear convergence rate independent of the time step size.

Watch

Extended reading notes

Core claim

The central claim is that the embedded ADMM iteration converges unconditionally for multiplier update step sizes α in the open interval (0, (√5 + 1)/2) and that the convergence is linear; the proof proceeds by minor adaptation of the prior Cahn-Hilliard ADMM analysis to the singular system produced by convex splitting of the Allen-Cahn equation with Flory-Huggins potential.

Load-bearing premise

The singular nonlinear system produced by convex splitting of the Allen-Cahn equation with Flory-Huggins potential is close enough in structure to the corresponding Cahn-Hilliard system that the existing ADMM convergence proof transfers after only minor changes.

Editorial extensions

If this is right

  • The solver can be applied at any time-step size without loss of convergence.
  • Bound preservation and discrete energy stability hold for the computed solutions.
  • No artificial separation condition away from the singularities at 0 and 1 is required.
  • The linear rate is independent of the spatial mesh size once α is fixed inside the stated interval.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same step-size restriction may apply to other singular potentials whose convex splitting yields algebraically similar systems.
  • In three-dimensional simulations the removal of time-step limits could reduce total wall-clock time by an order of magnitude for long-time coarsening studies.
  • An adaptive choice of α near the upper end of the interval might further accelerate practical runs while remaining provably convergent.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a convex-splitting time discretization of the Allen-Cahn equation with Flory-Huggins potential that yields a singular nonlinear algebraic system at each step, then introduces an ADMM solver for that system. Building on the ADMM framework of Li et al. (2026) for the Cahn-Hilliard equation, the authors claim two main theorems: unconditional convergence of the iteration for any multiplier step-size α ∈ (0, (√5 + 1)/2) and a rigorous linear convergence rate. The scheme is asserted to be bound-preserving and energy-stable; numerical experiments are presented to confirm the theoretical rates.

Significance. If the convergence statements are correctly established for the Allen-Cahn operator, the result supplies a practical, time-step-independent solver for a singular nonlinear system that arises in many phase-field applications. The explicit linear-rate guarantee and the removal of separation conditions would be a useful advance in the numerical analysis of singularly nonlinear phase-field models.

major comments (2)
  1. [Key theoretical contributions paragraph and the convergence analysis section (presumably §4)] The central claim that the unconditional convergence result for α ∈ (0, (√5 + 1)/2) carries over from the Cahn-Hilliard setting requires explicit verification that the contraction inequality used in Li et al. (2026) remains valid for the second-order Allen-Cahn operator. The CH discretization produces a fourth-order saddle-point structure whose proximal mapping and energy estimates exploit the bi-Laplacian; the AC convex splitting produces a second-order system whose Jacobian and monotonicity properties differ. The manuscript must either re-derive the key estimate (e.g., the bound on the multiplier update) for the AC case or demonstrate structural equivalence of the two saddle-point problems; without this step the α-range and linear-rate claims rest on an unverified transfer.
  2. [Linear convergence theorem (likely Theorem 4.2 or equivalent)] The linear convergence theorem likewise depends on the same contraction mapping. If the contraction constant derived for the CH operator does not hold verbatim for the AC discretization, the linear-rate statement may require a different α interval or additional assumptions on the spatial mesh; the manuscript should state the precise contraction factor obtained after any necessary adaptation.
minor comments (2)
  1. [Abstract] The abstract contains a stray '%' character ('convex splitting schemes to %preserve this bound'), which appears to be an artifact of LaTeX editing and should be removed.
  2. [Preliminaries and notation] Notation for the discrete inner product, the proximal mapping, and the multiplier update should be introduced with a single consistent symbol set in the preliminary section rather than redefined inline in the convergence proofs.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on the transfer of the ADMM convergence analysis. We agree that the differing operator orders require explicit verification and will revise the manuscript to strengthen this part of the presentation.

read point-by-point responses

  1. Referee: [Key theoretical contributions paragraph and the convergence analysis section (presumably §4)] The central claim that the unconditional convergence result for α ∈ (0, (√5 + 1)/2) carries over from the Cahn-Hilliard setting requires explicit verification that the contraction inequality used in Li et al. (2026) remains valid for the second-order Allen-Cahn operator. The CH discretization produces a fourth-order saddle-point structure whose proximal mapping and energy estimates exploit the bi-Laplacian; the AC convex splitting produces a second-order system whose Jacobian and monotonicity properties differ. The manuscript must either re-derive the key estimate (e.g., the bound on the multiplier update) for the AC case or demonstrate structural equivalence of the two saddle-point problems; without this step the α-range and linear-rate claims rest on an unverified transfer.

    Authors: We acknowledge the structural differences between the second-order Allen-Cahn and fourth-order Cahn-Hilliard saddle-point problems. Section 4 adapts the contraction argument by exploiting the shared convexity of the Flory-Huggins splitting and the monotonicity of the proximal mapping for the AC operator. However, the adaptation steps are not written out in full detail. We will revise by inserting a new subsection that re-derives the key multiplier-update bound directly for the AC discretization, confirming that the identical α interval (0, (√5 + 1)/2) is admissible without invoking bi-Laplacian estimates. revision: yes

  2. Referee: [Linear convergence theorem (likely Theorem 4.2 or equivalent)] The linear convergence theorem likewise depends on the same contraction mapping. If the contraction constant derived for the CH operator does not hold verbatim for the AC discretization, the linear-rate statement may require a different α interval or additional assumptions on the spatial mesh; the manuscript should state the precise contraction factor obtained after any necessary adaptation.

    Authors: We agree that the precise contraction factor must be stated after adaptation. In the revision we will compute the explicit contraction constant for the AC case (which turns out to be the same as in the CH setting because the relevant monotonicity constant is controlled by the convex splitting rather than the differential order) and will record it explicitly in the statement of the linear-convergence theorem, together with a short remark confirming that no additional mesh assumptions are needed. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; independent proofs claimed for Allen-Cahn case

full rationale

The paper cites Li et al. (2026) for the base ADMM framework on the related Cahn-Hilliard equation but explicitly presents its own key contributions as new proofs of unconditional convergence (for α in (0, (√5+1)/2)) and linear convergence specifically for the Allen-Cahn discretization with Flory-Huggins potential. The cited work involves different authors and a distinct fourth-order system, so the extension does not reduce to self-citation load-bearing or self-definitional equivalence. No quotes or equations in the provided material show a prediction reducing to a fitted input by construction, an ansatz smuggled via citation, or renaming of known results. The derivation is presented as self-contained new analysis building on an external reference.

Assumptions & free parameters 0 free parameters · 0 assumptions · 0 invented entities

Only the abstract is available; no explicit free parameters, axioms, or invented entities are stated in the provided text.

how reviews work

0 comments
Cite this review

Pith. "Pith review of An Unconditionally Linearly Convergent ADMM Approach for the Allen-Cahn Equation with Flory-Huggins Potential." pith.science (2026). https://pith.science/paper/R45KGALF

@misc{pith2026260526623,
  author       = {Pith},
  title        = {Pith review of: An Unconditionally Linearly Convergent ADMM Approach for the Allen-Cahn Equation with Flory-Huggins Potential},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/R45KGALF}},
  note         = {Machine review of arXiv:2605.26623}
}
abstract

The Allen-Cahn equation with Flory-Huggins potential is a fundamental and crucial model in phase field simulation for describing phase separation phenomena, which serves as a core tool in diverse branches of natural sciences. The numerical simulation of the Allen-Cahn equation is of great importance but poses significant challenges due to the strong nonlinearity and the presence of logarithmic singularities at $u=0,1$ in the Flory-Huggins potential. In this paper, we consider convex splitting schemes to %preserve this bound and guarantee unconditional unique solvability, which reduces the numerical simulation to solving a singular nonlinear system arising from spatial discretization at each time step. We propose an iterative solver that is specifically designed for such systems based on the alternating direction method of multipliers (ADMM) approach. The scheme possesses properties such as bound preserving and discrete energy stability. Building upon the recent unconditionally convergent ADMM framework for the Cahn-Hilliard equation (Li et al., 2026), our key theoretical contributions are twofold: (a) a proof of unconditional convergence when the multiplier update step size $\alpha \in (0,\frac{\sqrt{5}+1}{2})$; (b) a rigorous establishment of the linear convergence for the embedded ADMM solver. This effectively liberates the solver from time-step constraints or strict separation conditions. Comprehensive numerical experiments validate our proposed ADMM framework, where its theoretical predictions are fully substantiated in practice, showcasing efficiency and robustness.

Figures

Figures reproduced from arXiv: 2605.26623 by the authors.

Figure 1
Figure 1. Convergence test [PITH_FULL_IMAGE:figures/full_fig_p021_1.png] view at source ↗
Figure 2
Figure 2. The dynamical evolution over time t. The initial condition is (5.7) u0(x, y) = 0.01 + 0.98 rand(x, y). where rand(x, y) is the function producing the random numbers in (0, 1). We set the convergence criterion for Algorithm 1 as: (5.8) max{r, s} ≤ γ = 10−8 . For θ = 3.0, the number of iterations, the energy, 1 − max(u) and min(u) over time t are given in [PITH_FULL_IMAGE:figures/full_fig_p022_2.png] view at source ↗
Figure 3
Figure 3. 2D simulation The parameter ϵ is chosen as the variable parameter in this experiment, which is taken as ϵ = 0.05, 0.10, 0.15. The initial condition is (5.10) u0(x, y, z) = 0.45 + 0.10 rand(x, y, z). where rand(x, y, z) is the function producing the random numbers in (0, 1). The convergence criterion is: (5.11) max{r, s} ≤ γ = 10−8 . For ϵ = 0.05, the number of iterations, the energy, 1 − max(u) and min(u) over time … view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: 2D simulation: the dynamical evolution over time t. (b) θ = 4.0 [PITH_FULL_IMAGE:figures/full_fig_p024_4.png]
Figure 4
Figure 4. Figure 4: 2D simulation: the dynamical evolution over time t. 6. Conclusion In this paper, we present a novel iterative solver for the nonlinear system at each time step, which naturally arises from applying the convex splitting scheme to the Allen-Cahn equation with the Flory-H…
Figure 5
Figure 5. Figure 5: 3D simulation remains an open theoretical challenge in the field [17,43] that warrants future study. (c) Conducting comprehensive benchmark comparisons. Although our theoretical convergence is rigorously validated, evaluating the empirical computational effi￾ciency aga…
Figure 6
Figure 6. Figure 6: 3D simulation: the dynamical evolution over time t. (b) ϵ = 0.10 [PITH_FULL_IMAGE:figures/full_fig_p027_6.png]
Figure 6
Figure 6. Figure 6: 3D simulation: the dynamical evolution over time t. [7] W. Chen, C. Wang, X. Wang, and S. M. Wise, Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential, Journal of Computational Physics: X 3 (2019), 100031. [8…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

47 extracted references · 1 canonical work pages

  1. [1]

    S. M. Allen and J. W. Cahn, A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening , Acta Metallurgica 27 (1979), no. 6, 1085–1095

    1979

  2. [2]

    Baskaran, J

    A. Baskaran, J. Lowengrub, C. Wang, and S. M. Wise, Convergence analysis of a second order convex splitting scheme for the modified phase field crystal equation , SIAM Journal on Numerical Analysis 51 (2013), no. 5, 2851–2873

    2013

  3. [3]

    Boyd and L

    S. Boyd and L. Vandenberghe, Convex optimization , Cambridge University Press, 2004

    2004

  4. [4]

    Chen, Phase-field models for microstructure evolution , Annual Review of Materials Re- search 32 (2002), no

    L. Chen, Phase-field models for microstructure evolution , Annual Review of Materials Re- search 32 (2002), no. 1, 113–140

    2002

  5. [5]

    L. Chen, D. Sun, and K. Toh, A note on the convergence of ADMM for linearly con- strained convex optimization problems , Computational Optimization and Applications 66 (2017), no. 2, 327–343

    2017

  6. [6]

    W. Chen, S. Conde, C. Wang, X. Wang, and S. M. Wise, A linear energy stable scheme for a thin film model without slope selection , Journal of Scientific Computing 52 (2012), no. 3, 546–562. A LINEARLY CONVERGENT SOL VER FOR THE AC EQUATION WITH FH POTENTIAL 27 (a) ϵ = 0.05 Figure 6. 3D simulation: the dynamical evolution over time t. (b) ϵ = 0.10 Figure 6...

    2012

  7. [7]

    W. Chen, C. Wang, X. Wang, and S. M. Wise, Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential , Journal of Computational Physics: X 3 (2019), 100031

    2019

  8. [8]

    Chen, Generation and propagation of interfaces for reaction-diffusion equations , Journal of Differential Equations 96 (1992), no

    X. Chen, Generation and propagation of interfaces for reaction-diffusion equations , Journal of Differential Equations 96 (1992), no. 1, 116–141

    1992

Show all 47 references
  1. [9]

    Concus, G

    P. Concus, G. H. Golub, and D. P. O’Leary, A generalized conjugate gradient method for the numerical solution of elliptic partial differential equations , Sparse matrix computations, 1976, pp. 309–332

    1976

  2. [10]

    J. W. Cooley and J. W. Tukey, An algorithm for the machine calculation of complex Fourier series, Mathematics of Computation 19 (1965), no. 90, 297–301

    1965

  3. [11]

    J. E. Dennis and R. B. Schnabel, Numerical methods for unconstrained optimization and nonlinear equations, SIAM, 1996

    1996

  4. [12]

    Du and X

    Q. Du and X. Feng, The phase field method for geometric moving interfaces and their nu- merical approximations, Handbook of Numerical Analysis 21 (2020), 425–508

    2020

  5. [13]

    Q. Du, L. Ju, X. Li, and Z. Qiao, Maximum principle preserving exponential time differencing schemes for the nonlocal Allen–Cahn equation , SIAM Journal on Numerical Analysis 57 (2019), no. 2, 875–898

    2019

  6. [14]

    2, 317–359

    , Maximum bound principles for a class of semilinear parabolic equations and expo- nential time-differencing schemes , SIAM Review 63 (2021), no. 2, 317–359

    2021

  7. [15]

    D. J. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation , MRS Online Proceedings Library (OPL) 529 (1998), 39

    1998

  8. [16]

    Flory, Thermodynamics of high polymer solutions , Journal of Chemical Physics 10 (1942), no

    P. Flory, Thermodynamics of high polymer solutions , Journal of Chemical Physics 10 (1942), no. 1, 51–61

    1942

  9. [17]

    Gomez and T

    H. Gomez and T. J. R. Hughes, Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models , Journal of Computational Physics 230 (2011), no. 13, 5310–5327

    2011

  10. [18]

    Guillén-González and G

    F. Guillén-González and G. Tierra, On linear schemes for a Cahn–Hilliard diffuse interface model, Journal of Computational Physics 234 (2013), 140–171. A LINEARLY CONVERGENT SOL VER FOR THE AC EQUATION WITH FH POTENTIAL 29

    2013

  11. [19]

    J. Kou, S. Sun, and X. Wang, A novel energy factorization approach for the diffuse-interface model with Peng–Robinson equation of state , SIAM Journal on Scientific Computing 42 (2020), no. 1, B30–B56

    2020

  12. [20]

    Kroemer and T

    M. Kroemer and T. Laux, Quantitative convergence of the nonlocal Allen–Cahn equation to volume-preserving mean curvature flow , Mathematische Annalen 391 (2025), no. 3, 4455– 4472

    2025

  13. [21]

    D. Li, C. Quan, and J. Xu, Stability and convergence of Strang splitting. Part I: scalar Allen-Cahn equation , Journal of Computational Physics 458 (2022), 111087

    2022

  14. [22]

    R. Li, S. Liang, and Z. Qiao, Overcoming logarithmic singularities in the Cahn-Hilliard equa- tion with Flory-Huggins potential: An unconditionally convergent ADMM approach , arXiv preprint arXiv:2506.09361 (2026, accepted for publication by Mathematics of Computation)

    2026

  15. [23]

    MacNamara and G

    S. MacNamara and G. Strang, Operator splitting, Splitting methods in communication, imag- ing, science, and engineering, 2017, pp. 95–114

    2017

  16. [24]

    K. W. Morton and D. F. Mayers, Numerical solution of partial differential equations , Cam- bridge University Press, 1994

    1994

  17. [25]

    P. Neal, E. Chu, B. Peleato, and J. Eckstein, Distributed optimization and statistical learning via the alternating direction method of multipliers , Foundations and Trends® in Machine learning 3 (2011), no. 1, 1–122

    2011

  18. [26]

    Nishiura and H

    Y. Nishiura and H. Fujii, Stability of singularly perturbed solutions to systems of reaction- diffusion equations , SIAM journal on Mathematical Analysis 18 (1987), no. 6, 1726–1770

    1987

  19. [27]

    J. M. Ortega and W. C. Rheinboldt, Iterative solution of nonlinear equations in several variables, SIAM, 2000

    2000

  20. [28]

    Provatas and K

    N. Provatas and K. Elder, Phase-field methods in materials science and engineering , John Wiley & Sons, 2011

    2011

  21. [29]

    R. T. Rockafellar, The multiplier method of Hestenes and Powell applied to convex program- ming, Journal of Optimization Theory and Applications 12 (1973), no. 6, 555–562

    1973

  22. [30]

    W. R. Rudin, Principles of mathematical analysis , McGraw-Hill, 1976

    1976

  23. [31]

    J. Shen, T. Tang, and J. Yang, On the maximum principle preserving schemes for the gen- eralized Allen–Cahn equation , Communications in Mathematical Sciences 14 (2016), no. 6, 1517–1534

    2016

  24. [32]

    J. Shen, C. Wang, X. Wang, and S. M. Wise, Second-order convex splitting schemes for gradient flows with Ehrlich–Schwoebel type energy: application to thin film epitaxy , SIAM Journal on Numerical Analysis 50 (2012), no. 1, 105–125

    2012

  25. [33]

    J. Shen, J. Xu, and J. Yang, The scalar auxiliary variable (SA V) approach for gradient flows , Journal of Computational Physics 353 (2018), 407–416

    2018

  26. [34]

    Stehlík and J

    P. Stehlík and J. Volek, Maximum principles for discrete and semidiscrete reaction-diffusion equation, Discrete Dynamics in Nature and Society 2015 (2015), no. 1, 791304

    2015

  27. [35]

    Steinbach, Phase-field model for microstructure evolution at the mesoscopic scale , Annual Review of Materials Research 43 (2013), no

    I. Steinbach, Phase-field model for microstructure evolution at the mesoscopic scale , Annual Review of Materials Research 43 (2013), no. 1, 89–107

    2013

  28. [36]

    Strang, On the construction and comparison of difference schemes , SIAM Journal on Numerical Analysis 5 (1968), no

    G. Strang, On the construction and comparison of difference schemes , SIAM Journal on Numerical Analysis 5 (1968), no. 3, 506–517

    1968

  29. [37]

    Tang and J

    T. Tang and J. Yang, Implicit-explicit scheme for the Allen-Cahn equation preserves the maximum principle , Journal of Computational Mathematics (2016), 451–461

    2016

  30. [38]

    Tryggvason, R

    G. Tryggvason, R. Scardovelli, and S. Zaleski, Direct numerical simulations of gas–liquid multiphase flows , Cambridge university press, 2011

    2011

  31. [39]

    S. Wang, R. Sekerka, A. Wheeler, B. Murray, S. Coriell, R. Braun, and G. McFadden, Thermodynamically-consistent phase-field models for solidification , Physica D: Nonlinear Phenomena 69 (1993), no. 1-2, 189–200

    1993

  32. [40]

    X. Wang, J. Kou, and J. Cai, Stabilized energy factorization approach for Allen–Cahn equation with logarithmic Flory–Huggins potential , Journal of Scientific Computing 82 (2020), no. 2, 25

    2020

  33. [41]

    Weierstrass, Mathematische werke von Karl Weierstrass: bd

    K. Weierstrass, Mathematische werke von Karl Weierstrass: bd. Abhandlungen I. 1894 , Vol. 1, Mayer & Müller, 1894

  34. [42]

    Wheeler, W

    A. Wheeler, W. Boettinger, and G. McFadden, Phase-field model for isothermal phase tran- sitions in binary alloys , Physical Review A 45 (1992), no. 10, 7424

    1992

  35. [43]

    X. Yang, Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends , Journal of Computational Physics 327 (2016), 294–316. 30 PENG JIANG, SHENGTONG LIANG, AND TIAO LU

    2016

  36. [44]

    X. Yang, J. Zhao, and Q. Wang, Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method , Journal of Computational Physics 333 (2017), 104–127

    2017

  37. [45]

    M. Yuan, W. Chen, C. Wang, S. M. Wise, and Z. Zhang, A second order accurate in time, energy stable finite element scheme for the Flory-Huggins-Cahn-Hilliard equation , Advances in Applied Mathematics and Mechanics 14 (2022), no. 6

    2022

  38. [46]

    J. Zhao, Q. Wang, and X. Yang, Numerical approximations for a phase field dendritic crystal growth model based on the invariant energy quadratization approach , International Journal for Numerical Methods in Engineering 110 (2017), no. 3, 279–300

    2017

  39. [47]

    J. Zhao, X. Yang, Y. Gong, X. Zhao, X. Yang, J. Li, and Q. Wang, A general strategy for numerical approximations of non-equilibrium models–Part I: Thermodynamical systems , International Journal of Numerical Analysis and Modeling 15 (2018), no. 6, 884–918. Corresponding author...

    2018

Pith tools

Reviewed June 29, 2026 · model on record in the stance chip above.