心电仿真中电流偶极子求解的并行算法研究.pdfVIP

心电仿真中电流偶极子求解的并行算法研究 郭永刚 周荷琴 冯焕清 陈立群 (中国科学技术太学 自动化系，合肥 230026) (中国科学技术大学电子科技系，合肥 230026) E—mail：hqzhou@uste．edu．∞ 摘 要 该文舟绍了心电正问题研充的井秆环境．井以心电仿真中的电流偶极子莩解曲倒．讨论在曙光 1000井扦机上 实现井秆心脏模型的方法，比较井秆计算的性能 ．为在国产并秆计，}机上进行心电正问题研究打下了基础。 美■词 井行算法 MP[ 电流偶极子 心 电仿真 文章量号 1002—8331一(2002)03-0244-03 文麓标谀码 A 中围分类号 m 9l TheCaiculation ofCurrentDipolesin ECG Simulation for ParaileiPlaftorm GuoYonggang ZIIollHeqinj ⅡgHuanqing~ Chen Liqunj (Dcpt．ofAutomafi0n，USTC，Hefei230026) (Dept．ofElectronicScienceandTechnology，USTC．Herei23O026) 舢 shcI： ThepaperdescribestheparallelcircumstancesforthestIldyofECG forwardproblem．Bytaking the caleuhttion ofcunentdipolesas the 42X81~1pie，itdiscusses the method to impiementthe paralld heart model On Dawi~1000higIIperformanceparallelmachineandanalyzestheperformanceofparallela1印 thIII．Itbuildsupthe eo面 dence{rIrthesutdyofECG forwardproblem on thenationalproducedparalld computer Keywords：parallelalgorithm，Ma~agePassing Interface，currentdipole，ECG simulation 1 引言 算领域 已经广泛地使用 MPI作为消息传递接口的标准。笔者 心 电正问题研究需要求解大量心肌细胞 、复杂 的传导系 在曙光 1000上的MPl并行环境下，使用标准C语盲实现电流 统、不规则的人体躯干模型构成的大系坑的多个参数间的美 偶极子求解的并行算法 。 系，计算量大．心电仿真精度的提高直接以计算成本的提高为 前提。若在微型机上进行，必须简化算法，从而会降低仿真精 3 心电正问题巾电流偶极子的理论厦意义 度 心舷建筷要对大量的心肌单元进行枢似的计算 ．裉适合在 3．1 心电场方程 并行机上进行， 显著提高计算速度和精度。国产曙光 1000高 心电场是电流场。场中任一点的电流密度， 性能并行计算机．支持C、Fortran90等高级程序设计语言，具有 完善的并行程序开发环境和工具，支持分布式图形库．适用于 J：j js (” 心电正问题研究。该文 心电仿真中电流偶极子的计算为例 ． 式中，^是由电场产生的电流 ，^是由电场中存在的源产 讨论在并行机上实现心电仿真的方法。 生的电流，称为源电流．或非保守场致电流。由于心电坜可视为