Unit test Rd gamma function for high dimensions

[arvoelke]

Should be able to test known values for large dimensions. Can I replace this with a scipy function, or use a better number of iterations? Can also vectorize the first phase of Rd once the gamma value is computed.