def METHOD_NAME(ad_ag1, ad_ag2):
    expected = np.zeros((len(ad_ag1.universe.trajectory),
                         ad_ag1.atoms.n_atoms))
    # calculate distances with PBCs using dist()
    for i, ts in enumerate(ad_ag1.universe.trajectory):
        expected[i] = calc_bonds(ad_ag1, ad_ag2, box=ad_ag1.dimensions)
    return expected