The mps::solve() routine now operates on pairs of sites.