Abstract:Analyzing geometric properties of high-dimensional loss functions, such as local curvature and the existence of other optima around a certain point in loss space, can help provide a better understanding of the interplay between neural network structure, implementation attributes, and learning performance. In this work, we combine concepts from high-dimensional probability and differential geometry to study how curvature properties in lower-dimensional loss representations depend on those in the original loss space. We show that saddle points in the original space are rarely correctly identified as such in expected lower-dimensional representations if random projections are used. The principal curvature in the expected lower-dimensional representation is proportional to the mean curvature in the original loss space. Hence, the mean curvature in the original loss space determines if saddle points appear, on average, as either minima, maxima, or almost flat regions. We use the connection between expected curvature in random projections and mean curvature in the original space (i.e., the normalized Hessian trace) to compute Hutchinson-type trace estimates without calculating Hessian-vector products as in the original Hutchinson method. Because random projections are not suitable to correctly identify saddle information, we propose to study projections along dominant Hessian directions that are associated with the largest and smallest principal curvatures. We connect our findings to the ongoing debate on loss landscape flatness and generalizability. Finally, for different common image classifiers and a function approximator, we show and compare random and Hessian projections of loss landscapes with up to about $7\times 10^6$ parameters.

Estimating the Hessian by Back-propagating Curvature

Exact Calculation of the Hessian Matrix for the Multilayer Perceptron

The Hessian by blocks for neural network by backward propagation

Visualizing high-dimensional loss landscapes with Hessian directions

M-FAC: Efficient Matrix-Free Approximations of Second-Order Information

Detecting negative eigenvalues of exact and approximate Hessian matrices in optimization

Optimizing Neural Networks with Kronecker-factored Approximate Curvature

Approximating Hessian matrices using Bayesian inference: a new approach for quasi-Newton methods in stochastic optimization

Compressed Gradient Methods With Hessian-Aided Error Compensation

A matrix algebra approach to approximate Hessians

Estimating the Jacobian matrix of an unknown multivariate function from sample values by means of a neural network

Analytic Insights into Structure and Rank of Neural Network Hessian Maps

Enhance Curvature Information by Structured Stochastic Quasi-Newton Methods

Convergence of Hessian estimator from random samples on a manifold with boundary

On the Local Hessian in Back-Propagation

Revisiting Scalable Hessian Diagonal Approximations for Applications in Reinforcement Learning

On the Curvature Estimates for Hessian Equations

Series of Hessian-Vector Products for Tractable Saddle-Free Newton Optimisation of Neural Networks

First and zeroth-order implementations of the regularized Newton method with lazy approximated Hessians