The development of hierarchical matrices as a tool to accurately approximate dense matrices of large dimensions allows addressing problems of scales that would be difficult to address using the conventional representations. This opens up an important range of opportunities in various domains of application. However, the value of this mathematical tool could be diminished if it cannot make an efficient use of modern hardware architectures. In this work, we implement a recently proposed algorithm to solve Sylvester equations involving two types of hierarchical matrices, HODLR and HALR. We also perform an experimental evaluation on a multicore CPU system to detect the principal performance bottlenecks, and discuss the main challenges to develop high-performance implementation of this type of algorithms.

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Implementation of the Sylvester Equation for HODLR and HALR \(\mathcal {H}\) -Matrices

  • Gonzalo Berger,
  • Florencia Uslenghi,
  • Ernesto Dufrechou,
  • Pablo Ezzatti

摘要

The development of hierarchical matrices as a tool to accurately approximate dense matrices of large dimensions allows addressing problems of scales that would be difficult to address using the conventional representations. This opens up an important range of opportunities in various domains of application. However, the value of this mathematical tool could be diminished if it cannot make an efficient use of modern hardware architectures. In this work, we implement a recently proposed algorithm to solve Sylvester equations involving two types of hierarchical matrices, HODLR and HALR. We also perform an experimental evaluation on a multicore CPU system to detect the principal performance bottlenecks, and discuss the main challenges to develop high-performance implementation of this type of algorithms.