Abstract
In this paper, we describe tensor methods for large systems of nonlinear equations based on Krylov subspace techniques for approximately solving the linear systems that are required in each tensor iteration. We refer to a method in this class as a tensor-Krylov algorithm. We describe comparative testing for a tensor-Krylov implementation versus an analogous implementation based on a Newton-Krylov method. The test results show that tensor-Krylov methods are much more efficient and robust than Newton-Krylov methods on hard nonlinear equations problems.
Similar content being viewed by others
References
A. Bouaricha, “Solving large sparse systems of nonlinear equations and nonlinear least squares problems using tensor methods on sequential and parallel computers,” Ph.D. Thesis, Cmputer Science Department, University of Colorado at Boulder, 1992.
A. Bouaricha and R.B. Schnabel, “TENSOLVE: A software package for solving systems of nonlinear equations and nonlinear least squares problems using tensor methods,” Preprint MCS-P463-0894, Mathematics and Computer Science Division, Argonne National Laboratory, 1994.
A. Bouaricha and R.B. Schnabel, “Tensor methods for large, sparse systems of nonlinear equations,” Preprint MCS-P473-1094, Mathematics and Computer Science Division, Argonne National Laboratory, 1994.
P.N. Brown and Y. Saad, “Globally convergent techniques in nonlinear Newton-Krylov algorithms,” Technical Report, L.-316, Computing and Mathematics Research Division, Lawrence Livermore National Laboratory, 1989.
P.N.Brown and Y.Saad, “Hybrid Krylov methods for nonlinear systems of equations,” SIAM J. Stat. Comp., vol. 11, pp. 450–481, 1990.
J.E.Dennis and R.B.Schnabel, Numerical Methods for Unconstrained Optimization and Nonlinear Equations, Prentice-Hall: Englewood Cliffs, N.J., 1983.
D.Feng, P.Frank, and R.B.Schnabel, “Local convergence analysis of tensor methods for nonlinear equations,” Math. Prog., vol. 62, pp. 427–459, 1993.
Y.Saad, “Krylov subspace methods for solving unsymmetric linear systems,” Math. Comp., vol. 37, pp. 105–126, 1981.
Y.Saad and M.H.Schultz, “Conjugate gradient-like algorithms for solving nonsymmetric linear systems,” Math. Comp., vol. 44, pp. 417–424, 1985.
Y.Saad and M.H.Schultz, “GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems,” SIAM J. Sci. Stat. Comput., vol. 7, pp. 856–869, 1986.
R.B.Schnabel and P.D.Frank, “Tensor methods for nonlinear equations,” SIAM J. Numer. Anal., vol. 21, pp. 815–843, 1984.
R.B.Schnabel, J.E.Koontz and B.E.Weiss, “A modular system of Algorithms of unconstrained minimization,” ACM Trans. Math. Softw., vol. 11, pp. 419–440, 1985.
J.N.Shadid and R.S.Tuminaro, “Sparse iterative algorithm software for large-scale MIMD machines: An initial discussion and implementation,” Technical Report Sand91–0059, Sandia National Laboratories, Albuquerque, N.M., 1991.
V. Simonsini and E. Gallopoulos, “Convergence properties of block GMRES for solving systems with multiple right-sides,” Technical Report 1316, Center for Supercomputing Research and Development, University of Illinois at Urbana-Champaign, 1992.
V.Simonsini and E.Gallopoulos, “An iterative method for nonsymmetric systems with multiple right handhand sides,” SIAM J. Sci. Stat. Comput., vol. 16, pp. 917–933, 1995.
P.Sonneveld, “CGS: A fast Lanczos-type solver for nonsymmetric linear system,” SIAMJ. Sci. Stat. Comput., vol. 10, pp. 36–52, 1989.
H.Van DerVost, “Bi-CGSTAB: A fast and smoothly converging variant of Bi-CG for the solution of nonsymmetric linear systems,” SIAM J. Sci. Stat. Comput., vol. 13, pp. 631–644, 1992.
Author information
Authors and Affiliations
Additional information
Part of this work was performed while the author was research associate at CERFACS (Centre Européen de Recherche et de Formation Avancée en Calcul Scientifique).
Research supported in part by the Office of Scientific Computing, U.S. Department of Energy, under Contract W-31-109-Eng-38.
Rights and permissions
About this article
Cite this article
Bouaricha, A. Tensor-Krylov methods for large nonlinear equations. Comput Optim Applic 5, 207–232 (1996). https://doi.org/10.1007/BF00248265
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF00248265