HODLR2D: A new class of Hierarchical matrices

04/12/2022
by   V A Kandappan, et al.
0

This article introduces HODLR2D, a new hierarchical low-rank representation for a class of dense matrices arising out of N body problems in two dimensions. Using this new hierarchical framework, we propose a new fast matrix-vector product that scales almost linearly. We apply this fast matrix-vector product to accelerate the iterative solution of large dense linear systems arising out of radial basis function interpolation and discretized integral equation. The space and computational complexity of HODLR2D matrix-vector products scales as 𝒪(pN log(N)), where p is the maximum rank of the compressed matrix subblocks. We also prove that p ∈𝒪(log(N)log(log(N))), which ensures that the storage and computational complexity of HODLR2D matrix-vector products remain tractable for large N. Additionally, we also present the parallel scalability of HODLR2D as part of this article.

READ FULL TEXT

page 4

page 15

research
07/30/2023

HODLR3D: Hierarchical matrices for N-body problems in three dimensions

This article introduces HODLR3D, a class of hierarchical matrices arisin...
research
08/14/2019

Accuracy Controlled Structure-Preserving H^2-Matrix-Matrix Product in Linear Complexity with Change of Cluster Bases

H^2-matrix constitutes a general mathematical framework for efficient c...
research
12/19/2022

Matrix recovery from matrix-vector products

Can one recover a matrix efficiently from only matrix-vector products? I...
research
01/30/2023

Algebraic Inverse Fast Multipole Method: A fast direct solver that is better than HODLR based fast direct solver

This article presents a fast direct solver, termed Algebraic Inverse Fas...
research
09/12/2019

The Fast and Free Memory Method for the efficient computation of convolution kernels

We introduce the Fast Free Memory method (FFM), a new fast method for th...
research
10/11/2021

Sparse recovery of elliptic solvers from matrix-vector products

In this work, we show that solvers of elliptic boundary value problems i...
research
09/13/2022

Numerical rank of singular kernel functions

We study the rank of sub-matrices arising out of kernel functions, F(x,y...

Please sign up or login with your details

Forgot password? Click here to reset