added logic to incorporate the spin dependence of comparisons like 'mass <= Mmax(Mtov, spin)' based on Breu+Rezzolla 2016. This forces us to assume a maximum spin, which should be computed based on the EoS itself but is currently taken to be a fixed value that's characteristic of most EoS. We may be able to improve upon that in the future