func BetaBiasFunc(fTh RadialFuncType) RadialFuncType { return func(h *Halo, r float64) float64 { fThR := func(r float64) float64 { return fTh(h, r) } return num.Derivative(fThR, r)(r) * r / fThR(r) } }
func (h *Halo) DPThermalDr(ppt PressureProfileType, pt PressurePopulationType, r float64) float64 { p := func(r float64) float64 { return h.ThermalPressure(ppt, pt, r) } return num.Derivative(p, r)(r) / cosmo.MpcMks }