I have been reading Nualart's notes on Malliavin calculus and I am aware of his derivation of the integration by parts formula. We consider a Hilbert space $(H,\langle \cdot,\cdot\rangle )$. If we let $(B_{t})_{t\in [0,\infty)}$ be a $d$-dimensional Brownian motion and for $h\in H$ set $W(h)=\sum_{i=1}^{d}\int_{0}^{\infty}h_{i}(s)\text{d}B_{s}^{i}$. Then we have $$ \mathbb{E}\left[ \langle DF, h\rangle \right] = \mathbb{E}\left[ FW(h) \right] $$ where $DF$ is the weak derivative of $F$. The proof of this relies on choosing a suitable orthonormal system and using that things are normally distributed. However, on the wikipedia page (https://en.wikipedia.org/wiki/Malliavin_calculus) it is stated that it is possible to derive the same (perhaps in a less general setting) identity by Girsanov's theorem.
The details are omitted, so I have tried to fill them in myself, but I am not quite sure that I obtain the same formula as Nualart. The problem seems to me, to be that Girsanov's theorem concerns change of measure, whereas the above identity seems to be under the same probability measure $\mathbb{P}$. As on Wikipedia, we set $\varphi_{t} = \int_{0}^{t}h_{s}\text{d}\langle B\rangle_{s} = \int_{0}^{t}h_{s}\text{d}s$ for a suitable predictable process $h$.
I then use a Girsanov change of measure using the Doléans-Dade exponential $$ \mathcal{E}(\varepsilon h \star B)_{t} = e^{\varepsilon\int_{0}^{t}h_{s}\text{d}B_{s}-\varepsilon^{2}\frac{1}{2}\int_{0}^{t}h_{s}^{2}\text{d}s} $$ If we assume that the underlying filtration is the one generated by $B$, we then set (ignoring that the measures should be restricted to intervals) $$ \frac{\text{d}\mathbb{Q}}{\text{d}\mathbb{P}} = \mathcal{E}(\varepsilon h \star B)_{t} $$ then by Girsanov's theorem we have $$ \text{d}B_{t} = \varepsilon h_{t}\text{d}t + \text{d}B^{\mathbb{Q}}_{t} $$ where $B^{\mathbb{Q}}$ is a Brownian motion under the (equivalent) measure $\mathbb{Q}$. Then we would have \begin{align*} \mathbb{E^P}\left[ F(B)\mathcal{E}(\varepsilon h \star B) \right] &= \int_{\Omega}F(B)\mathcal{E}(\varepsilon h \star B)\text{d}\mathbb{P} \\ &= \int_{\Omega}F(B)\text{d}\mathbb{Q} \\ &= \int_{\Omega}F(B^{\mathbb{Q}} + \varepsilon\varphi)\text{d}\mathbb{Q} \\ &= \mathbb{E^Q}\left[ F(B^{\mathbb{Q}} + \varepsilon\varphi) \right] \end{align*} Differentiating this wrt. $\varepsilon$ and evaluating at $\varepsilon = 0$, we obtain the same kind of formula as in Nualart, but it seems to me that we are still taking expectations under different measures. Am I missing something, or is it even a problem?