It seems like the identities of curl of gradient, divergence of curl, and the simple derivations of electromagnetic waves from Maxwell equations all rely on the symmetry (interchangeability of their order) of partial derivatives.
But mathematically speaking, this tacitly assumes that the electric/magnetic fields in question are continuously double differentiable, something which needs justification.
Why do we assume that in the real world, the fields are smooth enough for the identities to work, atleast almost everywhere?