---
abstract: 'Darboux transformations for linear operators on regular two dimensional lattices are reviewed. The six point scheme is considered as the master linear problem, whose various specifications, reductions, and their sublattice combinations lead to other linear operators together with the corresponding Darboux transformations. The second part of the review deals with multidimensional aspects of (basic reductions of) the four point scheme, as well as the three point scheme.'
address:
- 'Adam Doliwa, Wydzia[ł]{} Matematyki i Informatyki, Uniwersytet Warmińsko-Mazurski w Olsztynie, ul. Żo[ł]{}nierska 14, 10-561 Olsztyn, Poland'
- 'Maciej Nieszporski, Katedra Metod Matematycznych Fizyki, Uniwersytet Warszawski, ul. Hoża 74, 00-682 Warszawa, Poland'
author:
- Adam Doliwa
- Maciej Nieszporski
title: |
    Darboux transformations for linear operators\
    on two dimensional regular lattices
---

Introduction
============

Darboux transformations are a well known tool in the theory of integrable systems [@ms; @RogersSchief; @GuHuZhou]. The classical Darboux transformation [@DarbouxIV] deals with a Sturm–Liouville problem (the one dimensional stationary Schrödinger equation) generating, at the same time, new potentials and new wave functions from given ones thus providing solutions to the Korteweg–de Vries hierarchy of $1+1$ dimensional integrable systems. However, as it is clearly stated in [@DarbouxIV], it was an earlier work of Moutard [@Moutard] which inspired Darboux. In that work one can find the proper transformation that applies to $2+1$ dimensional integrable systems. Note that the initial area of applications of the Darboux transformations, which preceded the theory of integrable systems, was the theory of conjugate and asymptotic nets where the large body of results on Darboux transformations was formulated [@Jonas; @DarbouxIV; @Eisenhart-TS; @Tzitzeica; @Bianchi; @Lane; @Finikov].

Most of the techniques that allow us to find solutions of integrable non-linear differential equations have been successfully applied to difference equations. These include mutually interrelated methods such as bilinearization method [@Hirota] and the Sato approach [@DJM-I; @DJM], direct linearization method [@FQC; @FWN-Capel], inverse scattering method [@Abl-Lad], the nonlocal $\bar{\partial}$ dressing method [@BoKo], the algebro-geometric techniques [@Krichever-4p].

Also the method of the Darboux transformation has been successfully applied to the discrete integrable systems. The present paper aims to review application of the Darboux transformation technique for the equations that can be regarded as discretizations of second order linear differential equations in two dimensions , their distinguished subclasses and systems of such equations.

While presenting the results we are trying to keep the relation with continuous case. However, we are aware of weak points of this way of exposition. The theory of discrete integrable systems [@Suris; @DIS] is reacher but also, in a sense, simpler then the corresponding theory of integrable partial differential equations. In the course of a limiting procedure, which gives differential systems from the discrete ones, various symmetries and relations between different discrete systems are lost. The classical example is provided (see, for example [@BoKo-KP]) by the hierarchy of the Kadomsev–Petviashvilii (KP) equations, which can be obtained from a single Hirota–Miwa equation — the opposite way, from differential to discrete, involves all equations of the hierarchy [@Miwa].

The structure of the paper is as follows. In section \[dodane\] we expose main ideas in the theory of Darboux transformations in multidimension. In Section \[sec:6p\] we present the construction of the Darboux transformation for the discrete second order linear problem — the 6-point scheme — which can be considered as a discretization of the general second order linear partial differential equation in two variables . Then we discuss various specifications and reductions of the 6-point scheme. We separately present in Section \[sec:s-a\] the Darboux transformations for discrete self-adjoint two dimensional linear systems on the square, triangular and the honeycomb lattices, and their relation to the discrete Moutard transformation. Section \[sec:4p\] is devoted to detailed presentation of the Darboux transformations for systems of the 4-point linear problems, their various specifications and the corresponding permutability theorems. Section \[sec:4p-red\] is dedicated to reductions of the fundamental transformation compatible with additional restrictions on the form of the four point scheme. Finally, in Section \[3p-sch\] we review the Darboux transformations for the 3-point linear problem, and the corresponding celebrated Hirota’s discrete KP nonlinear system.

Jonas fundamental transformation and its basic reductions {#dodane}
=========================================================

The main idea included in the paper has its origin in Jonas paper [@Jonas] where fundamental transformation for conjugate nets has been presented. Neglecting the geometrical context of Jonas paper we only would like to say that the fundamental transformation acts on solutions of compatible system of second order linear differential equations on function $\psi: {\mathbb R}^n \to {\mathbb R}^m$ $$\psi,_{x^ix^j}=A^{ij} \psi,_{x^i}+A^{ji} \psi,_{x^j}, \qquad i,j=1,...,n  \qquad i \ne j$$ where coefficients $A^{ij}$ and $A^{ji}$ are functions of $(x^i)$, and for $n>2$ they have to obey a nonlinear differential equation (compatibility conditions); we use notation that subscript preceded by coma denotes partial differentiation with respect to indicated variables.

Every fundamental transformation $\psi \mapsto \tilde{\psi}$ can be presented as composition of suitably chosen

1.  radial transformation $\psi \mapsto \psi_r$ $$\psi_r=\theta \psi$$

2.  Combescure transformation $\psi_r \mapsto \tilde{\psi_r}$ $${\tilde{\psi_r}},_{x^i}=\alpha_i {\psi_r},_{x^i}$$

3.  radial transformation $\tilde{\psi_r} \mapsto \tilde{\psi}$ $$\tilde{\psi}=\gamma \tilde{\psi_r}$$

We refer interested reader to Eisenhart book [@Eisenhart-TS] for further details. As the practice shows the fundamental transformation is indeed fundamental - in the sense that most (if not all) of the Darboux transformations are reductions of the fundamental transformation.

Here we follow this basic idea discovered almost one hundred years ago. We start from presenting two dimensional difference operators $L$ (and corresponding difference equations $L\Psi=0$, where $\Psi$ is function of discrete variables $m$ and $n$) together with the transformation $(\Psi,L) \mapsto(\tilde{\Psi},\tilde{L})$ which is composition of

1.  gauge transformation $(\Psi,L)\mapsto (\Psi_r,L_r)$ (see subsection \[GE\]) $$L_r =\Phi \circ L \circ \Theta, \qquad \Psi_r=\frac{\Psi}{\Theta}$$ where $\Theta$ is a solution of $L\Theta=0$ whereas $\Phi$ is a solution of $L^{\dagger}\Phi=0$ where $L^{\dagger}$ denotes operator formally adjoint to operator $L$. We emphasize that this particular choice of functions $\Theta$ and $\Phi$ is essential, for it guarantees existence of function $\tilde{\Psi}_r$ in the next transformation

2.  transformations $(\Psi_r, L_r)\mapsto(\tilde{\Psi}_r,\tilde{L}_r)$ that take either the form (details are given in the text below) $$\label{f1}
    \begin{array}{l}
    \left[ 
    \begin{matrix}
    {\Delta_1}\tilde{\Psi}_r  \cr {\Delta_2}\tilde{\Psi}_r 
    \end{matrix}
    \right] = 
     \left[ 
    \begin{matrix}
     \alpha & \beta\cr
    \gamma& \delta
    \end{matrix}
    \right] 
    \left[ 
    \begin{matrix}
    {\Delta_1}\Psi_r \cr  {\Delta_2}\Psi_r 
    \end{matrix}
    \right]
    \end{array}$$ or the form $$\label{f2}
    \begin{array}{l}
    \left[ 
    \begin{matrix}
    {\Delta_1}\tilde{\Psi}_r  \cr {\Delta_2}\tilde{\Psi}_r 
    \end{matrix}
    \right] = 
     \left[ 
    \begin{matrix}
     \alpha & \beta\cr
    \gamma& \delta
    \end{matrix}
    \right] 
    \left[ 
    \begin{matrix}
    {\Delta_{-1}}\Psi_r \cr  {\Delta_{-2}}\Psi_r 
    \end{matrix}
    \right]
    \end{array}$$ where ${\Delta_1}$, ${\Delta_2}$ denote forward difference operators ${\Delta_1}f(m,n):= f(m+1,n)-f(m,n)$, ${\Delta_2}f(m,n):= f(m,n+1)-f(m,n)$, whereas ${\Delta_{-1}}$, ${\Delta_{-2}}$ denote backward difference operators ${\Delta_{-1}}f(m,n):= f(m-1,n)-f(m,n)$, ${\Delta_{-2}}f(m,n):= f(m,n-1)-f(m,n)$.

3.  gauge transformation $(\Psi_r,L_r)\mapsto (\tilde{\Psi}_r,\tilde{L}_r)$ $$L_r =R \circ L \circ S, \qquad \Psi_r=\frac{\Psi}{S}$$ where in general $R$ and $S$ are arbitrary functions (which should be specified when reductions or specifications of the transformations are considered).

(620,400) (90,360)[(215,90)]{} (-10,375)[(40,10)\[l\][6-point]{}]{} (80,390) (73,380)(87,380) (66,370)(80,370) (94,370) (120,375)[(40,10)\[l\][sect. \[Gc\]]{}]{} (-10,325)[(50,10)\[l\][$a \psi,_{xx}+b \psi,_{yy}+2 c \psi,_{xy} + g \psi,_x+ h \psi,_y=f \psi$]{}]{} (-10,350)[(50,10)\[l\][$A \Psi \! {_{(11)}}\! + \! B \Psi \! {_{(22)}}\! + \!2 C \Psi \! {_{(12)}}\!+\! G \Psi \! {_{(1)}}\!+\! H \Psi \! {_{(2)}}\! =  \! F \Psi$]{}]{}

(90,310)[(0,-1)[55]{}]{} (15,285)[(50,10)\[l\][Specification]{}]{} (0,270)[(50,10)\[l\][$A=0=B$, $C=\frac{1}{2}$ ]{}]{}

(20,250)[(0,-1)[80]{}]{} (160,250)[(0,-1)[80]{}]{} (20,250)[(1,0)[140]{}]{} (20,170)[(1,0)[140]{}]{} (30,225)[(40,10)\[l\][4-point]{}]{} (78,235)(88,235) (78,225)(88,225) (100,225)[(40,10)\[l\][sect. \[sec:4p-6p\], \[sec:fund\]]{}]{} (30,200)[(50,10)\[l\][$\Psi \! {_{(12)}}\!+\! G \Psi \! {_{(1)}}\!+\! H \Psi \! {_{(2)}}\! =  \! F \Psi$]{}]{} (30,175)[(50,10)\[l\][$\psi,_{xy} + g \psi,_x+ h \psi,_y=f \psi$]{}]{}

(90,165)[(0,-1)[55]{}]{} (15,140)[(50,10)\[l\][Specification]{}]{} (30,125)[(50,10)\[l\][$F=0$]{}]{}

(30,105)[(0,-1)[70]{}]{} (150,105)[(0,-1)[70]{}]{} (30,105)[(1,0)[120]{}]{} (30,35)[(1,0)[120]{}]{} (40,85)[(40,10)\[l\][3-point]{}]{} (78,95)(92,95) (85,85) (100,85)[(40,10)\[l\][sect. \[3MN\], \[3p-sch\]]{}]{} (40,65)[(50,10)\[l\][$\Psi \! {_{(12)}}\!+\! G \Psi \! {_{(1)}}\!+\! H \Psi \! {_{(2)}}\! = 0$ ]{} ]{} (40,45)[(50,10)\[l\][$ g \psi,_x+ h \psi,_y= \psi$]{}]{}

(165,240)[(40,10)\[l\][Ribaucour reduction]{}]{} (165,230)[(1,0)[90]{}]{} (255,230)[(0,-1)[15]{}]{}

(210,210)[(0,-1)[60]{}]{} (300,210)[(0,-1)[60]{}]{} (210,210)[(1,0)[90]{}]{} (210,150)[(1,0)[90]{}]{} (220,195)[(40,10)\[l\][See section \[Qr\]]{}]{} (215,175)[(40,10)\[l\][$\Psi$ lies on a quadric]{}]{} (215,160)[(40,10)\[l\][$\psi$ lies on a quadric]{}]{}

(145,130)[(0,1)[35]{}]{} (425,130)[(-1,0)[280]{}]{} (425,130)[(0,1)[15]{}]{} (330,135)[(50,10)\[l\][Moutard reduction]{}]{}

(330,210)[(0,-1)[60]{}]{} (423,210)[(0,-1)[60]{}]{} (330,210)[(1,0)[93]{}]{} (330,150)[(1,0)[93]{}]{} (340,195)[(50,10)\[l\][discrete Moutard]{}]{} (345,185)[(40,10)\[l\][sect. \[Mr\], \[bkp\]]{}]{} (332,170)[(50,10)\[l\][$\Psi \! {_{(12)}}\!+\! \Psi \!= \! M(\Psi \! {_{(1)}}\!+\!  \Psi \! {_{(2)}}) $]{}]{} (350,155)[(50,10)\[l\][$\psi,_{xy} =f \psi$]{}]{}

(423,180)[(1,0)[4]{}]{} (427,210)[(0,-1)[60]{}]{} (527,210)[(0,-1)[60]{}]{} (427,210)[(1,0)[100]{}]{} (427,150)[(1,0)[100]{}]{} (440,195)[(50,10)\[l\][adjoint discrete]{}]{} (478,185)[(50,10)\[l\][sect. \[Mr\]]{}]{} (430,185)[(50,10)\[l\][Moutard,]{}]{} (428,170)[(50,10)\[l\][$\Psi \! {_{(12)}}\!\!+\!\! \Psi \!= \! M \!{_{(1)}}\! \Psi \! {_{(1)}}\! \!+ \! M  \!{_{(2)}}\!  \Psi \! {_{(2)}}$]{}]{} (447,155)[(50,10)\[l\][$\psi,_{xy} =f \psi$]{}]{}

(335,365)[(135,70)]{} (273,385)[ (60,10)\[l\][7-point]{}]{} (268,375)[ (60,10)\[l\][selfadjoint]{}]{} (328,395)(342,395) (321,385)(335,385) (349,385) (328,375)(342,375) (355,380)[(40,10)\[l\][sect. \[7ps\]]{}]{} (270,360)[(100,10)\[l\][${\mathcal A}\! {_{(1)}}\! \Psi\! {_{(1)}}\! +\! {\mathcal A}\!  \Psi{_{(-1)}}\! +\! 
   {\mathcal B}\! {_{(2)}}\! \Psi\! {_{(2)}}\! +\!{\mathcal B} \! \Psi{_{(-2)}}\! 
   +$]{}]{} (270,350)[(100,10)\[l\][$  {\mathcal C}\! {_{(1)}}\! \Psi{_{(1-2)}}\! +\! {\mathcal C}\! {_{(2)}}\! \Psi\! {_{(-12)}}\! =\! {\mathcal F} \! \Psi$]{}]{}

(272,335)[(100,10)\[l\][$ (a \psi,_x\!+c \psi,_y\!\!),_x \!+(b \psi,_y\! +c\psi,_x\!\!),_y\!=\!f \psi$]{}]{}

(475,365)[(135,70)]{} (415,380)[(50,10)\[l\][Honeycomb]{}]{} (483,395)(497,395) (476,385) (504,385) (483,375)(497,375)

(410,355)[(50,10)\[l\][See section \[hcb\] for details]{}]{}

(405,215)[(0,1)[115]{}]{} (405,330)[(-1,0)[10]{}]{} (260,315)[(50,10)\[l\][Specification ${\mathcal C}=0$]{}]{} (405,330)[(1,0)[10]{}]{} (350,325)[(0,-1)[10]{}]{} (410,275)[(50,10)\[l\][Sublattice approach]{}]{} (410,265)[(50,10)\[l\][on quasi-regular]{}]{} (410,255)[(50,10)\[l\][ rhombic lattice]{}]{}

(380,215)[(0,1)[20]{}]{} (290,225)[(50,10)\[l\][Sublattice approach]{}]{} (290,215)[(50,10)\[l\][on square lattice]{}]{}

(335,275)[(130,70)]{} (276,295)[ (60,10)\[l\][5-point]{}]{} (271,285)[ (60,10)\[l\][selfadjoint]{}]{} (334,305) (324,295) (334,295)(344,295) (334,285) (352,290)[(40,10)\[l\][sect. \[5r\]]{}]{} (290,270)[(100,10)\[l\][${\mathcal A}\!{_{(1)}}\!\Psi\!{_{(1)}}\!+\!{\mathcal A} \!\Psi\!{_{(-1)}}+$]{}]{} (290,260)[(100,10)\[l\][$
   {\mathcal B}\!{_{(2)}}\Psi{_{(2)}}\!+\!{\mathcal B}\! \Psi\!{_{(-2)}}\!
    =\!{\mathcal F}\! \Psi$]{}]{}

(290,245)[(100,10)\[l\][$(a \psi,_x),_x +(b \psi,_y),_y=f \psi$]{}]{} (125,120)[(0,1)[45]{}]{} (595,120)[(-1,0)[470]{}]{} (595,120)[(0,1)[25]{}]{} (510,125)[(50,10)\[l\][Goursat reduction]{}]{}

(550,235)[(0,-1)[85]{}]{} (640,235)[(0,-1)[85]{}]{} (550,235)[(1,0)[90]{}]{} (550,150)[(1,0)[90]{}]{} (557,225)[(50,10)\[l\][discrete Goursat]{}]{} (560,215)[(40,10)\[l\][sect. \[Gr\], \[ckp\]]{}]{}

(551,200)[(100,10)\[l\][ $\Psi\! {_{(12)}}\!=\!\frac{ {q}{_{(2)}}}{ {q}} \Psi \!{_{(1)}}\!+ \!\Psi \!{_{(2)}}- $]{}]{} (551,185)[(100,10)\[l\][ $\frac{\tau {_{(12)}}\tau}{\tau {_{(1)}}\tau {_{(2)}}} \frac{ {q}{_{(2)}}}{ {q}} \Psi$]{}]{} (551,170)[(100,10)\[l\][ ${q}^2:= \frac{\tau {_{(1)}}\tau{_{(2)}}- \tau {_{(12)}}\tau }{\tau^2 {_{(1)}}}$]{}]{} (551,150)[(100,10)\[l\][ $\psi,_{xy} = \frac{p,_y}{p} \psi,_x+p^2\psi$]{}]{}

We pay special attention to subclasses of operators $L$ that admit Darboux transformations. A trivial, from the point of view of integrable systems, examples are subclasses obtained by fixing a gauge (see section \[GE\]). Very important classes of operators are obtained by imposing conditions on the mentioned above functions $\Theta$ and $\Phi$. For example in the case of self-adjoint 7-point scheme we discuss in section \[7ps\], the relation is $\Phi=\Theta$. We refer to such procedures as to reductions. The further examples of reductions are given in sections \[Gr\], \[Mr\] and \[sec:4p-red\]. We reserve separate name “specification” to the cases when the operator is reduced, whereas the transformations remain essentially unaltered (i.e. no constraint on functions $\Theta$ and $\Phi$ is necessary). Examples of specifications are given in sections \[sec:4p-6p\], \[5r\] and \[3MN\].

With the notion of “specification” another important issue appears, since continuous counterparts of the specifications presented here can be viewed as a choice of particular gauge and independent variables. The interesting (because not well understood) aspects of discrete integrable systems come from the fact there is no theory how to change discrete independent variables so that not to destroy underlying integrable phenomena. Sublattice approach, widely used in theoretical physics, can be regarded, to some extent, as counterpart of change of independent variables. To some extent, because the only case studied in detail from the point of view of Darboux transformations is the case of self-adjoint equation , which in the discrete case consists of Moutard case (section \[Mr\]), self-adjoint case (section \[sec:s-a\]) and their mutual relations (c.f. [@DGNS-4-5; @DNS-4-7]). Interrelations between results presented here are summarized in Figure 1.

Large part of the paper is dedicated to the specification in which the matrix in is diagonal. In this case one can consider systems of four point operators defined on lattices with arbitrary number of independent variables. The multidimensional lattices are extensively discussed in sections \[sec:4p\] and \[sec:4p-red\], where compact elegant expressions for superpositions of fundamental transformations are presented either.

Two dimensional systems {#sec:6p}
=======================

In this section we present discretizations of equation and its subclasses covariant under Darboux transformations.

General case {#Gc}
------------

Out of the schemes that can serve as a discretization of the 2D equation $$\label{2Dgeneral} 
a \psi,_{xx}+b \psi,_{yy}+2 c \psi,_{xy}
+ g \psi,_x+ h \psi,_y=f \psi$$ the following 6-point scheme $$\label{6}
A \Psi {_{(11)}}+B \Psi {_{(22)}}+2 C \Psi {_{(12)}}+G \Psi {_{(1)}}+H \Psi {_{(2)}}=  F \Psi$$ deserve a special attention (coefficients $A$, $B$ etc. and dependent variable $\Psi$ are functions on ${\mathbb Z}^2$, subscript in brackets denote shift operators, $f(m,n){_{(1)}}=f(m+1,n)$, $f(m,n){_{(2)}}=f(m,n+1)$, $f(m,n){_{(11)}}=f(m+2,n)$, $f(m,n){_{(22)}}=f(m,n+2)$ and $f(m,n){_{(12)}}=f(m+1,n+1)$). The scheme admits decomposition $[(\alpha_1 T_1+\alpha_2 T_2+\alpha_3) 
(\beta_1 T_1+\beta_2 T_2+\beta_3)+\gamma]\psi=0$ (where $T_1$ and $T_2$ are forward shift operators in first and second direction respectively) and therefore its Laplace transformations can be constructed (see [@Laplace; @DarbouxIV; @Weiss; @DCN; @NovVes; @Adler-l; @Nieszporski-ll] and Section \[sec:Laplace-tr\] for the notion of the Laplace transformations of multidimensional linear operators). What more important from the point of view of this review the scheme is covariant under a fundamental Darboux transformation [@Nieszporski-6p] (the transformation is often referred to as binary Darboux transformation in soliton literature). Indeed,

\[th:6p-Dt\] Given a non-vanishing solution $\Theta$ of $$\begin{array}{l}
\label{6t}
 A \Theta{_{(11)}}+B \Theta{_{(22)}}+2 C \Theta{_{(12)}}+G \Theta{_{(1)}}+H \Theta{_{(2)}}=  F \Theta
\end{array}$$ and a non-vanishing solution $\Phi$ of the equation adjoint to equation $$\begin{array}{l}
\label{6ad}
A{_{(-1-1)}}\Phi{_{(-1-1)}}+B{_{(-2-2)}}\Phi{_{(-2-2)}}+2 C{_{(-1-2)}}\Phi{_{(-1-2)}}+ G{_{(-1)}}\Phi{_{(-1)}}+H{_{(-2)}}\Phi{_{(-2)}}=F \Phi
\end{array}$$ (negative integers in brackets in subscript denote backward shift e.g. $f(m,n){_{(-1)}}=f(m-1,n)$ etc.) the existence of auxiliary function $P$ is guaranteed $$\label{P}
\begin{array}{l}
{\Delta_{-1}}(\Phi\Theta {_{(12)}}P) = 
B{_{(-2)}}\Phi{_{(-2)}}\Theta{_{(2)}}+B \Phi \Theta{_{(22)}}+
C{_{(-1)}}\Phi{_{(-1)}}\Theta{_{(2)}}+C \Phi \Theta{_{(12)}}+H \Phi \Theta{_{(2)}}\\
{\Delta_{-2}}(\Phi\Theta {_{(12)}}P) =
-(A{_{(-1)}}\Phi{_{(-1)}}\Theta{_{(1)}}+A \Phi \Theta{_{(11)}}+
C{_{(-2)}}\Phi{_{(-2)}}\Theta{_{(1)}}+C \Phi \Theta{_{(12)}}+G \Phi \Theta{_{(1)}}).
\end{array}$$ Then equation can be rewritten as $$\begin{array}{l}  {\Delta_2}\left[ (C{_{(-2)}}+ P{_{(-2)}})\Phi{_{(-2)}}\Theta{_{(1)}}{\Delta_1}\left(\frac{\Psi}{\Theta} \right) +
B{_{(-2)}}\Phi{_{(-2)}}\Theta{_{(2)}}{\Delta_2}\left(\frac{\Psi}{\Theta} \right) \right]+\\
{\Delta_1}\left[ A{_{(-1)}}\Phi{_{(-1)}}\Theta{_{(1)}}{\Delta_1}\left(\frac{\Psi}{\Theta} \right) -(P{_{(-1)}}- C{_{(-1)}})\Phi{_{(-1)}}\Theta{_{(2)}})
{\Delta_2}\left(\frac{\Psi}{\Theta} \right) \right]=0
\end{array}$$ which in turn guarantees the existence of functions $\Psi'$ such that $$\label{DT6}
\begin{array}{l}
\left[ 
\begin{matrix}
{\Delta_1}(\Psi') \cr {\Delta_2}(\Psi')
\end{matrix}
\right] = 
 \left[ 
\begin{matrix}
 (C{_{(-2)}}+ P{_{(-2)}})\Phi{_{(-2)}}\Theta{_{(1)}}& B{_{(-2)}}\Phi{_{(-2)}}\Theta{_{(2)}}\cr
-A{_{(-1)}}\Phi{_{(-1)}}\Theta{_{(1)}}&(P{_{(-1)}}- C{_{(-1)}})\Phi{_{(-1)}}\Theta{_{(2)}}\end{matrix}
\right] 
\left[ 
\begin{matrix}
{\Delta_1}\left(\frac{\Psi}{\Theta} \right)\cr  {\Delta_2}\left(\frac{\Psi}{\Theta} \right)
\end{matrix}
\right]
\end{array}$$ Assuming that matrix in is invertible on the whole lattice, i.e. everywhere and finally on introducing $\tilde{\Psi}$ via $$\label{S}
\tilde{\Psi}=\frac{\Psi'}{S}$$ and taking the opportunity of multiplying resulting equation by non-vanishing function $R$, we arrive at the conclusion that the function $\tilde{\Psi}$ satisfies equation of the form but with new coefficients $$\label{dcoef}
\begin{array}{l}
\tilde{A}=\frac{R S{_{(11)}}\Phi \Theta{_{(11)}}}{D{_{(1)}}} A,
\qquad
\tilde{B}=\frac{R S{_{(22)}}\Phi \Theta{_{(22)}}}{D{_{(2)}}} B,
\qquad
\tilde{F}=\frac{R S \Phi \Theta }{D} F,
\\ 
\\
\tilde{C} =\frac{R S{_{(12)}}}{2}\left(
\frac{\Theta{_{(11)}}\Phi{_{(1-2)}}(C+P){_{(1-2)}}}{D{_{(1)}}}+\frac{\Theta{_{(22)}}\Phi{_{(-12)}}(C-P){_{(-12)}}}{D{_{(2)}}}\right), 
\\
\\
{\tilde{G}}=-{R S{_{(1)}}}\left( \frac{\Theta{_{(11)}}\Phi{_{(1-2)}}(C+P){_{(1-2)}}+\Theta{_{(11)}}\Phi A}{D{_{(1)}}}+
                                 \frac{\Theta{_{(2)}}\Phi{_{(-1)}}(C-P){_{(-1)}}+\Theta{_{(1)}}\Phi{_{(-1)}}A{_{(-1)}}}{D} \right),
\\
\\
{\tilde{H}}=-{R S{_{(2)}}}\left( \frac{\Theta{_{(22)}}\Phi{_{(-12)}}(C-P){_{(-12)}}+\Theta{_{(22)}}\Phi B}{D{_{(2)}}}
                                +\frac{\Theta{_{(1)}}\Phi{_{(-2)}}(C+P){_{(-2)}}+\Theta{_{(2)}}\Phi{_{(-2)}}B{_{(-2)}}}{D}\right),
\end{array}$$

The family of maps $\Psi \mapsto \tilde{\Psi}$ given by - we refer to as Darboux transformations of equation .

Gauge equivalence {#GE}
-----------------

We say that two linear operators $L$ and $L'$ are gauge equivalent if one can find functions say $\Phi$ and $\Theta$ such that $L'=\Phi \circ L \circ \Theta$ (where $\circ$ stands for composition of operators). The idea to consider equivalence classes of linear operators with respect to gauge rather than single operator itself, goes to Laplace and Darboux papers [@Laplace; @DarbouxIV]. In the continuous case it reflects in the fact that one can confine himself e.g. to equations (the so called affine gauge) $$a \psi,_{xx}+b \psi,_{yy}+ \psi,_{xy}
+ g \psi,_x+ h \psi,_y=0$$ or to (in this case we would like to introduce the name [*basic gauge*]{}) $$a \psi,_{xx}+b \psi,_{yy}+ 2c\psi,_{xy}
+ k,_y \psi,_x -k,_x \psi,_y=0$$ without loss of generality.

The Darboux transformation can be viewed as transformation acting on equivalence classes (with respect to gauge) of equation (compare [@Nieszporski-6p]). Therefore one can confine himself to particular elements of equivalence class. Commonly used choice is to confine oneself to the affine gauge i.e. to equations that obeys $$A + B  + 2C + G + H - F  = 0.$$ If one puts S = const in (this condition is not necessary) then the above constraint is preserved under the Darboux transformation. One can consider further specification of the gauge $$A + B  + 2C + G + H - F  = 0, \quad A{_{(-12)}}+ B{_{(1-2)}}+ 2C  + G {_{(2)}}+ H {_{(1)}}- F {_{(12)}}= 0$$ This choice of gauge we would like to refer as to basic gauge of equation . Note that if the equation is in a basic gauge its formal adjoint is in a basic gauge too.

Specification to 4-point scheme and its reductions {#sec:4p-6p}
--------------------------------------------------

In the continuous case due to possibility of changing independent variables one can reduce, provided equation is hyperbolic, to canonical form $$\psi,_{xy}+ g \psi,_x+ h \psi,_y=f \psi.$$ In the discrete case similar result can be obtained in a different way. Taking a glance at (\[dcoef\]) we can notice that coefficients $A$, $B$ and $F$ transform in a very simple manner. In particular, if any of these coefficients equals zero then its transform equals zero too. Let us stress that in this case if we do so, we do not impose any constraints on transformation data, transformations remains essentially the same. First we shall concentrate on the case when two out of three mentioned functions vanish.

If we put $$\label{4ps}
A=0, \quad B=0, \quad C=\frac{1}{2}$$ then from equation $$\tilde{A}=0, \quad \tilde{B}=0$$ and one can adjust the function $R$ so that $$\tilde{C}=\frac{1}{2}$$ so we arrive at the 4-point scheme $$\label{4pabc}
\Psi {_{(12)}}=\alpha \Psi {_{(1)}}+\beta \Psi {_{(2)}}+\gamma \Psi.$$ and its fundamental transformation c.f. [@BoKo; @MQL; @MDS; @TQL].

We observe that the form of equation is covariant under the gauge $L\mapsto \frac{1}{g {_{(12)}}} L g$ and to identify whether two equations are equivalent or not we use the invariants of the gauge $$\kappa=\frac{\alpha{_{(2)}}\beta }{\gamma{_{(2)}}}, \quad n=\frac{\alpha \beta {_{(1)}}}{\gamma{_{(1)}}}$$ Two equations are equivalent if their corresponding invariants $\kappa$ and $n$ are equal [@Nieszporski-ll].

### Goursat equation {#Gr}

In this subsection we discuss the discretization of class of equations $$\psi,_{xy} = \frac{p,_x}{p} \psi,_y+p^2\psi$$ which is referred to as Goursat equation. The discrete counterpart of Goursat equation arose from the surveys on Egorov lattices [@Schief-talk] and symmetric lattices [@DS-sym] and can be written in the form [@Nieszporski-ll] $$\label{dG}
\Psi {_{(12)}}=\frac{ {q}{_{(2)}}}{ {q}} \Psi {_{(1)}}+ \Psi {_{(2)}}-\frac{\tau {_{(12)}}\tau}{\tau {_{(1)}}\tau {_{(2)}}} \frac{ {q}{_{(2)}}}{ {q}} \Psi$$ where functions $ {q}$ and $\tau$ are related via $${q}^2= \frac{\tau {_{(1)}}\tau{_{(2)}}- \tau {_{(12)}}\tau }{\tau^2 {_{(1)}}}.$$ The gauge invariant characterization of the discrete Goursat equation is either $$n{_{(2)}}^2=\kappa \kappa{_{(12)}}\frac{(1+\kappa {_{(1)}})(1+\kappa {_{(2)}})}{(1+\kappa )(1+\kappa {_{(12)}})}$$ or $$\kappa{_{(1)}}^2=n n{_{(12)}}\frac{(1+n {_{(1)}})(1+n {_{(2)}})}{(1+n )(1+n {_{(12)}})}$$

The Goursat equation can be isolated from the others 4-point schemes in the similar way that Goursat did it over hundred years ago [@Goursat; @Nieszporski-ll] i.e. as the equation such that one of its Laplace transformations maps solutions of the equation to solutions of the adjoint equation. Therefore in this case if $\Theta$ obeys $$\label{dG1}
\Theta {_{(12)}}=\frac{ {q}{_{(2)}}}{ {q}} \Theta {_{(1)}}+ \Theta {_{(2)}}-\frac{\tau {_{(12)}}\tau}{\tau {_{(1)}}\tau {_{(2)}}} \frac{ {q}{_{(2)}}}{ {q}} \Theta$$ then its Laplace transformation $$\label{fGt}
\Phi=\frac{1}{{q}^2{_{(2)}}}\frac{\tau{_{(2)}}}{\tau {_{(12)}}} {\Delta_1}\Theta {_{(2)}}$$ is a solution of equation adjoint to equation [@Nieszporski-ll]. Equation is the constraint we impose on transformation data (i.e. functions $\Theta$ and $\Phi$) in the fundamental transformation , (we recall we have already put $A$ and $B$ equal to zero). In addition if $\Theta$ obeys then $${\Delta_1}\left(\frac{\tau}{\tau {_{(2)}}}\Theta^2 \right)={\Delta_2}\left[\frac{1}{{q}^2}\frac{\tau}{\tau {_{(1)}}}({\Delta_1}\Theta)^2 \right]$$ and as a result there exists function ${\vartheta}$ such that $${\Delta_1}{\vartheta}=\frac{1}{{q}^2}\frac{\tau}{\tau {_{(1)}}}({\Delta_1}\Theta)^2, \quad{\Delta_2}{\vartheta}=\frac{\tau}{\tau {_{(2)}}}\Theta^2$$ Now it can be shown that one can put $$\Phi \Theta {_{(12)}}(P-\frac{1}{2})=-{\vartheta}{_{(12)}}, \quad \Phi \Theta {_{(12)}}(P+\frac{1}{2})= 
\left[\sqrt{\frac{\tau {_{(2)}}}{\tau {_{(1)}}}}\frac{\sqrt{{\Delta_1}{\vartheta}{\Delta_2}{\vartheta}}}{q} \right] {_{(2)}}-{\vartheta}{_{(2)}}$$ and the transformation takes form $$\label{GDT}
\begin{array}{l}
\left[ 
\begin{matrix}
{\Delta_1}(\sqrt{\frac{{\tau }}{{\tau{_{(2)}}}}} \frac{{\vartheta}}{\sqrt{{\Delta_2}{\vartheta}}}\tilde{\Psi}) \cr 
{\Delta_2}(\sqrt{\frac{{\tau }}{{\tau {_{(2)}}}}} \frac{{\vartheta}}{\sqrt{{\Delta_2}{\vartheta}}}\tilde{\Psi})
\end{matrix}
\right] = 
 \left[ 
\begin{matrix}
  \sqrt{\frac{\tau {_{(2)}}}{\tau {_{(1)}}}}\frac{\sqrt{{\Delta_1}{\vartheta}{\Delta_2}{\vartheta}}}{q} -{\vartheta}& 0 \cr
0 & - {\vartheta}{_{(2)}}\end{matrix}
\right] 
\left[ 
\begin{matrix}
{\Delta_1}\left(\sqrt{\frac{{\tau}}{{\tau {_{(2)}}}}}  \frac{\Psi}{\sqrt{{\Delta_2}{\vartheta}}}\right)\cr  
{\Delta_2}\left(\sqrt{\frac{{\tau}}{{\tau {_{(2)}}}}}  \frac{\Psi}{\sqrt{{\Delta_2}{\vartheta}}} \right)
\end{matrix}
\right]
\end{array}$$ which is the discrete version of Goursat transformation [@Goursat]. The transformation rule for the field $\tau$ is $$\tilde\tau=\tau {\vartheta}.$$

### Moutard equation {#Mr}

In this subsection we discuss discretization Moutard equation and its (Moutard) transformation [@Moutard] $$\psi,_{xy} = f \psi.$$ Moutard equation is self-adjoint equation (which allowed us to impose reduction $\Phi=\Theta$ in the continuous analogue of transformation c.f. [@Nieszporski-6p]). The point is that opposite to the continuous case there is no appropriate self-adjoint 4-point scheme

We do have the reduction of fundamental transformation that can be regarded as discrete counterpart of Moutard transformation. Namely, class of equations that can be written in the form $$\label{dM}
\Psi {_{(12)}}+ \Psi=M(\Psi {_{(1)}}+ \Psi {_{(2)}})$$ we refer nowadays to as discrete Moutard equation. It appeared in the context of integrable systems in [@DJM] and then its Moutard transformations have been studied in detail in [@NiSchief]. Gauge invariant characterization of the class of discrete Moutard equations is [@Nieszporski-ll] $$n{_{(2)}}n=\kappa {_{(1)}}\kappa$$ Let us trace this reduction on the level of fundamental transformation. Putting $2C=-F$ and $G=H=:M F$ the equations , and take respectively form $$\begin{array}{l}
\label{6M}
F (\Psi {_{(12)}}+\Psi)= G (\Psi{_{(1)}}+ \Psi{_{(2)}})      
\end{array}$$

$$\begin{array}{l}
\label{6tM}
F(\Theta{_{(12)}}+\Theta)= G (\Theta{_{(1)}}+ \Theta{_{(2)}}) 
\end{array}$$

$$\begin{array}{l}
\label{6adM}
F{_{(-1-2)}}\Phi{_{(-1-2)}}+F \Phi= G{_{(-1)}}\Phi{_{(-1)}}+ G{_{(-2)}}\Phi {_{(-2)}}\end{array}$$

The crucial observation is: if the function $\Theta$ satisfies equation then the function $\Phi$ given by $$\label{PT}
\begin{array}{l}
\Phi=\frac{1}{F} (\Theta{_{(1)}}+ \Theta{_{(2)}})
\end{array}$$ satisfies equation [@Nieszporski-ll]. If we put $$2 P=\frac{\Theta{_{(1)}}- \Theta{_{(2)}}}{\Phi}$$ then equations will be automatically satisfied. If in addition we put in $S=\Theta$ then Darboux transformation takes form (c.f. [@NiSchief]) $$\label{MDT}
\left[ 
\begin{matrix}
{\Delta_1}(\Theta \tilde{\Psi}) \cr {\Delta_2}(\Theta \tilde{\Psi})
\end{matrix}
\right]
= \left[ 
\begin{matrix} \Theta \Theta {_{(1)}}& 0\cr
 0  & -\Theta \Theta {_{(2)}}\end{matrix}
\right]
\left[ 
\begin{matrix}  {\Delta_1}\frac{\Psi}{\Theta}\cr {\Delta_2}\frac{\Psi}{\Theta}
\end{matrix}
\right]$$ and it serves as transformation $\Psi \mapsto \tilde{\Psi}$ that maps solutions of equation into solutions of another discrete Moutard equation $$\tilde{\Psi}{_{(12)}}+ \tilde{\Psi}=\widetilde{M}( \tilde{\Psi}{_{(1)}}+ \tilde{\Psi} {_{(2)}}), \quad 
\widetilde{M}=\frac{\Theta {_{(1)}}\Theta {_{(2)}}}{\Theta {_{(12)}}\Theta}M
$$ We recall that $\Theta$ an arbitrary non-vanishing fixed solution of the equation $$\label{Ms}
 \Theta {_{(12)}}+ \Theta=M(\Theta {_{(1)}}+ \Theta {_{(2)}})$$

Every equation that is gauge equivalent to we refer to as an adjoint discrete equation Moutard equation, its gauge invariant characterization is [@Nieszporski-ll] $$\kappa {_{(12)}}\kappa {_{(1)}}= n{_{(12)}}n {_{(2)}}.$$ Equation can be regarded as potential version of equation

Self-adjoint case {#sec:s-a}
-----------------

In the continuous case there is direct reduction of the fundamental transformation for equation to the Moutard type transformation for self-adjoint equation [@Nieszporski-6p] $$\label{sag}
(a \psi,_x+c \psi,_y),_x +(b \psi,_y +c\psi,_x),_y=f \psi.$$ In this Section we present difference analogues of equation which allow for the Darboux transformation.Opposite to the continuous case the transformation will not be direct reduction of the fundamental transformation for the 6-point scheme . The self-adjoint discrete operators studied below are however intimately related to discrete Moutard equation which provides the link between their Darboux transformations and the fundamental transformation for equation .

### 7-point self-adjoint scheme {#7ps}

The following 7-point linear system $$\label{7eq}
{\mathcal A}{_{(1)}}\Psi{_{(1)}}+{\mathcal A} \Psi{_{(-1)}}+
   {\mathcal B}{_{(2)}}\Psi{_{(2)}}+{\mathcal B} \Psi{_{(-2)}}+  {\mathcal C}{_{(1)}}\Psi{_{(1-2)}}+{\mathcal C}{_{(2)}}\Psi{_{(-12)}}={\mathcal F} \Psi,$$ allows for the Darboux transformation [@NSD].

\[th:tp-D\] Given scalar solution $\varTheta$ of the linear equation , then $\tilde\Psi$ given as solution of the following system $$\begin{pmatrix}
{\Delta_1}(\tilde{\Psi}  \varTheta) \\ {\Delta_2}( \tilde{\Psi} \varTheta) 
\end{pmatrix} =
\begin{pmatrix}
 {\mathcal C} \varTheta{_{(-1)}}\varTheta{_{(-2)}}&
-\varTheta{_{(-2)}}({\mathcal B} \varTheta  +
{\mathcal C} \varTheta{_{(-1)}}) \\
\varTheta{_{(-1)}}\left( {\mathcal A} \varTheta  \!+\! {\mathcal C}
 \varTheta{_{(-2)}}\right) &
- {\mathcal C} \varTheta{_{(-1)}}\varTheta{_{(-2)}}\end{pmatrix}
\begin{pmatrix}
{\Delta_{-1}}\frac{\Psi}{\varTheta}
\\
{\Delta_{-2}}\frac{\Psi}{\varTheta} 
 \end{pmatrix}$$ satisfies the 7-point scheme with the new fields given by $$\tilde{{\mathcal A}}{_{(1)}}=\frac{\varTheta \varTheta{_{(1)}}{\mathcal A}}{\varTheta{_{(-2)}}{\mathcal P}},  
\qquad
\tilde{{\mathcal B}}{_{(2)}}=\frac{\varTheta \varTheta{_{(2)}}{\mathcal B}}{\varTheta{_{(-1)}}{\mathcal P}},
, \qquad
\tilde{{\mathcal C}}=\frac{{\mathcal C}{_{(-1-2)}}\varTheta{_{(-1)}}\varTheta{_{(-2)}}}
{\varTheta{_{(-1-2)}}{\mathcal P}{_{(-1-2)}}},$$ $$\tilde{{F}}=\varTheta \left(\frac{\tilde{{\mathcal A}}{_{(1)}}}{\varTheta{_{(1)}}}+   
\frac{\tilde{{\mathcal A}} }{\varTheta{_{(-1)}}}+  
\frac{\tilde{{\mathcal B}}{_{(2)}}}{\varTheta{_{(2)}}} + 
 \frac{\tilde{{\mathcal B}}}{\varTheta{_{(-2)}}}+
 \frac{\tilde{{\mathcal C}}{_{(1)}}}{\varTheta{_{(1-2)}}}+
\frac{\tilde{{\mathcal C}}{_{(2)}}}{\varTheta{_{(-12)}}}\right)$$ where $${\mathcal P}= \varTheta {\mathcal A}{\mathcal B} +
\varTheta{_{(-1)}}{\mathcal C} {\mathcal A}+\varTheta{_{(-2)}}{\mathcal C}  {\mathcal B}.$$

![The bipartite quasi-regular rhombic tiling and the triangular and honeycomb lattices[]{data-label="fig:QRRT"}](QRRT.eps "fig:"){width="7cm"} ![The bipartite quasi-regular rhombic tiling and the triangular and honeycomb lattices[]{data-label="fig:QRRT"}](T-H-lattices1.eps "fig:"){width="7cm"}

As it was shown in [@DNS-4-7] the self-adjoint 7-point scheme can be obtained from the system of Moutard equations imposed consistently on quadrilaterals of the bipartite quasi-regular rhombic tiling (see Fig. \[fig:QRRT\]), which is a particular case of the approach considered in [@BobMerSur]. Then the Moutard transformations can be also restricted to the triangular sublattice leading to Theorem \[th:tp-D\].

### Specification to 5-point scheme {#5r}

The 7-point scheme admits specification ${\mathcal C}=0$ (alternatively one can put ${\mathcal A}=0$ or ${\mathcal B}=0$) and as result we obtain specification to 5-point self-adjoint scheme [@NSD]. $$\label{5eq}
 {\mathcal A}{_{(1)}}\Psi{_{(1)}}+{\mathcal A} \Psi{_{(-1)}}+
   {\mathcal B}{_{(2)}}\Psi{_{(2)}}+{\mathcal B} \Psi{_{(-2)}}={\mathcal F} \Psi .$$ The self-adjoint 5-point scheme and its Darboux transformation can be also obtained form the Moutard equation on the (bipartite) square lattice [@DGNS-4-5].

### The honeycomb lattice {#hcb}

It is well known that the triangular and honeycomb grids are dual to each other (see Fig. \[fig:QRRT\]). Restriction of the system of the Moutard equations on the rhombic tiling to the honeycomb sublattice gives [@DNS-4-7] the following linear system $$\begin{aligned}
 \label{eq:hcl1}
\frac{1}{{\mathcal A}}(\Psi^+ - \Psi^-) + 
\frac{1}{{\mathcal B}}(\Psi^+_{(-1 2)} - \Psi^-) +
\frac{1}{{\mathcal C}_{(1 2)}}(\Psi^+_{(2)} - \Psi^-) & = 0,\\
\label{eq:hcl2}
\frac{1}{{\mathcal A}}(\Psi^- - \Psi^+) + 
\frac{1}{{\mathcal B}_{(1 -2)}}(\Psi^-_{(1 -2)} - \Psi^+) +
\frac{1}{{\mathcal C}_{(1)}}(\Psi^-_{(-2)} - \Psi^+) & = 0.\end{aligned}$$

Because $\Psi^-$ and $\Psi^+$ satisfy separately the self-adjoint 7-point schemes , but with different coefficients, then the linear problem - can be considered as a relation between two equations . This is the Laplace transformation between self-adjoint 7-point schemes studied in [@NovVes; @NovDyn].

The corresponding restriction of the Moutard transformation gives the Darboux transformation for the honeycomb linear problem.

\[th:DT-hc\] Given scalar solution $(\varTheta^+,\varTheta^-)$ of the honeycomb linear system - then the solution $(\tilde\Psi^+,\tilde\Psi^-)$ of the system $$\begin{aligned}
\label{eq:DT-hc1}
\tilde\Psi^+_{(2)} - \tilde\Psi^- & = \frac{{\mathcal C}}{{\mathcal R}}\left( 
\varTheta^-_{(-2)}\Psi^-_{(-1)} - \theta^-_{(-1)}\Psi^-_{(-2)} \right),\\
\label{eq:DT-hc2}
\tilde\Psi^+ - \tilde\Psi^- & = \frac{{\mathcal A}_{(-1)}}{{\mathcal R}}\left(
\varTheta^-_{(-1-2)}\Psi^-_{(-2)} -\varTheta^-_{(-2)}\Psi^-_{(-1-2)} \right), \\
\label{eq:DT-hc3}
\tilde\Psi^- - \tilde\Psi^+_{(-1 2)} & = \frac{{\mathcal B}_{(-2)}}{{\mathcal R}}\left(
\varTheta^-_{(-1-2)}\Psi^-_{(-1)} -\varTheta^-_{(-1)}\Psi^-_{(-1-2)} \right),\end{aligned}$$ with ${\mathcal R}$ given by $${\mathcal R}= {\mathcal A}_{(-1)}{\mathcal B}_{(-2)} +
{\mathcal C} {\mathcal A}_{(-1)}+ 
{\mathcal C}  {\mathcal B}_{(-2)},$$ satisfies a new honeycomb linear system with the coefficients $$\label{eq:coeff-DT-hc}
\tilde{\mathcal A} = \frac{{\mathcal A}_{(-1)}}{{\mathcal R}} 
\varTheta^{-}_{(-2)}\varTheta^{-}_{(-1-2)}, \qquad
\tilde{\mathcal B} = \frac{{\mathcal B}_{(-2)}}{{\mathcal R}} 
\varTheta^{-}_{(-1)}\varTheta^{-}_{(-1-2)}, \qquad
\tilde{\mathcal C}_{(12)} = \frac{{\mathcal C}}{{\mathcal R}} 
\varTheta^{-}_{(-1)}\varTheta^{-}_{(-2)}.$$

Specification to 3-point scheme {#3MN}
-------------------------------

We end the review of two-dimensional case with specification to the 3-point scheme i.e. to a discretization of first order differential equation $$g \psi,_x+ h \psi,_y=f\psi.$$ If we put $$A=0,\quad B=0, \quad F=0$$ then according to $$\tilde{A}=0,\quad \tilde{B}=0, \quad \tilde{F}=0$$ It means that the fundamental transformation is also the Darboux transformation for the 3-point scheme [@DJM-II; @Nimmo-KP] $$2C \Psi {_{(12)}}+ G\Psi {_{(1)}}+H\Psi {_{(2)}}=0 .$$ This elementary scheme is the simplest one from the class considered here, but it deserves a special attention, because it leads to one of the most studied integrable discrete equation [@Hirota]. We confine ourselves to recalling briefly in Section \[3p-sch\] main results in this field.

To the end let us rewrite the 3-point scheme in the basic gauge $$(u{_{(2)}}-u{_{(1)}})\Psi {_{(12)}}+(u{_{(1)}}-u)\Psi{_{(1)}}-(u{_{(2)}}-u)\Psi {_{(2)}}=0$$

The four point systems {#sec:4p}
======================

In this Section we present the Darboux transformations for the four point scheme (the discrete Laplace equation) from the point of view of systems of such equations, and the corresponding permutability theorems. To keep the paper of reasonable size and in order to present the results from a simple algebraic perspective we do not discuss important relations of the subject to incidence and difference geometry [@DCN; @MQL; @Bobenko-DO; @CDS; @q-red; @Dol-Koe; @BQL; @CQL; @TQL; @DS-sym; @KingSchief1; @KingSchief; @KoSchief2; @KoSchiefSBKP; @Schief-JNMP; @Schief-II; @DGNS-4-5; @DNS-4-7] (see also [@DS-EMP; @BobSur] and earlier works [@Sauer2; @Sauer]), application of analytic [@BoKo; @DMS; @TQL; @DMMMS; @WTLS; @WTS; @WOGTS; @DMM; @DS-sym; @MM-2; @AD-tauQL] and algebro-geometric [@Krichever-4p; @KWZ; @AKV; @AD-DQL; @AD-Chicago; @BD; @DGNS-4-5; @BQL; @GrKri; @CQL] techniques of the integrable systems theory to construct large classes of solutions of the linear systems in question and solutions of the corresponding nonlinear discrete equations.

In order to simplify discussion of the Darboux transformations for systems of the 4-point schemes we fix (without loss of generality [@MQL]) the gauge to the affine one $$\label{eq:4p-affine}
{\boldsymbol{\Psi}}_{(ij)} = A_{ij(i)}{\boldsymbol{\Psi}}_{(i)}  + 
A_{ji(j)}{\boldsymbol{\Psi}}_{(j)}  +(1 - A_{ij(i)} - A_{ji(j)}) {\boldsymbol{\Psi}}, \qquad i\ne j,$$ where $A_{ij}:{{\mathbb Z}}^N\to{{\mathbb R}}$ are some functions constrained by the compatibility of the system . It is also convenient [@BoKo] to replace the second order linear system by a first order system as follows. The compatibility of allows for definition of the potentials (the Lamé coefficients) $H_i:{{\mathbb Z}}^N\to{{\mathbb R}}$ such that $$A_{ij} = \frac{H_{i(j)}}{H_i}, \; \; \; i\ne j \; .$$ The new wave functions ${\boldsymbol{\psi}}_i$ given by the decomposition $$\label{eq:H-X}
{{\Delta}}_i {\boldsymbol{\Psi}}= H_{i(i)} {\boldsymbol{\psi}}_i\; ,$$ satisfy the first order linear system $$\label{eq:lin-X}
\Delta_j{\boldsymbol{\psi}}_i = Q_{ij(j)}{\boldsymbol{\psi}}_j,    \; \; \; i\ne j \; ,$$ where the functions $Q_{ij}$, called the rotation coefficients, are calculated from the equation $$\label{eq:lin-H}
{{\Delta}}_iH_j =  Q_{ij} H_{i(i)}, \; \; \; i\ne j \; ,$$ which is called adjoint to . Both the equations and are compatible provided the fields $Q_{ij}$ satisfy the discrete Darboux equations [@BoKo] $$\label{eq:MQL-Q}
{{\Delta}}_kQ_{ij} = Q_{ik(k)} Q_{kj}\;\; , \qquad  i\neq j\neq k\neq i.$$

\[cor:4p-psi-i\] Due to the function ${\boldsymbol{\psi}}_i$ satisfies itself the four point equation $$\label{eq:4p-psi-i}
{\boldsymbol{\psi}}_{i(ij)} = {\boldsymbol{\psi}}_{i(i)} + \frac{Q_{ij(ij)}}{Q_{ij(j)}}{\boldsymbol{\psi}}_{i(j)} -
\frac{Q_{ij(ij)}}{Q_{ij(j)}}\left( 1 - Q_{ij(j)}Q_{ji(i)} \right) {\boldsymbol{\psi}}_i,
\qquad i\ne j,$$ while the function $\left( \frac{H_j}{Q_{ij}} \right)_{(ij)}$ satifies the adjoint of in the sense of .

The discrete Darboux equations imply existence of the potentials $\rho_i$ given as solutions of the compatible system $$\label{eq:rho-constr}
\frac{\rho_{i(j)}}{\rho_i} = 1 - Q_{ji(i)} Q_{ij(j)}, \qquad  i\ne j \; ,$$ and yet another potential $\tau$ such that $$\label{eq:tau}
\rho_i = \frac{\tau_{(i)}}{\tau} .$$ In terms of the $\tau$-function and the functions $$\label{eq:tau-ij}
\tau_{ij} = \tau Q_{ij} \; ,$$ the meaning of which will be given in section \[sec:Laplace-tr\], equations and can be rewritten [@DMMMS; @DS-sym] in the bilinear form $$\begin{aligned}
 \label{eq:Hir-ij}
\tau_{(ij)}  \tau & = \tau_{(i)} \tau_{(j)} - \tau_{ji(i)} \tau_{ij(j)} \; ,
\qquad i\ne j\\
 \label{eq:Hir-ijk}
\tau_{ij(k)}  \tau  & = \tau_{(k)} \tau_{ij} + \tau_{ik(k)} \tau_{kj} \;,
\qquad i\ne j \ne k \ne i .\end{aligned}$$

The vectorial fundamental (bilinear Darboux) transformation {#sec:fund}
-----------------------------------------------------------

We start with a simple algebraic fact, whose consequences will be discussed throughout the remaining part of this Section.

\[th:vect-Darb\] Given the solution $\boldsymbol{\phi}_i:\mathbb{Z}^N\to{{\mathbb U}}$, of the linear system , and given the solution $\boldsymbol{\phi}^*_i:\mathbb{Z}^N\to{{\mathbb W}}^*$, of the adjoint linear system . These allow to construct the linear operator valued potential $\boldsymbol{\Omega}[\boldsymbol{\phi},\boldsymbol{\phi}^*]:
\mathbb{Z}^N\to \mathrm{L}({{\mathbb W}},{{\mathbb U}})$, defined by $$\label{eq:Omega-Y-Y}
\Delta_i \boldsymbol{\Omega}[\boldsymbol{\phi},\boldsymbol{\phi}^*] = 
\boldsymbol{\phi}_i \otimes \boldsymbol{\phi}^*_{i(i)}, 
\qquad i = 1,\dots , N.$$ If ${{\mathbb W}}={{\mathbb U}}$ and the potential $\boldsymbol{\Omega}$ is invertible, $\boldsymbol{\Omega}[\boldsymbol{\phi},\boldsymbol{\phi}^*]\in\mathrm{GL}({{\mathbb W}})$, then $$\begin{aligned}
\tilde{{\boldsymbol{\phi}}}_i &={\boldsymbol{\Omega}}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*]^{-1}{\boldsymbol{\phi}}_i, \;\; \; i=1,...,N,
\label{eq:hatY}\\
\tilde{{\boldsymbol{\phi}}}^*_i&={\boldsymbol{\phi}}^*_i{\boldsymbol{\Omega}}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*]^{-1},\;\;\; i=1,...,N,\label{eq:hatY*}\end{aligned}$$ satisfy the linear systems - correspondingly, with the fields $$\tilde{Q}_{ij} = Q_{ij}-\langle {\boldsymbol{\phi}}^*_j|{\boldsymbol{\Omega}}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*]^{-1}|{\boldsymbol{\phi}}_i \rangle,
\;\;\; i,j=1,\dots,N, \;\;i\neq j, \label{eq:hatQ}.$$ In addition, $$\label{eq:hatOm}
{\boldsymbol{\Omega}}[\tilde{{\boldsymbol{\phi}}},\tilde{{\boldsymbol{\phi}}}^*] = C - {\boldsymbol{\Omega}}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*]^{-1},$$ where $C$ is a constant operator.

Notice that because of we have $\Psi = {\boldsymbol{\Omega}}[{\boldsymbol{\psi}},H]$.

The potentials $\rho_i$ and the $\tau$-function transform according to $$\begin{aligned}
\label{eq:fund-vect-rho}
\tilde{\rho}_i & = \rho_i(1 + {\boldsymbol{\phi}}^*_{i(i)}
\boldsymbol{\Omega}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*]^{-1}
{\boldsymbol{\phi}}_i),\\
\tilde{\tau} & = \tau \det\boldsymbol{\Omega}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*].\end{aligned}$$

Applying the above transformation one can produce new compatible (affine) four point linear problems from the old ones.

To obtain conventional transformation formulas consider [@TQL] the following splitting of the vector space ${{\mathbb W}}$ of Theorem \[th:vect-Darb\]: $${{\mathbb W}}= {{\mathbb E}}\oplus{{\mathbb V}}\oplus{{\mathbb F}}\;, \; \; \;
 {{\mathbb W}}^* = {{\mathbb E}}^*\oplus{{\mathbb V}}^*\oplus{{\mathbb F}}^* \; \; \; \; ,$$ if $${\boldsymbol{\phi}}_i = \begin{pmatrix} {\boldsymbol{\psi}}_i \\ {\boldsymbol{\theta}}_i \\ 0 \end{pmatrix} , \; \; \;
{\boldsymbol{\phi}}^*_i = ( 0 , {\boldsymbol{\theta}}^*_i, {\boldsymbol{\psi}}^*_i ) \; ,$$ then, the corresponding potential matrix is of the form $${\boldsymbol{\Omega}}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*] = \begin{pmatrix}
{{\mathbb I}}_{{\mathbb E}}& {\boldsymbol{\Omega}}[{\boldsymbol{\psi}},{\boldsymbol{\theta}}^*] & {\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\psi}}^*] \\
0       & {\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*] & {\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, {\boldsymbol{\psi}}^*] \\
0       &      0          &   {{\mathbb I}}_{{\mathbb F}}\end{pmatrix}$$ and its inverse is $${\boldsymbol{\Omega}}[{\boldsymbol{\phi}},{\boldsymbol{\phi}}^*]^{-1} =
\begin{pmatrix}
{{\mathbb I}}_{{\mathbb E}}& -{\boldsymbol{\Omega}}[{\boldsymbol{\psi}},{\boldsymbol{\theta}}^*]  {\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]^{-1} &
-{\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\psi}}^*] +{\boldsymbol{\Omega}}[{\boldsymbol{\psi}},{\boldsymbol{\theta}}^*]{\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]^{-1}{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, {\boldsymbol{\psi}}^*]  \\
0       & {\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]^{-1} & -  {\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]^{-1}{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, {\boldsymbol{\psi}}^*] \\
0       &      0          &   {{\mathbb I}}_{{\mathbb F}}\end{pmatrix}.$$

Let us consider the case of $K$-dimensional transformation data space, ${{\mathbb V}}={{\mathbb R}}^K$, and ${{\mathbb F}}={{\mathbb R}}$, ${{\mathbb E}}= {{\mathbb R}}^M$, then the transformed solution $\tilde{{\boldsymbol{\Psi}}} = {\boldsymbol{\Omega}}[\tilde{{\boldsymbol{\psi}}}, 
\tilde{H}]$ of the four point scheme (recall that ${\boldsymbol{\Psi}}= {\boldsymbol{\Omega}}[{\boldsymbol{\psi}},H]$) up to a constant vector reads $$\label{eq:fund-x}
\tilde{{\boldsymbol{\Psi}}} = {\boldsymbol{\Psi}}- {\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]{\boldsymbol{\Omega}}[ {\boldsymbol{\theta}}, {\boldsymbol{\theta}}^*]^{-1}
{\boldsymbol{\Omega}}[ {\boldsymbol{\theta}}, H],$$ where the corresponding transformed solutions $\tilde{{\boldsymbol{\psi}}}_i$ of the linear problem and $\tilde{H}_i$ of the adjoint linear problem are given by equations $$\begin{aligned}
\tilde{{\boldsymbol{\psi}}}_i & = {\boldsymbol{\psi}}_i -{\boldsymbol{\Omega}}[{\boldsymbol{\psi}},{\boldsymbol{\theta}}^*]  
{\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]^{-1}{\boldsymbol{\theta}}_i ,\\
\tilde{H}_i & = H_i - {\boldsymbol{\theta}}^*_i {\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]^{-1}{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, H],\end{aligned}$$ and $$\label{eq:Q-fund}
\tilde{Q}_{ij} = Q_{ij} - {\boldsymbol{\theta}}^*_j {\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]^{-1}{\boldsymbol{\theta}}_i.$$ The scalar ($K=1$) fundamental transformation in the above form was given in [@KoSchief2].

To connect the above formalism to the results of Section \[sec:6p\] notice that given scalar solution $\theta_i$ of the linear system then the potential $\Theta = \Omega[\theta, H]$ is the scalar solution of the system of second order linear equations. Moreover, given the solution $\theta^*_i$, of the adjoint linear system then the functions $$\label{eq:sigma-i-def}
\sigma_i = \frac{\theta^*_i}{H_i}$$ satisfy [@TQL] the linear system $$\label{eq:sigma}
\Delta_j \sigma_i = (A_{ij}-1) \left( \sigma_{j(j)} - \sigma_{i(j)} \right),
\qquad i\ne j.$$ Equations imply that the functions $$\label{eq:Phi-ij}
\Phi_{ij} = \sigma_{i(ij)}
- \sigma_{j(ij)},
\qquad i\ne j,$$ satisfy the corresponding equations $$\label{eq:4p-affine-adj-ij}
\Phi_{ij(-i -j)} = A_{ij}\Phi_{ij(-i)} + 
A_{ji} \Phi_{ij(-j)} +( 1- A_{ij(i)} - A_{ji(j)})\Phi_{ij}  , \qquad i\ne j,$$ adjoint of equations . In the case $M=2$ we obtain therefore the data of the transformation being the solution $\Theta$ of the linear problem and the solution $\Phi=\Phi_{12}$ of its adjoint, thus we recover results of Theorem \[th:6p-Dt\] with the four point specification in the affine gauge; see [@Dol-Koe] for more detailed description. Notice that in order to describe fundamental transformation in the second order formalism for $M>2$ one should consider in addition the algebraic relations $$\Phi_{ij} = - \Phi_{ji}, \qquad \Phi_{ij(-i -j)} + \Phi_{jk(-j -k)} =
\Phi_{ik(-i -k)},$$ which are consequences of definition .

It is important to notice that the vectorial fundamental transformation can be obtained as a superposition of $K$ scalar transformations, which follows from the following observation.

\[th:perm-fund\] Assume the following splitting of the data of the vectorial fundamental transformation $${\boldsymbol{\theta}}_i = \left( \begin{array}{c} 
{\boldsymbol{\theta}}_i^a \\ {\boldsymbol{\theta}}_i^b \end{array} \right),\qquad
{\boldsymbol{\theta}}_i^* = \left( \begin{array}{cc} 
{\boldsymbol{\theta}}_{ai}^{*}  & {\boldsymbol{\theta}}_{b i}^{*} \end{array} \right),$$ associated with the partition ${{\mathbb R}}^K = {{\mathbb R}}^{K_a} \oplus {{\mathbb R}}^{K_b}$, which implies the following splitting of the potentials $$\label{eq:split-fund-1}
\boldsymbol{\Omega}[{\boldsymbol{\theta}},H] =  \left( \begin{array}{c} 
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,H] \\ 
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,H] \end{array} \right), \qquad
\boldsymbol{\Omega}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*] = \left( \begin{array}{cc}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}_a^*] &
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}_b^*] \\
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}_a^*] &
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}_b^*]\end{array} \right),$$ $$\label{eq:split-fund-2}
\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}^*] = \left( \begin{array}{cc} 
\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}_a^*]  &
\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}_b^*]\end{array} \right).$$ Then the vectorial fundamental transformation is equivalent to the following superposition of vectorial fundamental transformations:\
1\] Transformation ${\boldsymbol{\Psi}}\to{\boldsymbol{\Psi}}^{\{a\}}$ with the data ${\boldsymbol{\theta}}_i^a$, ${\boldsymbol{\theta}}_{ai}^*$ and the corresponding potentials $\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,H]$, $\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}_a^*]$, $\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}_a^*]$ $$\begin{aligned}
\label{eq:fund-vect-a}
{\boldsymbol{\Psi}}^{\{a\}}  & = {\boldsymbol{\Psi}}- 
\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}^*_a]
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,H],\\
\boldsymbol{\psi}_i^{\{a\}} & = \boldsymbol{\psi}_i -
\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}^*_a]
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
{\boldsymbol{\theta}}^a_i,
\\
H_i^{\{a\}} & = H_i - {\boldsymbol{\theta}}^*_{i a}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,H].\end{aligned}$$ 2\] Application on the result the vectorial fundamental transformation with the transformed data $$\begin{aligned}
{{\boldsymbol{\theta}}}_i^{b\{a\}} & = {\boldsymbol{\theta}}_i^b -
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}^*_a]
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
{\boldsymbol{\theta}}^a_i,
\\
{{\boldsymbol{\theta}}}_{i b}^{*\{a\}} & = {\boldsymbol{\theta}}_{i b}^* - 
{\boldsymbol{\theta}}^*_{i a}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a, {\boldsymbol{\theta}}_{b}^*],\end{aligned}$$ and potentials $$\begin{aligned}
{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}^b,H]^{\{a\}} & =
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,H] - 
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}^*_a]
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,H]=
\boldsymbol{\Omega}[{{\boldsymbol{\theta}}}^{b\{a\}},H^{\{a\}}],
\\
{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}^*_b]^{\{a\}} & =
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}^*_b] - 
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}^*_a]
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_b]=
\boldsymbol{\Omega}[{{\boldsymbol{\theta}}}^{b\{a\}},{{\boldsymbol{\theta}}}_b^{*\{a\}}],
\\
{\boldsymbol{\Omega}}[\boldsymbol{\psi},{\boldsymbol{\theta}}^*_b]^{\{a\}}  & =
\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}^*_b] - 
\boldsymbol{\Omega}[\boldsymbol{\psi},{\boldsymbol{\theta}}^*_a]
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_a]^{-1}
\boldsymbol{\Omega}[{\boldsymbol{\theta}}^a,{\boldsymbol{\theta}}^*_b]=
\boldsymbol{\Omega}[{\boldsymbol{\psi}}^{\{a\}},{{\boldsymbol{\theta}}}_b^{*\{a\}}],
\label{eq:fund-vect-potentials-slit}\end{aligned}$$ i.e., $$\label{eq:fund-vect-a-b}
\tilde{\boldsymbol{\Psi}}= {\boldsymbol{\Psi}}^{\{a,b\}}  = {\boldsymbol{\Psi}}^{\{a\}} - 
{\boldsymbol{\Omega}}[\boldsymbol{\psi},{\boldsymbol{\theta}}^*_b]^{\{a\}}
[{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}^b,{\boldsymbol{\theta}}^*_b]^{\{a\}}]^{-1}
{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}^b,H]^{\{a\}}.$$

The above formulas, apart from existence of the $\tau$-function, remain valid (eventually one needs the proper ordering of some factors) if, instead of the real field ${{\mathbb R}}$, we consider [@Dol-GAQL] arbitrary division ring. Notice, that because the structure of the transformation formulas is a consequence of equation then the formulas may be expressed in terms of quasi-determinants [@GGRW-quasideterminants] (recall, that roughly speaking, a quasi-determinant is the inverse of an element of the inverse of a matrix with entries in a division ring).

Reductions of the fundamental transformation
--------------------------------------------

Let us list basic reductions of the (scalar) fundamental transformation. We follow the nomenclature of [@TQL] which has origins in geometric terminology of transformations of conjugate nets [@Jonas; @DarbouxIV; @Eisenhart-TS; @Tzitzeica; @Bianchi; @Lane; @Finikov]. We provide also the terminology of modern theory of integrable systems [@ms; @OevelSchief], where the fundamental transformation is called the binary Darboux transformation. All the transformations presented in this section can be derived [@TQL] from the fundamental transformation through limiting procedures.

### The Lévy (elementary Darboux) transformation

Given a scalar solution $\theta_i$ of the linear problem the Lévy transform of ${\boldsymbol{\Psi}}={\boldsymbol{\Omega}}[{\boldsymbol{\psi}},H]$ is then given by $${{\mathcal L}}_i({\boldsymbol{\Psi}}) = {\boldsymbol{\Psi}}- \frac{{\boldsymbol{\Omega}}[\theta,H]}{\theta_i} {\boldsymbol{\psi}}_i.$$ The transformed Lamé coefficients and new wave functions are of the form $$\begin{aligned}
{{\mathcal L}}_i(H_i) & = \frac{1}{\theta_i}{\boldsymbol{\Omega}}[\theta,H] \; ,\\
{{\mathcal L}}_i(H_j) & = H_j - \frac{Q_{ij}}{\theta_i}{\boldsymbol{\Omega}}[\theta,H] \; , \\
{{\mathcal L}}_i({\boldsymbol{\psi}}_i) & = -{\boldsymbol{\psi}}_{i(i)} + \frac{\theta_{i(i)}}{\theta_i}{\boldsymbol{\psi}}_i \; ,\\
{{\mathcal L}}_i({\boldsymbol{\psi}}_j) & = {\boldsymbol{\psi}}_j - \frac{\theta_j}{\theta_i}{\boldsymbol{\psi}}_i \; .\end{aligned}$$ Within the nonlocal $\bar\partial$-dressing method the elementary Darboux transformation was introduced in [@BoKo].

### The adjoint Lévy (adjoint elementary Darboux) transformation

Given a scalar solution $\theta^*_i$ of the adjoint linear problem the adjoint Lévy transform ${{\mathcal L}}^*_i({\boldsymbol{\Psi}})$ of ${\boldsymbol{\Psi}}$ is given by $${{\mathcal L}}^*_i({\boldsymbol{\Psi}}) = {\boldsymbol{\Psi}}- \frac{H_i}{\theta^*_i}{\boldsymbol{\Omega}}[{\boldsymbol{\psi}},\theta^*].$$ The new Lamé coefficients and the wave functions are of the form $$\begin{aligned}
{{\mathcal L}}^*_i(H_i) & = H_{i(-i)}- \frac{\theta^*_{i(-i)}}{\theta_i^*}H_i\; ,\\
{{\mathcal L}}^*_i(H_j) & = H_j - \frac{\theta_j^*}{\theta_i^*}H_i \; , \\
{{\mathcal L}}^*_i({\boldsymbol{\psi}}_i) & =  \frac{1}{\theta_i^*}{\boldsymbol{\Omega}}[{\boldsymbol{\psi}},\theta^*] \; ,\\
{{\mathcal L}}^*_i({\boldsymbol{\psi}}_j) & = {\boldsymbol{\psi}}_j - \frac{Q_{ji}}{\theta_i^*}{\boldsymbol{\Omega}}[{\boldsymbol{\psi}},\theta^*] \; .\end{aligned}$$

As it was shown in [@TQL], the scalar fundamental transformation can be obtained as superposition of the Lévy transformation and its adjoint. The closed formulae for iterations of the Lévy transformations in terms of Casorati determinants, and analogous result for the adjoint Lévy transformation, was given in [@LiuManas].

The description of the Lévy transformation and its adjoint in the homogeneous formalism in the case of $M=2$ is given in [@Dol-Koe].

From analytic [@BoKo] and geometric [@TQL] point of view one can distinguish also the so called *Combescure transformation*, whose algebraic description is however very simple (the wave functions ${\boldsymbol{\psi}}_i$ are invariant). The Combescure transformation supplemented by the *projective (or radial* transformation, whose algebraic description is also trivial [@TQL]), generate the fundamental transformation. See also section \[dodane\] and \[Gc\] for generalization of the Combescure transformation.

### The Laplace (Schlesinger) transformation {#sec:Laplace-tr}

The following transformations does not involve any functional parameters, and can be considered as further degeneration of the Lévy (or its adjoint) reduction. The Laplace transformation of ${\boldsymbol{\Psi}}$ is given by $$\label{def:Lij}
{{\mathcal L}}_{ij}({\boldsymbol{\Psi}}) = {\boldsymbol{\Psi}}-\frac{H_j}{Q_{ij}}{\boldsymbol{\psi}}_i , \qquad i\ne j.$$ The Lamé coefficients of the transformed linear problems read $$\begin{aligned}
{{\mathcal L}}_{ij}(H_i) &=  \frac{H_j}{Q_{ij}}
\; , \label{eq:L-Hi} \\
{{\mathcal L}}_{ij}(H_j) &= 
\left( Q_{ij}{{\Delta}}_j \left(\frac{H_j}{Q_{ij}}  \right) \right)_{(-j)}
\; , \label{eq:L-Hj} \\
{{\mathcal L}}_{ij}(H_k) &= H_k - \frac{Q_{ik}}{Q_{ij}} H_j
\; , \; \; k\ne i,j, \label{eq:L-Hk}\end{aligned}$$ and the new wave functions read $$\begin{aligned}
{{\mathcal L}}_{ij}({\boldsymbol{\psi}}_i) & = - {{\Delta}}_i{\boldsymbol{\psi}}_i + \frac{{{\Delta}}_iQ_{ij}}{Q_{ij}} {\boldsymbol{\psi}}_i
\; ,\label{eq:L-Xi} \\
{{\mathcal L}}_{ij}({\boldsymbol{\psi}}_j) & = -\frac{1}{Q_{ij}}{\boldsymbol{\psi}}_i \; ,
\label{eq:L-Xj} \\
{{\mathcal L}}_{ij}({\boldsymbol{\psi}}_k) & = {\boldsymbol{\psi}}_k - \frac{Q_{kj}}{Q_{ij}} {\boldsymbol{\psi}}_i
\; ,\; \; k\ne i,j \label{eq:L-Xk} \; .\end{aligned}$$ The Laplace transformations satisfy generically the following identities $$\begin{aligned}
\label{eq:Lid-ij}
{{\mathcal L}}_{ij} \circ {{\mathcal L}}_{ji} & = \text{id} \; ,\\
{{\mathcal L}}_{jk} \circ {{\mathcal L}}_{ij} &= {{\mathcal L}}_{ik}, \\
{{\mathcal L}}_{ki} \circ {{\mathcal L}}_{ij} &= {{\mathcal L}}_{kj}. \label{eq:Lid-ijk}\end{aligned}$$ The Laplace transformation for the four point affine scheme was introduced in [@DCN] following the geometric ideas of [@Sauer] and independently in [@NovVes] using the factorization approach. The generalization for systems of four point schemes (quadrilateral lattices) was given in [@TQL].

As it was shown in [@Dol-Hir] the functions $\tau_{ij}$ defined by equation are $\tau$-functions of the transformed four point schemes ${{\mathcal L}}_{ij}({\boldsymbol{\Psi}})$ $$\tau_{ij} = {{\mathcal L}}_{ij}(\tau),$$ which, due to , leads [@DCN] to the discrete Toda system [@Hirota].

Distinguished reductions of the four point scheme {#sec:4p-red}
=================================================

In this section we study (systems of) four point linear equations subject to additional constraints, and we provide corresponding reductions of the fundamental transformation. The basic algebraic idea behind such reduced transformations lies in a relationships between solutions of the linear problem and its adjoint, which should be preserved by the fundamental transformation (see, for example, application of this technique in [@NimmoWillox] to reductions of the binary Darboux transformation for the Toda system). Some results presented here have been partially covered in Section \[sec:4p-6p\] but in a different setting.

The Moutard (discrete BKP) reduction {#bkp}
------------------------------------

Consider the system of discrete Moutard equations (the discrete BKP linear problem [@DJM; @NiSchief]) $$\label{eq:BKP-lin}
{\boldsymbol{\Psi}}_{(ij)} - {\boldsymbol{\Psi}}= F_{ij} ({\boldsymbol{\Psi}}_{(i)} - {\boldsymbol{\Psi}}_{(j)}) , \quad 1\leq i< j\leq N,$$ for suitable functions $F_{ij}:{{\mathbb Z}}^N\to{{\mathbb R}}$. Compatibility of the system implies existence of the potential $\tau^B:{{\mathbb Z}}^N\to{{\mathbb R}}$, in terms of which the functions $F_{ij}$ can be written as $$\label{eq:tau-B}
F_{ij} = \frac{\tau^B_{(i)}\tau^B_{(j)}}{\tau^B \, \tau^B_{(ij)}}, 
\qquad i\ne j ,$$ which satisfies system of Miwa’s discrete BKP equations [@Miwa] $$\label{eq:BKP-nlin}
\tau^B \tau^B_{(ijk)} = \tau^B_{(ij)}\tau^B_{(k)} - \tau^B_{(ik)}\tau^B_{(j)} + 
\tau^B_{(jk)}\tau^B_{(i)}, \quad 1\leq i< j < k \leq N.$$ The discrete Moutard system can be given [@BQL] the first order formulation - upon introducing the Lamé coefficients $$\label{eq:H-tau-B}
H_i = (-1)^{\sum_{k<i}m_k}\frac{\tau^B_{(-i)}}{\tau^B},$$ and the rotation coefficients (below we assume $i<j$) $$\begin{aligned}
 \label{eq:Q-BQL-ij}
Q_{ij} & = - (-1)^{\sum_{i\leq k<j}m_k}
\left( \frac{\tau^B_{(i-j)}}{\tau^B_{(-j)}} + \frac{\tau^B_{(i)}}{\tau^B}\right) 
\frac{\tau^B_{(-j)}}{\tau^B},\\
\label{eq:Q-BQL-ji}
Q_{ji} & = - (-1)^{\sum_{i\leq k<j}m_k}
\left( \frac{\tau^B_{(-i j)}}{\tau^B_{(-i)}} - \frac{\tau^B_{(j)}}{\tau^B}\right)
\frac{\tau^B_{(-i)}}{\tau^B},\end{aligned}$$ which in view of , gives the familiar relation between the $\tau$-functions of the KP and BKP hierarchies $$\tau = \left(\tau^B\right)^2.$$ The corresponding reduction of the fundamental transformation was given in [@BQL], where also a link with earlier work [@NiSchief] on the discrete Moutard transformation has been established.

\[prop:vect-fund-red-BQL\] Given solution $\boldsymbol{\theta}_i:{{\mathbb Z}}^N\to{{\mathbb R}}^K$ of the linear problem corresponding to the Moutard linear system and its first order form -. Denote by $\boldsymbol{\Theta} = \boldsymbol{\Omega}[\boldsymbol{\theta},H]$ the corresponding potential, which is also new vectorial solution of the linear problem .\
1) Then $$\label{eq:BQL-Y*}
\boldsymbol{\theta}^*_i = (-1)^{\sum_{k<i}m_k}\frac{\tau^B_{(-i)}}{\tau^B} 
(\boldsymbol{\Theta}_{(-i)}^T + \boldsymbol{\Theta}^T)$$ provides a vectorial solution of the adjoint linear problem, and the corresponding potential $\boldsymbol{\Omega}[\boldsymbol{\theta},\boldsymbol{\theta}^*]$ allows for the following constraint $$\label{eq:vect-Moutard-constr}
\boldsymbol{\Omega}[\boldsymbol{\theta},\boldsymbol{\theta}^*] + 
\boldsymbol{\Omega}[\boldsymbol{\theta},\boldsymbol{\theta}^*]^T =
2 \boldsymbol{\Theta}\otimes \boldsymbol{\Theta}^T.$$ 2) The fundamental vectorial transform $\tilde{\boldsymbol{\Psi}}$ of ${\boldsymbol{\Psi}}$, given by with the potentials $\boldsymbol{\Omega}$ restricted as above satisfies Moutard linear system and can be considered as the superposition of $K$ scalar reduced fundamental transforms.

In the scalar case and for $N=2$ inserting in equations and the functions $H_i$ and $\theta^*i$ given by and correspondingly, we obtain the formula in the gauge $F=1$ (modulo the corresponding change of sign).

Notice that given $\boldsymbol{\Theta}$ then, because of the constraint , to construct $\boldsymbol{\Omega}(\boldsymbol{\theta},\boldsymbol{\theta}^*)$ we need only its antisymmetric part $S(\boldsymbol{\Theta} | \boldsymbol{\Theta})$, which satisfies the system $$\Delta_i S(\boldsymbol{\Theta} | \boldsymbol{\Theta}) = 
\boldsymbol{\Theta}_{(i)}\otimes \boldsymbol{\Theta}^{\mathrm{t}} -
\boldsymbol{\Theta}\otimes \boldsymbol{\Theta}^{\mathrm{t}}_{(i)}.$$ This observation is the key element of the connection of the above reduction of the fundamental transformation with earlier results [@NiSchief] on the vectorial Moutard transformation for the system , where the formulas using Pfaffians were obtained (recall that determinant of a skew-symmetric matrix is a square of Pfaffian). In particular, the transformation rule for the $\tau^B$-function can be recovered $$\tilde\tau^B = \begin{cases}
\tau^B \: \mathrm{Pf} \: S(\boldsymbol{\Theta} | \boldsymbol{\Theta}) , & 
K \; \text{even} \\
\tau^B \: \mathrm{Pf} \left( \begin{array}{cc}
0 & - \boldsymbol{\Theta}^T , \\
\boldsymbol{\Theta} &  S(\boldsymbol{\Theta} | \boldsymbol{\Theta})
\end{array} \right), \qquad & K \; \text{odd}.
\end{cases}$$

The symmetric (discrete CKP) reduction {#ckp}
--------------------------------------

Consider the linear problem subject to the constraint [@DS-sym] which arose from studies on the Egorov lattices [@Schief-talk] $$\label{eq:symm}
\rho_i Q_{ji(i)}=\rho_j Q_{ij(j)} \; , \quad i\ne j \; .$$ Then the discrete Darboux equations can be rewritten [@Schief-JNMP] in the following quartic form $$\begin{split}
\left( \tau   \tau_{(ijk)} - \tau_{(i)}  \tau_{(jk)} -
\tau_{(j)}  \tau_{(ik)} -  \tau_{(k)}  \tau_{(ij)} \right)^2 
+ 4(\tau_{(i)}\tau_{(j)}\tau_{(k)}\tau_{(ijk)} + &
\tau  \tau_{(ij)}\tau_{(jk)} \tau_{(ik)} )= \\
4( \tau_{(i)}  \tau_{(jk)} \tau_{(j)}  \tau_{(ik)} + \tau_{(i)}  \tau_{(jk)} 
\tau_{(k)}  \tau_{(ij)}  + 
\tau_{(j)}  \tau_{(ik)}\tau_{(k)}  \tau_{(ij)}), &
\end{split}$$ which can be identified with equation derived in [@Kashaev] in connection with the star-triangle relation in the Ising model. According to [@Schief-JNMP], the above equation can be obtained from the CKP hierarchy via successive application of the corresponding reduction of the binary Darboux transformations.

Construction [@MM] of the reduction of the fundamental transformations which preserves the constraint makes use the following observation.

The following conditions are equivalent:\
1) The functions $Q_{ij}$, $\rho_i$ satisfy constraint ;\
2) Given a nontrivial solution ${\boldsymbol{\phi}}_i^*$ of the adjoint linear problem then $$\label{eq:C-ft-YY}
{\boldsymbol{\phi}}_i = \rho_i ({\boldsymbol{\phi}}^*_{i(i)})^T$$ provides a solution of the linear problem ;\
3) The corresponding potential $\boldsymbol{\Omega}[\boldsymbol{\phi},\boldsymbol{\phi}^*] $ allows for the constraint $$\label{eq:C-ft-O}
\boldsymbol{\Omega}[\boldsymbol{\phi},\boldsymbol{\phi}^*]^T=
\boldsymbol{\Omega}[\boldsymbol{\phi},\boldsymbol{\phi}^*].$$

When the data ${\boldsymbol{\theta}}_i$, ${\boldsymbol{\theta}}^*_i$, ${\boldsymbol{\Omega}}[{\boldsymbol{\theta}},{\boldsymbol{\theta}}^*]$ of the fundamental transformation satisfy conditions - then the new functions $\tilde{Q}_{ij}$, $\tilde{\rho}_i$ found by equations , are constrained by , i.e. $$\tilde{\rho}_i \tilde{Q}_{ji(i)}=\tilde{\rho}_j \tilde{Q}_{ij(j)} \; , \quad i\ne j \; .$$

The corresponding permutability principle has been proved in [@CQL].

To connect the symmetric reduction with the Goursat equation and the corresponding transformation (see Section \[Gr\]) notice that because of Corollary \[cor:4p-psi-i\] and equations - the function ${\boldsymbol{\psi}}_2$ satisfies equation with $q=Q_{21(1)}$. Then, due to the same Corollary and condition , given scalar solution $\theta_2$ of the four point equation of ${\boldsymbol{\psi}}_2$ then the function $$\left(\frac{\theta_1}{\rho_2
Q_{21(1)}}\right)_{(2)} = \left(\frac{\tau \Delta_1 \theta_2}{\tau_{(2)}
Q_{21(1)}^2}\right)_{(2)},$$ satisfies its adjoint; in the lat equation we have used the linear problem and equation .

Quadratic reduction {#Qr}
-------------------

Consider the system of four point equations which solution ${\boldsymbol{\Psi}}$ is subject to the following quadratic constraint $$\label{eq:quadric}
{\boldsymbol{\Psi}}^T  Q  {\boldsymbol{\Psi}}+  {\boldsymbol{a}}^T {\boldsymbol{\Psi}}+ c = 0  \; \; \; ;$$ here $Q$ is a non-degenerate symmetric matrix, ${\boldsymbol{a}}$ is a constant vector, $c$ is a scalar.

Notice that unlike in two previous reductions we fix (by giving the quadratic equation) the dimension of ${\boldsymbol{\Psi}}$.

Double discrete differentiation of equation in $i\ne j$ directions gives, after some algebra, the condition $$\label{eq:XQX}
{\boldsymbol{\psi}}_{i(j)}^T  Q  {\boldsymbol{\psi}}_j +  {\boldsymbol{\psi}}_{j(i)}^T  Q  {\boldsymbol{\psi}}_i  = 0,$$ analogous to that obtained in [@DMS] in order to characterize circular lattices [@Bobenko-DO; @CDS]. It implies [@DS-sym], in particular, that ${\boldsymbol{\psi}}_i^T  Q  {\boldsymbol{\psi}}_i$ satisfy the same equation as the potentials $\rho_i$.

As in two above reductions, the quadratic condition allows for a relation between solutions of the linear system and its adjoint . The following proposition can be easily derived from analogous results of [@q-red], where as the basic ingredient of the transformation was used the potential ${\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]$, but we present here its direct proof in the spirit of corresponding results found for the circular lattice [@KoSchief2; @LiuManas-SR].

Given a nontrivial solution ${\boldsymbol{\theta}}_i^*$ of the adjoint linear problem corresponding to the system of four point equations which solution ${\boldsymbol{\Psi}}$ is subject to the constraint then $$\label{eq:Q-ft-YY}
{\boldsymbol{\theta}}_i = ({\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]_{(i)} + {\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*])^T Q {\boldsymbol{\psi}}_i$$ provides a solution of the corresponding linear problem .

After some algebra using the equations satisfied by ${{\boldsymbol Y}}^*_i$ and ${\boldsymbol{\psi}}_i$ one gets $${{\Delta}}_j {\boldsymbol{\theta}}_i - Q_{ij(j)}{\boldsymbol{\theta}}_j = {\boldsymbol{\theta}}^{*T}_{j(ij)}\left(
{\boldsymbol{\psi}}_{i(j)}^T  Q  {\boldsymbol{\psi}}_j +  {\boldsymbol{\psi}}_{j(i)}^T  Q  {\boldsymbol{\psi}}_i \right),$$ which vanishes due to .

The following result gives the discrete Ribaucour reduction of the fundamental transformation.

Given solution $\boldsymbol{\theta}^*_i:{{\mathbb Z}}^N\to{{\mathbb R}}^K$ of the adjoint linear problem corresponding to the quadratic constraint .\
1) Then the potentials ${\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]$, ${\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, H]$ and ${\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, {\boldsymbol{\theta}}^*]$, where ${\boldsymbol{\theta}}_i$ is the solution of the linear problem constructed from ${\boldsymbol{\theta}}^*_i$ by means of formula , allow for the constraints $$\begin{aligned}
\label{eq:bphi-qv}
{\boldsymbol{\Omega}}[{\boldsymbol{\theta}},H]^T & = 2 {\boldsymbol{\psi}}^T Q  {\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]  +
 {\boldsymbol{a}}^T {\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]  \; \; , \\
\label{eq:bPhiC-qv}
{\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, {\boldsymbol{\theta}}^*] + {\boldsymbol{\Omega}}[{\boldsymbol{\theta}}, {\boldsymbol{\theta}}^*]^T & = 
 2 {\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]^T  Q {\boldsymbol{\Omega}}[{\boldsymbol{\psi}}, {\boldsymbol{\theta}}^*]  \; .\end{aligned}$$ 2) The Ribaucour reduction of the fundamental vectorial transform $\tilde{\boldsymbol{\Psi}}$ of ${\boldsymbol{\Psi}}$, given by with the potentials $\boldsymbol{\Omega}$ restricted as above satisfies the quadratic constraint and can be considered as the superposition of $K$ scalar Ribaucour transformations.

As it was explained in [@q-red], the Ribaucour transformations [@KoSchief2] of the circular lattice [@Bobenko-DO; @CDS] can be derived from the above approach after the stereographic projection from the Möbius sphere. The superposition principle for the Ribaucour transformation of circular lattices was derived also in [@LiuManas-SR].

The three point scheme {#3p-sch}
======================

In this final Section we present the vectorial Darboux transformations for the three point scheme . The corresponding nonlinear difference system , known as the Hirota–Miwa equation, is perhaps the most important and widely studied integrable discrete system. It was discovered by Hirota [@Hirota], who called it the discrete analogue of the two dimensional Toda lattice (see also [@LPS]), as a culmination of his studies on the bilinear form of nonlinear integrable equations. General feature of Hirota’s equation was uncovered by Miwa [@Miwa] who found a remarkable transformation which connects the equation to the KP hierarchy [@DKJM]. The Hirota-Miwa equation, called also the discrete KP equation, can be encountered in various branches of theoretical physics [@Saito; @KNS; @Zabrodin] and mathematics [@Shiota; @KWZ; @Octahedron].

Consider the linear system [@DJM-II] $$\label{eq:lin-dKP}
{\boldsymbol{\Psi}}_{(i)} - {\boldsymbol{\Psi}}_{(j)} = U_{ij} {\boldsymbol{\Psi}},  \qquad i \ne j \leq N,$$ whose compatibility leads to the following parametrization of the field $U_{ij}$ in terms of the potentials $r_i$ $$r_{i(j)} = r_i U_{ij},$$ and then to existence of the $\tau$ function $$r_i = (-1)^{\sum_{k<i}m_k}\frac{\tau^H_{(i)}}{\tau^H}$$ and, finally to the the discrete KP system [@Hirota; @Miwa] $$\label{eq:H-M}
\tau^H_{(i)}\tau^H_{(jk)} - \tau^H_{(j)}\tau^H_{(ik)} + \tau^H_{(k)}\tau^H_{(ij)},
\qquad i< j <k.$$ The same nonlinear systems arises from compatibility of $$\label{eq:lin-dKP*}
{\boldsymbol{\Psi}}^*_{(j)} - {\boldsymbol{\Psi}}^*_{(i)} = 
U_{ij} {\boldsymbol{\Psi}}^*_{(ij)} , 
\qquad i \ne j ,$$ called the adjoint of .

We present the Darboux transformation for the three point scheme in the way similar to that of Section \[sec:fund\] following the approach of [@Nimmo-KP], see however early works on the subject [@Saito-Saitoh; @Saitoh-Saito].

\[th:vect-Darb-KP\] Given the solution $\boldsymbol{\Phi}:\mathbb{Z}^N\to{{\mathbb U}}$, of the linear system , and given the solution $\boldsymbol{\Phi}^*:\mathbb{Z}^N\to{{\mathbb W}}^*$, of the adjoint linear system . These allow to construct the linear operator valued potential $\boldsymbol{\Omega}[\boldsymbol{\Phi},\boldsymbol{\Phi}^*]:
\mathbb{Z}^N\to \mathrm{L}({{\mathbb W}},{{\mathbb U}})$, defined by $$\label{eq:omega-p-p}
\Delta_i \boldsymbol{\Omega}[\boldsymbol{\Phi},\boldsymbol{\Phi}^*] = 
\boldsymbol{\Phi} \otimes \boldsymbol{\Phi}^*_{(i)}, 
\qquad i = 1,\dots , N.$$ If ${{\mathbb W}}={{\mathbb U}}$ and the potential $\boldsymbol{\Omega}$ is invertible, $\boldsymbol{\Omega}[\boldsymbol{\Phi},\boldsymbol{\Phi}^*]\in\mathrm{GL}({{\mathbb W}})$, then $$\begin{aligned}
\tilde{\boldsymbol{\Phi}} &=
{\boldsymbol{\Omega}}[{\boldsymbol{\Phi}},{\boldsymbol{\Phi}}^*]^{-1}\boldsymbol{\Phi}, 
\label{eq:hatp}\\
\tilde{\boldsymbol{\Phi}}^*&=\boldsymbol{\Phi}^*
{\boldsymbol{\Omega}}[{\boldsymbol{\Phi}},{\boldsymbol{\Phi}}^*]^{-1},\label{eq:hatp*}\end{aligned}$$ satisfy the linear systems and , correspondingly, with the fields $$\tilde{U}_{ij} = U_{ij} - 
(\boldsymbol{\Phi}^* {\boldsymbol{\Omega}}[\boldsymbol{\Phi},\boldsymbol{\Phi^*}]^{-1} 
\boldsymbol{\Phi})_{(i)}
+ (\boldsymbol{\Phi}^* {\boldsymbol{\Omega}}[\boldsymbol{\Phi},\boldsymbol{\Phi}^*]^{-1}
\boldsymbol{\Phi})_{(j)}.$$ In addition, $$\label{eq:hat-om}
{\boldsymbol{\Omega}}[\tilde{\boldsymbol{\Phi}},\tilde{\boldsymbol{\Phi}}^*] = C - 
{\boldsymbol{\Omega}}[\boldsymbol{\Phi},\boldsymbol{\Phi}^*]^{-1},$$ where $C$ is a constant operator.

The transformation rule for the potentials $r_i$ reads $$\tilde{r}_i = r_i (1 - \boldsymbol{\Phi}^*_{(i)} 
{\boldsymbol{\Omega}}[\boldsymbol{\Phi},\boldsymbol{\Phi}^*]^{-1} \boldsymbol{\Phi} ),$$ while using the technique of the bordered determinants [@Hirota-book] one can show that [@Nimmo-KP] $$\tilde{\tau}^H = \tau^H \det \boldsymbol{\Omega}
[\boldsymbol{\Phi},\boldsymbol{\Phi}^*]  \label{eq:hat-tau-dKP}.$$

More standard transformation formulas arise, when one splits, like in Section \[sec:fund\], the vector space ${{\mathbb W}}$ of Theorem \[th:vect-Darb-KP\] as follows $${{\mathbb W}}= {{\mathbb E}}\oplus{{\mathbb V}}\oplus{{\mathbb F}}\;, \; \; \;
 {{\mathbb W}}^* = {{\mathbb E}}^*\oplus{{\mathbb V}}^*\oplus{{\mathbb F}}^* \; \; \; \; ;$$ if $$\boldsymbol{\Phi} = \begin{pmatrix} 
 \boldsymbol{\Psi} \\ {\boldsymbol{\Theta}}\\ 0 \end{pmatrix} , \; \; \;
\boldsymbol{\Phi}^* = ( 0 , {\boldsymbol{\Theta}}^*, \boldsymbol{\Psi}^* ) \; .$$ Then the corresponding potential matrix and its inverse have the structure like those in Section \[sec:fund\], which gives [@Nimmo-KP] $$\begin{aligned}
\tilde{\boldsymbol{\Psi}} & = \boldsymbol{\Psi} -
{\boldsymbol{\Omega}}[\boldsymbol{\Psi},{\boldsymbol{\Theta}}^*]  {\boldsymbol{\omega}}[{\boldsymbol{\Theta}},{\boldsymbol{\Theta}}^*]^{-1}{\boldsymbol{\Theta}},\\
\tilde{\boldsymbol{\Psi}}^* & = \boldsymbol{\Psi}^* - 
{\boldsymbol{\Theta}}^* {\boldsymbol{\Omega}}[{\boldsymbol{\Theta}},{\boldsymbol{\Theta}}^*]^{-1}{\boldsymbol{\Omega}}[{\boldsymbol{\Theta}}, \boldsymbol{\Psi}^*],\\
\tilde{U}_{ij}&  = U_{ij} - ({\boldsymbol{\Theta}}^* {\boldsymbol{\Omega}}[{\boldsymbol{\Theta}},{\boldsymbol{\Theta}}^*]^{-1} {\boldsymbol{\Theta}})_{(i)}
+ ({\boldsymbol{\Theta}}^* {\boldsymbol{\Omega}}[{\boldsymbol{\Theta}},{\boldsymbol{\Theta}}^*]^{-1} {\boldsymbol{\Theta}})_{(j)}.\end{aligned}$$ Notice that one can consider [@FWN; @FWN-Capel; @Nimmo-NCKP] the three point linear problem in associative algebras. Then the structure of the transformation formulas implies the quasi-determinant interpretation [@GNO] of the above equations.

The permutability property for the Darboux transformations of the three point scheme [@Nimmo-NCKP] can be formulated exactly (cancel subscripts $i$) like Theorem \[th:perm-fund\].

Finally we remark that the binary Darboux transformation for the three point linear problem can be decomposed [@Nimmo-KP] into superposition of the elementary Darboux transformation and its adjoint, which can be described as follows. Given a scalar solution $\Theta$ of the linear problem then the elementary Darboux transformation $$\begin{aligned}
{{\mathcal D}}(\boldsymbol{\Psi}) &= \boldsymbol{\Psi}_{(i)} - \frac{\Theta_{(i)}}{\Theta}
\boldsymbol{\Psi}, \\
{{\mathcal D}}(\boldsymbol{\Psi}^*) &= \frac{1}{\Theta}{\boldsymbol{\Omega}}[\Theta,\boldsymbol{\Psi}^*]\; ,\\
{{\mathcal D}}(\tau) &= \Theta \tau\end{aligned}$$ leaves equations and invariant.

Analogously, given a scalar solution $\Theta^*$ of the linear problem then the adjoint elementary Darboux $$\begin{aligned}
{{\mathcal D}}^*(\boldsymbol{\Psi}) &=\frac{1}{\Theta^*}{\boldsymbol{\Omega}}[\Theta^*,\boldsymbol{\Psi}],
 \\
{{\mathcal D}}^*(\boldsymbol{\Psi}^*) &=  - \boldsymbol{\Psi}^*_{(-i)} + 
\frac{\Theta^*_{(-i)}}{\Theta^*}
\boldsymbol{\Psi}^*, \\
{{\mathcal D}}^*(\tau) &= \Theta^* \tau\end{aligned}$$ leaves equations and invariant.

The above transformations allow for vectorial forms, which can be conveniently written [@Nimmo-KP; @OHTI] in terms of Casorati determinants.

Summary and open problems
=========================

In the paper we aimed to present results on Darboux transformations of linear operators on two dimensional regular lattices. To put some order into the review we started from the six point scheme (and its Darboux transformation) as the master linear problem. The path between its various specifications and reductions has been visualized in Figure \[5p\]. We considered also in more detail the corresponding theory of the Darboux transformations of the systems of the four point schemes and their reductions. Finally, we briefly discussed the multidimensional aspects of the three point schemes. It is worth to mention that for systems of the three or four point schemes the Darboux transformations can be interpreted as a way to generate new dimensions. This is very much connected to the permutability of the transformations, which is a core of integrability of the corresponding nonlinear systems.

Separate issue touched here is such extention of four point schemes that can be regarded as an analogue of discretization of a diferential equation in arbitrary parametrization. More precisely we discussed here such an extension of general (AKP) case (section \[Gc\]) and Moutard–selfadjoint (BKP) case (sections \[Mr\] and \[sec:s-a\]). Goursat and Ribaucour reductions have not been investigated from this point of view. The multidimensional schemes that mimic equations governing conjugate nets (and their reductions) with arbitrary change of independend variables have been not exploited either.

We also mentioned another approach to the problem of construction of the unified theory of the Darboux transformations. It consists in isolating basic “bricks” in order to use their combinations to construct more involved linear operators together with their Darboux transformations. Such an idea has been applied to derive, starting from the Moutard reduction of the four point scheme, the self-adjoint operators on the square (the five point scheme), triangular (the seven point scheme) and the honeycomb grids. Recently it was shown in [@DODO] that the theory of systems of four point linear equations and their Laplace transformations follows from the theory of the three point systems. This means, that also the transformations of the four point scheme can be, in principle, derived from the three point scheme. An open question is if the six point scheme and its transformations can be decomposed in a similar way.

Finally, we would like to mention possibility of considering (hierarchies of) continuous deformations of the above lattice linear problems, which would lead to (hierarchies of) discrete-differential integrable — by construction — equations. Some aspects of deformations of the self-adjoint seven and five linear systems have been elaborated in [@SND; @Gilson; @SDN-Toda].

Acknowledgements {#acknowledgements .unnumbered}
================

The authors thank the Isaac Newton Institute for Mathematical Sciences (Cambridge, UK) for hospitality during the programme *Discrete Integrable Systems*.

[100]{}

M. J. Ablowitz, J. F. Ladik, *Nonlinear differential–difference equations* J. Math. Phys. **16** (1975) 598–603.

V. E. Adler, S. Ya. Startsev, *Discrete analogues of the Liouville equation*, Theor. Math. Phys. **121**(1999)Pages 271–284.

V. E. Adler, *Discrete equations on planar graphs*, J. Phys. A: Math. Gen. **34** (2001) 10453-60.

A. A. Akhmetshin, I. M. Krichever and Y. S. Volvovski, *Discrete analogues of the Darboux–Egoroff metrics*, Proc. Steklov Inst. Math. **225** (1999) 16–39.

M. Bia[ł]{}ecki, A. Doliwa, *Algebro-geometric solution of the discrete KP equation over a finite field out of a hyperelliptic curve*, Comm. Math. Phys. **253** (2005), 157–170.

L. Bianchi, *Lezioni di geometria differenziale*, Zanichelli, Bologna, 1924.

A. Bobenko, [*Discrete conformal maps and surfaces*]{}, \[in:\] [*Symmetries and Integrability of Difference Equations II*]{}, P. Clarkson and F. Nijhoff (eds.), Cambridge University Press.

A. I. Bobenko and Yu. B. Suris, *Discrete differential geometry: integrable structure*, AMS, Providence, 2009.

A. I. Bobenko, Ch. Mercat, Yu. B. Suris, *Linear and nonlinear theories of discrete analytic functions. Integrable structure and isomonodromic Green’s function*, J. Reine Angew. Math. **583** (2005) 117–161.

L. V. Bogdanov and B. G. Konopelchenko, *Lattice and $q$-difference Darboux–Zakharov–Manakov systems via $\bar\partial$ method*, J. Phys. A: Math. Gen. **28** L173–L178.

L. V. Bogdanov and B. G. Konopelchenko, *Analytic-bilinear approach to integrable hiererchies I. Generalizded KP hierarchy*, J. Math. Phys. **39** (1998) 4683–.

L. V. Bogdanov and B. G. Konopelchenko, *Analytic-bilinear approach to integrable hiererchies II. Multicomponent KP and 2D Toda hiererchies*, J. Math. Phys. **39** (1998) 4701–4728.

J. Cieśliński, A. Doliwa and P. M. Santini, [*The integrable discrete analogues of orthogonal coordinate systems are multidimensional circular lattices*]{}, Phys. Lett. A [**235**]{} (1997) 480–488.

G. Darboux, *Leçons sur les systémes orthogonaux et les coordonnées curvilignes*, Gauthier-Villars, Paris, 1910.

G. Darboux, *Leçons sur la théorie générale des surfaces. [I–IV]{}*, Gauthier – Villars, Paris, 1887–1896.

E. Date, M. Kashiwara, M. Jimbo and T. Miwa, *Transformation groups for soliton equations*, \[in:\] Proceedings of RIMS Symposium on Non-Linear Integrable Systems — Classical Theory and Quantum Theory (M. Jimbo and T. Miwa, eds.) World Science Publishing Co., Singapore, 1983, pp. 39–119.

E. Date, M. Jimbo and T. Miwa, *Method for generating discrete soliton equations. I*, J. Phys. Soc. Japan **51** (1982) 4116–24.

E. Date, M. Jimbo and T. Miwa, *Method for generating discrete soliton equations. II*, J. Phys. Soc. Japan **51** (1982) 4125–31.

E. Date, M. Jimbo and T. Miwa, *Method for generating discrete soliton equations. V*, J. Phys. Soc. Japan **52** (1983) 766–771.

A. Doliwa, [*Geometric discretisation of the Toda system*]{}, Phys. Lett. A [**234**]{} (1997) 187–192.

A. Doliwa, *Quadratic reductions of quadrilateral lattices*, J. Geom. Phys. **30** (1999) 169–186.

A. Doliwa, [*The Darboux-type transformations of integrable lattices*]{} , Rep. Math. Phys. [**48**]{} (2001) 59–66.

A. Doliwa, *Integrable multidimensional discrete geometry: quadrilateral lattices, their transformations and reductions*, \[in:\] Integrable Hierarchies and Modern Physical Theories, (H. Aratyn and A. S. Sorin, eds.) Kluwer, Dordrecht, 2001 pp. 355–389.

A. Doliwa, *Lattice geometry of the [H]{}irota equation*, \[in:\] SIDE III – Symmetries and Integrability of Difference Equations, D. Levi and O. Ragnisco (eds.), pp. 93–100, CMR Proceedings and Lecture Notes, vol. 25, AMS, Providence, 2000.

A. Doliwa, *Geometric discretization of the [K]{}oenigs nets*, J. Math. Phys. **44** (2003), 2234–2249.

A. Doliwa, *The B-quadrilateral lattice, its transformations and the algebro-geometric construction*, J. Geom. Phys. **57** (2007) 1171–1192.

A. Doliwa, [*The C-(symmetric) quadrilateral lattice, its transformations and the algebro-geometric construction*]{}, [arXiv:nlin.SI/0710.5820]{}.

A. Doliwa, *Geometric algebra and quadrilateral lattices*, `arXiv:nlin.SI/0801.0512`.

A. Doliwa, [*On $\tau$-function of the quadrilateral lattice*]{}, J. Phys. A (to appear), [arXiv:nlin.SI/0901.0112]{}.

A. Doliwa, [*Desargues maps and the Hirota-Miwa equations*]{}, [arXiv:nlin.SI/0906.1000]{}.

A. Doliwa and P. M. Santini, *Multidimensional quadrilateral lattices are integrable*, Phys. Lett. A **233** (1997), 365–372.

A. Doliwa and P. M. Santini, *The symmetric, [D]{}-invariant and [E]{}gorov reductions of the quadrilateral lattice*, J. Geom. Phys. **36** (2000) 60–102.

A. Doliwa and P. M. Santini, *Integrable systems and discrete geometry*, \[in:\] Encyclopedia of Mathematical Physics, J. P. François, G. Naber and T. S. Tsun (eds.), Elsevier, 2006, Vol. III, pp. 78-87.

A. Doliwa, S. V. Manakov and P. M. Santini, *$\bar\partial$-reductions of the multidimensional quadrilateral lattice: the multidimensional circular lattice*, Comm. Math. Phys. **196** (1998), 1–18.

A. Doliwa, P. M. Santini and M. Ma[ñ]{}as, *Transformations of quadrilateral lattices*, J. Math. Phys. **41** (2000) 944–990.

A. Doliwa, M. Mañas, L. Martínez Alonso, [*Generating quadrilateral and circular lattices in KP theory*]{}, Phys. Lett. A [**262**]{} (1999) 330–343.

A. Doliwa, M. Mañas, L. Martínez Alonso, E. Medina and P. M. Santini, [*Charged free fermions, vertex operators and transformation theory of conjugate nets*]{}, J. Phys. A [**32**]{} (1999) 1197–1216.

A. Doliwa, M. Nieszporski and P. M. Santini, [*Integrable lattices and their sub-lattices II. From the B-quadrilateral lattice to the self-adjoint schemes the triangular and the honeycomb lattices*]{}, J. Math. Phys. [**48**]{} (2007) art. no. 113056.

A. Doliwa, P. Grinevich, M. Nieszporski & P. M. Santini, [*Integrable lattices and their sub-lattices: from the discrete Moutard (discrete Cauchy-Riemann) 4-point equation to the self-adjoint 5-point scheme*]{}, J. Math. Phys. [**48**]{} (2007) art. no. 013513.

L. P. Eisenhart, *Transformations of surfaces*, Princeton University Press, Princeton, 1923.

S. P. Finikov, *Theorie der [K]{}ongruenzen*, Akademie-Verlag, Berlin, 1959.

I. Gelfand, S. Gelfand, V. Retakh and R. L. Wilson, *Quasideterminants*, Adv. Math. **193** (2005), 56–141.

C. R. Gilson, J. J. C. Nimmo *The relation between a 2D Lotka-Volterra equation and a 2D Toda Lattice*, Journal of Nonlinear Mathematical Physics, [**12**]{} Supplement (2005) 169-179.

C. R. Gilson, J. J. C. Nimmo, Y. Ohta, *Quasideterminant solutions of a non-Abelian Hirota–Miwa equation*, J. Phys. A: Math. Theor. **40** (2007) 12607-12617.

E. Goursat, *Sur une transformation de l’équation $s^2=4\lambda(x,y) pq$*, Bulletin de la S.M.F. 28 (1900) 1-6.

B. Grammaticos, Y. Kosmann-Schwarzbach, T. Tamizhmani (Eds.), *Discrete integrable systems*, Lect. Notes Phys. **644**, Springer, Berlin, 2004.

S. Grushevsky, I. Krichever, *Integrable discrete Schrödinger equations and characterization of Prym varieties by a pair of quadrisecants*, `arXiv:math.AG/0705.2829`.

C. H. Gu, H. S. Hu and Z. X. Zhou, *Darboux transformations in integrable systems. Theory and their applications to geometry*, Springer, Dordrecht, 2005

R. Hirota, *Discrete analogue of a generalized Toda equation*, J. Phys. Soc. Jpn. [**50**]{} (1981) 3785-3791.

R. Hirota, *The direct method in soliton theory*, Cambridge University Press, 2004.

H. Jonas, *[Ü]{}ber die [Transformation]{} der konjugierten [Systeme]{} and [ü]{}ber den gemeinsamen [Ursprung]{} der [Bianchischen]{} [P]{}ermutabilit[ä]{}tstheoreme*, Sitzungsberichte Berl. Math. Ges., **14** (1915) 96–118.

V. G. Kac and J. van de Leur, *The $n$-component KP hierarchy and representation theory*, \[in:\] Important developments in soliton theory, (A. S. Fokas and V. E. Zakharov, eds.) Springer, Berlin, 1993, pp. 302–343.

R. M. Kashaev, *On discrete three-dimensional equation associated with the local Yang–Baxter relation*, Lett. Math. Phys. **38** (1996), 389–397.

A. D. King and W. K. Schief, *Tetrahedra, octahedra and cubo-octahedra: integrable geometry of multi-ratios*, J. Phys. A **36** (2003) 785–802.

A. D. King and W. K. Schief, *Application of an incidence theorem for conics: Cauchy problem and integrability of the dCKP equation*, J. Phys. A **39** (2006) 1899-1913.

A. Knutson, T. Tao, C. Woodward, *A positive proof of the Littewood–Richardson rule using the octahedron recurrence*, Electron. J. Combin. **11** (2004) Research Paper 61, `arXive:math.CO/0306274`.

B. G. Konopelchenko and W. K. Schief, *Three-dimensional integrable lattices in [Euclidean]{} spaces: Conjugacy and orthogonality*, Proc. Roy. Soc. London A **454** (1998), 3075–3104.

B. G. Konopelchenko and W. K. Schief, *Reciprocal figures, graphical statics and inversive geometry of the Schwarzian BKP hierarchy*, Stud. Appl. Math. **109** (2002) 89–124.

I. M. Krichever, *Two-dimensional periodic difference operators and algebraic geometry*, Dokl. Akad. Nauk SSSR [**285**]{} (1985) 31–36.

I. M. Krichever, P. Wiegmann, and A. Zabrodin, *Elliptic solutions to difference non-linear equations and related many body problems*, Commun. Math. Phys. **193** (1998), 373–396.

A. Kuniba, T. Nakanishi and J. Suzuki, [*Functional relations in solvable lattice models, I: Functional relations and representation theory, II: Applications*]{}, Int. Journ. Mod. Phys. A [**9**]{} (1994) 5215.

G. Lamé, [*Leçons sur les coordonnées curvilignes et leurs diverses applications*]{}, Mallet–Bachalier, Paris, 1859.

E. P. Lane, *Projective differential geometry of curves and surfaces*, Univ. Chicago Press, Chicago, 1932.

P. S. Laplace, *Recherches sur le Calcul intégral aux différences partielles*, Mémoires de Mathématique et de Physique de l’Académie des Sciences 341–403 (1773).

D. Levi, R. Benguria, *Bäcklund transformations and nonlinear differential-difference equations*, Proc. Nat. Acad. Sci. USA **77** (1980) 5025–5027.

D. Levi, L. Pilloni and P. M. Santini, [*Integrable three-dimensional lattices*]{}, J. Phys. A [**14**]{} (1981), 1567.

Q. P. Liu, M. Ma[ñ]{}as, *Discrete Levy transformations and Casorati determinant solutions for quadrilateral lattices*, Phys. Lett. A **239** (1998) 159–166.

Q. P. Liu and M. Mañas, [*Superposition Formulae for the Discrete Ribaucour Transformations of Circular Lattices*]{}, Phys. Lett. A **249** (1998) 424–430.

M. Mañas, *Fundamental transformation for quadrilateral lattices: first potentials and $\tau$-functions, symmetric and pseudo-Egorov reductions*, J. Phys. A **34** (2001) 10413–10421.

M. Mañas, *From integrable nets to integrable lattices*, J. Math. Phys. **43** (2002) 2523–2546.

M. Mañas, A. Doliwa and P.M. Santini, *Darboux transformations for multidimensional quadrilateral lattices. I*, Phys. Lett. A **232** (1997) 99–105.

V. B. Matveev, M. A. Salle, [*Darboux transformations and solitons*]{}, Springer, Berlin, 1991.

T. Miwa, *On Hirota’s difference equations*, Proc. Japan Acad. **58** (1982) 9–12.

Th-F. Moutard, [*Sur la construction des équations de la forme $\frac{1}{z} \frac{\partial^2 z}{\partial x\partial y }=\lambda (x,y)$, qui admettent une integrale général explicite*]{} J. Ec. Pol. 45 (1878) 1.

M. Nieszporski, *Laplace ladder of discrete Laplace equations*, Theor. Math. Phys. **133** (2002) 1576-1584.

M. Nieszporski, *Darboux transformations for 6-point scheme*, J. Phys. A: Math. Theor. **40** (2007) 4193–4207.

P. Ma[ł]{}kiewicz, M. Nieszporski, *Darboux transformations for $q$-discretizations of second order difference equations*, J. Nonl. Math. Phys. **12** Suppl. 2 (2005) 231–239.

M. Nieszporski, P. M. Santini and A. Doliwa, [*Darboux transformations for 5-point and 7-point self-adjoint schemes and an integrable discretization of the 2D Schrödinger operator*]{}, Phys. Lett. A **323** (2004), 241–250.

F. W. Nijhoff, *Theory of integrable three-dimension nonlinear lattice equations*, Lett. Math. Phys. **9** (1985) 235–241.

F. W. Nijhoff, H. W. Capel, *The direct linearization approach to hierarchies of integrable PDEs in $2+1$ dimensions: I. Lattice equations and the differential-difference hierarchies*, Inverse Problems **6** (1990) 567–590.

F.W. Nijhoff, G.R.W. Quispel and H.W. Capel, Direct Linearization of Nonlinear Difference-Difference Equations, Phys. Lett. [**97A**]{} (1983) 125–128.

J. J. C. Nimmo, *Darboux transformations and the discrete KP equation*, J. Phys. A: Math. Gen. **30** (1997) 8693–8704.

J. J. C. Nimmo, *Darboux transformations for discrete systems*, Chaos, Solitons and Fractals **11** (2000) 115–120.

J. J. C. Nimmo, *On a non-Abelian Hirota-Miwa equation*, J. Phys. A: Math. Gen. **39** (2006) 5053–5065.

J. J. C. Nimmo and W. K. Schief, *Superposition principles associated with the [Moutard]{} transformation. [An]{} integrable discretisation of a (2+1)-dimensional sine-[Gordon]{} system*, Proc. R. Soc. London A **453** (1997), 255–279.

J. J. C. Nimmo, R. Willox, *Darboux transformations for the two dimensional Toda system*, Proc. R. Soc. London A **453** (1997) 2497–2525.

S. P. Novikov and I. A. Dynnikov, *Discrete spectral symmetries of low dimensional differential operators and difference operators on regular lattices and two-dimensional manifolds*, Russian Math. Surveys **52** (1997), 1057–1116. S. P. Novikov and A. P. Veselov, *Exactly solvable two-dimensional Schrödinger operators and Laplace transformations*, \[in:\] *Solitons, geometry, and topology: on a crossroad*, V. M. Buchstaber and S. P. Novikov (eds.) AMS Transl. **179** (1997), 109–132.

W. Oevel, W. Schief, *Darboux theorems and the KP hierarchy*, \[in:\] P. A. Clarkson (ed.), *Application of Analytic and Geometric Methods to Nonlinear Differential Equations*, Kluwer Academic Publishers, 1993, pp. 193–206.

Y. Ohta, R. Hirota, S. Tsujimoto, T. Imai, *Casorati and discrete Gram type determinant representations of solutions to the discrete KP hierarchy*, J. Phys. Soc. Japan **62** (1993) 1872–86.

C. Rogers and W. K. Schief, *Bäcklund and Darboux transformations. Geometry and modern applications in soliton theory*, Cambridge University Press, Cambridge, 2002.

S. Saito, *String theories and Hirota’s bilinear diference equation*, Pys. Rev. Lett. **59** (1987) 1798–1801.

S. Saito, N. Saitoh, *Gauge and dual symmetries and linearization of Hirota’s bilinear equations*, J. Math. Phys. **28** (1987) 1052–1055.

N. Saitoh, S. Saito, *General solutions to the Bäcklund transformation of Hirota’s bilinear difference equation*, J.Phys. Soc. Japan **56** (1987) 1664–1674.

P. M. Santini, M. Nieszporski and A. Doliwa, *An integrable generalization of the Toda law to the square lattice*, Physical Review E [**70**]{} (2004) art. no. 056615

P. M. Santini, A. Doliwa and M. Nieszporski, [*Integrable dynamics of Toda-type on the square and triangular lattices*]{}, Phys. Rev. E [**77**]{} (2008) art. no. 056601.

R. Sauer, *Projective Liniengeometrie*, de Gruyter, Berlin–Leipzig, 1937.

R. Sauer, *Differenzengeometrie*, Springer, Berlin, 1970.

W. K. Schief, *Finite analogues of Egorov-type coordinates*, talk given at *Nonlinear Systems, Solitons and Geometry*, Mathematisches Forschungsinstitut Oberwolfach, Germany, 1997.

W. K. Schief, *Lattice geometry of the discrete Darboux, KP, BKP and CKP equations. Menelaus’ and Carnot’s theorems*, J. Nonl. Math. Phys. **10** Supplement 2 (2003) 194–208.

W. K. Schief, *On the unification of classical and novel integrable surfaces. II. Difference geometry*, R. Soc. Lond. Proc. Ser. A **459** (2003) 373–391.

T. Shiota, *Characterization of Jacobian varieties in terms of soliton equations*, Invent. Math. **83** (1986) 333–382.

Yu. B. Suris, *The problem of integrable discretization: Hamiltonian approach*, Birkhäuser, Basel, 2003.

G. Tzitzéica, *Géométrie différentielle projective des réseaux*, Cultura Naţionala, Bucarest, 1923.

J. Weiss, [*Bäcklund transformations, focal surfaces and the two-dimensional Toda lattice*]{}, Phys. Lett. A 137 (1989), 365.

R. Willox, T. Tokihiro, J. Satsuma, *Nonautonomous discrete integrable systems*, Chaos, Solitons and Fractals **11** (2000), 121–135.

R. Willox, T. Tokihiro, I. Loris, J. Satsuma, *The fermionic approach to Darboux transformations*, Inverse Problems **14** (1998), 745–762.

R. Willox, Y. Ohta, C. R. Gilson, T. Tokihiro, J. Satsuma, *Quadrilateral lattices and eigenfunction potentials for $N$-component KP hierarchies*, Phys. Lett. A **252** (1999), 163–172.

A. Zabrodin, *A survey of Hirota’s difference equations*, Theor. Math. Phys. **113** (1997) 1347–1392.
