I have the following decomposition for the electric component of light:
$$\renewcommand{\vec}[1]{\mathbf{#1}}\vec{E}(\vec r)=\frac1{4\pi^2} \iint_\Omega \vec A(k_x, k_y) \mathrm{e}^{i \vec{k} \cdot \vec{r}} \, \mathrm{d}k_x \mathrm{d}k_y.$$
Similarly, the magnetic field is: $$\renewcommand{\vec}[1]{\mathbf{#1}}\vec{H}(\vec r)=\frac1{4\pi^2} \iint_\Omega \frac{\vec k}{\omega \mu_0} \times \vec A(k_x, k_y) \mathrm{e}^{i \vec{k} \cdot \vec{r}} \, \mathrm{d}k_x \mathrm{d}k_y.$$
Fine, this is the setting. Now I wish to compute the average Poynting vector $\langle \vec S \rangle$:
$$\langle \vec S \rangle = \frac12 \operatorname{Re} [\vec E(\vec r) \times \vec H(\vec r)^*].$$ Is there a way to express $\langle \vec S \rangle$ in a nice form? Integrals of functions of $\vec A$ for example. I get horrendous expressions with convolutions that don't give me a nice compact formula.
I know that this might actually be more mathematics than physics, but it might be possible that the physical boundary conditions give a better solution.