Repository logo

Infoscience

  • English
  • French
Log In
Logo EPFL, École polytechnique fédérale de Lausanne

Infoscience

  • English
  • French
Log In
  1. Home
  2. Academic and Research Output
  3. Journal articles
  4. A block algorithm for computing antitriangular factorizations of symmetric matrices
 
research article

A block algorithm for computing antitriangular factorizations of symmetric matrices

Bujanovic, Zvonimir  
•
Kressner, Daniel  
2016
Numerical Algorithms

Any symmetric matrix can be reduced to antitriangular form in finitely many steps by orthogonal similarity transformations. This form reveals the inertia of the matrix and has found applications in, e.g., model predictive control and constraint preconditioning. Originally proposed by Mastronardi and Van Dooren, the existing algorithm for performing the reduction to antitriangular form is primarily based on Householder reflectors and Givens rotations. The poor memory access pattern of these operations implies that the performance of the algorithm is bound by the memory bandwidth. In this work, we develop a block algorithm that performs all operations almost entirely in terms of level 3 BLAS operations, which feature a more favorable memory access pattern and lead to better performance. These performance gains are confirmed by numerical experiments that cover a wide range of different inertia.

  • Details
  • Metrics
Type
research article
DOI
10.1007/s11075-015-9983-8
Web of Science ID

WOS:000369063900003

Author(s)
Bujanovic, Zvonimir  
Kressner, Daniel  
Date Issued

2016

Publisher

Springer

Published in
Numerical Algorithms
Volume

71

Issue

1

Start page

41

End page

57

Subjects

Antitriangular factorization

•

Matrix inertia

•

Block algorithm

•

Symmetric matrix

Editorial or Peer reviewed

REVIEWED

Written at

EPFL

EPFL units
ANCHP  
Available on Infoscience
April 1, 2016
Use this identifier to reference this record
https://infoscience.epfl.ch/handle/20.500.14299/125499
Logo EPFL, École polytechnique fédérale de Lausanne
  • Contact
  • infoscience@epfl.ch

  • Follow us on Facebook
  • Follow us on Instagram
  • Follow us on LinkedIn
  • Follow us on X
  • Follow us on Youtube
AccessibilityLegal noticePrivacy policyCookie settingsEnd User AgreementGet helpFeedback

Infoscience is a service managed and provided by the Library and IT Services of EPFL. © EPFL, tous droits réservés