Skip to Main content Skip to Navigation
New interface
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
Contributor : Jean-Guillaume Dumas Connect in order to contact the contributor
Submitted on : Monday, January 4, 2021 - 4:06:57 PM
Last modification on : Tuesday, November 22, 2022 - 2:26:16 PM
Long-term archiving on: : Monday, April 5, 2021 - 8:53:03 PM


Files produced by the author(s)



Jean-Guillaume Dumas, Clément Pernet, Alexandre Sedoglavic. Some fast algorithms multiplying a matrix by its adjoint. Journal of Symbolic Computation, 2023, 115 (March–April), pp.285-315. ⟨10.1016/j.jsc.2022.08.009⟩. ⟨hal-03095393⟩



Record views


Files downloads