计算机科学 ›› 2026, Vol. 53 ›› Issue (6A): 250900126-6.doi: 10.11896/jsjkx.250900126
李琴1, 吴思原2, 杨浩原2, 杜沁2, 凌旭1, 肖国庆3
LI Qin1, WU Siyuan2, YANG Haoyuan2, DU Qin2, LING Xu1, XIAO Guoqing3
摘要: 预处理共轭梯度(PCG)法作为一种大规模稀疏矩阵迭代求解算法,广泛应用于科学工程计算、人工智能等领域。现有研究聚焦于使用深度学习生成预条件算子,以加快求解速度。然而,由于稀疏矩阵的空间复杂性,固定预条件算子不具备通用性,难以适用于所有的稀疏矩阵。为了解决这个问题,提出了一种基于深度学习的预条件算子自适应选择算法及其优化方法。首先,设计了一种卷积神经网络PCNN来捕获稀疏矩阵的空间结构特性;其次,基于此构建了一种结合多层感知机的自适应分类预测模型以选择最优预条件算子。在佛罗里达公开数据集上的实验结果表明,所提出方法的分类准确率优于MLP和SVM等深度学习方法,可达到70.49%;与基于Jacobi,ICCG和SSOR的PCG算法相比,所提出方法的性能分别提高了5.5,4.3和6.2倍。
中图分类号:
| [1] GAIHRE A,LI X S,LIU H.Gsofa:Scalable sparse symbolic lu factorization on gpus[J].IEEE Transactions on Parallel and Distributed Systems,2021,33(4):1015-1026. [2] KARSAVURAN M O,NG E G,PEYTON B W.GPU Accelerated Sparse Cholesky Factorization[C]//SC24-W:Workshops of the International Conference for High Performance Computing,Networking,Storage and Analysis.IEEE,2024:703-707. [3] LIU R Y,XU Z Q,SHANG H H,et al.Study of Iterative Solution Algorithm of Response Density Matrix in Density Functio-nal Perturbation Theory[J].Computer Science,2023,50(6):81-85. [4] HESTENES M R,STIEFEL E.Methods of conjugate gradients for solving linear systems[J].Journal of research of the National Bureau of Standards,1952,49(6):409-436. [5] YE Y,GUO H,WANG B,et al.Coupled Incomplete Cholesky and Jacobi Preconditioned Conjugate Gradient on the New Generation of Sunway Many-Core Architecture[J].IEEE Transactions on Computers,2023,72(11):3326-3339. [6] NAPOV A.An incomplete Cholesky preconditioner based on orthogonal approximations[J].SIAM Journal on Scientific Computing,2023,45(2):A729-A752. [7] HUANG G H,XU Y Z,YI X W,et al.Highly efficient iterative methods for solving linear equations of three-dimensional sphere discontinuous deformation analysis[J].International Journal for Numerical and Analytical Methods in Geomechanics,2020,44(9):1301-1314. [8] MENG X F,HAO X L,MA C H,et al.Research on Machine Learning for Scientific Discovery[J].Chinese Journal of Computers,2023,46(5):877-895. [9] MA H J.Recommendation Algorithm Based on Convolutional Neural Network and Constrained Probability Matrix Factorization[J].Computer Science,2020,47(S1):540-545. [10] LOCKHART S,BIENZ A,GROPP W,et al.Performance analysis and optimal node-aware communication for enlarged conjugate gradient methods[J].ACM Transactions on Parallel Computing,2023,10(1):1-25. [11] TIWARI M,VADHIYAR S.Pipelined Preconditioned Conjugate Gradient Methods for real and complex linear systems for distributed memory architectures[J].Journal of Parallel and Distributed Computing,2022,163:147-155. [12] YUAN F,YANG X,LI S,et al.Optimizing multi-grid preconditioned conjugate gradient method on multi-cores[J].IEEE Transactions on Parallel and Distributed Systems,2024,35(5):768-779. [13] GAO J,ZHOU Y,HE G,et al.A multi-GPU parallel optimization model for the preconditioned conjugate gradient algorithm[J].Parallel Computing,2017,63:1-16. [14] ZHANG T,XUE F.A New Preconditioned Nonlinear Conjugate Gradient Method in Real Arithmetic for Computing the Ground States of Rotational Bose-Einstein Condensate[J].SIAM Journal on Scientific Computing,2024,46(3):A1764-A1792. [15] LI Y,CHEN P Y,DU T,et al.Learning preconditioners for conjugate gradient PDE solvers[C]//International Conference on Machine Learning.PMLR,2023:19425-19439. [16] TRIFONOV V,RUDIKOV A,ILIEV O,et al.Learning from linear algebra:A graph neural network approach to preconditioner design for conjugate gradient solvers[J].arXiv:2405.15557,2024. [17] LIN S,YANG W,WANG H,et al.STM-multifrontal QR:streaming task mapping multifrontal QR factorization empowered by GCN[C]//Proceedings of the International Conference for High Performance Computing,Networking,Storage and Analysis.2021:1-14. [18] WANG H,YANG W,OUYANG R,et al.A heterogeneous parallel computing approach optimizing SpTTM on CPU-GPU via GCN[J].ACM Transactions on Parallel Computing,2023,10(2):1-23. [19] SUN Q,LIU Y,DUN M,et al.SpTFS:Sparse tensor format selection for MTTKRP via deep learning[C]//SC20:International Conference for High Performance Computing,Networking,Stora-ge and Analysis.IEEE,2020:1-14. [20] HESTENES M R,STIEFEL E.Methods of conjugate gradients for solving linear systems[J].Journal of research of the National Bureau of Standards,1952,49(6):409-436. [21] CHEN R S,YANG L,ZHAO L,et al.SSOR preconditioned CG method for the linear interference cancellation of asynchronous CDMA systems[J].International Journal of Electronics,2004,91(12):745-757. [22] ZHENG C,ZHANG J H.SSOR-PCG method used in simulation of geotechnical engineering with finite element method [J].Chinese Journal of Rock Mechanics and Engineering,2007,26(S1):2820-2826. [23] MUKUNOKI D,OZAKI K,OGITA T,et al.Conjugate gradient solvers with high accuracy and bit-wise reproducibility between CPU and GPU using Ozaki scheme[C]//The International Conference on High Performance Computing in Asia-Pacific Region.2021:100-109 [24] VAPNIK V,IZMAILOV R.Reinforced SVM method and memorization mechanisms[J].Pattern Recognition,2021,119:108018. [25] ZHAO Y,LI J,LIAO C,et al.Bridging the gap between deep learning and sparse matrix format selection[C]//Proceedings of the 23rd ACM SIGPLAN Symposium on Principles and Practice of Parallel Programming.2018:94-108. |
|
||