Exploiting Fast Matrix Arithmetic in Block Low-Rank Factorizations

Jeannerod, Claude-Pierre and Mary, Theo and Pernet, Clément and Roche, Daniel (2019) Exploiting Fast Matrix Arithmetic in Block Low-Rank Factorizations. [MIMS Preprint]

Text main.pdf Download (279kB)

Abstract

We consider the LU factorization of an $n\times n$ matrix $A$ represented as a block low-rank (BLR) matrix: most of its off-diagonal blocks are approximated by matrices of small rank $r$, which reduces the asymptoticcomplexity of computing the LU factorization of $A$ down to $\O(n^2r)$. In this article, our aim is to further reduce this complexity by exploiting fast matrix arithmetic, that is, the ability to multiply two $n\times n$ full-rank matrices together for $\O(n^\w)$ flops, where $\w<3$. This is not straightforward: simply accelerating the intermediate operations performed in the standard BLR factorization algorithm does not suffice to reduce the quadratic complexity in $n$, because these operations are performed on matrices whose size is too small. To overcome this obstacle, we devise a new BLR factorization algorithm that, by recasting the operations so as to work on intermediate matrices of larger size, can exploit more efficiently fast matrix arithmetic. This new algorithm achieves an asymptotic complexity of $\O(n^{(\w+1)/2}r^{(\w-1)/2})$, which represents an asymptotic improvement compared to the standard BLR factorization as soon as $\w<3$. In particular, for Strassen's algorithm, $\w\approx2.81$ yields a complexity of $\O(n^{1.904}r^{0.904})$. Our numerical experiments are in good agreement with this theoretical result.

Item Type:	MIMS Preprint
Subjects:	MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User:	Dr Theo Mary
Date Deposited:	03 Feb 2019 09:46
Last Modified:	03 Feb 2019 09:46
URI:	https://eprints.maths.manchester.ac.uk/id/eprint/2682

Available Versions of this Item

• Exploiting Fast Matrix Arithmetic in Block Low-Rank Factorizations. (deposited 03 Feb 2019 09:46) [Currently Displayed]
  • Exploiting Fast Matrix Arithmetic in Block Low-Rank Factorizations. (deposited 12 Apr 2019 07:26)
    • Improving the Complexity of Block Low-Rank Factorizations with Fast Matrix Arithmetic. (deposited 16 Oct 2019 23:23)
      • Improving the Complexity of Block Low-Rank Factorizations with Fast Matrix Arithmetic. (deposited 01 Dec 2019 11:31)

Actions (login required)

View Item