多处理机系统中微分方程的分块并行求解.pdfVIP

铢 t 996年3月 西 安 电 子 科 技 大 学 学 报 Mar．1996 第 23卷 第 1期 JOURNAL OF XIDIAN UNIVERSITY V01．23 No．1 2、b 多处理机系统中微分方程的分块并行求解 圭量堡 熊里旦 余晓勇 C=．』 ＼ (西安电子科技大学计 机学院 西安 7l007]) 、 卜＼ 摘要 着重讨论j一类微分方程——齐次Na~4er—Stokes方程(NSE)在多处理 机 系统中的并行求解 ．即分块并行求解．同时讨论 了该算法的网络实现及在实现过程 中进程 问通信的同步与传输控制的设计． 关键词 ；多处理机系统；并行；同步；进程通信 中图分类号 ：TP301．6 设所要求解的微分方程为 J，4-A一0 ∈(0·]· ∈。 r1) l(，。)一 () 0，z∈口 … 其中 ，)可以是向量函数， 为微分算子，z为空问坐标． 方程 (1)是 Navier—Stokes方程的特殊形式 ．即齐次NavierStokes问题．对这一非线睦方程 的求解若采用传统的序列算法 。：是不可行的，但 由于求解这一方程具有很重要的实际意义．所 以开发其并行性将显得尤其重要．文献 [1～3]给出了几种行之有效的方法，文中采用分块技 术．提出求解方程(1)的分块并行算法． 1 该方程的分块并行算法 假设方程(1)中 可以分解成 A — A一十 A2十 … + A 其 中 A．(i 1，2，…，)是新的微分算子． 令 A一 岫0，z)AL+ (￡，f) 2+ … + (，T) ， 一 ∑ (￡．f)．一∑ 。 其中 (·)一 L_-1)c≤(7f+号)r f1 = ， 一10 ≠， 军事电子预研基金资鼬项 目．收稿 日期 ：1995—0106 第 1期 李长虹等：多处理机糸统 中微分方程的分缺并行求睥 93 显然 弱近似算子 ，且保持原问题 (1)的初值．用方程 ()的解来近似方程 (1)的解 +帆 一 0 ∈ 【0]， ∈口 ㈤ lu (z，0)一 “o0 ，0) 一 0， ∈ 0 对于方程(4)可以用以下公式计算 一 ¨一 一一 ：南 + ( 一 )一一0 2A - [- + - -[( A2) + T一 0 A p+ 1 … (5) - -[( ， )一+ ： 一0 = l南+一南+…+：一] 由式 (