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