Efficient solution of symmetric eigenvalue problems using multigrid preconditioners in the locally optimal block conjugate gradient method

Electronic Transactions on Numerical Analysis

Volumn 15, Issue , 2003, Pages 38-55

  Knyazev, Andrew V ^a   Neymeyr, Klaus ^b  

^a UNIVERSITY OF COLORADO   (United States)

^b UNIVERSITY OF TÜBINGEN   (Germany)

Author keywords

Multigrid preconditioning; Preconditioned conjugate gradient iterative method; Symmetric eigenvalue problems

Indexed keywords

EID: 2942668468     PISSN: 10689613     EISSN: None     Source Type: Journal    
DOI: None     Document Type: Article

References (62)

1. G. ASTRAKHANTSEV, The iterative improvement of eigenvalues, USSR Comput. Math. and Math. Physics, 16,1 (1976), pp. 123-132.
2. Z. BAI, J. DEMMEL, J. DONGARRA, A. RUHE, AND H. VAN DER VORST, eds., Templates for the solution of algebraic eigenvalue problems, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2000.
3. R. E. BANK, Analysis of a multilevel inverse iteration procedure for eigenvalue problems, SIAM J. Numer. Anal., 19 (1982), pp. 886-898.
4. _, PLTMG: a software package for solving elliptic partial differential equations, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 1998. Users' guide 8.0.
5. R. E. BANK AND T. F. CHAN, PLTMGC: a multigrid continuation program for parameterized nonlinear elliptic systems, SIAM J. Sci. Statist. Comput., 7 (1986), pp. 540-559.
6. J. H. BRAMBLE, J. E. PASCIAK, AND A. V. KNYAZEV, A subspace preconditioning algorithm for eigenvector/eigenvalue computation, Adv. Comput. Math., 6 (1996), pp. 159-189.
7. A. BRANDT, S. MCCORMICK, AND J. RUGE, Multigrid methods for differential eigenproblems, SIAM J. Sci. Statist. Comput., 4 (1983), pp. 244-260.
8. V. E. BULGAKOV, M. V. BELYI, AND K. M. MATHISEN, Multilevel aggregation method for solving large-scale generalized eigenvalue problems in structural dynamics, Internat. J. Numer. Methods Engrg., 40 (1997), pp. 453-471.
9. Z. CAI, J. MANDEL, AND S. MCCORMICK, Multigrid methods for nearly singular linear equations and eigenvalue problems, SIAM J. Numer. Anal., 34 (1997), pp. 178-200.
10. S. COSTINER AND S. TA'ASAN, Adaptive multigrid techniques for large-scale eigenvalue problems: solutions of the Schrödinger problem in two and three dimensions, Phys. Rev. E (3), 51 (1995), pp. 3704-3717.
11. _, Simultaneous multigrid techniques for nonlinear eigenvalue problems: solutions of the nonlinear Schrödinger-Poisson eigenvalue problem in two and three dimensions, Phys. Rev. E (3), 52 (1995), pp. 1181-1192.
12. P. DEUFLHARD, T. FRIESE, AND F. SCHMIDT, A nonlinear multigrid eigenproblem solver for the complex Helmholtz equation, Scientific report 97-55, Konrad-Zuse-Zentrum für Informationstechnik Berlin, 1997.
13. P. DEUFLHARD, T. FRIESE, F. SCHMIDT, R. MÄRZ, AND H.-P. NOLTING, Effiziente Eigenmodenberechnung für den Entwurf integriert-optischer Chips. (Efficient eigenmodes computation for designing integrated optical chips)., in Hoffmann, Karl-Heinz (ed.) et al., Mathematik: Schlüsseltechnologie für die Zukunft. Verbundprojekte zwischen Universität und Industrie. Berlin: Springer. 267-279 , 1997.
14. P. DEUFLHARD, P. LEINEN, AND H. YSERENTANT, Concepts of an adaptive hierarchical finite element code, Impact Comput. Sci. Engrg., 1 (1989), pp. 3-35.
15. D. C. DOBSON, An efficient method for band structure calculations in 2D photonic crystals, J. Comput. Phys., 149 (1999), pp. 363-376.
16. D. C. DOBSON, J. GOPALAKRISHNAN, AND J. E. PASCIAK, An efficient method for band structure calculations in 3D photonic crystals, J. Comput. Phys., 161 (2000), pp. 668-679.
17. E. G. D'YAKONOV, Iteration methods in eigenvalue problems, Math. Notes, 34 (1983), pp. 945-953.
18. E. G. D'YAKONOV, Optimization in solving elliptic problems, CRC Press, Boca Raton, FL, 1996. Translated from the 1989 Russian original, Translation edited and with a preface by Steve McCormick.
19. E. G. D'YAKONOV AND A. V. KNYAZEV, Group iterative method for finding lower-order eigenvalues, Moscow University, Ser. XV, Computational Math. and Cybernetics, (1982), pp. 32-40.
20. _, On an iterative method for finding lower eigenvalues, Russian J. of Numerical Analysis and Math. Modelling, 7 (1992), pp. 473-486.
21. E. G. D'YAKONOV AND M. Y. OREKHOV, Minimization of the computational labor in determining the first eigenvalues of differential operators, Mathematical Notes, 27 (1980), pp. 382-391.
22. J. FATTEBERT AND J. BERNHOLC, Towards grid-based O(N) density-functional theory methods: Optimized nonorthogonal orbitals and multigrid acceleration, PHYSICAL REVIEW B (Condensed Matter and Materials Physics), 62 (2000), pp. 1713-1722.
23. S. K. GODUNOV, V. V. OGNEVA, AND G. K. PROKOPOV, On the convergence of the modified steepest descent method in application to eigenvalue problems, in American Mathematical Society Translations. Ser. 2. Vol. 105. (English) Partial differential equations (Proceedings of a symposium dedicated to Academician S. L. Sobolev)., Providence, R. I, 1976, American Mathematical Society, pp. 77-80.
24. W. HACKBUSCH, On the computation of approximate eigenvalues and eigenfunctions of elliptic operators by means of a multi-grid method, SIAM J. Numer. Anal., 16 (1979), pp. 201-215.
25. _, Multigrid eigenvalue computation, in Advances in multigrid methods (Oberwolfach, 1984), Vieweg, Braunschweig, 1985, pp. 24-32.
26. W. HACKBUSCH, Multigrid methods and applications, Springer-Verlag, Berlin, 1985.
27. V. HEUVELINE AND C. BERTSCH, On multigrid methods for the eigenvalue computation of nonselfadjoint elliptic operators, East-West J. Numer. Math., 8 (2000), pp. 275-297.
28. T. HWANG AND I. D. PARSONS, A multigrid method for the generalized symmetric eigenvalue problem: Part I algorithm and implementation, International Journal for Numerical Methods in Engineering, 35 (1992), p. 1663.
29. _, A multigrid method for the generalized symmetric eigenvalue problem: Part II performance evaluation, International Journal for Numerical Methods in Engineering, 35 (1992), p. 1677.
30. _, Multigrid solution procedures for structural dynamics eigenvalue problems, Computational Mechanics, 10 (1992), p. 247.
31. A. V. KNYAZEV, Computation of eigenvalues and eigenvectors for mesh problems: algorithms and error estimates, Dept. Numerical Math. USSR Academy of Sciences, Moscow, 1986. (In Russian).
32. _, Convergence rate estimates for iterative methods for mesh symmetric eigenvalue problem, Soviet J. Numerical Analysis and Math. Modelling, 2 (1987), pp. 371-396.
33. _, A preconditioned conjugate gradient method for eigenvalue problems and its implementation in a subspace, in International Ser. Numerical Mathematics, v. 96, Eigenwertaufgaben in Natur- und Ingenieurwissenschaften und ihre numerische Behandlung, Oberwolfach, 1990., Basel, 1991, Birkhauser, pp. 143-154.
34. _, Preconditioned eigensolvers-an oxymoron?, Electron. Trans. Numer. Anal., 7 (1998), pp. 104-123 (electronic). Large scale eigenvalue problems (Argonne, IL, 1997).
35. _, Preconditioned eigensolvers-an oxymoron?, Electron. Trans. Numer. Anal., 7 (1998), pp. 104-123 (electronic). Large scale eigenvalue problems (Argonne, IL, 1997).
36. _, Preconditioned eigensolvers: practical algorithms, in Templates for the Solution of Algebraic Eigenvalue Problems: A Practical Guide, Z. Bai, J. Demmel, J. Dongarra, A. Ruhe, and H. van der Vorst, eds., SIAM, Philadelphia, 2000, pp. 352-368. Section 11.3. An extended version published as a technical report UCD-CCM 143, 1999, at the Center for Computational Mathematics, University of Colorado at Denver, http://wwwmath. cudenver.edu/ccmreports/rep143.ps.gz.
37. A. V. KNYAZEV, Toward the optimal preconditioned eigensolver: locally optimal block preconditioned conjugate gradient method, SIAM J. Sci. Comput., 23 (2001), pp. 517-541 (electronic). Copper Mountain Conference (2000).
38. A. V. KNYAZEV, Toward the optimal preconditioned eigensolver: locally optimal block preconditioned conjugate gradient method, SIAM J. Sci. Comput., 23 (2001), pp. 517-541 (electronic). Copper Mountain Conference (2000).
39. A. V. KNYAZEV AND K. NEYMEYR, A geometric theory for preconditioned inverse iteration. III: A short and sharp convergence estimate for generalized eigenvalue problems, Linear Algebra Appl., 358 (2003), pp. 95-114.
40. A. V. KNYAZEV AND A. L. SKOROKHODOV, On exact estimates of the convergence rate of the steepest ascent method in the symmetric eigenvalue problem, Linear Algebra and Applications, 154-156 (1991), pp. 245-257.
41. P. LEINEN, W. LEMBACH, AND K. NEYMEYR, An adaptive subspace method for elliptic eigenproblems with hierarchical basis preconditioning, scientific report, Sonderforschungsbereich 382, Universitäten Tübingen und Stuttgart, 1997.
42. S. H. LUI, H. B. KELLER, AND T. W. C. KWOK, Homotopy method for the large, sparse, real nonsymmetric eigenvalue problem, SIAM J. Matrix Anal. Appl., 18 (1997), pp. 312-333.
43. J. MANDEL AND S. MCCORMICK, A multilevel variational method for Au = lBu on composite grids, J. Comput. Phys., 80 (1989), pp. 442-452.
44. J. MARTIKAINEN, T. ROSSI, AND J. TOIVANEN, Computation of a few smallest eigenvalues of elliptic operators using fast elliptic solvers, Communications in Numerical Methods in Engineering, 17 (2001), pp. 521-527.
45. S. F. MCCORMICK, A mesh refinement method for Ax = λBx, Math. Comput., 36 (1981), pp. 485-498.
46. K. NEYMEYR, A geometric theory for preconditioned inverse iteration. I: Extrema of the Rayleigh quotient, Linear Algebra Appl., 322 (2001), pp. 61-85.
47. _, A geometric theory for preconditioned inverse iteration, II: Sharp convergence estimates, Linear Algebra Appl., 322 (2001), pp. 87-104.
48. K. NEYMEYR, A geometric theory for preconditioned inverse iteration applied to a subspace, Math. Comp., 71 (2002), pp. 197-216.
49. _, A posteriori error estimation for elliptic eigenproblems, Numer. Linear Algebra Appl., 9 (2002), pp. 263-279.
50. Y. NOTAY, Combination of Jacobi-Davidson and conjugate gradients for the partial symmetric eigenproblem, Numer. Lin. Alg. Appl., 9 (2002), pp. 21-44.
51. S. OLIVEIRA, On the convergence rate of a preconditioned subspace eigensolver, Computing, 63 (1999), pp. 219-231.
52. E. E. OVTCHINNIKOV AND L. S. XANTHIS, Successive eigenvalue relaxation: a new method for generalized eigenvalue problems and convergence estimates, Proc. R. Soc. Lond. A, 457 (2001), pp. 441-451.
53. B. SAMOKISH, The steepest descent method for an eigenvalue problem with semibounded operators, Izvestiya Vuzov, Math., (1958), pp. 105-114. (In Russian).
54. G. L. G. SLEIJPEN AND H. A. VAN DER VORST, A Jacobi-Davidson iteration method for linear eigenvalue problems, SIAM J. Matrix Anal. Appl., 17 (1996), pp. 401-425.
55. P. SMIT AND M. H. C. PAARDEKOOPER, The effects of inexact solvers in algorithms for symmetric eigenvalue problems, Linear Algebra Appl., 287 (1999), pp. 337-357. Special issue celebrating the 60th birthday of Ludwig Elsner.
56. S. I. SOLOV'EV, Convergence of the modified subspace iteration method for nonlinear eigenvalue problems, Preprint SFB393/99-35, Sonderforschungsbereich 393 an der Technischen Universität Chemnitzs, Technische Universität D-09107 Chemnitz, Germany, 1999. Available at http://www.tu-chemnitz.de/sfb393/Files/PS/sfb99-35.ps.gz.
57. _, Preconditioned gradient iterative methods for nonlinear eigenvalue problems, Preprint SFB393/00-28, Sonderforschungsbereich 393 an der Technischen Universität Chemnitz, Technische Universität D-09107 Chemnitz, Germany, 2000. 17 pp. Available at http://www.tu-chemnitz.de/sfb393/Files/PS/sfb00-28.ps.gz.
58. L. STRAKHOVSKAYA, An iterative method for evaluating the first eigenvalue of an elliptic operator, USSR Comput. Math. and Math. Physics, 17,3 (1977), pp. 88-101.
59. L. STRAKHOVSKAYA AND R. P. FEDORENKO, Solution of the principal spectral problem and mathematical modeling of nuclear reactors, Comput. Math. Math. Phys., 40 (2000), pp. 880-888.
60. _, Solution of the principal spectral problem for a system of multigroup diffusion equations, Comput. Math. Math. Phys., 40 (2000), pp. 1312-1321.
61. H. YSERENTANT, Hierarchical bases, in ICIAM 91, O'Malley and R.E., eds., Washington, DC, 1992, SIAM, Philadelphia, pp. 256-276.
62. T. ZHANG, K. H. LAW, AND G. H. GOLUB, On the homotopy method for perturbed symmetric generalized eigenvalue problems, SIAM J. Sci. Comput., 19 (1998), pp. 1625-1645 (electronic).