The solution is based on the observation that rows in transfer matrices commute, for a certain parametrisation of these four Boltzmann weights. This came about as a modification of an alternate solution for the six-vertex model ; it makes use of elliptic theta functions.