fn interface_matrix(blocks: &Blocks) -> SMatrix<Complex, 4, 4>
Blocks [[b0, b1], [b2, b3]] of 2-by-2 polarization blocks as one 4-by-4 matrix.
[[b0, b1], [b2, b3]]