Speedup for angular_diameter_distance_z1z2, due to @nhmc
Update funcs docs to give Quantity return type
Change distmod to use units.mag, remove c_kms

* Following suggestion by @nhmc, use constants.c directly instead
  of c_kms
* For consistency, distmod also returns a quantity object.