Got the contr_hrr part of the integral code working.

Still a factor of 2 slower than the old PyQuante code.