def METHOD_NAME(self):
    fftpack = pytest.importorskip("scipy.fftpack")
    fft = fftpack.fft
    fftshift = fftpack.fftshift
    width = 20
    r_max = 10.0
    sigma = 0.1
    num_points = 10000
    box_size = r_max * 3.1
    box, points = freud.data.make_random_system(box_size, num_points, is2D=True)
    for w in (width, (width, width), [width, width]):
        gd = freud.density.GaussianDensity(w, r_max, sigma)
        # Test access
        with pytest.raises(AttributeError):
            gd.box
        with pytest.raises(AttributeError):
            gd.density
        gd.compute((box, points))
        # Test access
        gd.box
        gd.density
        # Verify the output dimensions are correct
        assert gd.density.shape == (width, width)
        assert np.prod(gd.density.shape) == np.prod(gd.width)
        myDiff = gd.density
        myFFT = fft(fft(myDiff[:, :], axis=1), axis=0)
        myDiff = (myFFT * np.conj(myFFT)).real
        myDiff = fftshift(myDiff)[:, :]
        npt.assert_equal(
            np.where(myDiff == np.max(myDiff)),
            (np.array([width // 2]), np.array([width // 2])),
        )