def METHOD_NAME(filename, read_labels=False):
    # Returns an ndarray of shape (n,49) where n is number of samples
    f = open(filename)
    digits_lst = []
    labels_lst = []
    for line in f:
        tok = line.split(",")
        digits = unpackbits(tok[0])
        digits_lst.append(digits)
        if read_labels:
            labels_lst.append(onehot(int(tok[1]), 10))
    if read_labels:
        return np.stack(digits_lst), np.stack(labels_lst)
    else:
        return np.stack(digits_lst)