The Polar Express: Optimal Matrix Sign Methods and Their Application to the Muon Algorithm
Not on the map. Why?
Abstract (arXiv)
Computing the polar decomposition and the related matrix sign function has been a well-studied problem in numerical analysis for decades. Recently, it has emerged as an important subroutine within the Muon optimizer for training deep neural networks. However, the requirements of this application differ sharply from classical settings: deep learning demands GPU-friendly algorithms that prioritize high throughput over high precision. We introduce Polar Express, a new method for computing the polar decomposition. Like Newton-Schulz and other classical polynomial methods, our approach uses only matrix-matrix multiplications, making it very efficient on GPUs. Inspired by earlier work of Chen & Chow and Nakatsukasa & Freund, Polar Express adapts the update rule at each iteration by solving a minimax optimization problem. We prove that this strategy minimizes error in a worst-case sense, allowing Polar Express to converge as rapidly as possible both in the early iterations and asymptotically. We also address finite-precision issues, making it practical to use in bfloat16. When integrated into Muon, our method yields consistent improvements in validation loss for a GPT-2 model trained on one to ten billion tokens from the FineWeb dataset, outperforming recent alternatives across a range of learning rates.
Key facts. Cited by 110 papers in the map (data through 14 Oct 2026). Builds on 11 of its 59 references that are in the map. First posted on arXiv 22 May 2025.
Builds on
11 of 59 references.
References · 11
- 2015Deep Residual Learning for Image Recognitioncited by 52,096
- 2024The FineWeb Datasets: Decanting the Web for the Finest Text Data at Scalecited by 1,160
- 2025Muon is Scalable for LLM Trainingcited by 428
- 2018Shampoo: Preconditioned Stochastic Tensor Optimizationcited by 426
- 2024Old Optimizer, New Norm: An Anthologycited by 193
- 2025Training Deep Learning Models with Norm-Constrained LMOscited by 192
- 2023A Spectral Condition for Feature Learningcited by 109
- 2025Practical Efficiency of Muon for Pretrainingcited by 80
- 2023A Distributed Data-Parallel PyTorch Implementation of the Distributed Shampoo Optimizer for Training Neural Networks At-Scalecited by 77
- 2024Modular Duality in Deep Learningcited by 77
- 2025Gluon: Making Muon & Scion Great Again! (Bridging Theory and Practice of LMO-based Optimizers for LLMs)cited by 72
Built on by
Citing papers per year: 2025: 13; 2026: 97.
110 papers.
Sort:
By area:
- 2026Parcae: Scaling Laws For Stable Looped Language Modelscited by 55
- 2025Isotropic Curvature Model for Understanding Deep Learning Optimization: Is Gradient Orthogonalization Optimal?cited by 37
- 2025MuonBP: Faster Muon via Block-Periodic Orthogonalizationcited by 35
- 2026On the Design of Qwen3.8-Next Architecture: Evaluation, Efficiency, and Training Stabilitycited by 35
- 2026Controlled LLM Training on Spectral Spherecited by 31
- 2026Toto 2.0: Time Series Forecasting Enters the Scaling Eracited by 30
- 2025An Exploration of Non-Euclidean Gradient Descent: Muon and its Many Variantscited by 28
- 2025Beyond the Ideal: Analyzing the Inexact Muon Updatecited by 27
- 2026The Newton-Muon Optimizercited by 18
- 2025Generalized Gradient Norm Clipping & Non-Euclidean (L0,L1)-Smoothnesscited by 17
- 2026MuonEq: Balancing Before Orthogonalization with Lightweight Equilibrationcited by 17
- 2026Muon2: Boosting Muon via Adaptive Second-Moment Preconditioningcited by 17
38 more
- 2026Adaptive Matrix Online Learning through Smoothing with Guarantees for Nonsmooth Nonconvex Optimizationcited by 14
- 2026Muon is Not That Special: Random or Inverted Spectra Work Just as Wellcited by 14
- 2026Insights on Muon from Simple Quadraticscited by 13
- 2026Powering Up Zeroth-Order Training via Subspace Gradient Orthogonalizationcited by 12
- 2025Preconditioned Norms: A Unified Framework for Steepest Descent, Quasi-Newton and Adaptive Methodscited by 11
- 2026Enhancing LLM Training via Spectral Clippingcited by 10
- 2026Non-Euclidean Gradient Descent Operates at the Edge of Stabilitycited by 9
- 2025FedMuon: Federated Learning with Bias-corrected LMO-based Optimizationcited by 8
- 2025Optimal Scaling Needs Optimal Normcited by 8
- 2025Optimizing Rank for High-Fidelity Implicit Neural Representationscited by 8
- 2026Manifold constrained steepest descent for smooth and closed-set optimizationcited by 8
- 2026Clarifying Shampoo: Adapting Spectral Descent to Stochasticity and the Parameter Trajectorycited by 8
- 2026Accelerating LLM Pre-Training through Flat-Direction Dynamics Enhancementcited by 8
- 2026Muon with Nesterov Momentum: Heavy-Tailed Noise and (Randomized) Inexact Polar Decompositioncited by 8
- 2026Rethinking Muon Beyond Pretraining: Spectral Failures and High-Pass Remedies for VLA and RLVRcited by 8
- 2026Can Muon Fine-tune Adam-Pretrained Models?cited by 7
- 2026TEON: Tensorized Orthonormalization Beyond Layer-Wise Muon for Large Language Model Pre-Trainingcited by 6
- 2026MUON+: Towards More Effective Muon via One Additional Normalization Step for LLM Pre-trainingcited by 6
- 2026Intrinsic Muon: Spectral Optimization on Riemannian Matrix Manifoldscited by 6
- 2026Pion: A Spectrum-Preserving Optimizer via Orthogonal Equivalence Transformationcited by 6
- 2026SOAP, Muon, and Beyond: Pushing LLM Pretraining Scalescited by 6
- 2026Where Does Warm-Up Come From? Adaptive Scheduling for Norm-Constrained Optimizerscited by 5
- 2026Muon with Finite Newton-Schulz: The Smoothing Benefit in Nonsmooth Nonconvex Optimizationcited by 5
- 2025Ky Fan Norms and Beyond: Dual Norms and Combinations for Matrix Optimizationcited by 4
- 2026PRISM: Distribution-free Adaptive Computation of Matrix Functions for Accelerating Neural Network Trainingcited by 4
- 2026Beyond Muon: MUD (MomentUm Decorrelation) for Faster Transformer Trainingcited by 4
- 2026Softsign: Smooth Sign in Your Optimizer For Better Parameter Heterogeneity Handlingcited by 4
- 2026LoRA-Muon: Spectral Steepest Descent on the Low-Rank Manifoldcited by 4
- 2026Dion3: Full-Stack Orthogonal Updatescited by 4
- 2026Anytime Training with Schedule-Free Spectral Optimizationcited by 3
- 2026CacheMuon: Using Temporal Preconditioning To Approximate Polar Factorcited by 3
- 2026Muown Implicitly Performs Angular Step-size Decaycited by 3
- 2026DMuon: Efficient Distributed Muon Training with Near-Adam Overheadcited by 3
- 2026CMuon: Accelerating and Stabilizing Diffusion Transformer Training via Chunked Momentum Orthogonalizationcited by 3
- 2026SCORE: Replacing Layer Stacking with Contractive Recurrent Depthcited by 2
- 2026PolarAdamW: Disentangling Spectral Control and Schur Gauge-Equivariance in Matrix Optimisationcited by 2
- 2026Gradient Clipping Beyond Vector Norms: A Spectral Approach for Matrix-Valued Parameterscited by 2
- 2026Tight Complexity of Muon under Heavy-Tailed Noisecited by 2
50 most cited of 110 All 110 (JSON)
- 2026From Spectral to Frobenius: When Muon Loses Its Spectral Implicit Bias14 Oct 2026
- 2026Monarch-Muon: Structure-Aware Muon for Fast and Memory-Efficient Pretraining12 Oct 2026
- 2026Spectrally Targeted Muon7 Oct 2026
- 2026MuonIO: Principled Norm-Aware Descent for Embedding Tables and Language Model Heads2 Oct 2026
- 2026SoftServe: A Scalable Quasi-Newton Method for Deep Learning1 Oct 2026
- 2026The Row Normalization Puzzle in Muon30 Sep 2026
- 2026Physical Muon: Orthogonalization as an Equilibrium Computation29 Sep 2026
- 2026LionMuon: Alternating Spectral and Sign Descent for Efficient Training28 Sep 2026
- 2026Convergence of Practical Muon27 Sep 2026
- 2026Cost-free Spectral Estimation for Adaptive Newton--Schulz in Matrix Optimizers27 Sep 2026
- 2026Muon Under Gradient Noise and the Limits of Orthogonalization Near Optima26 Sep 2026
- 2026Fast Differentiable SVD on GPU via Polar Decomposition26 Sep 2026
38 more
- 2026Convergence guarantees for Muon: New parameter regimes and generalizations24 Sep 2026
- 2026Design of the IBM Granite 5.0 TurboCTC ASR Model17 Sep 2026
- 2026Fenchel-Young Duality Gaps: Certified Early Stopping for Regularized Inverse Problems15 Sep 2026
- 2026Nonsmooth Optimization via Orthogonalized Momentum12 Sep 2026
- 2026Beyond the Matrix Sign: Quadratic Spectral Descent7 Sep 2026
- 2026On the Design of Qwen3.8-Next Architecture: Evaluation, Efficiency, and Training Stability31 Aug 2026
- 2026Curvature-Conditioned Multiscale Momentum with Sphere Constraints for LLM Pretraining28 Aug 2026
- 2026Muon with Finite Newton-Schulz: The Smoothing Benefit in Nonsmooth Nonconvex Optimization26 Aug 2026
- 2026A Physical Response-and-Memory Model for Muon Optimization24 Aug 2026
- 2026Scaling Muon for Diffusion Transformers21 Aug 2026
- 2026SingularClip: Preventing Spectral Collapse to Maintain Plasticity in Continual and Reinforcement Learning18 Aug 2026
- 2026Approximate Muon with low-rank adapters14 Aug 2026
- 2026MOON: Multi-Objective OrthoNormalized Updates for Multitask Learning12 Aug 2026
- 2026Dion3: Full-Stack Orthogonal Updates12 Aug 2026
- 2026Muon on the Stiefel Manifold Admits an Exact Closed-Form Update6 Aug 2026
- 2026MALT: Lightweight Curvature-Aware Muon via Diagonal Preconditioning5 Aug 2026
- 2026CMuon: Accelerating and Stabilizing Diffusion Transformer Training via Chunked Momentum Orthogonalization3 Aug 2026
- 2026Sharpness-Aware Minimization and Muon: Robustness under the Spectral Norm28 Jul 2026
- 2026Scale Weight Decay and Train Better26 Jul 2026
- 2026Singular value soft-thresholding via the polar decomposition24 Jul 2026
- 2026PoLoRA: A Preconditioned Orthogonalized LoRA Optimizer20 Jul 2026
- 2026Muse: Representation Geometry of Muon Beyond Normalized Momentum16 Jul 2026
- 2026Reassessing Muon for Matrix Factorization14 Jul 2026
- 2026SOAP, Muon, and Beyond: Pushing LLM Pretraining Scales13 Jul 2026
- 2026SLORR: Simple and Efficient In-Training Low-Rank Regularization9 Jul 2026
- 2026DMuon: Efficient Distributed Muon Training with Near-Adam Overhead25 Jun 2026
- 2026Muown Implicitly Performs Angular Step-size Decay22 Jun 2026
- 2026Towards Understanding the Power and Limits of the Muon Optimizer: A River-Valley Perspective19 Jun 2026
- 2026CacheMuon: Using Temporal Preconditioning To Approximate Polar Factor15 Jun 2026
- 2026Beyond a Single Explanation of the Adam--SGD Gap12 Jun 2026
- 2026Muonp: Muon with Fractional Spectral Powers11 Jun 2026
- 2026LoRA-Muon: Spectral Steepest Descent on the Low-Rank Manifold11 Jun 2026
- 2026PC Layer: Polynomial Weight Preconditioning for Improving LLM Pre-Training4 Jun 2026
- 2026Learned Subspace Compression for Communication-Efficient Pipeline Parallelism3 Jun 2026
- 2026A Note on Stability for Orthogonalized Matrix Momentum with Client Sampling1 Jun 2026
- 2026Exploiting weight-space symmetries for approximating curvature30 May 2026
- 2026How Much Orthogonalization Does Muon Need?29 May 2026
- 2026Softsign: Smooth Sign in Your Optimizer For Better Parameter Heterogeneity Handling29 May 2026
50 newest of 110 All 110 (JSON)
- 2026Laguna M.1/XS.2 Technical Reportcited by 1
- 2026Laguna M.1/XS.2 Technical Report26 May 2026
- 2026Design of the IBM Granite 5.0 TurboCTC ASR Modelcited by 0
- 2026Design of the IBM Granite 5.0 TurboCTC ASR Model17 Sep 2026
- 2026Parcae: Scaling Laws For Stable Looped Language Modelscited by 55
- 2025Isotropic Curvature Model for Understanding Deep Learning Optimization: Is Gradient Orthogonalization Optimal?cited by 37
- 2025MuonBP: Faster Muon via Block-Periodic Orthogonalizationcited by 35
- 2026Toto 2.0: Time Series Forecasting Enters the Scaling Eracited by 30
- 2025An Exploration of Non-Euclidean Gradient Descent: Muon and its Many Variantscited by 28
- 2025Beyond the Ideal: Analyzing the Inexact Muon Updatecited by 27
- 2026The Newton-Muon Optimizercited by 18
- 2025Generalized Gradient Norm Clipping & Non-Euclidean (L0,L1)-Smoothnesscited by 17
- 2026MuonEq: Balancing Before Orthogonalization with Lightweight Equilibrationcited by 17
- 2026Muon2: Boosting Muon via Adaptive Second-Moment Preconditioningcited by 17
- 2026Adaptive Matrix Online Learning through Smoothing with Guarantees for Nonsmooth Nonconvex Optimizationcited by 14
- 2026Muon is Not That Special: Random or Inverted Spectra Work Just as Wellcited by 14
12 most cited of 101 All 101 (JSON)
- 2026From Spectral to Frobenius: When Muon Loses Its Spectral Implicit Bias14 Oct 2026
- 2026Spectrally Targeted Muon7 Oct 2026
- 2026MuonIO: Principled Norm-Aware Descent for Embedding Tables and Language Model Heads2 Oct 2026
- 2026SoftServe: A Scalable Quasi-Newton Method for Deep Learning1 Oct 2026
- 2026The Row Normalization Puzzle in Muon30 Sep 2026
- 2026Physical Muon: Orthogonalization as an Equilibrium Computation29 Sep 2026
- 2026LionMuon: Alternating Spectral and Sign Descent for Efficient Training28 Sep 2026
- 2026Convergence of Practical Muon27 Sep 2026
- 2026Cost-free Spectral Estimation for Adaptive Newton--Schulz in Matrix Optimizers27 Sep 2026
- 2026Muon Under Gradient Noise and the Limits of Orthogonalization Near Optima26 Sep 2026
- 2026Fast Differentiable SVD on GPU via Polar Decomposition26 Sep 2026
- 2026Convergence guarantees for Muon: New parameter regimes and generalizations24 Sep 2026
12 newest of 101 All 101 (JSON)