Superfast Fourier Transform Using QTT Approximation

Sergey Dolgov, Boris Khoromskij, Dmitry Savostyanov

Research output: Contribution to journalArticle

Abstract

We propose Fourier transform algorithms using QTT format for data-sparse approximate representation of one- and multi-dimensional vectors (m-tensors). Although the Fourier matrix itself does not have a low-rank QTT representation, it can be efficiently applied to a vector in the QTT format exploiting the multilevel structure of the Cooley-Tukey algorithm. The m-dimensional Fourier transform of an n×⋯×n vector with n=2 d has O(md 2R 3) complexity, where R is the maximum QTT-rank of input, output and all intermediate vectors in the procedure. For the vectors with moderate R and large n and m the proposed algorithm outperforms the O(n mlog n) fast Fourier transform (FFT) algorithm and has asymptotically the same log-squared complexity as the superfast quantum Fourier transform (QFT) algorithm. By numerical experiments we demonstrate the examples of problems for which the use of QTT format relaxes the grid size constrains and allows the high-resolution computations of Fourier images and convolutions in higher dimensions without the 'curse of dimensionality'. We compare the proposed method with Sparse Fourier transform algorithms and show that our approach is competitive for signals with small number of randomly distributed frequencies and signals with limited bandwidth.

Original languageEnglish
Pages (from-to)915-953
Number of pages39
JournalJournal of Fourier Analysis and Applications
Volume18
Issue number5
DOIs
Publication statusPublished - 1 Oct 2012

Fingerprint

Fourier transform
Fourier transforms
Approximation
Sparse Data
Curse of Dimensionality
Fast Fourier transform
Convolution
Fast Fourier transforms
Higher Dimensions
High Resolution
Bandwidth
Numerical Experiment
Grid
Output
Demonstrate
Experiments

Keywords

  • Convolution
  • Fourier transform
  • High-dimensional problems
  • QTT
  • Quantum Fourier transform
  • Sparse Fourier transform
  • Tensor train format

Cite this

Dolgov, Sergey ; Khoromskij, Boris ; Savostyanov, Dmitry. / Superfast Fourier Transform Using QTT Approximation. In: Journal of Fourier Analysis and Applications. 2012 ; Vol. 18, No. 5. pp. 915-953.
@article{9d3eb9d478554009bd0f9c4b088dd639,
title = "Superfast Fourier Transform Using QTT Approximation",
abstract = "We propose Fourier transform algorithms using QTT format for data-sparse approximate representation of one- and multi-dimensional vectors (m-tensors). Although the Fourier matrix itself does not have a low-rank QTT representation, it can be efficiently applied to a vector in the QTT format exploiting the multilevel structure of the Cooley-Tukey algorithm. The m-dimensional Fourier transform of an n×⋯×n vector with n=2 d has O(md 2R 3) complexity, where R is the maximum QTT-rank of input, output and all intermediate vectors in the procedure. For the vectors with moderate R and large n and m the proposed algorithm outperforms the O(n mlog n) fast Fourier transform (FFT) algorithm and has asymptotically the same log-squared complexity as the superfast quantum Fourier transform (QFT) algorithm. By numerical experiments we demonstrate the examples of problems for which the use of QTT format relaxes the grid size constrains and allows the high-resolution computations of Fourier images and convolutions in higher dimensions without the 'curse of dimensionality'. We compare the proposed method with Sparse Fourier transform algorithms and show that our approach is competitive for signals with small number of randomly distributed frequencies and signals with limited bandwidth.",
keywords = "Convolution, Fourier transform, High-dimensional problems, QTT, Quantum Fourier transform, Sparse Fourier transform, Tensor train format",
author = "Sergey Dolgov and Boris Khoromskij and Dmitry Savostyanov",
year = "2012",
month = "10",
day = "1",
doi = "10.1007/s00041-012-9227-4",
language = "English",
volume = "18",
pages = "915--953",
journal = "Journal of Fourier Analysis and Applications",
issn = "1069-5869",
number = "5",

}

Superfast Fourier Transform Using QTT Approximation. / Dolgov, Sergey; Khoromskij, Boris; Savostyanov, Dmitry.

In: Journal of Fourier Analysis and Applications, Vol. 18, No. 5, 01.10.2012, p. 915-953.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Superfast Fourier Transform Using QTT Approximation

AU - Dolgov, Sergey

AU - Khoromskij, Boris

AU - Savostyanov, Dmitry

PY - 2012/10/1

Y1 - 2012/10/1

N2 - We propose Fourier transform algorithms using QTT format for data-sparse approximate representation of one- and multi-dimensional vectors (m-tensors). Although the Fourier matrix itself does not have a low-rank QTT representation, it can be efficiently applied to a vector in the QTT format exploiting the multilevel structure of the Cooley-Tukey algorithm. The m-dimensional Fourier transform of an n×⋯×n vector with n=2 d has O(md 2R 3) complexity, where R is the maximum QTT-rank of input, output and all intermediate vectors in the procedure. For the vectors with moderate R and large n and m the proposed algorithm outperforms the O(n mlog n) fast Fourier transform (FFT) algorithm and has asymptotically the same log-squared complexity as the superfast quantum Fourier transform (QFT) algorithm. By numerical experiments we demonstrate the examples of problems for which the use of QTT format relaxes the grid size constrains and allows the high-resolution computations of Fourier images and convolutions in higher dimensions without the 'curse of dimensionality'. We compare the proposed method with Sparse Fourier transform algorithms and show that our approach is competitive for signals with small number of randomly distributed frequencies and signals with limited bandwidth.

AB - We propose Fourier transform algorithms using QTT format for data-sparse approximate representation of one- and multi-dimensional vectors (m-tensors). Although the Fourier matrix itself does not have a low-rank QTT representation, it can be efficiently applied to a vector in the QTT format exploiting the multilevel structure of the Cooley-Tukey algorithm. The m-dimensional Fourier transform of an n×⋯×n vector with n=2 d has O(md 2R 3) complexity, where R is the maximum QTT-rank of input, output and all intermediate vectors in the procedure. For the vectors with moderate R and large n and m the proposed algorithm outperforms the O(n mlog n) fast Fourier transform (FFT) algorithm and has asymptotically the same log-squared complexity as the superfast quantum Fourier transform (QFT) algorithm. By numerical experiments we demonstrate the examples of problems for which the use of QTT format relaxes the grid size constrains and allows the high-resolution computations of Fourier images and convolutions in higher dimensions without the 'curse of dimensionality'. We compare the proposed method with Sparse Fourier transform algorithms and show that our approach is competitive for signals with small number of randomly distributed frequencies and signals with limited bandwidth.

KW - Convolution

KW - Fourier transform

KW - High-dimensional problems

KW - QTT

KW - Quantum Fourier transform

KW - Sparse Fourier transform

KW - Tensor train format

UR - http://www.scopus.com/inward/record.url?scp=84866519225&partnerID=8YFLogxK

U2 - 10.1007/s00041-012-9227-4

DO - 10.1007/s00041-012-9227-4

M3 - Article

AN - SCOPUS:84866519225

VL - 18

SP - 915

EP - 953

JO - Journal of Fourier Analysis and Applications

JF - Journal of Fourier Analysis and Applications

SN - 1069-5869

IS - 5

ER -