I am trying to apply the propagation of singularities theorem to a distribution $u \in D’(M \times M)$ that verifies $Pu = f$, with $P$ a linear differential operator and $f \in D’(M \times M)$, as usual. But I do no see how to do it properly.
My first idea is to apply the theorem to $P_xu_y(x)=f_y(x)$ for $x \in M$ and for $y$ in some small enough compact set (are there results that could help this kind of proof ?)
And my second idea : $(P_x \otimes 1_y)u(x,y) = f(x,y) $.
At the end I want to know how to pass from the wavefront set of $u_y$ to the one of $u$. Is there something as the pullback of a wavefront set to do so ?
Any ideas, suggestions or references would be nice, thanks.
Edit (after Deane Yang and Bazin’s comments) : Thank you very much for your answers, here are more details about the context.
$M = \mathbb{R} \times \Sigma$ is a Lorentzian manifold equipped with a lorentzian metric $g$, that is globally hyperbolic (i.e. there exists a (Cauchy) hypersurface $\Sigma$). We have $g = -dt^2 +h$, $t$ being the variable in $\mathbb{R} $ and $h$ the induced Riemannian metric on $\Sigma$.
$P$ is a second order linear differential operator acting on sections of $S \rightarrow M$, complex vector bundle over $M$, that is normally hyperbolic, meaning that its principal symbol is given by the metric : $\sigma_P(\xi)= g(\xi^{\#},\xi ^{\#})\cdot id_S $. Moreover, $P$ has product structure over $M$, i.e. $P$ takes the form $$ P = \frac{\partial^2}{\partial t^2} + \Delta $$ with $\Delta$ a Laplace-type operator on $\Sigma$, independent of $t$.