np.where(mode, 0, A * m * (p_tilde ** (m - 1)) * (p_ref ** m)),