×

改进的ST算法和数值实验。 (英语) Zbl 1031.65041号

摘要:最近G.H.戈卢布J.-Y.袁[BIT 42,No.4,814-822(2002;Zbl 1017.65020号)]提出了矩阵的ST分解。然而,到目前为止,其数值稳定性尚未讨论。在这里,我们对ST分解的数值行为进行了初步研究。我们还提出了改进(改进算法)以提高算法的数值稳定性。针对一些著名的测试矩阵,对Golub-Yuan算法和我们的改进算法进行了数值测试。
所有测试都包括与LU(或Cholesky)分解的比较,而不进行旋转。这些数值测试表明,Golub Yuan算法及其修改版本具有合理的数值稳定性。特别是,改进算法对于稀疏矩阵是稳定的。此外,在稠密矩阵的情况下,它比Golub-Yuan算法更稳定。

MSC公司:

65平方英尺 线性系统和矩阵反演的直接数值方法
65层50 稀疏矩阵的计算方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Datta,B.N.,《数值线性代数及其应用》(1994年),布鲁克斯和科尔:布鲁克斯和科勒太平洋格罗夫,加利福尼亚州
[2] Demmel,J.W。;Higham,N.J.,具有快速3级BLAS的块算法的稳定性,ACM Trans。数学。软质。,18274-291(1992年)·Zbl 0892.65016号
[3] Demmel,J.W。;新泽西州海姆。;Schreiber,R.S.,块LU分解的稳定性,数值。线性代数应用。,2, 173-190 (1995) ·Zbl 0834.65010号
[4] Fausett,L.V.,《使用Matlab进行应用数值分析》(1999),普伦蒂斯·霍尔:普伦蒂斯霍尔·恩格尔伍德克利夫斯,新泽西州·Zbl 0943.65002号
[5] Golub,G.H。;Van Loan,C.F.,《矩阵计算》(1996),约翰霍普金斯大学出版社:约翰霍普金大学出版社,马萨诸塞州巴尔的摩
[6] Golub,G.H。;Yuan,J.-Y.,ST:对称三角分解及其应用。第一部分:定理和算法,BIT,42,814-822(2002)·Zbl 1017.65020号
[7] 格雷戈里,R.T。;Karney,D.L.,《测试计算算法的矩阵集合》(1969),威利跨科学:威利跨学科纽约·Zbl 0195.44803号
[8] Higham,N.J.,《694算法:MATLAB中的测试矩阵集合》,ACM Trans。数学。软质。,17, 289-305 (1991) ·Zbl 0900.65120号
[9] Higham,N.J.,《MATLAB测试矩阵工具箱》,第237号数值分析报告(1993年),曼彻斯特计算数学中心:曼彻斯特计算机数学中心
[10] Higham,N.J.,《数值算法的准确性和稳定性》(1996年),工业和应用数学学会:工业与应用数学学会,宾夕法尼亚州费城·Zbl 0847.65010号
[11] 凯勒,C。;新墨西哥州古尔德。;Wathen,A.J.,《不定线性系统的约束预处理》(1999年),卢瑟福·阿普尔顿实验室:卢瑟福阿普尔顿实验室Chilton,Didcot
[12] 达夫,I.S。;格里姆斯,R.G。;Lewis,J.G.,《Harwell-Being稀疏矩阵集合用户指南》(第1版),波音计算机服务研究与技术部,报告RAL-92-086(1992),地图集中心,卢瑟福-阿普尔顿实验室:阿特拉斯中心,卢瑟福-阿普尔顿实验室,奇尔顿,迪科特
[13] Todd,J.,Hilbert矩阵有限段的条件,(Taussky,O.,在确定特征值时对线性方程组解的贡献。在确定特征值时对线性方程组解的贡献,应用数学系列,39(1954),国家标准局,美国商务部:美国商务部国家标准局,华盛顿特区),106-116
[14] Varah,J.M.,长矩阵,线性代数应用。,187, 269-278 (1993) ·Zbl 0782.15014号
[15] Wathen,A.J.,Galerkin质量矩阵的实际特征值界限,IMA J.Numer。分析。,7, 449-457 (1987) ·Zbl 0648.65076号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。