The member function compute pair(...) of the integral object can then be called repeatedly to iterate over all combinations of angular components (n xi , n yi , n zi ) and (n x j , n y j , n z j ) belonging to the shells with total angular momenta l i and l j :