標題: Parallel solution of large-scale elgenvalue. problem for master equation in protein folding dynamics
作者: Li, Yiming
Yu, Shao-Ming
Li, Yih-Lang
資訊工程學系
電信工程研究所
Department of Computer Science
Institute of Communications Engineering
關鍵字: master equations;eigenvalue problem;implicitly restarted Arnoldi method;Jacobi-Davidson method;QR method;parallelization;PC-based Linux cluster;MPI
公開日期: 1-May-2008
摘要: It is known that a master equation characterizes time evolution of trajectories and transition of states in protein folding dynamics. Solution of the master equation may require calculating eigenvalues for the corresponding eigenvalue problem. In this paper, we numerically study the folding rate for a dynamic problem of protein folding by solving a large-scale eigenvalue problem. Three methods, the implicitly restarted Arnoldi. Jacobi-Davidson. and QR methods are employed in solving the corresponding large-scale eigenvalue problem for the transition matrix of master equation. Comparison shows that the QR method demands tremendous computing resource when the length of sequence L > 10 due to extremely large size of matrix and CPU time limitation. The Jacobi-Davidson method may encounter convergence issue, for cases of L > 9. The implicitly restarted Arnoldi method is suitable for solving problems among them. Parallelization of the implicitly restarted Arnoldi method is successfully implemented on a PC-based Linux cluster. The parallelization scheme mainly partitions the operation of matrix. For the Arnoldi factorization. we replicate the upper Hessenberg matrix H-m for each processor, and distribute the set of Arnoldi vectors V-m among processors. Each processor performs its own operation. The algorithm is implemented on a PC-based Linux cluster with message passing interface (MPI) libraries. Numerical experiment performing on our 32-nodes PC-based Linux cluster shows that the maximum difference among processors is within 10%. A 23-times speedup and 72% parallel efficiency are achieved when the matrix size is greater than 2 x 106 on the 32-nodes PC-based Linux cluster. This parallel approach enables us to explore large-scale dynamics of protein folding. (C) 2007 Elsevier Inc. All rights reserved.
URI: http://dx.doi.org/10.1016/j.jpdc.2007.09.002
http://hdl.handle.net/11536/9426
ISSN: 0743-7315
DOI: 10.1016/j.jpdc.2007.09.002
期刊: JOURNAL OF PARALLEL AND DISTRIBUTED COMPUTING
Volume: 68
Issue: 5
起始頁: 678
結束頁: 685
Appears in Collections:Articles


Files in This Item:

  1. 000255765100009.pdf

If it is a zip file, please download the file and unzip it, then open index.html in a browser to view the full text content.