JPSOL: A Parallel Numerical Algebraic Solver Driven by Application Features

doi:10.19596/j.cnki.1001-246x.8771

Abstract

Abstract:

JPSOL(J Parallel Solver Library for Numerical Algebra Problems) is introduced, including the software architecture, matrix vector data structure, three kinds of algorithm libraries (linear, nonlinear and eigenvalue) and domain specific solvers. Then, the high parallel scalability of JPSOL are demonstrated by the testing results of basic iterative methods. Finally, the effect and robustness of JPSOL are demonstrated by several typical practical applications.

Key words: numerical simulation, parallel computing, solver, Krylov subspace method, preconditioned algorithm

CLC Number:

O246

Shaoliang HU, Xiaowen XU, Hengbin AN, Ran XU, Ronghong FAN. JPSOL: A Parallel Numerical Algebraic Solver Driven by Application Features[J]. Chinese Journal of Computational Physics, 2024, 41(1): 110-121.

Figures/Tables 23

References 42

1	FALGOUT R D, YANG U M. Hypre: A library of high performance preconditioners[C]//International Conference on Computational Science: Computational Science-ICCS 2002. Amsterdam, The Netherlands: Springer, 2002: 632-641.
2	JOLIVET P , ROMAN J E , ZAMPINI S .KSPHPDDM and PCHPDDM: Extending PETSc with advanced Krylov methods and robust multilevel overlapping Schwarz preconditioners[J].Computers & Mathematics with Applications,2021,84,277-295.
3	HERNANDEZ V , ROMAN J E , VIDAL V .SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems[J].ACM Transactions on Mathematical Software,2005,31(3):351-362. DOI
4	AMESTOY P R, DUFF I S, L'EXCELLENT J Y, et al. MUMPS: A general purpose distributed memory sparse solver[C]//International Workshop on Applied Parallel Computing: PARA 2000: Applied Parallel Computing. New Paradigms for HPC in Industry and Academia. Bergen, Norway: Springer, 2000: 121-130.
5	LI Xiaoye .An overview of SuperLU: Algorithms, implementation, and user interface[J].ACM Transactions on Mathematical Software,2005,31(3):302-325. DOI
6	SCHENK O , GÄRTNER K , FICHTNER W , et al.PARDISO: A high-performance serial and parallel sparse linear solver in semiconductor device simulation[J].Future Generation Computer Systems,2001,18(1):69-78. DOI
7	DAVIS T A .Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method[J].ACM Transactions on Mathematical Software,2004,30(2):196-199. DOI
8	HEROUX M A , BARTLETT R A , HOWLE V E , et al.An overview of the Trilinos project[J].ACM Transactions on Mathematical Software,2005,31(3):397-423. DOI
9	RUGE J W, STVBEN K. Algebraic multigrid[M]//MCCORMICK S F. Multigrid Methods. Philadelphia, PA: Society for Industrial and Applied Mathematics, 1987: 73-130.
10	胡少亮, 徐小文, 郑宇腾, 等.系统级封装应用中时谐Maxwell方程大规模计算的求解算法: 现状与挑战[J].计算物理,2021,38(2):131-145.
11	XU Xiaowen , MO Zeyao .Algebraic interface-based coarsening AMG preconditioner for multi-scale sparse matrices with applications to radiation hydrodynamics computation[J].Numerical Linear Algebra with Applications,2017,24(2):e2078. DOI
12	丁琪, 尚月强.非定常Navier-Stokes方程基于两重网格离散的有限元并行算法[J].计算物理,2020,37(1):10-18.
13	李凌霄, 翟传磊, 谢辉, 等.一种求解三维热辐射输运方程的整体预处理迭代方法及并行计算[J].计算物理,2021,38(3):269-279.
14	丁永龙, 胡琳萍, 张瑞勤.一种基于迭代子空间直接求逆算法的高效子空间混合算法[J].计算物理,2021,38(4):418-422.
15	MO Zeyao , ZHANG Aiqing , CAO Xiaolin , et al.JASMIN: A parallel software infrastructure for scientific computing[J].Frontiers of Computer Science in China,2010,4(4):480-488. DOI
16	LIU Qingkai , MO Zeyao , ZHANG Aiqing , et al.JAUMIN: A programming framework for large-scale numerical simulation on unstructured meshes[J].CCF Transactions on High Performance Computing,2019,1(1):35-48. DOI
17	XU Xiaowen , YUE Xiaoqiang , MAO Runzhang , et al.JXPAMG: A parallel algebraic multigrid solver for extreme-scale numerical simulations[J].CCF Transactions on High Performance Computing,2023,5(1):72-83. DOI
18	SAAD Y. Iterative methods for sparse linear systems[M]. 2nd ed. Philadelphia, PA: Society for Industrial and Applied Mathematics, 2003.
19	MORGAN R B .GMRES with deflated restarting[J].SIAM Journal on Scientific Computing,2002,24(1):20-37. DOI
20	PARKS M L , DE STURLER E , MACKEY G , et al.Recycling Krylov subspaces for sequences of linear systems[J].SIAM Journal on Scientific Computing,2006,28(5):1651-1674. DOI
21	CAI Xiaochuan , SARKIS M .A restricted additive Schwarz preconditioner for general sparse linear systems[J].SIAM Journal on Scientific Computing,1999,21(2):792-797. DOI
22	DE STERCK H , YANG U M , HEYS J J , et al.Reducing complexity in parallel algebraic multigrid preconditioners[J].SIAM Journal on Matrix Analysis and Applications,2006,27(4):1019-1039. DOI
23	VANěK P , MANDEL J , BREZINA M .Algebraic multigrid by smoothed aggregation for second and fourth order elliptic problems[J].Computing,1996,56(3):179-196. DOI
24	VAN\VEK P , BREZINA M , MANDEL J .Convergence of algebraic multigrid based on smoothed aggregation[J].Numerische Mathematik,2001,88(3):559-579. DOI
25	WALKER H F , NI Peng .Anderson acceleration for fixed-point iterations[J].SIAM Journal on Numerical Analysis,2011,49(4):1715-1735. DOI
26	TOTH A , KELLEY C T .Convergence analysis for Anderson acceleration[J].SIAM Journal on Numerical Analysis,2015,53(2):805-819. DOI
27	SAAD Y .Numerical methods for large eigenvalue problems[M].Philadelphia, PA: Society for Industrial and Applied Mathematics,2011.
28	STEWART G W .A Krylov-Schur algorithm for large eigenproblems[J].SIAM Journal on Matrix Analysis and Applications,2002,23(3):601-614. DOI
29	KNYAZEV A V .Toward the optimal preconditioned eigensolver: Locally optimal block preconditioned conjugate gradient method[J].SIAM Journal on Scientific Computing,2001,23(2):517-541. DOI
30	KNYAZEV A V , ARGENTATI M E , LASHUK I , et al.Block locally optimal preconditioned eigenvalue Xolvers (BLOPEX) in Hypre and PETSc[J].SIAM Journal on Scientific Computing,2007,29(5):2224-2239. DOI
31	SLEIJPEN G L G , VAN DER VORST H A .A Jacobi-Davidson iteration method for linear eigenvalue problems[J].SIAM Journal on Matrix Analysis and Applications,1996,17(2):401-425. DOI
32	ROMERO E , ROMAN J E .A parallel implementation of Davidson methods for large-scale eigenvalue problems in SLEPc[J].ACM Transactions on Mathematical Software,2014,40(2):1-29.
33	张宁, 李瑜, 谢和虎, 等.一种求解特征值问题的广义共轭梯度算法[J].中国科学(数学),2021,51(8):1297-1320.
34	赵振国, 李光荣, 童杰, 等.封装结构强电磁脉冲多物理效应并行计算程序研制[J].强激光与粒子束,2018,30(8):083001. DOI
35	王卫杰, 赵振国, 胡少亮, 等.芯片-系统电磁脉冲耦合的高性能全波电磁模拟[J].强激光与粒子束,2021,33(12):107-114.
36	王卫杰, 胡少亮, 郑宇腾, 等.并行预处理有限元方法及其在系统级封装结构电磁模拟中的应用[J].电子学报,2021,49(1):58-63.
37	HIPTMAIR R , XU Jinchao .Nodal auxiliary space preconditioning in H(curl) and H(div) spaces[J].SIAM Journal on Numerical Analysis,2007,45(6):2483-2509. DOI
38	BENZI M , GOLUB G H , LIESEN J .Numerical solution of saddle point problems[J].Acta Numerica,2005,14,1-137. DOI
39	TIAN Rong , ZHOU Mozhen , WANG Jingtao , et al.A challenging dam structural analysis: large-scale implicit thermo-mechanical coupled contact simulation on Tianhe-Ⅱ[J].Computational Mechanics,2019,63(1):99-119. DOI
40	AN Hengbin , MO Zeyao , WANG Jingtao , et al.Shear decoupled parallel scalable preconditioners for nonlinear thermo-mechanical coupled contact applications[J].Journal of Scientific Computing,2021,90(1):4.
41	PATANKAR S V , SPALDING D B .A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows[J].International Journal of Heat and Mass Transfer,1972,15(10):1787-1806. DOI
42	胡少亮, 许开龙, 徐然, 等.求解压力Poisson方程的混合粗化代数多重网格算法[J].计算物理,2023,40(5):527-534.

核组数	总核数	运行时间/s	并行效率/%
512	33 280	46.20	100.00
1 024	66 560	46.44	99.50
2 048	133 120	46.68	98.99
4 096	266 240	47.06	98.19
8 162	530 530	47.38	97.53
16 384	1 064 960	48.08	96.09
32 768	2 129 920	63.76	72.46

核组数	总核数	运行时间/s	并行效率/%
512	33 280	46.20	100.00
1 024	66 560	46.44	99.50
2 048	133 120	46.68	98.99
4 096	266 240	47.06	98.19
8 162	530 530	47.38	97.53
16 384	1 064 960	48.08	96.09
32 768	2 129 920	63.76	72.46

核数	AMG_CG	AMG_BiCG	AMG_GMRES
512	16.23 (40)	19.40 (26)	16.60 (40)
1 024	8.18 (40)	9.18 (24)	8.01 (39)
2 048	4.44 (41)	5.06 (25)	4.38 (39)
4 096	2.69 (42)	3.07 (25)	2.69 (40)
8 192	1.68 (41)	2.02 (27)	1.74 (40)
16 384	1.40 (41)	1.51 (26)	1.46 (41)

核数	AMG_CG	AMG_BiCG	AMG_GMRES
512	16.23 (40)	19.40 (26)	16.60 (40)
1 024	8.18 (40)	9.18 (24)	8.01 (39)
2 048	4.44 (41)	5.06 (25)	4.38 (39)
4 096	2.69 (42)	3.07 (25)	2.69 (40)
8 192	1.68 (41)	2.02 (27)	1.74 (40)
16 384	1.40 (41)	1.51 (26)	1.46 (41)

核数	AMG_CG	AMG_BiCG	AMG_GMRES
32	2.03 (23)	2.54 (14)	2.21 (23)
256	3.14 (30)	3.76 (20)	3.32 (30)
2 048	4.73 (41)	5.39 (25)	4.68 (39)
16 384	7.62 (58)	8.04 (34)	7.83 (56)