def METHOD_NAME():
    """ Test if the first argument allows to use multiple centers"""
    stack = load_structure(TEST_FILE)
    # calculate oxygen RDF for water with and without a selection
    oxygen = stack[:, stack.atom_name == 'OW']
    interval = np.array([0, 10])
    n_bins = 100
    # averaging individual calculations
    bins1, g_r1 = rdf(oxygen[:, 1].coord, oxygen[:, 2:], interval=interval,
                    bins=n_bins, periodic=False)
    bins2, g_r2 = rdf(oxygen[:, 0].coord, oxygen[:, 2:], interval=interval,
                    bins=n_bins, periodic=False)
    mean = np.mean([g_r1, g_r2], axis=0)
    # this should give the same result as averaging for oxygen 0 and 1
    bins, g_r = rdf(oxygen[:, 0:2].coord, oxygen[:, 2:], interval=interval,
                    bins=n_bins, periodic=False)
    assert np.allclose(g_r, mean, rtol=0.0001)