当解由椭圆边值问题离散而得到的大型线性方程组的时候,迭代法和基于不完全LU分解的预处理子经常结合起来使用。我们研究了一类依赖于松弛参数的不完全LU分解,给出了松弛不完全LU分解(RILU)的存在性条件(Axelsson and Lindskog,1986)。基于Elman(1986,1989)的分析结果,我们也给出了在预处理迭代过程中RILU预处理操作中涉及到的三角因子计算的稳定性条件。在顾桂定的文章中(2001),他的分析是基于一类五对角矩阵。我们的分析包括对一类由七点差分得到的七对角阵的理论分析与数值实验。这类七对角矩阵常常由三维椭圆边值问题用七点差分离散产生。数值实验表明这些条件与松弛参数在得到一个有效的RILU预处理子时起了重要的作用。
展开▼