Radiative transfer with reciprocal transactions: Numerical method and its implementation