Skip to Main content Skip to Navigation
Journal articles

Some fast algorithms multiplying a matrix by its adjoint

Abstract : We present a non-commutative algorithm for the multiplication of a 2 × 2 block-matrix by its adjoint, defined by a matrix ring anti-homomorphism. This algorithm uses 5 block products (3 recursive calls and 2 general products)over C or in positive characteristic. The resulting algorithm for arbitrary dimensions is a reduction of multiplication of a matrix by its adjoint to general matrix product, improving by a constant factor previously known reductions. We prove also that there is no algorithm derived from bilinear forms using only four products and the adjoint of one of them. Second we give novel dedicated algorithms for the complex field and the quaternions to alternatively compute the multiplication taking advantage of the structure of the matrix-polynomial arithmetic involved. We then analyze the respective ranges of predominance of the two strategies. Finally we propose schedules with low memory footprint that support a fast and memory efficient practical implementation over a prime field.
Document type :
Journal articles
Complete list of metadata

https://hal.archives-ouvertes.fr/hal-03095393
Contributor : Jean-Guillaume Dumas <>
Submitted on : Monday, January 4, 2021 - 4:06:57 PM
Last modification on : Wednesday, September 22, 2021 - 11:01:33 AM
Long-term archiving on: : Monday, April 5, 2021 - 8:53:03 PM

Files

aat_adjoint.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-03095393, version 1
  • ARXIV : 2101.01025

Citation

Jean-Guillaume Dumas, Clément Pernet, Alexandre Sedoglavic. Some fast algorithms multiplying a matrix by its adjoint. Journal of Symbolic Computation, Elsevier, In press. ⟨hal-03095393⟩

Share

Metrics

Record views

133

Files downloads

155