Gradient Flow Equations for Deep Linear Neural Networks: A Survey from a Network Perspective
Joel Wendin, Claudio Altafini
TL;DR
The paper analyzes gradient flow dynamics of deep linear networks with quadratic loss by reformulating the training dynamics in terms of a block-shift adjacency matrix $A$, revealing a polynomial, nilpotent, and isospectral ODE with conservation laws. It characterizes a loss landscape with infinitely many global minima and saddle points, none of which are strict local maxima or minima, and introduces a quotient-space viewpoint to classify critical points by invariant loss values. The work further develops simplified dynamics under special initializations (alignment, decoupling, 0-balance, and $\delta$-balanced) and provides detailed analyses of stability, convergence rates, and examples, along with extensions to discrete updates and non-quadratic losses. Collectively, these results offer a global, nonlocal perspective on learning dynamics in deep linear networks, illuminating both the convergence mechanisms and the rich geometry of the loss surface, with potential implications for understanding implicit regularization and initialization effects in broader deep learning contexts.
Abstract
The paper surveys recent progresses in understanding the dynamics and loss landscape of the gradient flow equations associated to deep linear neural networks, i.e., the gradient descent training dynamics (in the limit when the step size goes to 0) of deep neural networks missing the activation functions and subject to quadratic loss functions. When formulated in terms of the adjacency matrix of the neural network, as we do in the paper, these gradient flow equations form a class of converging matrix ODEs which is nilpotent, polynomial, isospectral, and with conservation laws. The loss landscape is described in detail. It is characterized by infinitely many global minima and saddle points, both strict and nonstrict, but lacks local minima and maxima. The loss function itself is a positive semidefinite Lyapunov function for the gradient flow, and its level sets are unbounded invariant sets of critical points, with critical values that correspond to the amount of singular values of the input-output data learnt by the gradient along a certain trajectory. The adjacency matrix representation we use in the paper allows to highlight the existence of a quotient space structure in which each critical value of the loss function is represented only once, while all other critical points with the same critical value belong to the fiber associated to the quotient space. It also allows to easily determine stable and unstable submanifolds at the saddle points, even when the Hessian fails to obtain them.
