Vectorization of Direct Fock Matrix Construction in DIRAC-DHF Calculations

Iterations of the Fock matrix construction from two-electron integrals completely dominates the total cost of direct self-consistent-field (DSCF) calculations, especially for the four-component Dirac-Hartree-Fock (DHF) relativistic molecular-orbital method. We revive here the notion of vectorization...

Full description

Saved in:
Bibliographic Details
Published in:Journal of nuclear science and technology 2002-02, Vol.39 (2), p.195
Main Authors: MOCHIZUKI, Yuji, MATSUMURA, Masayuki, YOKURA, Tetsui-chi, HIRAHARA, Yukio, IMAMURA, Toshiyuki
Format: Article
Language:English
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Iterations of the Fock matrix construction from two-electron integrals completely dominates the total cost of direct self-consistent-field (DSCF) calculations, especially for the four-component Dirac-Hartree-Fock (DHF) relativistic molecular-orbital method. We revive here the notion of vectorization as helpful to accelerate DSCF-DHF calculations. DIRAC is a representative state-of-the-art DHF code on which we consider the vectorized Fock matrix construction. While the integral generation is straightforwardly vectorizable, the contraction of generated integrals with associated densities into the Fock matrix elements is unable to vectorize simply due to frequent address conflicts. In this report, we present a vectorization of the direct Fock matrix construction by introducing a working-array technique. Total vectorization efficiency for DHF runs is improved from 50-60% by the original DIRAC to more than 90% by the modified version.
ISSN:0022-3131
1881-1248
DOI:10.3327/jnst.39.195