High-dimensional Optimization with Low Rank Tensor Sampling and Local Search
Konstantin Sozykin, Andrei Chertkov, Anh-Huy Phan, Ivan Oseledets, Gleb Ryzhakov
TL;DR
TESALOCS tackles $d$-dimensional, non-convex optimization by fusing a TT-format low-rank discrete surrogate for global exploration with a local optimizer for refinement, updating the TT cores via SGD on a log-likelihood loss. The method achieves scalable exploration with memory that grows linearly in $d$, and it demonstrates improvements over a suite of gradient-based and gradient-free baselines under a fixed budget $M$. Through 20 non-convex, 100-dimensional benchmarks, TESALOCS attains orders-of-magnitude better results, illustrating its effectiveness for challenging high-dimensional problems. The work offers a practical, memory-efficient framework for global-to-local optimization with potential impacts on hyperparameter tuning and materials/design applications.
Abstract
We present a novel method called TESALOCS (TEnsor SAmpling and LOCal Search) for multidimensional optimization, combining the strengths of gradient-free discrete methods and gradient-based approaches. The discrete optimization in our method is based on low-rank tensor techniques, which, thanks to their low-parameter representation, enable efficient optimization of high-dimensional problems. For the second part, i.e., local search, any effective gradient-based method can be used, whether existing (such as quasi-Newton methods) or any other developed in the future. Our approach addresses the limitations of gradient-based methods, such as getting stuck in local optima; the limitations of discrete methods, which cannot be directly applied to continuous functions; and limitations of gradient-free methods that require large computational budgets. Note that we are not limited to a single type of low-rank tensor decomposition for discrete optimization, but for illustrative purposes, we consider a specific efficient low-rank tensor train decomposition. For 20 challenging 100-dimensional functions, we demonstrate that our method can significantly outperform results obtained with gradient-based methods like Conjugate Gradient, BFGS, SLSQP, and other methods, improving them by orders of magnitude with the same computing budget.
