def METHOD_NAME():
    """
    Checks if score_genes output agrees with pre-computed reference values.
    The reference values had been generated using the same code
    and stored as a pickle object in ./data
    """
    adata = paul15()
    sc.pp.normalize_per_cell(adata, counts_per_cell_after=10000)
    sc.pp.scale(adata)
    sc.tl.score_genes(adata, gene_list=adata.var_names[:100], score_name='Test')
    with Path(HERE, 'score_genes_reference_paul2015.pkl').open('rb') as file:
        reference = pickle.load(file)
    # assert np.allclose(reference, adata.obs.Test.values)
    assert np.array_equal(reference, adata.obs.Test.values)