Training linear SVMs in linear time
Presents a cutting-plane algorithm based on an equivalent structural optimization formulation that provably trains linear Support Vector Machines in linear time, dramatically outperforming traditional decomposition methods on large, high-dimensional datasets.
Modern data applications in domains such as text categorization, search ranking, and computational drug design frequently process datasets containing hundreds of thousands to millions of examples with tens of thousands of features. Linear Support Vector Machines provide state-of-the-art accuracy for these high-dimensional, sparse data problems. However, traditional training algorithms exhibit super-linear computational scaling relative to the number of examples, making large-scale training slow or computationally intractable. For ordinal regression tasks, existing solvers experience quadratic growth in effective problem size, typically stalling on datasets with more than a few thousand records.
The article develops and evaluates a cutting-plane optimization algorithm, implemented as SVM-Perf, designed to train linear Support Vector Machines in provably linear time for binary classification and near-linear time for ordinal regression. The objective is to eliminate scalability bottlenecks while maintaining exact equivalence to standard soft-margin models.
To evaluate this approach, the author reformulates the standard training problem into an equivalent structural formulation that shares a single slack variable across exponentially many subset constraints. This structure allows a cutting-plane method to greedily identify the most violated constraint in each iteration, guaranteeing convergence within a constant number of steps that does not grow with the training set size. The method was evaluated across five diverse, large-scale benchmarks spanning text classification, physics, and forestry, comprising up to 804,414 examples and 99,757 features. Performance was directly benchmarked against the standard decomposition solver SVM-Light.
The experimental and theoretical results demonstrate substantial performance gains. First, the algorithm achieves provable computational complexity of O(sn) for classification and O(sn log n) for ordinal regression, where s is the average number of non-zero features and n is the sample size. Second, SVM-Perf delivers dramatic speedups over SVM-Light, executing approximately 100 times faster on large text datasets (150 CPU-seconds versus 20,075 CPU-seconds on Reuters CCAT) and roughly 150 times faster on forest cover data. Third, empirical training times grew sub-linearly with the number of examples (roughly proportional to n^0.8), driven by an iteration count that remains bounded or even decreases on larger samples. Fourth, SVM-Perf successfully solved ordinal regression problems on over 800,000 examples in minutes, whereas SVM-Light became computationally intractable beyond approximately 4,000 examples. Finally, these efficiency gains caused no loss in statistical accuracy; predictive performance remained virtually identical between solvers across regularized parameters.
These findings mean that organizations can train linear models on massive datasets at a fraction of the computational time and infrastructure cost previously required. The algorithm also introduces an intuitive stopping tolerance directly tied to training loss, preventing wasteful over-optimization. For operations reliant on search ranking or ranking metrics like the area under the receiver operating characteristic curve, the approach turns previously intractable ranking tasks into routine computations.
Decision-makers and engineering teams managing large-scale text or sparse tabular data should adopt cutting-plane methods like SVM-Perf in place of legacy decomposition solvers for linear classification and ranking. Future development should focus on implementing parallel, out-of-core memory variants of the algorithm to process datasets that exceed system RAM, as well as exploring sampling and matrix approximation techniques to adapt the approach to non-linear kernel models.
These performance advantages apply specifically to linear models trained on sparse, high-dimensional data. For dense data with few features, other specialized solvers remain competitive, and extending the cutting-plane method directly to standard non-linear kernels introduces substantial computational overhead. Within the defined scope of linear models on sparse data, the theoretical proofs and consistent experimental benchmarks provide high confidence in the algorithm's speed and reliability.
- Paper: Pegasos: primal estimated sub-gradient solver for SVM, Shai Shalev-Shwartz et al. (2007). Pegasos introduces a stochastic sub-gradient primal solver that builds upon the foundational cutting-plane objective formulations established in the source paper.
- Paper: LIBLINEAR: A Library for Large Linear Classification, Rong-En Fan et al. (2008). LIBLINEAR extends the linear SVM cutting-plane principles of the source paper into a comprehensive, high-performance open-source library for massive linear classification tasks.
