it seems a bug in align_two_reconstruction

It seems a bug here:

def align_two_reconstruction(r1, r2, common_tracks, threshold):
"""Estimate similarity transform between two reconstructions."""
t1, t2 = r1.points, r2.points

if len(common_tracks) > 6:

<<<<<<<<<<<<<
p1 = np.array([t1[t[0]].coordinates for t in common_tracks])
p2 = np.array([t2[t[1]].coordinates for t in common_tracks])
<<<<<<<<<<<<<

I think it should be:
p1 = np.array([t1[t].coordinates for t in common_tracks])
p2 = np.array([t2[t].coordinates for t in common_tracks])
here, t stands for the common shot id

    # 3 samples / 100 trials / 50% outliers = 0.99 probability
    # with probability = 1-(1-(1-outlier)^model)^trial
    T, inliers = multiview.fit_similarity_transform(
        p1, p2, max_iterations=100, threshold=threshold)
    if len(inliers) > 0:
        return True, T, inliers
return False, None, None

Start at align_two_reconstruction and inspect how common_tracks is represented and consumed by multiview.fit_similarity_transform. Verify whether each entry is a common shot ID or a pair of IDs, then confirm the corrected indexing produces valid similarity-transform inputs and preserves the function's success and inlier behavior.