High-performance generation of the Hamiltonian and Overlap matrices in FLAPW methods

详细信息    查看全文

• 作者：Edoardo Di Napoli^a ; ^b ; ^c ; ^{e.di.napoli@fz-juelich.de} ; Elmar Peise^c ; ^{peise@aices.rwth-aachen.de} ; Markus Hrywniak^d ; ^{markus.hrywniak@rwth-aachen.de} ; Paolo Bientinesi^c ; ^{pauldj@aices.rwth-aachen.de}
• 关键词：Density functional theory ; High-performance computing ; Dense linear algebra ; Matrix generation ; Performance portability ; FLAPW ; FLEUR
• 刊名：Computer Physics Communications
• 出版年：2017
• 出版时间：February 2017
• 年：2017
• 卷：211
• 期：Complete
• 页码：61-72
• 全文大小：618 K
• 卷排序：211

文摘

One of the greatest efforts of computational scientists is to translate the mathematical model describing a class of physical phenomena into large and complex codes. Many of these codes face the difficulty of implementing the mathematical operations in the model in terms of low level optimized kernels offering both performance and portability. Legacy codes suffer from the additional curse of rigid design choices based on outdated performance metrics (e.g. minimization of memory footprint). Using a representative code from the Materials Science community, we propose a methodology to restructure the most expensive operations in terms of an optimized combination of dense linear algebra (BLAS3) kernels. The resulting algorithm guarantees an increased performance and an extended life span of this code, enabling larger scale simulations.