Pseudo-differential representation of the metaplectic transform and its application to fast algorithms

Abstract

The metaplectic transform (MT), also known as the linear canonical transform, is a unitary integral mapping which is widely used in signal processing and can be viewed as a generalization of the Fourier transform. For a given function $\psi$ on an $N$-dimensional continuous space $\text{q}$, the MT of $\psi$ is parameterized by a rotation (or more generally, a linear symplectic transformation) of the $2N$-dimensional phase space $(\text{q},\text{p})$, where $\text{p}$ is the wavevector space dual to $\text{q}$. Here, we derive a pseudo-differential form of the MT. For small-angle rotations, or near-identity transformations of the phase space, it readily yields asymptotic \textit{differential} representations of the MT, which are easy to compute numerically. Rotations by larger angles are implemented as successive applications of $K\gg 1$ small-angle MTs. The algorithm complexity scales as $O(K{N}^3{N}_p)$, where $N_p$ is the number of grid points. We present a numerical implementation of this algorithm and discuss how to mitigate the associated numerical instabilities.

arXiv comment16 pages, 9 figures, 4 appendices ©2019 Optical Society of America. One print or electronic copy may be made for personal use only. Systematic reproduction and distribution, duplication of any material in this paper for a fee or for commercial purposes, or modifications of the content of this paper are prohibited

References in corpus (3)

• Quasioptical modeling of wave beams with and without mode conversion: I. Basic theory
• Synchrotron radiation representation in phase space
• Two-dimensional nonseparable discrete linear canonical transform based on CM-CC-CM-CC decomposition

Cited by in corpus (10)

• Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics
• Restoring geometrical optics near caustics using sequenced metaplectic transforms
• Metaplectic geometrical optics for modeling caustics in uniform and nonuniform media
• Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms
• Quasioptical modeling of wave beams with and without mode conversion: IV. Numerical simulations of waves in dissipative media
• Modeling circulating cavity fields using the discrete linear canonical transform
• Exactly unitary discrete representations of the metaplectic transform for linear-time algorithms
• Demonstration of Metaplectic Geometrical Optics for Reduced Modeling of Plasma Waves
• Regarding the extension of metaplectic geometrical optics to modelling evanescent waves in ray-tracing codes
• Gaussian beams and caustic avoidance in gravitational optics