| 1 | Jun 1, '23 | Eigenvalues of random matrices with iid entries converge to the Wigner circle law. This is universal (does not depend on the law of the entries), as proved by Tao and Vu in 2010. F | 16.5× | 330k | 1.5k | 239 | 327 |
| 2 | May 24, '23 | Taylor series expansion of sine. | 15.3× | 306k | 2.9k | 493 | 282 |
| 3 | May 27, '24 | Markov chain Monte Carlo methods are used to sample from a Gibbs distribution without knowing the normalizing constant. As the temperature epsilon gets small, samples cluster close | 14.1× | 281k | 1.4k | 207 | 620 |
| 4 | Dec 26, '22 | Dijkstra’s algorithm computes the geodesic distance on a graph in O(n*log(n)) operations using a front propagation algorithm. | 12.8× | 255k | 1.9k | 265 | 265 |
| 5 | Jun 12, '24 | Monte Carlo integration approximates integrals at a rate of 1/sqrt(n), independent of the dimension. | 12.5× | 249k | 1.9k | 340 | 638 |
| 6 | Jan 3, '23 | Shepard interpolation is a surprisingly simple multi-dimensional interpolation method. It uses singular radial basis functions inversely proportional to the distances to the sample | 12.2× | 243k | 1.8k | 247 | 304 |
| 7 | Dec 28, '23 | There are 3 types of non-degenerated quadrics in 3-D: Ellipsoid, Hyperboloid of one sheet, Hyperboloid of two sheets. | 12.1× | 241k | 2.4k | 459 | 468 |
| 8 | Jul 10, '23 | Heat, Wave and Schrodinger equations are fundamental linear partial differential equations of physics. | 11.6× | 232k | 2.2k | 453 | 393 |
| 9 | Jul 3, '24 | Mathematical modeling: part 1/2. | 11.6× | 231k | 3.2k | 419 | 1.3k |
| 10 | Apr 20, '23 | The mean-shift algorithm is a clustering method that progressively evolves points toward local maximum of a kernel density estimator. | 11.3× | 226k | 1.8k | 297 | 395 |
| 11 | Apr 3, '24 | Fokker-Planck equation equivalently describes the movement of a random particle with a drift (as a stochastic ODE) and the evolution of its density (as a PDE). | 11.3× | 225k | 1.8k | 246 | 661 |
| 12 | Jun 11, '23 | A 1 hidden layer perceptron can approximate arbitrary continuous functions using a large enough number of neutrons (Cybenko's theorem). | 11.3× | 225k | 1.9k | 332 | 477 |
| 13 | Jan 12, '24 | The central limit theorem equivalently reads as the convergence of iterated convolutions. | 10.9× | 217k | 1.6k | 237 | 534 |
| 14 | Jan 2, '23 | A joint distribution encodes the dependencies between its two marginal distributions (here in the simple case of a density with respect to a product of measures). | 10.8× | 215k | 1.4k | 185 | 259 |
| 15 | Apr 9, '23 | The harmonic oscillator is the canonical linear second order ODE. Parametrized by a single damping parameter. Solutions are (possibly complex) exponentials. Interpolates between pe | 10.7× | 214k | 1.1k | 228 | 171 |
| 16 | May 18, '23 | The fractional Laplacian is a differential operator generalizing the Laplacian. They are the only rotation/translation/scale invariant linear operator. | 10.7× | 213k | 1.8k | 294 | 293 |
| 17 | Dec 29, '22 | My periodic reminder that my book on the discrete Fourier transforms (incl. FFT, finite fields valued FT, representation theory, and more) is available for free. | 10.4× | 207k | 1.8k | 356 | 649 |
| 18 | Dec 16, '22 | Heat diffusion vs. wave equation on a surface. | 10.4× | 207k | 2.2k | 326 | 210 |
| 19 | Jun 4, '23 | Comparison of the Wasserstein, Hellinger, Kullback-Leibler and reverse KL on the space of Gaussian distributions. | 10.4× | 206k | 1.4k | 252 | 350 |
| 20 | Jun 10, '23 | Fourier is to convolution what Legendre is to inf-convolution. | 10.2× | 203k | 1.4k | 236 | 247 |
| 21 | Dec 29, '22 | Eigenfaces are singular vectors of an image dataset. A principal components analysis is applied to reduce dimensionality and perform face recognition. | 9.9× | 197k | 1.1k | 179 | 188 |
| 22 | Feb 23, '23 | Training a multilayer perceptron with a single hidden layer corresponds to the evolution of a sparse probability measure (sum of Dirac masses) over the neuron's parameter domain (h | 9.1× | 182k | 1.4k | 225 | 336 |
| 23 | Jan 22, '25 | The Mathematics of Artificial Intelligence: In this introductory and highly subjective survey, aimed at a general mathematical audience, I showcase some key theoretical concepts un | 9.1× | 181k | 1.9k | 401 | 1.9k |
| 24 | Dec 22, '22 | Stein’s lemma is a surprisingly simple yet useful characterization of the Gaussian distribution. Enables to derive sharp bounds for convergence rates in the central limit theorem i | 8.9× | 178k | 1.2k | 150 | 183 |
| 25 | Apr 19, '23 | The Helmholtz-Hodge decomposition split vectors fields in two orthogonal subspaces: incompressible and irrotational. In 2-D, corresponds to the sum of a gradient field and a rotate | 8.6× | 171k | 1.3k | 231 | 214 |
| 26 | May 12, '23 | Lorenz attractor is the set of limit trajectories for a simplified mathematical model for atmospheric convection. It has a fractal geometry. | 8.5× | 169k | 1.5k | 281 | 157 |
| 27 | May 28, '23 | Optimization can be carried over on manifolds using gradient descent, where the manifold gradient is the projection of an extrinsic gradient on the tangent plane. | 8.1× | 162k | 1.4k | 231 | 293 |
| 28 | May 29, '24 | Convolutional neural networks are shift-invariant representations obtained by iterating convolutions and pointwise non-linearities. Introduced by Le Cun. | 8.1× | 161k | 1.2k | 132 | 411 |
| 29 | Oct 28, '23 | Inverse of Laplacian matrix can be used to define covariance kernels with positive pairwise correlations between nodes. | 8.0× | 159k | 1.2k | 197 | 379 |
| 30 | Feb 18, '23 | Schrodinger’s problem is an approximation (regularization using diffusion) of the Optimal Transport problem. | 7.9× | 157k | 1.6k | 273 | 265 |
| 31 | Jul 30, '23 | The Shrödinger problem defines a dynamic interpolation between two distributions using Brownian bridges. The 0 temperature limit is Optimal Transport. | 7.8× | 155k | 1.3k | 249 | 310 |
| 32 | Apr 12, '23 | Only harmonic spring and gravity central forces produce periodic motions. | 7.8× | 155k | 1.3k | 223 | 163 |
| 33 | Jan 26, '23 | Sampling a signal is equivalent to periodizing its Fourier transform. For non-smooth signals (with non-compact Fourier transform), it creates aliazing. | 7.5× | 149k | 1.3k | 198 | 252 |
| 34 | Mar 13, '23 | Parabolic PDEs (e.g. heat) smooth out singularities. Hyperbolic PDEs (e.g. wave) displace singularities. | 7.5× | 149k | 1.2k | 220 | 188 |
| 35 | May 9, '23 | Behavior of dynamical system depends on the eigenvalues. In 2D it can be classified according to the trace and determinant. | 7.4× | 148k | 1.3k | 216 | 229 |
| 36 | Jun 18, '23 | The mean curvature motion is the most fundamental curve evolution. It is the equivalent for curves of the heat equation. It is independent of parameterisation (intrinsic) and minim | 7.4× | 147k | 1.3k | 221 | 253 |
| 37 | Feb 13, '23 | Burgers’ equation is one of the simplest nonlinear partial differential equations, with solutions exhibiting shocks. | 7.3× | 146k | 1.4k | 253 | 207 |
| 38 | Jun 7, '23 | A Brownian motion is approximated by a discrete random walk with Gaussian increments. | 7.2× | 144k | 1.6k | 255 | 252 |
| 39 | Mar 19, '23 | Harmonic functions are obtained by solving Laplace equation and define interpolating functions when using Dirichlet boundary conditions. | 7.1× | 143k | 1.3k | 198 | 164 |
| 40 | Jun 26, '23 | The Levenberg-Marquardt is a standard method for non-linear least-squares, combining the best of gradient descent and Newton methods by replacing the Hessian by the Jacobian term o | 7.1× | 141k | 1.1k | 220 | 318 |
| 41 | May 8, '23 | The convex hull is the smallest enclosing convex set. It is also the intersection of all half-spaces containing the points, and the set of all possible convex combinations. | 6.9× | 139k | 1.1k | 180 | 182 |
| 42 | May 5, '23 | Taylor series of exp(it) up to degree 15. | 6.9× | 138k | 1.3k | 208 | 125 |
| 43 | Mar 21, '23 | The simple gravity pendulum is one of the simplest nonlinear second order ODE. Exhibits elliptic-shape orbits in phase space for small enough initial speeds. | 6.9× | 137k | 1.4k | 269 | 173 |
| 44 | Mar 6, '24 | Heat equation vs wave equation inside a planar domain. Boundary conditions matter for PDEs. | 6.6× | 131k | 2.5k | 452 | 603 |
| 45 | Feb 21, '23 | Adding an independent Gaussian variable is the same as doing a heat diffusion on the density. | 6.6× | 131k | 1.1k | 159 | 119 |
| 46 | Apr 27, '23 | The damped harmonic oscillator (linear) vs simple gravity pendulum (non-linear). The two most celebrated second order dynamical systems. | 6.5× | 129k | 1.2k | 234 | 170 |
| 47 | Aug 5, '24 | "Transformers are Universal In-context Learners": in this paper, we show that deep transformers with a fixed embedding dimension are universal approximators for an arbitrarily larg | 6.4× | 128k | 1.5k | 302 | 849 |
| 48 | Apr 24, '23 | Impact of the interpolation method on the level set of the resulting function. | 6.3× | 125k | 1.1k | 152 | 180 |
| 49 | Nov 5, '23 | Generalized means interpolate between the min, the geometric mean, the arithmetic mean and the max. | 6.3× | 125k | 1.1k | 173 | 219 |
| 50 | Aug 18, '23 | Positive semi-definite matrices come in many flavors! | 6.2× | 125k | 1.1k | 235 | 271 |
| 51 | Mar 28, '23 | Gauss-Luca theorem shows that the roots of the derivative of a polynomial belong to the convex hull of the roots of the polynomial. | 6.2× | 123k | 1.2k | 207 | 132 |
| 52 | Jun 20, '23 | Fourier transforms comes with lots of flavors. Sampling and periodization are dual operations which interleave these transforms. | 6.1× | 122k | 1.5k | 284 | 258 |
| 53 | Nov 17, '24 | The discrete Fourier basis is a sampling of the continuous Fourier basis. They are both orthogonal. This is (almost) magic! | 6.1× | 121k | 1.2k | 161 | 469 |
| 54 | Mar 22, '23 | Taylor series approximate within its radius of convergence a function by a polynomial. | 6.0× | 120k | 1.4k | 230 | 153 |
| 55 | Dec 26, '23 | Heavy-ball dynamics (second-order ODEs) induce an advection equation (Liouville equation) of a point cloud (in the mean-field limit of an infinite number of particles: a density) o | 6.0× | 120k | 1.3k | 214 | 319 |
| 56 | Nov 14, '23 | Oldies but goldies: F. Chung, Spectral Graph Theory. Everything you want to know about the Laplacian of a graph! | 6.0× | 119k | 1.1k | 179 | 332 |
| 57 | Aug 1, '23 | I have written a very short mathematical summary of the main ingredients of denoising diffusion models. Most of this is directly borrowed from the amazingly clear slides of @Valent | 5.8× | 116k | 1.3k | 232 | 745 |
| 58 | May 30, '23 | A sub-group of Möbius transforms are those preserving the Poincaré upper-half (equivalently the disk). It corresponds to PSL(2) the 2x2 matrices with determinant 1. | 5.7× | 114k | 1.2k | 221 | 150 |
| 59 | Jun 17, '23 | A Riemannian manifold is locally an Euclidean space. An embedded surface (e.g. the Earth) is an example. It defines a Riemannian tensors field (first fundamental form) in parametri | 5.7× | 113k | 1.3k | 266 | 196 |
| 60 | May 10, '23 | Self-organizing maps of Kohonen defines a map between a 2D space to the data, by minimizing a quantization error and a regularity penalty. | 5.6× | 111k | 1.2k | 211 | 228 |
| 61 | Apr 15, '24 | Solving linear systems comes in different flavors. | 5.5× | 111k | 1.7k | 238 | 778 |
| 62 | May 5, '23 | Slides of my talk of my tutorial talk "the mathematics of neural networks" at the annual meeting of Société Française d'Optique next Thursday at @ENS_ULM | 5.5× | 110k | 1.1k | 199 | 512 |
| 63 | May 31, '24 | The gradient field defines the steepest descent direction. The gradient flow dynamic defines a segmentation of the space into attraction basins of the local minimizers. | 5.5× | 109k | 1.2k | 204 | 352 |
| 64 | Dec 18, '22 | The Delaunay triangulation is the dual of the Voronoi diagram. The most fundamental tool of computational geometry. | 5.4× | 107k | 1.1k | 200 | 176 |
| 65 | Mar 28, '24 | In the static regime, the electric/magnetic fields generated by point sources are the superposition of radial monopoles/dipoles. | 5.2× | 104k | 2.1k | 332 | 599 |
| 66 | Mar 4, '23 | Reaction-diffusions are nonlinear PDEs which describe the formation of a surprisingly rich family of patterns. Introduced by Alan Turing in 1952. | 5.0× | 99k | 1.1k | 202 | 160 |
| 67 | Feb 12, '23 | Fourier approximation of a cat. | 4.8× | 96k | 1.2k | 212 | 93 |
| 68 | Feb 15, '23 | Comparing 2D linear interpolation methods of order 0, 1 and 3. | 4.8× | 95k | 1.1k | 139 | 124 |
| 69 | Apr 21, '24 | The optimal computation of gradients for the composition of functions is an optimal parenthesis problem. Forward and backward (backpropagation) are two extreme cases. Backward is o | 4.7× | 94k | 1.4k | 240 | 667 |
| 70 | Jan 26, '24 | Gradient descent does not always converge. Trajectories can have infinite length. Picture & insights from @GuillaumeG_ | 4.6× | 91k | 1.3k | 202 | 349 |
| 71 | Nov 30, '23 | Many parametric probability distribution families can be written as maximum entropy distributions. | 4.5× | 90k | 1.4k | 255 | 365 |
| 72 | Jan 8, '24 | Spherical harmonics is an orthogonal basis of eigenvectors of the Laplacian on the sphere. Eigenspaces of fixed frequency are rotation invariant. The basic tool to perform linear a | 4.4× | 88k | 1.4k | 228 | 374 |
| 73 | Jul 6, '24 | Lloyd's algorithm is the continuous counterpart of k-means. Optimizes the optimal quantization energy. | 4.4× | 88k | 1.9k | 302 | 661 |
| 74 | Jun 16, '23 | The integral curves of a vector field tangent to a manifold stays on the manifold. | 4.4× | 87k | 1k | 163 | 113 |
| 75 | Jan 10, '24 | The three ways to solve Fokker-Planck equations: PDE evolution, stochastic (independent) or deterministic (coupled) particle systems. | 4.3× | 85k | 1.2k | 236 | 382 |
| 76 | Mar 12, '24 | Conformal maps (differentiable over the complex plane) can be defined infinitesimally as circle packings. | 4.0× | 80k | 1.5k | 255 | 421 |
| 77 | Aug 1, '24 | Large-scale N-body simulations require the evaluation of sums involving long-range interaction kernels. Fast multipoles and Barnes–Hut methods approximate these sums in n*log(n) op | 3.7× | 75k | 1.1k | 195 | 367 |
| 78 | Mar 20, '24 | Heat vs wave equations. | 3.3× | 66k | 1.3k | 208 | 300 |
| 79 | Apr 5, '24 | The Runge-Kutta formulas are the workhorse for the numerical integration of ordinary differential equations. shttps://en.wikipedia.org/wiki/Runge%E2%80%93Kutta_methods | 3.1× | 62k | 1.2k | 189 | 295 |
| 80 | Jan 28, '24 | The Brownian motion (aka Wiener process) is the scaling limit of a random walk. The Brownian bridge forces the path to interpolate between two points. | 2.9× | 59k | 1k | 152 | 269 |
| 81 | Dec 2, '24 | Optimal transport computes an interpolation between two distributions using an optimal coupling. Flow matching, on the other hand, uses a simpler “independent” coupling, which is t | 2.9× | 59k | 1.4k | 197 | 666 |
| 82 | Jul 14, '24 | The Minkowski sum is commutative, associative, and distributive with respect to union. Defines an interpolation between sets. | 2.9× | 57k | 1.1k | 157 | 338 |
| 83 | Mar 24, '20 | I have put online the course notes for my master course on optimization for ML. Covers convex analysis, gradient descent, SGD, Lasso, autodiff, MLP. | 2.8× | 0 | 2.4k | 624 | 643 |
| 84 | Mar 19, '24 | Oldies but goldies: Claude Shannon, A Mathematical Theory of Communication, 1948. Defines, among many things, the notion of entropy as a fundamental limit for compression. | 2.7× | 53k | 1.1k | 211 | 268 |
| 85 | Jul 24, '24 | Non-circular gears are ubiquitous and generate non-constant rotation speeds. The movement of the second gear is obtained from the first one from the no-slip condition, after solvin | 2.4× | 49k | 1.3k | 226 | 322 |
| 86 | Dec 5, '17 | Computational Optimal Transport: the book. All you ever wanted to know about OT: theoretical insights, algorithms and applications! | 2.1× | 0 | 1.9k | 534 | 362 |
| 87 | Sep 24, '24 | Approximation of a function on the sphere using an increasing number of spherical harmonics. | 2.1× | 42k | 1.1k | 163 | 295 |
| 88 | Jul 5, '22 | I wrote a Python tutorial on Conformal Prediction, which is a simple yet powerful method to compute confidence intervals for any prediction method. | 2.1× | 0 | 1.7k | 272 | 664 |
| 89 | Jun 16 | The alpha version of my new book "Optimal Transport for Machine Learners" is out, with in particular an online version with interactive figures | 1.9× | 39k | 490 | 104 | 443 |
| 90 | Oct 30, '22 | Lloyd’s algorithm is the continuous counterpart of k-means. Minimizes the optimal quantization energy, which is an optimal transport distance to free Dirac masses. | 1.8× | 0 | 1.8k | 277 | 399 |
| 91 | Nov 5, '22 | The uncertainty principle states that a function cannot be well concentrated both in space and in frequency. | 1.8× | 0 | 1.8k | 327 | 332 |
| 92 | Aug 5, '21 | Heat equation vs wave equation inside a planar domain. Boundary conditions matter for PDEs. | 1.6× | 0 | 1.9k | 310 | 153 |
| 93 | Aug 8, '20 | Taylor series expansion of sine. | 1.5× | 0 | 1.8k | 318 | 86 |
| 94 | Sep 16, '22 | Linearly interpolating matrices produces intriguing interpolations of the eigenvalues. | 1.3× | 0 | 1.4k | 194 | 225 |
| 95 | Jul 1, '22 | Shepard interpolation generalizes nearest neighbor interpolation. | 1.2× | 0 | 1.3k | 166 | 211 |
| 96 | Jan 5, '20 | Lloyd’s algorithm is the continuous counterpart of k-means. Optimizes the optimal quantization energy, which is an optimal transport distance to free Dirac masses. | 1.1× | 0 | 1.2k | 235 | 159 |
| 97 | Jul 23, '22 | Page rank is the leading eigenvector of a stochastic matrix. Can be computed efficiently for sparse graphs using power iterations. | 1.1× | 0 | 1.1k | 166 | 265 |
| 98 | Sep 14, '21 | Solving linear systems comes with different flavors. | 1.1× | 0 | 1.1k | 175 | 208 |
| 99 | Sep 8, '19 | Heat diffusion vs. wave equation on a surface. | 1.1× | 0 | 1.1k | 296 | 53 |
| 100 | Aug 11, '22 | Spherical harmonics are the equivalent on the sphere of the Fourier basis. Useful to perform low-frequency approximation of functions. | 1.0× | 0 | 1.1k | 161 | 159 |
| 101 | Apr 19, '20 | Gauss-Luca theorem shows that the roots of derivative of a polynomial belongs to the convex hull of the roots of the polynomial. | 1.0× | 0 | 1k | 255 | 74 |