* Add gamma family functions from SciPy with C code for GPU

The C code is MIT-licensed.

* Remove commented code and bring in missing functions

* Indicate scalar definitions for Elementwise extension

* Upgrade to float64 as previously

* Add docstrings and remove chi2sf duplicate