| 伍召清,张定宗,郭文峰,马骏,陈坤杰.SCELT程序功能拓展及并行开发[J].南华大学学报(自然科学版),2025,(1):88~96.[WU Zhaoqin,ZHANG Dinzong,GUO Wenfeng,MA Jun,CHEN Kunjie.SCELT Program Function Expansion and Parallel Development[J].Journal of University of South China(Science and Technology),2025,(1):88~96.] |
| SCELT程序功能拓展及并行开发 |
| SCELT Program Function Expansion and Parallel Development |
| 投稿时间：2024-06-20 |
| DOI： |
| 中文关键词: MHD 并行 MPI 特征值 |
| 英文关键词:MHD Parallel MPI Eigenvalue |
| 基金项目:国家自然科学基金项目(11805239;12075282;11775268);湖南省教育厅科学研究项目青年项目(21B0648) |
|
| 摘要点击次数: 0 |
| 全文下载次数: 1 |
| 中文摘要: |
| 本文对托卡马克的符号辅助计算特征值和线性求解程序(symbolic computation aided eigenvalue and linear code for Tokamaks,SCELT)进行了并行开发及功能拓展。具体为优化了非均匀网格加密的功能,提高了程序的收敛速度和收敛精度;增加了程序对稀疏矩阵功能的支持,减少了程序对内存的消耗;使用并行思想和信息传递接口(message passing interface,MPI)技术对程序数值离散模块并行开发,通过多进程并行提高了程序数值离散的效率。同时,基于特征值问题计算的可扩展库(the scalable library for eigenvalue problem computations,SLEPc)特征值求解库为程序添加了大规模稀疏矩阵特征值并行求解功能,提高了程序研究问题的速度以及规模。通过计算内扭曲模(m=1,n=1)和电阻撕裂模(m=2,n=1)的线性增长率和模结构验证了并行及拓展后的程序的可靠性。此外,通过与初始版本的运行结果进行比较,展示了优化后的程序在运行速度、运行规模上的巨大优势,以及程序在研究复杂的磁流体动力学(magnetohydrodynamic,MHD)特征值问题的潜力。 |
| 英文摘要: |
| This paper describes the parallel development and functional expansion of SCELT (symbolic computation aided eigenvalue and linear code for Tokamaks). Specifically, it optimizes the function of non-uniform grid encryption, improving the convergence speed and accuracy of the program; adds support for sparse matrix functions to the program, reducing its memory consumption; uses a parallel mindset and MPI (message passing interface) technology to parallelize the program's numerical discrete module, enhancing the efficiency of the program's numerical discretization through multi-process parallelism. Additionally, based on the SLEPc eigenvalue solver library, it incorporates large-scale sparse matrix eigenvalue parallel solving capabilities into the program, increasing the speed and scale at which research problems can be addressed. The reliability of the parallelized and expanded program is verified by computing the linear growth rates and mode structures of the internal kink mode (m=1, n=1) and the resistive tearing mode (m=2, n=1). Furthermore, by comparing with the results of the original version, the paper demonstrates the significant advantages of the optimized program in terms of running speed and scale, as well as its potential for studying complex MHD (magnetohydrodynamic) eigenvalue problems. |
| 查看全文 查看/发表评论 下载PDF阅读器 |
| 关闭 |
|
|
|