  The physical diagnosis of the solar atmosphere is achieved by solving the
polarized radiative transfer problem for plasmas in Non-Local Thermodynamical
Equilibrium (NLTE). This scenario poses theoretical challenges for integrating
the radiative transfer equation (RTE) efficiently. Namely, current methods are
limited to constant propagation matrices, thus imposing local solutions. To
spark significant advances, this paper lays the foundations of a formalism that
achieves an efficient non-local integration of the RTE based on the Magnus
expansion. First, we revisit the problem and its solutions in Jones and Stokes
formalisms. Looking at them as equivalent representations of the
Lorentz/Poincar\'e group of rotations, we interpret the RTE in terms of Lie
group theory to show the suitability of the Magnus expansion for obtaining
non-local solutions. We then present a detailed algebraic characterization of
the propagation matrix and combine it with the Magnus expansion to reformulate
the homogenous solution to the RTE in Stokes formalism. Thus, we obtain a
compact evolution operator supporting arbitrary variations of the propagation
matrix to first order in the Magnus expansion. Finally, we reformulate the
corresponding inhomogeneous solution as an equivalent homogeneous system, which
is then solved with the Magnus expansion again. This gives the first efficient
and consistent formal solution of the RTE that furthermore is non-local,
natively accurate, and that separates the integration from the formal solution.
Such disruptive formulation leads to a new whole family of numerical radiative
transfer methods and suggests accelerating NLTE syntheses and inversions with
non-local radiative transfer. With minor cosmetic changes, our results are
valid for other universal physical problems sharing the Lorentz/Poincar\'e
algebra of the RTE and special relativity
