Out-of-Core Solution of Eigenproblems for Macromolecular Simulations

We consider the solution of large-scale eigenvalue problems that appear in the motion simulation of complex macromolecules on desktop platforms. To tackle the dimension of the matrices that are involved in these problems, we formulate out-of-core (OOC) variants of the two selected eigensolvers, that basically decouple the performance of the solver from the storage capacity. Furthermore, we contend with the high computational complexity of the solvers by off-loading the arithmetically-intensive parts of the algorithms to a hardware graphics accelerator.

[1]  I. Bahar,et al.  Normal mode analysis of biomolecular structures: functional mechanisms of membrane proteins. , 2010, Chemical reviews.

[2]  Sivan Toledo,et al.  A survey of out-of-core algorithms in numerical linear algebra , 1999, External Memory Algorithms.

[3]  Enrique S. Quintana-Ortí,et al.  Exploring large macromolecular functional motions on clusters of multicore processors , 2013, J. Comput. Phys..

[4]  Robert A. van de Geijn,et al.  A Parallel Eigensolver for Dense Symmetric Matrices Based on Multiple Relatively Robust Representations , 2005, SIAM J. Sci. Comput..

[5]  B. Parlett,et al.  Multiple representations to compute orthogonal eigenvectors of symmetric tridiagonal matrices , 2004 .

[6]  Francisco D. Igual,et al.  A Run-Time System for Programming Out-of-Core Matrix Algorithms-by-Tiles on Multithreaded Architectures FLAME Working Note # 43 , 2010 .

[7]  Edoardo Di Napoli,et al.  Solving dense generalized eigenproblems on multi-threaded architectures , 2011, Appl. Math. Comput..

[8]  Zhaojun Bai,et al.  Optimizing Halley's Iteration for Computing the Matrix Polar Decomposition , 2010, SIAM J. Matrix Anal. Appl..

[9]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[10]  Pablo Chacón,et al.  iMod: multipurpose normal mode analysis in internal coordinates , 2011, Bioinform..

[11]  Christian H. Bischof,et al.  Algorithm 807: The SBR Toolbox—software for successive band reduction , 2000, TOMS.

[12]  Rolf Dach,et al.  Technical Report 2012 , 2013 .

[13]  Lars Skjærven,et al.  Normal mode analysis for proteins , 2009 .

[14]  Gregory A Voth,et al.  Systematic multiscale simulation of membrane protein systems. , 2009, Current opinion in structural biology.

[15]  Gene H. Golub,et al.  Matrix computations , 1983 .

[16]  Nicholas J. Higham,et al.  Stable and Efficient Spectral Divide and Conquer Algorithms for the Symmetric Eigenvalue Decomposition and the SVD , 2013, SIAM J. Sci. Comput..

[17]  Enrique S. Quintana-Ortí,et al.  Applying OOC Techniques in the Reduction to Condensed Form for Very Large Symmetric Eigenproblems on GPUs , 2012, 2012 20th Euromicro International Conference on Parallel, Distributed and Network-based Processing.