Gradient Method with Retards and Generalizations

From MaRDI portal





Let \({\mathbf A}\) be a real symmetric positive definite \(n \times n\) matrix and \(b\) be a real \(n\)-dimensional vector. The unique solution of the equation \({\mathbf A}x= b\) is the argument of the global minimum of the function \(f(x)\equiv \frac{1}{2} x^{T} {\mathbf A}x - b^{T}x\). Let us denote \(g(x)= {\mathbf A}x - b\) and let \(\lambda (x)\) be the minimizer of \(\varphi(\lambda) \equiv f(x - \lambda g(x))\). The minimum of \(f\) may be calculated using the classical steepest-descent iteration \(x_{k+1}=x_{k}-\lambda (x_{k})g(x_{k})\). The authors propose to use the iteration of the form \(x_{k+1}=x_{k}-\lambda (x_{\nu})g(x_{k})\) where \(\nu = \nu (k)\) is arbitrarily chosen in the set \(\{ k,k-1,\dots, \max \{ 0, k-m\}\}\). A practical motivation for the introduction of this iteration is that it allows to compute \(\lambda (x_{k})\) simultaneously in parallel with \(x_{k+1}\). For \(\nu (k)=k\) we obtain the steepest descent method, for \(\nu (k)= k-1 \) -- the iteration of \textit{J. Barzilai} and \textit{J. M. Borwein} [IMA J. Numer. Anal. 8, No. 1, 141-148 (1988; Zbl 0638.65055)]. The authors prove the global convergence of such iterations and discuss some strategies for the choice of \(\nu (k)\). In the convergence analysis, the authors consider the proposed iteration as a particular case of \(x_{k+1}=x_{k} - \frac {1}{\alpha _{k}} g(x_{k})\) where the scalar \(\alpha _{k}\) is the Rayleigh quotient of the matrix \({\mathbf A}\) at an \(n\)-dimensional vector, whose choice depends on the index \(k\). The proposed iteration was tested on two numerical problems, using 10 different strategies to choose the \(\nu (k)\).




Cited in
(82)








This page was built for publication: Gradient Method with Retards and Generalizations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4243561)