  The numerically stable evaluation of scattering matrix elements near the
infrared limit of gauge theories is of great importance for the success of
collider physics experiments. We present a novel algorithm that utilizes double
precision arithmetic and reaches higher precision than a naive quadruple
precision implementation at smaller computational cost. The method is based on
physics-driven modifications to propagators, vertices and external
polarizations.
