Abstract:The solution of (generalized) eigenvalue problems for symmetric or Hermitian matrices is a common subtask of many numerical calculations in electronic structure theory or materials science. Solving the eigenvalue problem can easily amount to a sizeable fraction of the whole numerical calculation. For researchers in the field of computational materials science, an efficient and scalable solution of the eigenvalue problem is thus of major importance. The ELPA-library is a well-established dense direct eigenvalue solver library, which has proven to be very efficient and scalable up to very large core counts. In this paper, we describe the latest optimizations of the ELPA-library for new HPC architectures of the Intel Skylake processor family with an AVX-512 SIMD instruction set, or for HPC systems accelerated with recent GPUs. We also describe a complete redesign of the API in a modern modular way, which, apart from a much simpler and more flexible usability, leads to a new path to access system-specific performance optimizations. In order to ensure optimal performance for a particular scientific setting or a specific HPC system, the new API allows the user to influence in straightforward way the internal details of the algorithms and of performance-critical parameters used in the ELPA-library. On top of that, we introduced an autotuning functionality, which allows for finding the best settings in a self-contained automated way. In situations where many eigenvalue problems with similar settings have to be solved consecutively, the autotuning process of the ELPA-library can be done "on-the-fly". Practical applications from materials science which rely on so-called self-consistency iterations can profit from the autotuning. On some examples of scientific interest, simulated with the FHI-aims application, the advantages of the latest optimizations of the ELPA-library are demonstrated.

NEP: a module for the parallel solution of nonlinear eigenvalue problems in SLEPc

Introduction to StarNEig -- A Task-based Library for Solving Nonsymmetric Eigenvalue Problems

Nonlinearizing Two-parameter Eigenvalue Problems

Solving large-scale nonlinear eigenvalue problems by rational interpolation approach and resolvent sampling based Rayleigh-Ritz method

Task-based, GPU-accelerated and Robust Library for Solving Dense Nonsymmetric Eigenvalue Problems

Riemannian Linearized Proximal Algorithms for Nonnegative Inverse Eigenvalue Problem.

A Nonlinear Generalization of the Bauer-Fike Theorem and Novel Iterative Methods for Solving Nonlinear Eigenvalue Problems

An inexact Matrix-Newton method for solving NEPv

Resolvent sampling based Rayleigh-Ritz method for large-scale nonlinear eigenvalue problems

Optimizations of the Eigensolvers in the ELPA Library

NE/SQP: A robust algorithm for the nonlinear complementarity problem

GenELPA:A Parallel General Eigenvalue Solver

Automatic rational approximation and linearization of nonlinear eigenvalue problems

Solving a Class of Nonlinear Eigenvalue Problems by Newton's Method

A Parallel Function Evaluation Approach for Solution to Large-Scale Equation-Oriented Models.

Avoiding discretization issues for nonlinear eigenvalue problems

Disguised and new Quasi-Newton methods for nonlinear eigenvalue problems

NonlinearSolve.jl: High-Performance and Robust Solvers for Systems of Nonlinear Equations in Julia

Locally Unitarily Invariantizable NEPv and Convergence Analysis of SCF

Sylvester-based preconditioning for the waveguide eigenvalue problem

Hierarchical subspace evolution method for super large parallel computing: A linear solver and an eigensolver as examples