您的位置:山东大学 -> 科技期刊社 -> 《山东大学学报(理学版)》

J4 ›› 2010, Vol. 45 ›› Issue (9): 65-69.

• 第十届中国Rough集与软计算学术会议专栏 • 上一篇    下一篇

解大型线性方程组的轮换重新开始Krylov子空间方法

陆峰   

  1. 江苏广播电视大学公共管理系, 江苏 南京 210036
  • 收稿日期:2010-01-06 出版日期:2010-09-16 发布日期:2010-10-12

Alternately restarted Krylov subspace methods for large linear systems of equations

LU Feng   

  1. Department of Public Management, Jiangsu Radio and TV University, Nanjing 210036, Jiangsu, China
  • Received:2010-01-06 Online:2010-09-16 Published:2010-10-12
  • About author:LU Feng (1975-), Female, Lecture, her research mainly focuses on computational mathematics. Email: lufeng@jstvu.edu.cn
  • Supported by:

    Supported by the Rearch Project on High Professional Education of Jiangsu Province(09GZQ041)

摘要:

重新开始Krylov子空间方法(包括Galerkin法和最小二乘法)是求解大型线性方程组的一类流行和重要的方法。然而,这类方法容易在收敛过程中发生中断或停滞现象。为了解决这一问题,本文提出一种新的重新开始格式,称之为轮换重新开始格式。该格式的基本思想是通过轮流使用方程组系数矩阵与其转置矩阵来生成Krylov子空间。轮换重新开始Krylov方法的迭代残量容易在各个特征向量方向上取得大致相等的收敛量,从而使得收敛得到改善。数值实验结果表明轮换重新开始Krylov子空间方法能够有效解决收敛失败的问题。

关键词: 线性方程组;迭代法;收敛;Krylov子空间方法;重新开始

Abstract:

The restarted Krylov subspace methods, including the Galerkin method and the leastsquares method, are popular and important for solving large linear systems of equations. However, the Galerkin method may suffer from serious breakdown, and the leastsquares method may encounter complete stagnation. To overcome the problems, a new restarting scheme, called the alternately restarting scheme, is proposed in this paper. The underlying idea is to use the Krylov subspaces generated by the coefficient matrix and its transpose alternately. We show that for an alternately restarted Krylov method, its residual tends to get the same reduction in every eigenvector direction, and therefore its convergence can be significantly improved. Numerical experiments are conducted, which indicate that the alternately restarted Krylov subspace methods are efficient and robust.

Key words:  linear systems of equations; iterative methods; convergence; Krylov subspace methods; restarting

No related articles found!
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!