I am reading Normal Approximations with Malliavin Calculus: From Stein's Method to Universality by Stein and Peccati and got stuck on problem 1.7.4.
We define the OU-semigroup of operators to act on $f\in \mathcal S$ (Schwartz space) by the following:
$$P_t f(x)=\int_{-\infty}^{\infty} f(e^{-t}x+\sqrt{1-e^{-2t}}y)d\gamma(y)$$
where $\gamma$ is the standard normal measure on $\Bbb R$.
We define the OU-process as the solution to
$$dX_t=\sqrt 2 dB_t-X_t dt, \ X_0=x$$
I've shown the solution to this is
$$X_t=e^{-t}x+\sqrt{2}\int_0^te^{-(t-s)}dB_s$$
Prove that $P_tf(x)=E(f(X_t))$ for $f\in \mathcal S$
I applied Ito's fomula to $f(X_t)$ and I got:
\begin{align*}df(X_t)&=f'(X_t)dX_t+\frac{1}{2}f''(X_t)2dt\\ &=f'(X_t)(\sqrt 2 dB_t-X_t dt)+f''(X_t)dt\\ &=(-f'(X_t)X_t+f''(X_t))dt+\sqrt{2}f'(X_t) dB_t\end{align*}
Noting that the generator $L$ of $P_t$ acts on $f$ by $Lf(x)=-xf'(x)+f''(x)$ and that $f(X_0)=x$ we have
$$f(X_t)=f(x)+\int_0^t L f(X_s) ds+\int_0^t \sqrt{2} f'(X_s) dB_s$$
Taking expectations of both sides and using Fubini/martingale property of Ito integral gives:
$$E(f(X_t))=f(x)+\int_0^t E(Lf(X_s)) ds$$
Here I am stuck. We also know that $L=\frac{d}{ds}\Bigg|_{s=0} P_s$ and $L=-\delta D$ if those are of use.