Let me give an example of a rather roundabout derivation which uses this fact.
Suppose that instead of doing the usual theorist’s $c=1$ trick, I wish to fully work out a strange unit system where
$$
\begin{align}
\nabla \cdot E&=c\rho,&\nabla\times E &=-\mathring B\\
\nabla\cdot B&=0,&\nabla\times B&=J + \mathring E
\end{align}
$$
where $\mathring A = \dot A/c.$ One can check that the units align properly with $\nabla\cdot J+ c\mathring\rho=0$being a valid continuity equation and curling the curl one finds e.g. $$\overset{\,\scriptsize\circ\circ}B -\nabla^2 B=\square B=\nabla\times J$$ which has the right d'Alembert operator with the right wave velocity, so it looks very promising, is just not Gaussian/CGS or SI.
Well, there is a lot of work and cross-checking in rebuilding all of your knowledge from the ground up, and relativity would be a good guide. We start with the standard definition of vector potential, $\nabla\cdot B=0$ implies $B =\nabla\times A$ for some $A$, which means $\nabla\times(E+\mathring A)=0$ which we use to say $E = -\mathring A - \nabla\varphi$ for some $\varphi$.
Defining $\lambda = \nabla\cdot A + \mathring\varphi$, the other two Maxwell equations say$$\square\varphi = c\rho +\mathring\lambda,\\
\square A= J -\nabla\lambda,$$and our gauge freedom means that mapping $A\mapsto A +\nabla \psi$ while $\varphi\mapsto \varphi-\mathring \psi$ preserves $E, B$ while mapping $\lambda\mapsto \lambda -\square\psi$, which we can solve for zero to force $\lambda \mapsto 0$, the Lorenz gauge. Then since $(c\rho, J) = J^\bullet$ is a 4-vector and $\square$ is covariant, we find that the appropriate 4-potential is $A^\bullet=(\varphi, A)$, no division or multiplication by $c$. Use the $({+}\,{–}\,{–}\,{–})$ metric, $A_\bullet = (\varphi, -A),$ we are ready to form the field tensor.
We have $F_{\mu\nu}=\partial_\mu A_\nu -\partial_\nu A_\mu$ which means $$F_{\bullet\bullet}=\begin{bmatrix}0&-E_x&-E_y&-E_z\\
E_x&0&B_z&-B_y\\
E_y&-B_z&0&B_x\\
E_z&B_y&-B_x&0\end{bmatrix}.$$
Now, you ask in what sense this tensor is a bilinear function, well, there are two answers to that, one is just, it is a matrix, so obviously it is a bilinear function. That is if we take two 4-vectors and write them as raw linear algebra column vectors $\mathbf u,\mathbf v$, and regard this matrix as $\mathbf M$, then this implements the bilinear function $$F(u, v) = \mathbf u^T \mathbf M \mathbf v$$ which is also a Lorentz-invariant scalar. Note that this is a raw transpose of the column, no components are negated here because those transforms are already absorbed into the above matrix.
The second answer is a bit more physical, this should have the function of connecting a 4-velocity to a 4-force. The 4-force can be regarded as a covector in the sense that it takes a small 4-displacement and tells you how much work is done on that displacement. This leads to a slightly redundant situation because of course the displacement that we would want is in the direction of the four velocity of the charged particle that we are tracking, so we find that we actually want $F(v, v)$ for the same 4-velocity when all is said and done, but since F is antisymmetric that will inevitably be zero! In 4D the “length” of the 4-velocity is actually fixed, so no “4-work” is ever truly done. Nevertheless we might want to take the dual of the 4-force and think of it as a change in 4-momentum per unit proper time, obviously the change has to be “perpendicular” to the 4-momentum per the above but that doesn't make it zero.
So let's do that. We find that we want a force
$${\mathrm dp^\bullet\over\mathrm d\tau} \propto
\begin{bmatrix}1&&&\\
&-1&&\\
&&-1&\\
&&&-1\end{bmatrix}
\begin{bmatrix}0&-E_x&-E_y&-E_z\\
E_x&0&B_z&-B_y\\
E_y&-B_z&0&B_x\\
E_z&B_y&-B_x&0\end{bmatrix}
\begin{bmatrix}\gamma c\\
\gamma v_x\\
\gamma v_y\\
\gamma v_z\end{bmatrix},$$
and we know that we want for example ${\mathrm dp^x\over\mathrm dt}=\gamma^{-1} {\mathrm dp^x\over\mathrm d\tau}=qE_x$ for a purely electrostatic case. Thus we have,
$$
{\mathrm dp^\mu\over\mathrm d\tau}=-\frac{q}{c}\eta^{\mu\nu}F_{\nu\sigma} v^\sigma,$$
and we thus come to the final discovery that the proper version of the Lorentz force law in these units is exactly the same as in CGS,
$$
F = q \left( E +\frac vc\times B\right).$$
There's probably an easier way to see that, but it's a nice application of the bilinear properties of the electromagnetic force tensor.
Another application of this fact is that $E^2-B^2=\frac12 F_{\mu\nu}F^{\nu\mu}$ is a Lorentz invariant scalar field, and exploiting a hidden symmetry in such antisymmetric tensors (Hodge star?), so is $E\cdot B \propto F_{\mu\nu} (\star F)^{\nu\mu}.$ Both of these results come from the coordinate invariance of trace, which can roughly be stated as the slightly more geometric, “any [m, n]-tensor can be expressed as a finite sum of outer products of vectors and [m–1, n]-tensors, or a finite sum of products of covectors and [m, n–1]-tensors.” Then the procedure to “contract” an [m, n]-tensor is to choose two indices to contract, express it as a finite sum of outer products of [m–1,n–1]-tensors times a vector times a covector, feed all of those vectors to covectors to create invariant scalars, a scalar times a tensor is a tensor and a sum of tensors of the same shape is a tensor, so there you go.
Either way, can you not first Post your idea of each of those and only then ask about their differences?
– Robbie Goodwin Feb 13 '23 at 21:59