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)$.