  We summarize our efforts to create a numerical implementation of the Local
Analytical Subtraction Scheme (LASS) for obtaining NNLO QCD predictions in
electron-positron collisions. We focus on the regularization of double-real
radiation contributions, which are typically the most computationally intensive
part of NNLO calculations. We describe the structure of the subtraction terms
in LASS and our approach to numerically validate them. Using arbitrary
precision arithmetic, we demonstrate the proper convergence of individual
subtraction terms, overlap removal terms, and spurious singularity
cancellations in various singular limits. We show results for the specific case
of three-jet production in $e^+e^-$ collisions. Our results validate the
correctness of both the LASS formalism and our numerical implementation,
setting the stage for full NNLO calculations using this method.
