English

Approximating Matrices with Multiple Symmetries

Numerical Analysis 2014-12-01 v2 Numerical Analysis

Abstract

If a tensor with various symmetries is properly unfolded, then the resulting matrix inherits those symmetries. As tensor computations become increasingly important it is imperative that we develop efficient structure preserving methods for matrices with multiple symmetries. In this paper we consider how to exploit and preserve structure in the pivoted Cholesky factorization when approximating a matrix AA that is both symmetric (A=ATA=A^T) and what we call {\em perfect shuffle symmetric}, or {\em perf-symmetric}. The latter property means that A=ΠAΠA = \Pi A\Pi where Π\Pi is a permutation with the property that Πv=v\Pi v = v if vv is the vec of a symmetric matrix and Πv=v\Pi v = -v if vv is the vec of a skew-symmetric matrix. Matrices with this structure can arise when an order-4 tensor A\cal A is unfolded and its elements satisfy A(i1,i2,i3,i4)=A(i2,i1,i3,i4)=A(i1,i2,i4,i3)=A(i3,i4,i1,i2).{\cal A}(i_{1},i_{2},i_{3},i_{4}) = {\cal A}(i_{2},i_{1},i_{3},i_{4}) ={\cal A}(i_{1},i_{2},i_{4},i_{3}) ={\cal A}(i_{3},i_{4},i_{1},i_{2}). This is the case in certain quantum chemistry applications where the tensor entries are electronic repulsion integrals. Our technique involves a closed-form block diagonalization followed by one or two half-sized pivoted Cholesky factorizations. This framework allows for a lazy evaluation feature that is important if the entries in A\cal A are expensive to compute. In addition to being a structure preserving rank reduction technique, we find that this approach for obtaining the Cholesky factorization reduces the work by up to a factor of 4.

Keywords

Cite

@article{arxiv.1411.6296,
  title  = {Approximating Matrices with Multiple Symmetries},
  author = {Charles Van Loan and Joseph Vokt},
  journal= {arXiv preprint arXiv:1411.6296},
  year   = {2014}
}

Comments

20 pages

R2 v1 2026-06-22T07:09:10.313Z