Skip to main content
1 of 2
Carlo Beenakker
  • 188.1k
  • 18
  • 448
  • 651

Q1: Can one define the operator $p(x,T)$ (with $T=d/dx$) using the Cauchy integral formula?

That will not work, for example, take $p(x,T)=Tx-xT=1$, but $$p(x,T) = \oint_{C} p(x,z)(zI-T)^{-1}dz= \oint_{C} (zx-xz)(zI-T)^{-1}dz=0\neq 1.$$

Q2: Is there some way to compute the action of $p(x,T)$ on a (smooth) function $f(x)$?

This problem is a common one in physics, for $p$ of the form $p(x,T)=\exp[A(x)+B(T)]$. One then starts from the Trotter formula $$e^{A+B} = \lim_{n \rightarrow \infty} (e^{A/n}e^{B/n})^n$$ and inserts Fourier transform operations ${\cal F}$, $$e^{A(x)+B(T)}f=\lim_{n \rightarrow \infty} (e^{A/n}{\cal F}^{-1}e^{\tilde{B}/n}{\cal F})^n f,\;\;\tilde{B}={\cal F}B{\cal F}^{-1}.$$ The operators $A$ and $\tilde{B}$ are both diagonal, so they act on a function simply by multiplication. The time consuming step is the Fourier transform, if $f$ is discretized in $N$ points that operation takes a time that scales as $N\log N$.

Carlo Beenakker
  • 188.1k
  • 18
  • 448
  • 651