Fastest known method for computing Potts model partition function (Bedini and Jacobsen's "A tree-decomposed transfer matrix for computing exact Potts model partition functions for arbitrary graphs, with applications to planar graph colourings") uses a "tree decomposition" of the square lattice. How can this method be extended to computing two-point correlation function? An obvious approach is to re-run it for every pair of vertices fixed to a particular value, but that ends up with a lot of redundant computation, is there a more efficient way?
Example splitting scheme for grid
