ILU(k) 预处理器：为一般稀疏矩阵计算填充级别为 k 的 ILU 分解。-matlab开发

此函数使用方形稀疏矩阵 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 年。