I wonder whether there exists a straightforward extension of the Ito isometry to multidimensional processes.
In the one-dimensional case the Ito isometry can be written as
$\mathbb{E}[ (\int_0^T X_t \; \mathrm{d}W_t)^2 ] = \mathbb{E}[ (\int_0^T X_t^2 \;\mathrm{d}_t) ]$.
If now $X_t$ is a vector of random variables instead, do I get something along these lines:
$\mathbb{E}[ (\int_0^T X_t \; \mathrm{d}W_t) (\int_0^T X_t^\top \; \mathrm{d}W_t^\top) ] = \mathbb{E}[ (\int_0^T X_t X_t^\top \;\mathrm{d}_t) ]$
????
Short answer: Yes, what you wrote is basically correct.
Longer anser:
${\bf Statement.}$ Let ${\bf B}: [0, T]\times \Omega \rightarrow \mathbf{R}^p$ be the standard $p$-dimensional Brownian motion and ${\bf M}: [0, T]\times \Omega \rightarrow \mathbf{R}^{p\times p}$ be a matrix-valued stochastic process adapted to the natural filtration of the Brownian motion. Then $$\mathbf{E} \left[ \left\Vert\int_0^T {\bf M}_t d{\bf B}_t\right\Vert_2^2\right] = \mathbf{E} \left[ \int_0^T \left\Vert{\bf M}_t\right\Vert_F^2 dt\right],$$ where the norm in the right hand side is the Frobenius norm of the matrix-valued process.