def METHOD_NAME(self, num_points):
    # Ensure that each point's Gaussian sums to 1
    width = 20
    r_max = 9.9
    sigma = 2
    box_size = width
    gd = freud.density.GaussianDensity(width, r_max, sigma)
    box, points = freud.data.make_random_system(box_size, num_points, is2D=True)
    gd.compute(system=(box, points))
    # This has discretization error as well as single-precision error
    assert np.isclose(np.sum(gd.density), num_points, rtol=1e-4)