mrpy.base.core.log_mass_mode¶
-
mrpy.base.core.
log_mass_mode
(logHs, alpha, beta)¶ The mode of the log-space MRP weighted by mass.
Parameters: - logHs, alpha, beta: array_like
Shape parameters of the MRP distribution.
Returns: - lmm : array_like
The log-space mass mode of the MRP.
Examples
This function:
>>> log_mass_mode(14.0,-1.8,0.7) 1.67016714698e+13
yields the same result as generating the mode via differentiation:
>>> from mrpy.base import stats >>> from scipy.interpolate import InterpolatedUnivariateSpline as spline >>> m = np.linspace(13.0,14.0,200) >>> # Add 1 to alpha to generate mass weighting >>> mrp = stats.TGGDlog(14.0,-1.8+1,0.7,m[0]).pdf(m,log=True) >>> s = spline(m,mrp,k=4) >>> 10**s.derivative().roots()[0] 1.67016715e+13