def METHOD_NAME(aidx1, aidx2, hidden_nids=hidden_nids, toy_params=toy_params):
    if aidx1 == aidx2:
        return 0
    rng = np.random.RandomState(int(aidx1 + aidx2))
    same = hidden_nids[int(aidx1)] == hidden_nids[int(aidx2)]
    mu, sigma = ut.dict_take(toy_params[same], ['mu', 'sigma'])
    return np.clip(rng.normal(mu, sigma), 0, np.inf)