跳到主要内容
研究论文

HSL_MI28:一种高效、鲁棒的有限内存不完全Cholesky因子分解码

出版:2014年7月8日出版历史
跳过抽象节

摘要

本文的重点是设计和开发一个新的健壮高效的通用不完全Cholesky因式分解包HSL_MI28,该包可在HSL数学软件库中获得。它实现了一种有限内存方法,利用了正半定Tismenetsky-Kaporin修改方案的思想,通过合并中间内存,是Lin和Moré广泛使用的ICFS算法的推广。不完全因子的密度及其计算中使用的内存量都由用户控制。HSL_MI28的性能是通过大量的数值实验来证明的,这些实验涉及大量的实际应用中产生的测试问题。数值实验用于隔离缩放、排序和删除策略的影响,以评估它们在开发鲁棒代数不完全因子分解预条件中的有用性,并为HSL_MI28选择默认设置。它们还说明了使用适量中间内存的显著优势。此外,结果表明,在内存有限的情况下,可以获得高质量但稀疏的通用预处理器。与ICFS、基于级别的不完全因子分解代码以及最先进的直接求解器进行了比较。

工具书类

  1. Patrick R.Amestoy、Timothy A.Davis和Iain S.Duff。2004.算法837:AMD,一种近似最小度排序算法。ACM事务处理。数学。柔和。30, 3, 381--388.谷歌学者谷歌学者数字图书馆数字图书馆
  2. Owe Axelsson、Igor Kaporin、Igor-Konshin、Andrey Kucherov、Maya Neytcheva、Ben Polman和Alex Yeremin。2000.线性弹性中一组基准问题的代数求解方法的比较。STW项目NNS.4683最终报告。奈梅亨大学数学系。谷歌学者谷歌学者
  3. Owe Axelsson和Niels Munksgaard。1983.固定存储分配的不完全因子分解分析。预处理方法:分析与应用。计算数学专题,第1卷,Gordon&Breach,纽约,219-241。谷歌学者谷歌学者
  4. 劳伦斯科·贝朗·达维加、维塔利·吉里亚、康斯坦丁·利普尼科夫和吉安马科·曼齐尼。2009.多边形网格上Stokes问题的模拟有限差分方法。J.计算。物理。228, 19, 7215--7232.谷歌学者谷歌学者数字图书馆数字图书馆
  5. 米歇尔·本齐。2002.大型线性系统的预处理技术:一项调查。J.计算。物理。182, 2, 418--477.谷歌学者谷歌学者数字图书馆数字图书馆
  6. Michele Benzi、John C.Haws和Miroslav Tůma。2000.预处理高度不确定和非对称矩阵。SIAM J.科学。计算。22, 4, 1333--1353.谷歌学者谷歌学者数字图书馆数字图书馆
  7. 米歇尔·本齐和M.Tůma。2003.正定矩阵的鲁棒不完全因子分解预处理器。数字。线性代数应用。10, 5--6, 385--400.谷歌学者谷歌学者交叉引用交叉引用
  8. 伊丽莎白·H·卡希尔和詹姆斯·麦基。1969.减少稀疏对称矩阵的带宽。ACM第24届全国会议记录。ACM出版社,157--172。谷歌学者谷歌学者数字图书馆数字图书馆
  9. 蒂莫西·戴维斯(Timothy A.Davis)和胡一凡(Yifan Hu),2011年。佛罗里达大学稀疏矩阵收藏。ACM事务处理。数学。柔和。38,第1条。谷歌学者谷歌学者数字图书馆数字图书馆
  10. 伊丽莎白·多兰(Elizabeth D.Dolan)和豪尔赫·莫雷(Jorge J.Moré)。2002.使用性能档案对优化软件进行基准测试。数学。掠夺。91, 2, 201--213.谷歌学者谷歌学者交叉引用交叉引用
  11. 伊恩。S.Duff和Jacko Koster。2001.关于将大型条目排列到稀疏矩阵对角线的算法。SIAM J.矩阵分析。申请。1973年4月22日至996年。谷歌学者谷歌学者数字图书馆数字图书馆
  12. 伊恩·达夫(Iain S.Duff)和圣普莱特(Stéphane Pralet)。2005.稀疏对称不定问题的缩放和旋转策略。SIAM J.矩阵分析。申请。27, 2, 313--340.谷歌学者谷歌学者数字图书馆数字图书馆
  13. 罗兰·弗伦德(Roland W.Freund)和诺埃尔·纳希蒂加尔(Nöel M.Nachtigal)。非厄米矩阵的look-ahead Lanczos算法的实现,第二部分。技术代表,TR 90-46。RIACS,NASA艾姆斯研究中心。谷歌学者谷歌学者
  14. 托马斯·乔治、安舒尔·古普塔和维韦克·萨林。2009年,2011年修订。SPD系统预处理器性能的实证分析。Res.Rep.RC 24737,IBM公司。谷歌学者谷歌学者
  15. 托马斯·乔治(Thomas George)、安舒尔·古普塔(Anshul Gupta)和维维克·萨林(Vivek Sarin),2012年。SPD系统预处理器性能的实证分析。ACM事务处理。数学。柔和。38,第24条。谷歌学者谷歌学者数字图书馆数字图书馆
  16. 尼古拉斯·古尔德(Nicholas I.M.Gould)、胡一凡(Yifan F.Hu)和珍妮弗·斯科特(Jennifer A.Scott)。2007.大型稀疏对称线性方程组的稀疏直接对称解算器的数值评估。ACM事务处理。数学。柔和。33,第10条。谷歌学者谷歌学者数字图书馆数字图书馆
  17. 尼古拉斯·古尔德(Nicholas I.M.Gould)和珍妮弗·斯科特(Jennifer A.Scott)。2004.大型稀疏对称线性方程组直接解的HSL包的数值评估。ACM事务处理。数学。柔和。30, 3, 300--325.谷歌学者谷歌学者数字图书馆数字图书馆
  18. Anshul Gupta和Thomas George,2010年。改进不完全因式分解预处理性能的自适应技术。SIAM J.科学。计算。32, 1, 84--110.谷歌学者谷歌学者数字图书馆数字图书馆
  19. 乔纳森·霍格和珍妮弗·斯科特。2008.缩放对稀疏对称不定解算器性能的影响。技术代表,RAL-TR-2008-007,卢瑟福阿普尔顿实验室。谷歌学者谷歌学者
  20. 乔纳森·霍格(Jonathan D.Hogg)和詹妮弗·斯科特(Jennifer A.Scott)。2011年,HSL_MA97:用于稀疏对称系统的位兼容多波前代码。技术代表,RAL-TR-2011-024,卢瑟福阿普尔顿实验室。谷歌学者谷歌学者
  21. 乔纳森·霍格(Jonathan D.Hogg)和詹妮弗·斯科特(Jennifer A.Scott)。2014.严格稀疏不定系统的旋转策略。ACM事务处理。数学。柔和。40,1,第4条。内政部:http://dx.doi.org/10.1145/2513109.2513113谷歌学者谷歌学者数字图书馆数字图书馆
  22. 高铁。2013年,用于大规模科学计算的Fortran代码集。http://www.hsl.rl.ac.uk。谷歌学者谷歌学者
  23. David Hysom和Alex Pothen。2001.不完全因子预处理的可扩展并行算法。SIAM J.科学。计算。22, 6, 2194--2215.谷歌学者谷歌学者数字图书馆数字图书馆
  24. 艾伦·詹宁斯和乔治·M·马利克。1977年,部分消除。J.Inst.数学。申请。20, 3, 307--316.谷歌学者谷歌学者交叉引用交叉引用
  25. Alan Jennings和George M.Malik。1978。用共轭梯度法求解稀疏线性方程组。国际期刊数字。方法工程12,1,141--158。谷歌学者谷歌学者交叉引用交叉引用
  26. 马克·琼斯和保罗·普拉斯曼。1995年。算法740:计算改进的不完全Cholesky因子分解的Fortran子程序。ACM事务处理。数学。柔和。21, 1, 18--19.谷歌学者谷歌学者数字图书馆数字图书馆
  27. 马克·琼斯和保罗·普拉斯曼。1995年b。一种改进的不完全Cholesky因子分解。ACM事务处理。数学。柔和。21, 1, 5--17.谷歌学者谷歌学者数字图书馆数字图书馆
  28. 伊戈尔·卡波林(Igor E.Kaporin),1998年。基于广义对称正定矩阵的高质量预处理单位T型U加号(&P);单位T型R+R(右)T型单位分解。数字。线性代数应用。5, 6, 483--509.谷歌学者谷歌学者交叉引用交叉引用
  29. 乔治·卡里皮斯和维平·库马尔。1997年10月。METIS:用于划分非结构化图、划分网格和计算稀疏矩阵的填充-约简顺序的软件包(3.0版)。技术代表,明尼苏达大学计算机科学系和陆军HPC研究中心。谷歌学者谷歌学者
  30. 林志仁(Chih-Jen Lin)和豪尔赫·莫雷(Jorge J.Moré)。1999.记忆有限的不完全Cholesky因子分解。SIAM J.科学。计算。21, 1, 24--45.谷歌学者谷歌学者数字图书馆数字图书馆
  31. 康斯坦丁·利普尼科夫、米哈伊尔·沙什科夫、丹尼尔·斯维亚茨基和尤里·瓦西列夫斯基。2007.非结构化三角形和形状规则多边形网格上扩散方程的单调有限体积格式。J.计算。物理。227, 1, 492--512.谷歌学者谷歌学者数字图书馆数字图书馆
  32. 康斯坦丁·利普尼科夫、丹尼尔·斯维亚茨基和尤里·瓦西列夫斯基。2009.多边形网格上扩散方程的无插值单调有限体积法。J.计算。物理。228, 3, 703--716.谷歌学者谷歌学者数字图书馆数字图书馆
  33. 托马斯·曼特乌费尔(Thomas A.Manteuffel)。正定线性系统的不完全因子分解技术。数学。公司。34, 150, 473--497.谷歌学者谷歌学者
  34. 约翰·雷德(John K.Reid)和詹妮弗·斯科特(Jennifer A.Scott)。1999.为小剖面和波前订购对称稀疏矩阵。国际期刊数字。方法工程45,1737-1755。谷歌学者谷歌学者交叉引用交叉引用
  35. 丹尼尔·鲁伊斯。2001.一种平衡矩阵中的行和列范数的缩放算法。技术代表,RAL-TR-2001-034,英国牛津郡奇尔顿卢瑟福阿普尔顿实验室。谷歌学者谷歌学者
  36. 丹尼尔·鲁伊斯和博拉·乌萨尔。2011.矩阵缩放的对称保持算法。法国格勒诺布尔INRIA RR-7552技术代表。谷歌学者谷歌学者
  37. 尤瑟夫·萨阿德。1994.ILUT:双阈值不完整因式分解。数字。线性代数应用。1, 4, 387--402.谷歌学者谷歌学者交叉引用交叉引用
  38. 优素福·萨阿德。1996.稀疏线性系统的迭代方法。PWS出版公司,波士顿。谷歌学者谷歌学者数字图书馆数字图书馆
  39. Jennifer A.Scott和Miroslav Tůma。2011.结构在不完全因子分解预条件中的重要性。位51、2、385--404。谷歌学者谷歌学者交叉引用交叉引用
  40. Jennifer A.Scott和Miroslav Tůma。2014.关于不完全Cholesky因子分解的半正定修正方案。SIAM J.科学。计算。36、2、A634--A667。谷歌学者谷歌学者交叉引用交叉引用
  41. 斯科特·斯隆。1986年。稀疏矩阵的剖面和波前缩减算法。国际期刊数字。方法工程23,239--251。谷歌学者谷歌学者交叉引用交叉引用
  42. 斯科特·斯隆。1989年。用于轮廓和波前缩减的Fortran程序。国际期刊数字。方法工程28,2651-2679。谷歌学者谷歌学者交叉引用交叉引用
  43. 制作了Suarjana和Kincho H.Law。1995。基于值和空间约束的稳健不完全因子分解。国际期刊数字。方法工程38,1703-1719。谷歌学者谷歌学者交叉引用交叉引用
  44. 米伦·蒂斯梅内茨基(Miron Tismenetsky)。1991年。求解大型稀疏线性系统的新预处理技术。线性代数应用。154--156, 331--353.谷歌学者谷歌学者
  45. 山崎一郎、白昭君、陈文斌和理查德·斯卡利塔。2009年。用于多长度尺度对称正定线性系统的高质量预处理技术。数字。数学。理论、方法和应用。2, 4, 469--484.谷歌学者谷歌学者交叉引用交叉引用

索引术语

  1. HSL_MI28:一种高效稳健的有限内存不完全Cholesky分解码

      建议

      评论

      登录选项

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

      登录

      完全访问权限

      • 发布于

        数学软件上的封面图像ACM事务
        ACM数学软件汇刊 第40卷第4期
        2014年6月
        154页
        国际标准编号:0098-3500
        EISSN公司:1557-7295
        内政部:10.1145/2639949
        期刊目录

        版权所有©2014 ACM

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

        出版商

        计算机协会

        美国纽约州纽约市

        出版历史

        • 出版:2014年7月8日
        • 修订过的:2013年10月1日
        • 认可的:2013年10月1日
        • 收到:2013年4月1日
        发布于汤姆斯第40卷第4期

        权限

        请求有关此文章的权限。

        请求权限

        检查更新

        限定符

        • 研究论文
        • 研究
        • 参考

      PDF格式

      以PDF文件查看或下载。

      PDF格式

      电子阅读器

      使用eReader联机查看。

      电子阅读器