跳到主要内容

算法894:关于基于sep-inverse估计的函数的块Schur-Parlett算法

出版:2009年4月7日出版历史
跳过抽象节

摘要

FORTRAN 95软件用于计算指数积分器所需φ-函数的矩阵值。库中的子程序接受一个完整、对角或上拟三角矩阵作为其参数,该矩阵具有四种精度中的一种精度的实项或复项。实现了两种不同的算法,一种是缩放和平方方法,另一种是改进的块Schur-Parlett算法。在后一种算法中,基于sep-inverse函数的估计,将递归的三乘三分块应用于参数。sep-inverse函数的估计由Hager-Higham估计器执行,该估计器作为LAPACK中的子程序xLACON实现。我们对块Schur-Parlett算法的修改与数值实验结果一起进行了描述。

跳过补充材料部分

补充材料

参考文献

  1. 安德森·E.、白·Z、毕晓夫·C.H.、布莱克福德·S、德梅尔·J.·W、东加拉·J.J.、杜·克罗兹·J.、格林鲍姆·A、哈马林·S.J.,麦肯尼·A和索伦森·D·C·1999。LAPACK用户指南,第三版,工业和应用数学学会,宾夕法尼亚州费城。谷歌学者谷歌学者数字图书馆数字图书馆
  2. Bai,Z.和Demmel,J.W.,1993年。关于交换实Schur形式的对角块。Lin.Alg。申请。186、73-95之间。谷歌学者谷歌学者交叉引用交叉引用
  3. Bartels,R.H.和Stewart,G.W.1972年。算法432:矩阵方程AX+XB&equals的解;C.社区。ACM 15、9、820--826。谷歌学者谷歌学者数字图书馆数字图书馆
  4. Berland,H.、Skaflestad,B.和Wright,W.M.,2007年。EXPINT——指数积分器的MATLAB包。ACM事务处理。数学。柔和。33,1,第4条。谷歌学者谷歌学者数字图书馆数字图书馆
  5. Blackford,L.S.、Demmel,J.、Dongarra,J.,Duff,I.、Hammarling,S.、Henry,G.、Heroux,M.、Kaufman,L.、Lumsdaine,A.、Petitet,A.、Pozo,R.、Remington,K.和Whaley,R.C.,2002年。一组更新的基本线性代数子程序(BLAS)。ACM事务处理。数学。柔和。第28页,第2135-151页。谷歌学者谷歌学者数字图书馆数字图书馆
  6. Davies,P.I.和Higham,N.J.,2003年。计算矩阵函数的Schur-Parlett算法。SIAM J.矩阵分析。申请。25, 2, 464--485.谷歌学者谷歌学者数字图书馆数字图书馆
  7. Frank,W.L.1958年。通过行列式求值以及Danielwski和Wielandt的方法计算复矩阵的特征值。《社会工业杂志》。申请。数学。6, 4, 378--392.谷歌学者谷歌学者交叉引用交叉引用
  8. 新泽西州海姆,1988年。FORTRAN码,用于估计实矩阵或复矩阵的一个范数,并应用于条件估计。ACM事务处理。数学。柔和。14, 4, 381--396.谷歌学者谷歌学者数字图书馆数字图书馆
  9. 新泽西州海姆,2002年。数值算法的准确性和稳定性,第二版,工业和应用数学学会,宾夕法尼亚州费城。谷歌学者谷歌学者数字图书馆数字图书馆
  10. 新泽西州海姆,2005年。重新讨论了矩阵指数的缩放和平方方法。SIAM J.矩阵分析。申请。26, 4, 1179--1193.谷歌学者谷歌学者数字图书馆数字图书馆
  11. 新泽西州海姆和蒂瑟,F.2000。矩阵1-范数估计的块算法及其在1-范数伪谱中的应用。SIAM J.矩阵分析。申请。21, 4, 1185--1201.谷歌学者谷歌学者数字图书馆数字图书馆
  12. 英特尔公司。2008。英特尔Fortran编译器。http://developer.intel.com/software/products/。谷歌学者谷歌学者
  13. Jonsson,I.和Kágström,B.2002年。求解三角形系统的递归分块算法。一、单侧耦合Sylvester型矩阵方程。ACM事务处理。数学。柔和。28, (4), 392--415.谷歌学者谷歌学者数字图书馆数字图书馆
  14. Kágström,B.和Poromaa,P.1992。具有sep的三角Sylvester方程的分布式共享存储块算法&负;1估计值。SIAM J.矩阵分析。申请。13, 1, 90--101.谷歌学者谷歌学者数字图书馆数字图书馆
  15. Kágström,B.和Westin,L.1989年。求解广义Sylvester方程的条件估计的广义Schur方法。IEEE传输。自动。续34、745--751。谷歌学者谷歌学者交叉引用交叉引用
  16. Kassam,A.K.和Trefethen,L.N.2005年。刚性PDE的四阶时间步进。SIAM J.科学。计算。26, 4, 1214--1233.谷歌学者谷歌学者数字图书馆数字图书馆
  17. Koikari,S.,2007年。改进的缩放和平方方法的误差分析。计算。数学。申请。53, 1293--1305.谷歌学者谷歌学者数字图书馆数字图书馆
  18. Minchev,B.和W.M.W.,2005年。一阶半线性问题指数积分器综述。技术报告编号:2/2005。挪威科技大学。挪威特隆赫因。谷歌学者谷歌学者
  19. Moler,C.B.和Van Loan,C.F.2003。二十五年后,十九种计算矩阵指数的可疑方法。SIAM修订版45,1,3--49。谷歌学者谷歌学者交叉引用交叉引用
  20. Parlett,B.N.1976年。三角矩阵函数元素之间的递归。Lin.阿尔及利亚。申请。14, 117--121.谷歌学者谷歌学者交叉引用交叉引用
  21. Sidje,R.B.1998年。EXPOKIT:计算矩阵指数的软件包。ACM事务处理。数学。柔和。24, 1, 130--156.谷歌学者谷歌学者数字图书馆数字图书馆
  22. 太阳微系统公司。2008。Sun Fortran 95编译器。http://developers.sun.com/sunstudio/。谷歌学者谷歌学者
  23. Tommila,M.2005年。apfloat,一个C++高性能任意精度算法包,2.41版。http://www.apfloat.org/。谷歌学者谷歌学者
  24. Torczon,V.1991年。关于多向搜索算法的收敛性。SIAM J.Optimiz公司。1, 1, 123--145.谷歌学者谷歌学者交叉引用交叉引用
  25. 沃特,A.2006。G95手册。网址:http://www.g95.org/。谷歌学者谷歌学者
  26. Ward,R.C.1977年。矩阵指数的数值计算及精度估计。SIAM J.数字。分析。14, 4, 600--610.谷歌学者谷歌学者数字图书馆数字图书馆

索引术语

  1. 算法894:在块上,基于sep逆估计的Γ-函数的Schur-Parlett算法

      建议

      评论

      登录选项

      检查您是否可以通过登录凭据或您的机构访问本文。

      登录

      完全访问权限

      • 发布于

        数学软件上的封面图像ACM事务
        ACM数学软件汇刊 第36卷第2期
        2009年3月
        149页
        国际标准编号:0098-3500
        EISSN公司:1557-7295
        内政部:10.1145/1499096
        期刊目录

        版权所有©2009 ACM

        如果复制品不是为了盈利或商业利益而制作或分发的,并且复制品的第一页载有本通知和完整引文,则允许免费制作本作品的全部或部分数字或硬拷贝以供个人或课堂使用。必须尊重ACM以外的其他人对本作品组成部分的版权。允许用信用证进行摘要。要以其他方式复制或重新发布,在服务器上发布或重新发布到列表,需要事先获得特定许可和/或付费。从请求权限[电子邮件保护]

        出版商

        计算机协会

        美国纽约州纽约市

        出版历史

        • 出版:2009年4月7日
        • 认可的:2008年9月1日
        • 修订过的:2008年3月1日
        • 收到:2007年7月1日
        发布于汤姆斯第36卷第2期

        权限

        请求有关此文章的权限。

        请求权限

        检查更新

        限定符

        • 研究论文
        • 研究
        • 推荐

      PDF格式

      以PDF文件查看或下载。

      PDF格式

      电子阅读器

      使用eReader联机查看。

      电子阅读器