ENH: Subhalo abundance (SHAM) code

[Fox-Davidson]

Description

This PR adds code to compute the subhalo abundance matching of internally generated halo/subhalo and galaxy catalogues based on a monotonic relation between their masses and the quenching model proposed in Peng et al 2010. Galaxy parameters given as examples are from the Schechter functions found in Birrer et al 2015.

Checklist

• [x] Follow the Contributor Guidelines
• [x] Write unit tests
• [x] Write documentation strings
• [x] Assign someone from your working team to review this pull request
• [x] Assign someone from the infrastructure team to review this pull request