Fast, Efficient Calculations of the Two-Body Matrix Elements of the Transition Operators for Neutrinoless Double Beta Decay