def METHOD_NAME(self):
    box = freud.box.Box.square(10)
    num_points = 100
    points = np.random.rand(num_points, 3) * box.L - box.L / 2
    orientations = rowan.random.rand(num_points)
    self.obj.METHOD_NAME((box, points), orientations, neighbors={"r_max": 3})