def METHOD_NAME(self):
    """Ensures that GaussianDensity can call compute
    multiple times with different data"""
    width = 20
    r_max = 9.9
    sigma = 2
    gd = freud.density.GaussianDensity(width, r_max, sigma)
    for num_points in [1, 10, 100]:
        system = freud.data.make_random_system(width, num_points, is2D=False)
        values = np.random.rand(num_points)
        gd.compute(system, values)
        assert np.isclose(np.sum(gd.density), np.sum(values), rtol=1e-4)