Transition matrix for a two-dimensional random walk in a torus: compute matrix in R

Question

I would like to build a transition matrix for a torus, in order to model a two dimensional random walk.

I know that for a circular one dimensional random walk, the co