@decorator_util.jit()<EOL>def total_sparse_pixels_from_mask(mask, unmasked_sparse_grid_pixel_centres):
total_sparse_pixels = <NUM_LIT:0><EOL>for unmasked_sparse_pixel_index in range(unmasked_sparse_grid_pixel_centres.shape[<NUM_LIT:0>]):<EOL><INDENT>y = unmasked_sparse_grid_pixel_centres[unmasked_sparse_pixel_index, <NUM_LIT:0>]<EOL>x = unmasked_sparse_grid_pixel_centres[unmasked_sparse_pixel_index, <NUM_LIT:1>]<EOL>if not mask[y,x]:<EOL><INDENT>total_sparse_pixels += <NUM_LIT:1><EOL><DEDENT><DEDENT>return total_sparse_pixels<EOL>
Given the full (i.e. without removing pixels which are outside the regular-masks) pixelization grid's pixel centers
    and the regular-masks, compute the total number of pixels which are within the regular-masks and thus used by the
    pixelization grid.

    Parameters
    -----------
    mask : ccd.masks.Mask
        The regular-masks within which pixelization pixels must be inside
    unmasked_sparse_grid_pixel_centres : ndarray
        The centres of the unmasked pixelization grid pixels.