def METHOD_NAME(sample_shape, batch_shape):
    alphas = torch.ones(batch_shape)
    dist = RejectionStandardGamma(alphas)
    x = dist.rsample(sample_shape)
    assert x.shape == sample_shape + batch_shape