... bond.¹

This is reminiscent of the method of [3] for Ising cellular automata, which conserves a quantity measured along bonds, rather than at sites.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

....²

One way around this, however, is to consider the amplitudes to be no longer in a single dimension; amplitudes that are consistent with the l_ij's being differences between them could then be found.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... site.³

Various other methods were tried: choosing some $\Delta p$ based on l^ij and solving for a new l^ij that conserves energy; and vice versa. These were found to have no solution for $\Delta p$ for certain values of l^ij.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... repeated.⁴

Processor i also has a thread running that updates pair (0,i), so each pair updates, on average, twice in time T. This is arbitrary.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

... conditions⁵

In the model, the boundary conditions are imposed by giving the sites at the boundary infinite mass.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.