skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Accelerating nuclear configuration interaction calculations through a preconditioned block iterative eigensolver

Journal Article · · Computer Physics Communications
ORCiD logo [1];  [2];  [1];  [1];  [3];  [3]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Computational Research Division
  2. Michigan State Univ., East Lansing, MI (United States). Dept. of Computer Science and Engineering
  3. Iowa State Univ., Ames, IA (United States). Dept. of Physics and Astronomy

In this paper, we describe a number of recently developed techniques for improving the performance of large-scale nuclear configuration interaction calculations on high performance parallel computers. We show the benefit of using a preconditioned block iterative method to replace the Lanczos algorithm that has traditionally been used to perform this type of computation. The rapid convergence of the block iterative method is achieved by a proper choice of starting guesses of the eigenvectors and the construction of an effective preconditioner. These acceleration techniques take advantage of special structure of the nuclear configuration interaction problem which we discuss in detail. The use of a block method also allows us to improve the concurrency of the computation, and take advantage of the memory hierarchy of modern microprocessors to increase the arithmetic intensity of the computation relative to data movement. Finally, we also discuss the implementation details that are critical to achieving high performance on massively parallel multi-core supercomputers, and demonstrate that the new block iterative solver is two to three times faster than the Lanczos based algorithm for problems of moderate sizes on a Cray XC30 system.

Research Organization:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE Office of Science (SC), Nuclear Physics (NP); Michigan State Univ. (United States)
Grant/Contract Number:
AC02-05CH11231; SC0008485; FG02-87ER40371; GE100082; DESC0008485
OSTI ID:
1439235
Alternate ID(s):
OSTI ID: 1550011
Journal Information:
Computer Physics Communications, Vol. 222; ISSN 0010-4655
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 32 works
Citation information provided by
Web of Science

References (19)

Ab initio no core shell model journal March 2013
An iteration method for the solution of the eigenvalue problem of linear differential and integral operators journal October 1950
Improving the scalability of a symmetric iterative eigensolver for multi-core platforms: IMPROVING THE SCALABILITY OF A SYMMETRIC ITERATIVE EIGENSOLVER journal September 2013
Minimization Principles for the Linear Response Eigenvalue Problem II: Computation journal January 2013
Adaptively Compressed Exchange Operator for Large-Scale Hybrid Density Functional Calculations with Applications to the Adsorption of Water on Silicene journal February 2017
An indefinite variant of LOBPCG for definite matrix pencils journal August 2013
Factorization in large-scale many-body calculations journal December 2013
The Shell-Model Code NuShellX@MSU journal June 2014
A High Performance Block Eigensolver for Nuclear Configuration Interaction Calculations journal June 2017
Methods of conjugate gradients for solving linear systems journal December 1952
A Jacobi–Davidson Iteration Method for Linear Eigenvalue Problems journal April 1996
Convergence Theory for Preconditioned Eigenvalue Solvers in a Nutshell journal November 2015
A geometric theory for preconditioned inverse iteration I: Extrema of the Rayleigh quotient journal January 2001
A geometric theory forpreconditioned inverse iterationII: Convergence estimates journal January 2001
Krylov subspace methods for solving large unsymmetric linear systems journal September 1981
GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems journal July 1986
Reducing the profile of sparse symmetric matrices journal December 1976
Basis selection in LOBPCG journal October 2006
Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method journal January 2001

Cited By (7)

Few- and many-nucleon systems with semilocal coordinate-space regularized chiral two- and three-body forces journal February 2019
Ab initio calculations of p -shell nuclei up to N 2 LO in chiral Effective Field Theory journal July 2019
Experimental study of the low-lying negative-parity states in Be 11 using the B 12 ( d , He 3 ) Be 11 reaction journal December 2019
Effective interactions in the s d shell journal November 2019
Few- and many-nucleon systems with semilocal coordinate-space regularized chiral two- and three-body forces text January 2018
Effective interactions in the sd shell text January 2019
Model Order Reduction Algorithm for Estimating the Absorption Spectrum journal September 2017

Similar Records

Accelerating nuclear configuration interaction calculations through a preconditioned block iterative eigensolver
Journal Article · Mon Jan 01 00:00:00 EST 2018 · Computer Physics Communications · OSTI ID:1439235

A comparison of eigensolvers for large-scale 3D modal analysis using AMG-preconditioned iterative methods.
Journal Article · Sat Jan 01 00:00:00 EST 2005 · Proposed for publication in International Journal for Numerical Methods in Engineering. · OSTI ID:1439235

A High Performance Block Eigensolver for Nuclear Configuration Interaction Calculations
Journal Article · Thu Jun 01 00:00:00 EDT 2017 · IEEE Transactions on Parallel and Distributed Systems · OSTI ID:1439235