Large-scale optimization problems arising from the discretization of problems involving PDEs sometimes admit solutions that can be well approximated by low-rank matrices. In this paper, we will exploit this low-rank approximation property by solving the optimization problem directly over the set of low-rank matrices. In particular, we introduce a new multilevel algorithm, where the optimization variable is constrained to the Riemannian manifold of fixed-rank matrices. In contrast to most other multilevel algorithms where the rank is chosen adaptively on each level in order to control the perturbation due to the low-rank truncation, we can keep the ranks (and thus the computational complexity) fixed throughout the iterations. Furthermore, classical implementations of line searches based on Wolfe conditions enable computing a solution where the numerical accuracy is limited to about the square root of the machine epsilon. Here, we propose an extension to Riemannian manifolds of the line search of Hager and Zhang, which uses approximate Wolfe conditions that enable computing a solution on the order of the machine epsilon. Numerical experiments demonstrate the computational efficiency of the proposed framework.
翻译:由PDE问题分解引起的大规模优化问题有时会承认低级矩阵可以非常接近于低级矩阵的解决方案。 在本文中,我们将利用这种低级近似属性,直接解决低级矩阵的优化问题。 特别是,我们引入了一个新的多级算法, 优化变量受Riemannian固定级矩阵的方块限制。 与大多数其他多级算法不同, 后者的级别是适应性地选择的, 以控制低级拖线造成的扰动, 我们可以在整个迭代中保持等级( 因而是计算复杂性 ) 固定。 此外, 典型的基于 Wolfe 条件的线搜索方法可以计算出一种解决方案, 其数字精确性仅限于机器 Epsilon 的平方根。 在此, 我们建议扩展Hager 和 Zhang 的线搜索的里曼式的方块, 后者使用近乎Wolfe 的条件, 从而可以计算出机器 epslon 的顺序上的解决方案。 Nummerical 实验显示拟议框架的计算效率 。