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
}