Matrix exponentials and parallel prefix computation in a quantum control problem

T. Auckenthaler, M. Bader, T. Huckle, A. Spörl, K. Waldherr

Research output: Contribution to journalArticlepeer-review

9 Scopus citations

Abstract

Quantum control plays a key role in quantum technology, in particular for steering quantum systems. As problem size grows exponentially with the system size, it is necessary to deal with fast numerical algorithms and implementations. We improved an existing code for quantum control concerning two linear algebra tasks: the computation of the matrix exponential and efficient parallelisation of prefix matrix multiplication. For the matrix exponential we compare three methods: the eigendecomposition method, the Padé method and a polynomial expansion based on Chebyshev polynomials. We show that the Chebyshev method outperforms the other methods both in terms of computation time and accuracy. For the prefix problem we compare the tree-based parallel prefix scheme, which is based on a recursive approach, with a sequential multiplication scheme where only the individual matrix multiplications are parallelised. We show that this fine-grain approach outperforms the parallel prefix scheme by a factor of 2-3, depending on parallel hardware and problem size, and also leads to lesser memory requirements. Overall, the improved linear algebra implementations not only led to a considerable runtime reduction, but also allowed us to tackle problems of larger size on the same parallel compute cluster.

Original languageEnglish
Pages (from-to)359-369
Number of pages11
JournalParallel Computing
Volume36
Issue number5-6
DOIs
StatePublished - May 2010

Keywords

  • Chebyshev polynomials
  • Matrix exponential
  • Parallel matrix multiplication
  • Parallel prefix problem
  • Quantum control algorithm

Fingerprint

Dive into the research topics of 'Matrix exponentials and parallel prefix computation in a quantum control problem'. Together they form a unique fingerprint.

Cite this