  As the new-generation precision experiments such as MOLLER and P2 look for
physics beyond Standard Model, it is becoming increasingly important to
evaluate the higher-order electroweak radiative corrections to a sub-percent
level of uncertainty. However, due to propagators with different masses and
higher-order tensor Feynman integrals, the two-loop calculations involving
thousands of Feynman graphs become a demanding task requiring novel
computational approaches. In this paper, we describe our dispersive sub-loop
insertion approach and develop two-loop integrals using two-point functions
basis which is applicable to wide range of processes.
