code: def find_inouts_st(voxel_grid, datapts, **kwargs):
      
      tol = kwargs.get('tol', 10e-8)
      filled = [0 for _ in range(len(voxel_grid))]
      for idx, bb in enumerate(voxel_grid):
          pts_inside = is_point_inside_voxel(bb, datapts, tol=tol)
          if pts_inside:
              filled[idx] = 1
      return filled
docs: Single-threaded ins and outs finding (default)
  
      :param voxel_grid: voxel grid
      :param datapts: data points
      :return: in-outs