此函数使用方形稀疏矩阵 A 的填充级别 k 计算不完全 LU 分解的 LU 因子。 Matlab 的 ilu 函数提供基于阈值的丢弃策略的 ILU(0) 和 ILU,但不提供具有一般填充级别的 ILU(k)克。 该代码是作为 C++ MEX 文件实现的,并针对速度进行了优化,同时尝试最大限度地减少内存使用。 请注意,对于 ilu(0) 分解,应该使用 Matlab 的内置 ilu 函数,因为它已经针对这种情况进行了优化,并且速度会更快。 参考: 1. Saad, Y.,稀疏线性系统的迭代方法,SIAM,宾夕法尼亚州费城,第 2 版,2003 年。 2. Dong, X. 和 Cooperman, G.,计算机集群和多核机器上面向任务的并行 ILU(k) 预处理,预印本,arXiv:0803.0048 [cs.DC],2011 年。
2021-09-23 21:44:58
9KB
matlab
1