def METHOD_NAME(self):
    # HIDL> query_disc, 8, [ 0.17101007,  0.03015369,  0.98480775],6,listpix,/DEG,NESTED=0
    # HIDL> print,listpix
    #           4
    np.testing.assert_array_equal(
        query_disc(self.NSIDE, self.vec, self.radius, inclusive=False),
        np.array([4]),
    )