def METHOD_NAME(self, ctx: np.recarray, imts, mean, sig, tau, phi):
    """
    See :meth:`superclass method
    <.base.GroundShakingIntensityModel.compute>`
    for spec of input and result values.
    """
    R = ctx.rrup
    M = ctx.mag
    # get constants
    Ssr = get_Ssr_term(ctx.vs30)
    Shr = get_Shr_term(ctx.vs30)
    rake = ctx.rake
    F = get_fault_term(rake)
    # compute mean
    mean[:] = -3.512 + (0.904 * M) - (
        1.328 * np.log(np.sqrt(R**2 + (0.149 * np.exp(0.647 * M))**2))) \
        + (1.125 - 0.112 * np.log(R) - 0.0957 * M) * F \
        + (0.440 - 0.171 * np.log(R)) * Ssr \
        + (0.405 - 0.222 * np.log(R)) * Shr
    # standard deviations from mean (pg 164; more robust than
    # estimate using magnitude
    mean = np.exp(mean)
    sig[:] = 0.39
    sig[mean < 0.068] = 0.55
    idx = np.logical_and(mean >= 0.068, mean <= 0.21)
    sig[idx] = 0.173 - 0.140 * np.log(mean[idx])