Solving the trust-region subproblem (TRS) plays a key role in numerical optimization and many other applications. Based on a fundamental result that the solution of TRS of size $n$ is mathematically equivalent to finding the rightmost eigenpair of a certain matrix pair of size $2n$, eigenvalue-based methods are promising due to their simplicity. For $n$ large, the implicitly restarted Arnoldi (IRA) and refined Arnoldi (IRRA) algorithms are well suited for this eigenproblem. For a reasonable comparison of overall efficiency of the algorithms for solving TRS directly and eigenvalue-based algorithms, a vital premise is that the two kinds of algorithms must compute the approximate solutions of TRS with (almost) the same accuracy, but such premise has been ignored in the literature. To this end, we establish close relationships between the two kinds of residual norms, so that, given a stopping tolerance for IRA and IRRA, we are able to determine a reliable one that GLTR should use so as to ensure that GLTR and IRA, IRRA deliver the converged approximate solutions with similar accuracy. We also make a convergence analysis on the residual norms by the Generalized Lanczos Trust-Region (GLTR) algorithm for solving TRS directly, the Arnoldi method and the refined Arnoldi method for the equivalent eigenproblem. A number of numerical experiments are reported to illustrate that IRA and IRRA are competitive with GLTR and IRRA outperforms IRA.
翻译:解决信任区域的子问题(TRS)在数字优化和许多其他应用中发挥着关键作用。基于一个基本结果,即对于直接解决TRS和基于电子价值的算法的总体效率进行合理的比较,一个关键的前提是,两种算法必须用(最接近)相同的精确度来计算TRS的大致解决办法,但在文献中却忽略了这种前提。为此,我们建立了两种剩余规范之间的密切关系,因此,鉴于IRA和IRRA的容忍度已经停止,我们可以确定一种可靠的方法,使GLTR用来直接解决TRS和基于电子价值的算法的总体效率作比较,以确保TRI和IRA的等值的算法能够直接报告TRA的近似方法,而IRRRR和RRA的精确度也能够直接将IRAR的ITR和IRA的精确度方法与IRARA的精确度进行对比。