共辘梯度法求解点源三维地电场 - 北京玛格泰克科技发展有限公司.pdf

第 卷 第 期 地 球 物 理 学 报 , 年 月 了 , 利 用 不 完 全 共 辘梯 度 法 求解 点源 三 维地 电场 吴 小平 徐果 明 李 时灿 中国科 技 大学 地 球 与空 间科 学 系 , 合肥 摘 要 点源 三 维地 电场 的求解 是 一 大 型 数值计算 问题 本 文用 有 限差 分方 法 求解 , 最 后 形成 一 个 线 性 方 程 组 月 , 这 里 月 是 大 型 稀 疏 的 带 状 对 称 矩 阵 解 大 型 稀 疏 方 程 组 的完 全 琦 分解 直接算法 , 一 般要 求 巨 大 的机器 内存来存储 系数矩 阵 , 而 且计算速度极慢 因 此 引入 不 完全 助 共 扼梯 度 算法 及 按行 索 引 的稀 疏存储模 式 , 充分利用 月 的稀 疏性 , 使得计算速度 大大提 高, 而 内存 要 求则大大减少 , 因此 算法 是地 电三 维正 演 的强 有力工 具 词 , , , 正 关键 共扼梯度 法 不 完全 共扼梯度法 按行索引稀疏存储模式 地 电三 维 演 引 言 、 , ’ 点源 三 维地 电场 的数值 计 算 方 法 如有 限差 分 有 限元 法 等 在理 论 上 比较成 熟 , 但 求解 最后 形成 的大 型线性 方程 组 月 则不 易实现 主要 是 以下 两个方 面 的原 因 机器 内存要 求 巨大 对于