Journal of Applied Mathematics
- J. Appl. Math.
- Volume 2014 (2014), Article ID 507175, 7 pages.
Riemannian Gradient Algorithm for the Numerical Solution of Linear Matrix Equations
A Riemannian gradient algorithm based on geometric structures of a manifold consisting of all positive definite matrices is proposed to calculate the numerical solution of the linear matrix equation . In this algorithm, the geodesic distance on the curved Riemannian manifold is taken as an objective function and the geodesic curve is treated as the convergence path. Also the optimal variable step sizes corresponding to the minimum value of the objective function are provided in order to improve the convergence speed. Furthermore, the convergence speed of the Riemannian gradient algorithm is compared with that of the traditional conjugate gradient method in two simulation examples. It is found that the convergence speed of the provided algorithm is faster than that of the conjugate gradient method.
J. Appl. Math., Volume 2014 (2014), Article ID 507175, 7 pages.
First available in Project Euclid: 26 March 2014
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
Duan, Xiaomin; Sun, Huafei; Zhao, Xinyu. Riemannian Gradient Algorithm for the Numerical Solution of Linear Matrix Equations. J. Appl. Math. 2014 (2014), Article ID 507175, 7 pages. doi:10.1155/2014/507175. https://projecteuclid.org/euclid.jam/1395855275