def METHOD_NAME(self, other, **kwargs):
    """Fit the landmarks in `self` to those in `other` in a least squares
    sense, i.e., the transformation matrix moves points from `self` to
    `other`.
    kwargs to pass to `fit_matched_points_analytical`
    """
    assert isinstance(other, Montage)
    assert self.landmarks and other.landmarks
    names = set(self.get_landmark_names()).intersection(
        set(other.get_landmark_names())
    )
    return fit_matched_points_analytical(
        self.get_landmark_pos(names), other.get_landmark_pos(names), **kwargs
    )