一种用于计算矢量有限元方程组的不完全分解预处理方法.pdfVIP

第 卷第 期 重庆邮电大学学报 自然科学版 年 月 一种用于计算矢量有限元方程组的不完全分解预处理方法 田瑾 龚利 史小卫 徐 乐 西安电子科技大学 天线与微波技术国家重点实验室陕西 西安 华东师范大学 科技创新与发展战略研究中心上海 摘 要 提出一种不完全分解预处理方法 并结合迭代法计算矢量有限元方程组 预处理方法采用基于拓展乔里 斯基分解的多波前法对有限元方程组的系数矩阵进行分解和更新并采用基本线性代数系统库函数计算稠密矩阵 乘来保证算法内层循环的高效率 该预处理算法在对系数矩阵进行数值分解前引入缩放矩阵以改善矩阵条件数 针对有限元方程组系数矩阵稀疏或部分稀疏的特性提出一种新的舍弃策略以保证不完全分解的精度和提高预条 件子的构造时间 通过与直接法对比 从时间花费与内存占用两方面 分析了该算法的计算性能 理论和数值实 验表明 提出的预处理方法能大大减少计算时间与分解过程所占用的内存 同时保证了计算的准确性和有效性 关键词矢量有限元 拓展乔里斯基 多波前 不完全分解预处理 中图分类号 文献标识码 文章编号 问题 金建铭 对该方法进行了较为全面 系 引 言 统的介绍和说明包括闭域和开域问题的几种求解 矢量有限元 方法 由于现实中研究的问题多是无限区域的因 是近 年来受到众多学者关注的一种求解边值问 此需要加入边界条件将开域问题转化为闭域问题来 题的电磁计算方法 已被广泛应用于求解各类电磁 求解 通常采用的边界条件有 种 吸收边界条件 收稿日期 修订日期 基金项目 国家自然科学基金 重 庆 邮 电大 学 学 报 自然科学版