Implicity Restarted Block Arnoldi (IRBA) Method
Implicity Restarted Block Arnoldi (IRBA) Method
批准号:
0311786
负责人:
James Baglama
金额:
$6.5万
依托单位国家:
美国
项目类别:
Standard Grant
财政年份:
2003
资助国家:
美国
项目状态:
已结题
起止时间:
2003-08-15 至 2005-07-31
中文摘要
大型稀疏非对称特征值问题最常用的近似技术是基于Arnoldi过程。1992年,Sorensen开发了Implantable Restarted Arnoldi(伊拉)方法。这导致了名为ARPACK(1996)的FORTRAN软件包。在1997年,我们创建了一个块版本,隐式重启块Lanczos(IRBL)方法求解对称特征值问题。当有多个或非常接近的特征值时,我们的方法优于ARPACK。2002年,我们开发了一个MATLAB代码irbleigs.m,实现了IRBL方法。分组码具有许多优点,特别是,它允许使用3级BLAS矩阵-矩阵乘积来提高性能,并且可以定位多个或集群特征值,而不需要任何紧缩技术。虽然,Lehoucq和Maschhoff在1997年创建了一个块Arnoldi方法,这是一个简单的推广(伊拉)的方法,有涉及的实施问题,需要进一步调查之前,一个可行的计算机代码可以创建。我们建议发展隐式重新开始块Arnoldi(IRBA)方法来解决非对称特征值问题,并创建一个公共领域的计算机代码。我们的方法将使用Householder变换而不是Gram-Schimdt过程来保持强正交性,并且将基于等效的Wu和Simon(2001)方法。解决特征值问题是我们生活中常见的一部分,无论我们是否意识到这一点。我将开发的软件包将非常有助于解决许多应用,包括:飞机降噪,建筑物和桥梁的振动分析,量子力学,图像去模糊和人脸识别。
英文摘要
The most commonly used approximation techniques for large sparse nonsymmetric eigenvalue problems are based on the Arnoldi process. In 1992, Sorensen developed the Implicitly Restarted Arnoldi (IRA) method. This resulted in the FORTRAN software package called ARPACK (1996). In 1997, we created a block-version, the Implicitly Restarted Block Lanczos (IRBL) method for solving symmetric eigenvalue problems. When there are multiple or very close eigenvalues our method outperforms ARPACK. In 2002, we developed a MATLAB code irbleigs.m that implements the IRBL method. A block code has many advantages, in particular, it allows the use of Level 3 BLAS matrix-matrix products to increase performance and can locate multiple or clustered eigenvalues without the need of any deflation techniques. Although, Lehoucq and Maschhoff created a block Arnoldi method in 1997 which is a straightforward generalization of the (IRA) method there are implementation issues involved that warrant further investigation before a viable computer code can be created. We propose to develop the implicitly restarted block Arnoldi (IRBA) method to solve the nonsymmetric eigenvalue problem and create a public domain computer code. Our method will use Householder transformations instead of the Gram-Schimdt process to maintain strong orthogonality and will be based on the equivalent Wu and Simon's (2001) approach. Solving eigenvalue problems are a common part of our life whether we realize it or not. The software package that I will develop will be very helpful in solving numerous applications, including; noise reduction in airplanes, vibrational analysis in buildings and bridges, quantum mechanics, image deblurring, and face recognition.
期刊论文(0)
专著(0)
科研奖励(0)
会议论文
海外基金