tensor calculus - Calculating electromagnetic invariant in matrix form

I'm kind of confused. I want to calculate the electromagnetic invariant $I := F^{\mu\nu}F_{\mu\nu}$, but I'm not sure what is the easiest way to do so. So, I was trying to do it in matrix form, i.e. defining $$\mathbf{F}:=\begin{pmatrix}0 & -E_{1} & -E_{2} & -E_{3}\\ E_{1} & 0 & -B_{3} & B_{2}\\ E_{2} & B_{3} & 0 & -B_{1}\\ E_{3} & -B_{2} & B_{1} & 0 \end{pmatrix}$$ and then calculate the quantity $I$, but I'm not sure how to obtain the matrix form for $I$ (By "matrix form" I mean expressed in terms of the matrix $\mathbf{F}$).

So, I have two questions, 1) what is the easiest way to calculate $I$?, and 2) how to obtain the matrix form for $I$ starting with $I := F^{\mu\nu}F_{\mu\nu}$?. I'm using the metric $\eta:=diag(1,-1,-1,-1)$.

Answer

Notice that if we define matrices $F = (F^{\mu\nu})$ and $\eta = (\eta_{\mu\nu})$, then notice that $$I = F_{\mu\nu}F^{\mu\nu} = \eta_{\mu\alpha}\eta_{\nu\beta}F^{\alpha\beta}F^{\mu\nu} = \eta_{\mu\alpha}F^{\alpha\beta}\eta_{\nu\beta}F^{\mu\nu} = -\eta_{\mu\alpha}F^{\alpha\beta}\eta_{\beta\nu}F^{\nu\mu} =-\mathrm {tr}(\eta F\eta F)$$ where $\mathrm{tr}$ denotes the trace and I have used antisymmetry and symmetry of $F$ and $\eta$ respectively. I'm guessing this is the type of matrix expression you're after?

Cheers!