def METHOD_NAME(
    reference: types.Reconstruction, candidate: types.Reconstruction
) -> np.ndarray:
    common_points = set(reference.points.keys()).intersection(
        set(candidate.points.keys())
    )
    return np.array(
        [
            reference.points[p].coordinates - candidate.points[p].coordinates
            for p in common_points
        ]
    )