(Before marking this question as a duplicate, please consider I've read this post but it I didn't find the answers to it quite satisfactory regarding my doubt).
I'm trying to derive the expression of Poynting's vector (averaged over time) for plane harmonic waves. I start from the well-known expression for harmonic waves and then substitute the fields by the ones corresponding to plane harmonic waves: $$\langle \vec{S}\rangle = \frac{1}{2}Re(\vec{E}_0\times \vec{H}_0^\star) \\ =\frac{1}{2\mu_0\omega}Re(\vec{E}_0\times (\vec{k}^\star\times \vec{E}_0^\star)) \\ =\frac{1}{2\mu_0\omega}Re(\vec{k}^\star(\vec{E}_0\cdot \vec{E}_0^\star))$$ where $\cdot$ is the dot product defined in $\mathbb{C}^n$. The next step would be to admit $\vec{E}_0\cdot \vec{E}_0^\star=|\vec{E}_0|^2$, but I don't find that correct since a vector's norm is given by the dot product of said vector with itself. That is to say, $|\vec{E}_0|^2\equiv \langle\vec{E}_0,\vec{E}_0\rangle = \vec{E}_0\cdot \vec{E}_0\neq \vec{E}_0\cdot \vec{E}_0^\star$.
Of course, I know my issue has everything to do with the fact that I'm not understanding how to use the dot product correctly in this operation, and I'm also aware that, for a complex vector, the dot product with any other complex vector requires to transpose and conjugate the first one, such that: $\langle\vec{u},\vec{v}\rangle=\vec{u}^\star\cdot\vec{v}$ where $\cdot$ is now just notation to represent the product of a row vector with a column vector. It would all make sense if the dot product in my calculations were just the product of these row and column vectors, but I don't think it's that way since the identity I've used to get from the second to the third line is: $$\vec a\times(\vec b\times \vec c) = \vec b(\vec a\cdot \vec c)-\vec c(\vec a\cdot \vec b)$$
EDIT:
I think my problem arises from the identity:
$$\vec a\times(\vec b\times \vec c) = \vec b(\vec a\cdot \vec c)-\vec c(\vec a\cdot \vec b)$$
Is this dot an actual dot product defined in $\mathbb{C}^n$ or is it just notation for the product of a row and a column vector such that: $\vec{E}_0^\star\cdot\vec{E}_0=<\vec{E}_0,\vec{E}_0>=|\vec{E}_0|^2$?