REVIEW 2 major objections 7 minor 2 cited by
The geometry of the deep linear network
T0 review · 2 major / 7 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Deep linear network training is exactly Riemannian gradient descent combined with a Boltzmann entropy that selects among degenerate minimizers.
desk verdict A valuable expository synthesis of deep linear network geometry, but its central thermodynamic theorems are deferred to an unavailable companion paper, so the manuscript is not self-contained. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is a chain of geometric identifications. The G-balanced varieties are the invariant sets of the gradient flow; the SVD parametrization $W_p = Q_p \Lambda Q_{p-1}^T$, with $\Lambda=\Sigma^{1/N}$, exhibits each balanced manifold as orbits of the orthogonal group and makes tangent directions explicit. The metric $g_N$ is carried from the Frobenius inner product upstairs by Riemannian submersion, and its inverse operator $A_{N,W}$ is diagonalized by the singular-vector matrices $u_k v_l^T$ with eigenvalue $(\sigma_k^2-\sigma_l^2)/(\sigma_k^{2/N}-\sigma_l^{2/N})$ for $k\neq l$ and $N \sigma_k^{2-2/N}$ on the diagonal blocks. The entropy $S(W)=\log \mathrm{vol}(O_W)$ is the logarithm of the group-orbit volume, computed from the pullback metric through an orthonormal basis built from Chebyshev polynomials; its closed form is a ratio of Vandermonde determinants. The stochastic side is carried by Brownian motion on the quotient manifold and an explicit Itô SDE whose drift contains the mean-curvature correction that becomes Coulomb repulsion in singular-value coordinates.
What would settle it
Run the balanced gradient flow with $E(W)=\frac12((W_{11}-1)^2+(W_{22}-1)^2)$ for $d=2$ from random full-rank balanced initial conditions and compare the empirical low-rank limits with the minimizers of $F_\beta(W)=E(W)-\beta^{-1}S(W)$. If the selected limit is not among the free-energy minimizers as $\beta$ and $N$ vary, or if rank does not drop, the entropic selection principle fails.
Extended reading notes
Core claim
The central discovery is that overparametrization in a DLN does not merely add redundant coordinates; it organizes phase space into invariant algebraic varieties $M_G$ defined by the quadratic balance equations $W_{p+1}^T W_{p+1} = W_p W_p^T - G_p$, with the case $G=0$ giving balanced manifolds $M_r$ of rank $r$. On $M_r$ the map from weights to the end-to-end matrix $W$ is a Riemannian submersion, so the flow of $W$ is exactly the Riemannian gradient flow of $E$ with respect to the metric $g_N(Z,Z)=\mathrm{Tr}(Z^T A_{N,W}^{-1} Z)$, where $A_{N,W}$ is a positive operator diagonalized in singular-vector coordinates by explicit eigenvalues. The same submersion structure makes the volume of a group orbit $O_W = \varphi^{-1}(W) \cap M$ computable in closed form as $c_d^{N-1} \sqrt{\mathrm{van}(\Sigma^2)/\mathrm{van}(\Sigma^{2/N})}$, and this volume defines the Boltzmann entropy $S(W)=\log \mathrm{vol}(O_W)$. With Brownian motion on $(M_d,g_N)$, these ingredients produce the Riemannian Langevin equation for $W$, whose drift is $-\mathrm{grad}_{g_N}F_\beta$ with $F_\beta = E - \beta^{-1} S$; the entropy gradient is expressed through a diagonal matrix $\Sigma'$ whose off-diagonal terms act like Coulombic repulsion between singular values. This gives a complete thermodynamic description of training in the full-rank, distinct-singular-value regime.
Load-bearing premise
The exact entropy formula, the metric identification, and the Riemannian Langevin equation are proved on the full-rank balanced manifold with distinct singular values; the rank-deficient limit that training actually approaches is not covered by these proofs.
Editorial extensions
If this is right
- On balanced manifolds, the end-to-end matrix follows a Riemannian gradient flow, so overparametrization does not change the energy-decay identity $dE/dt = -\|\mathrm{grad}_{g_N}E\|^2$; convergence to critical points happens at a rate set by the geometry, not by the number of parameters.
- For degenerate losses, minimizers of the free energy $F_\beta = E - \beta^{-1} S$ supply an entropic selection principle; the paper conjectures that the minimizing set is a single point $A_\beta$ and that training with round-off noise selects its $\beta\to\infty$ limit.
- Stochastic noise in the gauge directions is not a nuisance but the microscopic source of entropy: the Riemannian Langevin equation shows that tangential fluctuations produce a deterministic motion by mean curvature of the group orbits, exactly as in Dyson Brownian motion.
- The explicit formulas extend to an infinite-depth limit $N\to\infty$, where the metric and entropy survive through the operator $A_{\infty,W}=\int_0^1 (WW^T)^{1-s} Z (W^T W)^s ds$ even though the upstairs group-orbit picture disappears.
- The Bures-Wasserstein geometry of positive definite matrices arises as the $N=2$, $W_2=W_1^T$ specialization, connecting DLN training to optimal transport and to known gradient flows on matrix manifolds.
Reading between the lines
- The paper leaves implicit that if the thermodynamic description extends past full rank, implicit regularization is an entropic force: the volume of weight symmetries biases training toward macrostates with larger orbit volume, and this bias should be visible in the distribution of low-rank limits even before any explicit penalty is added.
- One testable extension the paper does not run is to compare the empirical limits of balanced gradient flow with tiny gauge noise against the minimizers of $F_\beta$ for the $d=2$ matrix-completion energy $E(W)=\frac12((W_{11}-1)^2+(W_{22}-1)^2)$; the theory predicts that the selected rank-one matrix $\begin{pmatrix}1&a\\1/a&1\end{pmatrix}$ is the free-energy minimizer for each $\beta$ and $N$.
- The anisotropic splitting of noise into gauge and observable directions predicts a concrete experimental signature: adding isotropic weight noise should change the selected low-rank limit in a computable way, whereas noiseless gradient flow should approach the entropic selection only through round-off errors.
- The Simons-cone connection suggests that the rank-deficient balanced variety may carry singular minimal-surface structure; one could test whether tangential noise on low-rank varieties drives them toward the balanced manifold by motion by curvature, which would explain rank collapse as a geometric phenomenon rather than a pure optimization effect.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This expository paper develops a geometric theory of training dynamics in deep linear networks (DLNs). The proofs in Section 7 establish Theorems 1–3: the phase space M_d^N is foliated by invariant G-balanced varieties; on balanced manifolds the end-to-end matrix evolves by the Riemannian gradient flow of the loss E with respect to the exactly computable metric g_N in (4.16), diagonalized explicitly in Lemma 1. The paper then constructs a thermodynamic picture: the fiber of the projection phi over W in the balanced manifold is an O_d^{N-1} orbit whose volume is claimed to be the Vandermonde-type formula (6.2) (Theorem 10); the metric g_N arises by Riemannian submersion (Theorem 13); an explicit orthonormal basis for the balanced manifold is given in (8.20)–(8.24) (Theorem 14); and a Riemannian Langevin equation with explicit free-energy gradient (12.12)–(12.13) is proposed as a stochastic extension, modeled on Dyson Brownian motion (Section 11). Theorems 10, 13, 14, and 18 are attributed to the author's unpublished companion manuscript [34] and are presented only as sketches.
Significance. If correct, this framework is a remarkably complete exact description of DLN training: an invariant foliation, a parameter-free Riemannian metric, entropy as the logarithm of a group-orbit volume, and an explicit noise model with a Coulomb-like repulsion term (12.8). The proofs that are actually present (Section 7, Lemmas 1–2) are direct, reproducible matrix computations, and the paper is candid about its modeling axiom (the RLE is phenomenological, Section 10.1) and about the rank-deficient gap (Section 14.4). I independently verified several central formulas: (7.13)–(7.14) follow from a geometric sum; (7.15) is the correct reciprocal metric; the interior (N−1)×(N−1) blocks of the pullback metric (8.10) have determinant 2^{N−1}(σ_k^2−σ_l^2)/(σ_k^{2/N}−σ_l^{2/N}), which yields (6.2); and (12.13) matches the g_N-gradient of −β^{−1}S computed from (6.2). The main weakness is verifiability: the headline quantitative claims are deferred to an unpublished, unlocatable companion [34], so the manuscript cannot currently serve as a standalone rigorous reference.
major comments (2)
- [§6.2, §8, §8.3, §12.2; ref. [34]] The paper's central new quantitative claims—Theorem 10 (entropy formula (6.2)), Theorem 14 (orthonormal basis (8.20)–(8.24)), and Theorem 18 (Brownian motion SDE (12.10)), together with the projected RLE (12.12)–(12.13)—are quoted from the author's companion manuscript [34], which has no venue or preprint identifier. The manuscript itself locates the omissions: Section 8 calls Theorem 14 'the technical core of [34]' and Section 12.2 says 'Then we show in [34] that…'. The determinant evaluation converting the block-tridiagonal metric (8.9)–(8.10) into (6.2), the Chebyshev-polynomial diagonalization behind (8.20)–(8.24), and the Itô/generator calculation behind (12.10) are not supplied. Because the entropy formula and the RLE are advertised in the abstract as the paper's contribution, this is a load-bearing verification gap rather than a presentation detail. I verified that the deferred claims are consistent with what is shown: the interior (N−1)×(N−1) blocks of (8.10) have determinant 2^{N−1}(σ_k^2−σ_l^2)/(σ_k^{2/N}−σ_l^{2/N}), which yields (6.2), and (12.13) agrees with the g_N-gradient of the entropy from (6.2). The concern is therefore completeness, not evident error. Please include the proofs (an appendix or supplement would suffice) or update [34] to a publicly verifiable form with results keyed to specific statements.
- [§5.3, §6.1, §13.3–§13.4, §14.1, §14.4; abstract] The entropy and RLE results are established only on the full-rank stratum with distinct singular values, where the SVD parametrization z is a local bijection (Section 5.3), and Section 14.4 explicitly concedes that singularities at repeated and repeated-zero singular values are not understood. Yet the motivating phenomena—rank drops during training (Section 13.3), convergence to low-rank energy minimizers (Section 14.1), and the suggested 'entropic origin for implicit regularization' (Sections 1 and 13.4)—live in the excluded rank-deficient regime, and the entropic selection principle (6.9) is explicitly conjectural (Section 6.3). I credit the paper for flagging the limitation in Sections 6.1 and 14.4; however, the abstract and introduction claim 'exact formulas for a Boltzmann entropy' and a thermodynamic framework without this scoping caveat. Please state in the abstract and introduction that the exact entropy and RLE formulas hold on the full-rank, distinct-singular-value stratum and that the rank-deficient and repeated-spectrum cases are open.
minor comments (7)
- [§6.2, Theorem 10] The theorem as stated ('For each W ∈ (M_d, g_N)') omits the hypothesis that the singular values of W be distinct; formula (6.2) and the parametrization (5.15)–(5.17) degenerate at repeated singular values (cf. Section 5.3). State the hypothesis explicitly.
- [§4.2, Eq. (4.8)] The range '1 ≤ p ≤ n' should be '1 ≤ p ≤ N − 1'; the integer n is not defined in this section.
- [§10.1, §12.1, abstract] The title 'RLE for stochastic gradient descent' and the abstract's phrase 'stochastic gradient descent of free energy' should be qualified: Section 10.1 states that the RLE is a phenomenological model that is 'not stochastic gradient descent (SGD) corresponding to batch processing'. Suggest rewording to avoid the implication of equivalence to practical SGD.
- [§14.1, Eqs. (14.3)–(14.4)] The symbol 'a ∈ R\{0}' is printed after the limiting matrix (14.4), but it parameterizes the family of rank-one matrices (14.3); move it to (14.3).
- [Ref. [34]] Reference [34] should be updated to include a preprint identifier or a 'to appear' venue; as printed it cannot be located, which compounds the verification gap described in Major Comment 1.
- [§3.3, §4.2, Eq. (4.20)] The symbol M_r is used both for the rank-r leaf of the balanced variety (Section 3.3) and for the rank-r matrix manifold (4.20); the relationship is explained, but the double use is a readability hazard.
- [§6.2, Eq. (6.3)] The phrase 'van(A) denoted' should read 'van(A) denotes'.
Circularity Check
No significant circularity: the balanced-variety and Riemannian-gradient results are derived from the group-action geometry; the entropy/RLE theorems are deferred to a same-author companion, which is a verification gap rather than a circular reduction.
full rationale
The derivations actually carried out in the manuscript are self-contained matrix computations: Theorem 1 (invariance of G-balanced varieties), Theorem 2 (reduced flow on balanced manifolds), Lemma 1 (spectral decomposition of A_{N,W}), and the explicit diagonal metric (7.15). Equation (4.15) is identified with a Riemannian gradient flow under the metric defined in (4.16); this is a definitional identification via the duality (4.18)-(4.19), not a fitted or circular step. The genuinely novel thermodynamic objects are not derived in this paper: Theorem 10 (entropy formula (6.2)), Theorem 14 (orthonormal basis (8.20)-(8.24)), and Theorem 18 (Brownian motion SDE (12.10)) are quoted from the same-author companion [34], with the paper stating 'This is the technical core of [34]' and 'Then we show in [34] that...'. That is an omitted-proof/verification gap: if [34] is unavailable or contains an error, the central formulas lack support in this manuscript. However, deferring a theorem to a companion paper is not equivalent to making the result identical to its inputs. No parameter is fitted to data, no prediction is a renamed fit, and no uniqueness theorem is imported to exclude alternatives. Section 14.4 honestly states that repeated and zero singular values are not understood; that is a scope limitation, not circularity. The self-citation is heavy and load-bearing, but the manuscript does not exhibit the equation-level reduction needed to call the derivation circular.
Assumptions & free parameters
free parameters (2)
- inverse temperature beta
- anisotropy parameter kappa
assumptions (3)
- standard math Existence and uniqueness of ODE and Ito SDE solutions, and Brownian motion on Riemannian manifolds.
- domain assumption Theorems 10, 13, 14, and 18 from Menon and Yu are accepted as established results.
- ad hoc to paper The Riemannian Langevin equation is a faithful model for noise in the null directions of the DLN.
Cite this review
Pith. "Pith review of The geometry of the deep linear network." pith.science (2026). https://pith.science/paper/ZWRBBPI5
@misc{pith2026241109004,
author = {Pith},
title = {Pith review of: The geometry of the deep linear network},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZWRBBPI5}},
note = {Machine review of arXiv:2411.09004}
}
read the original abstract
This article provides an expository account of training dynamics in the Deep Linear Network (DLN) from the perspective of the geometric theory of dynamical systems. Rigorous results by several authors are unified into a thermodynamic framework for deep learning. The analysis begins with a characterization of the invariant manifolds and Riemannian geometry in the DLN. This is followed by exact formulas for a Boltzmann entropy, as well as stochastic gradient descent of free energy using a Riemannian Langevin Equation. Several links between the DLN and other areas of mathematics are discussed, along with some open questions.
Figures
Figures from the paper (2 more)
Forward citations
Cited by 2 Pith papers
-
Implicit Bias of SGD in Multivariate ReLU Networks: Effective Width Collapse
Noisy SGD in the mean-field regime forces wide multivariate ReLU networks to an effective width of at most 2P-1, yielding a continuous piecewise-affine predictor whose hyperplanes are non-redundant with respect to the...
-
The Effect of Depth on the Expressivity of Deep Linear State-Space Models
Deeper linear state-space models can match the expressive power of shallower wider ones, and under norm constraints they do so with exponentially smaller parameter norms.
Reference graph
Works this paper leans on
-
[34]
G. Menon and T. Yu , An entropy formula for the deep linear network , 2024
work page 2024
- [1]
- [2]
- [3]
-
[4]
S. Arora, N. Cohen, W. Hu, and Y. Luo , Implicit regularization in deep matrix factoriza- tion, in Advances in Neural Information Processing Systems 32: Annual Conference on Neural Information Processing Systems 2019, NeurIPS 2019, December 8-14, 2019, Vancouver, BC, Canada, 2019, pp. 7411–7422
work page 2019
-
[5]
B. Bah, H. Rauhut, U. Terstiege, and M. Westdickenberg , Learning deep linear neural networks: Riemannian gradient flows and convergence to global minimizers , Information and Inference: A Journal of the IMA, 11 (2022), pp. 307–353
2022
-
[6]
J. Baik, P. Deift, and T. Suidan , Combinatorics and random matrix theory , vol. 172 of Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, 2016
work page 2016
-
[7]
P. Baldi and K. Hornik, Neural networks and principal component analysis: Learning from examples without local minima , Neural networks, 2 (1989), pp. 53–58
work page 1989
Show all 40 references
-
[8]
P. F. Baldi and K. Hornik , Learning in linear neural networks: A survey , IEEE Transac- tions on neural networks, 6 (1995), pp. 837–858
1995
-
[9]
E. A. Barbaˇ sin and N. N. Krasovski˘ i , On stability of motion in the large , Doklady Akad. Nauk SSSR (N.S.), 86 (1952), pp. 453–456
1952
-
[10]
Bhatia, T
R. Bhatia, T. Jain, and Y. Lim , On the Bures–Wasserstein distance between positive defi- nite matrices, Expositiones Mathematicae, 37 (2019), pp. 165–191
2019
-
[11]
A. M. Bloch, R. W. Brockett, and T. S. Ratiu , Completely integrable gradient flows , Communications in Mathematical Physics, 147 (1992), pp. 57–74
1992
-
[12]
Br ´echet, K
P. Br ´echet, K. Papagiannouli, J. An, and G. Mont ´ufar, Critical points and conver- gence analysis of generative deep linear networks trained with Bures-Wasserstein loss , in International Conference on Machine Learning, PMLR, 2023, pp. 3106–3147
2023
-
[13]
H. T. M. Chu, S. Ghosh, C. T. Lam, and S. S. Mukherjee , Implicit regularization via spectral neural networks and non-linear matrix sensing , arXiv:2402.17595, (2024)
2024 arXiv
-
[14]
Cohen, G
N. Cohen, G. Menon, and Z. Veraszto , Deep linear networks for matrix completion—an infinite depth limit , SIAM Journal on Applied Dynamical Systems, 22 (2023), pp. 3208–3232. THE GEOMETRY OF THE DEEP LINEAR NETWORK 37
2023
-
[15]
DeVore, B
R. DeVore, B. Hanin, and G. Petrova , Neural network approximation , Acta Numerica, 30 (2021), pp. 327–444
2021
-
[16]
C. C. J. Domin ´e, N. Anguita, A. M. Proca, L. Braun, D. Kunin, P. A. M. Mediano, and A. M. Saxe , From lazy to rich: Exact learning dynamics in deep linear networks , arXiv:2409.14623, (2024)
2024 arXiv
-
[17]
S. S. Du, W. Hu, and J. D. Lee , Algorithmic regularization in learning deep homogeneous models: Layers are automatically balanced , Advances in Neural Information Processing Sys- tems, 31 (2018)
2018
-
[18]
R. Ge, C. Jin, and Y. Zheng, No spurious local minima in nonconvex low rank problems: A unified geometric analysis , in International Conference on Machine Learning, PMLR, 2017, pp. 1233–1242
2017
-
[19]
Gunasekar, J
S. Gunasekar, J. D. Lee, D. Soudry, and N. Srebro , Implicit bias of gradient descent on linear convolutional networks , Advances in Neural Information Processing Systems, 31 (2018)
2018
-
[20]
Gunasekar, B
S. Gunasekar, B. E. Woodworth, S. Bhojanapalli, B. Neyshabur, and N. Srebro , Implicit regularization in matrix factorization , in Advances in Neural Information Processing Systems 30: Annual Conference on Neural Information Processing Systems 2017, December 4-9, 2017, Long Be...
2017
-
[21]
E. P. Hsu , Stochastic analysis on manifolds , vol. 38 of Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, 2002
2002
-
[22]
Huang, D
C.-P. Huang, D. Inauen, and G. Menon , Motion by mean curvature and Dyson Brownian motion, Electronic Communications in Probability, 28 (2023), pp. 1–10
2023
-
[23]
Ikeda and S
N. Ikeda and S. W atanabe, Stochastic differential equations and diffusion processes, vol. 24 of North-Holland Mathematical Library, North-Holland Publishing Co., Amsterdam; Kodan- sha, Ltd., Tokyo, second ed., 1989
1989
-
[24]
Inauen and G
D. Inauen and G. Menon , Stochastic Nash evolution , arXiv preprint arXiv:2312.06541, (2023)
2023 arXiv
-
[25]
K. Kohn, T. Merkh, G. Mont ´ufar, and M. Trager , Geometry of linear convolutional networks, SIAM J. Appl. Algebra Geom., 6 (2022), pp. 368–406
2022
-
[26]
K. Kohn, G. Mont´ufar, V. Shahverdi, and M. Trager, Function space and critical points of linear convolutional networks , SIAM J. Appl. Algebra Geom., 8 (2024), pp. 333–362
2024
-
[27]
L ojasiewicz, Sur les trajectoires du gradient d’une fonction analytique , in Geometry seminars, 1982–1983 (Bologna, 1982/1983), Univ
S. L ojasiewicz, Sur les trajectoires du gradient d’une fonction analytique , in Geometry seminars, 1982–1983 (Bologna, 1982/1983), Univ. Stud. Bologna, Bologna, 1984, pp. 115– 117
1982
-
[28]
Lapedes and R
A. Lapedes and R. F arber , How neural nets work , in Evolution, learning and cognition, World Sci. Publ., Teaneck, NJ, 1988, pp. 331–346
1988
-
[29]
J. P. LaSalle , Stability theory and invariance principles , in Dynamical systems (Proc. In- ternat. Sympos., Brown Univ., Providence, R.I., 1974), Vol. I, Academic Press, New York- London, 1976, pp. 211–222
1974
-
[30]
J. M. Lee, Introduction to Riemannian manifolds , vol. 176 of Graduate Texts in Mathemat- ics, Springer, Cham, second ed., 2018
2018
-
[31]
Marcotte, R
S. Marcotte, R. Gribonval, and G. Peyr ´e, Abide by the law and follow the flow: Con- servation laws for gradient flows , Advances in Neural Information Processing Systems, 36 (2024)
2024
-
[32]
, Keep the momentum: Conservation laws beyond euclidean gradient flows , arXiv preprint arXiv:2405.12888, (2024)
2024 arXiv
-
[33]
Menon, Pattern theory: old and new , Lecture notes, Brown University, June 2023
G. Menon, Pattern theory: old and new , Lecture notes, Brown University, June 2023
2023
-
[35]
Milnor, Morse theory, vol
J. Milnor, Morse theory, vol. No. 51 of Annals of Mathematics Studies, Princeton University Press, Princeton, NJ, 1963. Based on lecture notes by M. Spivak and R. Wells
1963
-
[36]
G. M. Nguegnang, H. Rauhut, and U. Terstiege , Convergence of gradient descent for learning linear neural networks , 2021
2021
-
[37]
D. A. Roberts, S. Yaida, and B. Hanin , The principles of deep learning theory , vol. 46, Cambridge University Press, Cambridge, MA, USA, 2022
2022
-
[38]
A. M. Saxe, J. L. McClelland, and S. Ganguli, Exact solutions to the nonlinear dynamics of learning in deep linear neural networks , arXiv preprint arXiv:1312.6120, (2013)
2013 arXiv
-
[39]
11537–11546
, A mathematical theory of semantic development in deep neural networks, Proceedings of the National Academy of Sciences, 116 (2019), pp. 11537–11546. 38 GOVIND MENON
2019
-
[40]
Simon , Theorems on regularity and singularity of energy minimizing maps , Lectures in Mathematics ETH Z¨ urich, Birkh¨ auser Verlag, Basel, 1996
L. Simon , Theorems on regularity and singularity of energy minimizing maps , Lectures in Mathematics ETH Z¨ urich, Birkh¨ auser Verlag, Basel, 1996. Based on lecture notes by Norbert Hungerb¨ uhler. Division of Applied Mathematics, Brown University, 182 George St., Providence...
1996
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.