Repository logo

Infoscience

  • English
  • French
Log In
Logo EPFL, École polytechnique fédérale de Lausanne

Infoscience

  • English
  • French
Log In
  1. Home
  2. Academic and Research Output
  3. Journal articles
  4. A mixed precision LOBPCG algorithm
 
research article

A mixed precision LOBPCG algorithm

Kressner, Daniel  
•
Ma, Yuxin
•
Shao, Meiyue  
May 11, 2023
Numerical Algorithms

The locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm is a popular approach for computing a few smallest eigenvalues and the corresponding eigenvectors of a large Hermitian positive definite matrix A. In this work, we propose a mixed precision variant of LOBPCG that uses a (sparse) Cholesky factorization of A computed in lower precision as the preconditioner. To further enhance performance, a mixed precision orthogonalization strategy is proposed. To analyze the impact of reducing precision in the preconditioner on performance, we carry out a rounding error and convergence analysis of PINVIT, a simplified variant of LOBPCG. Our theoretical results predict and our numerical experiments confirm that the impact on convergence remains marginal. In practice, our mixed precision LOBPCG algorithm typically reduces the computation time by a factor of 1.4-2.0 on both CPUs and GPUs.

  • Details
  • Metrics
Type
research article
DOI
10.1007/s11075-023-01550-9
Web of Science ID

WOS:000986804900004

Author(s)
Kressner, Daniel  
Ma, Yuxin
Shao, Meiyue  
Date Issued

2023-05-11

Publisher

SPRINGER

Published in
Numerical Algorithms
Subjects

Mathematics, Applied

•

Mathematics

•

symmetric eigenvalue problem

•

lobpcg algorithm

•

mixed precision algorithm

•

iterative refinement

Editorial or Peer reviewed

REVIEWED

Written at

EPFL

EPFL units
ANCHP  
Available on Infoscience
June 5, 2023
Use this identifier to reference this record
https://infoscience.epfl.ch/handle/20.500.14299/197958
Logo EPFL, École polytechnique fédérale de Lausanne
  • Contact
  • infoscience@epfl.ch

  • Follow us on Facebook
  • Follow us on Instagram
  • Follow us on LinkedIn
  • Follow us on X
  • Follow us on Youtube
AccessibilityLegal noticePrivacy policyCookie settingsEnd User AgreementGet helpFeedback

Infoscience is a service managed and provided by the Library and IT Services of EPFL. © EPFL, tous droits réservés