Lax equations for relativistic ${\rm GL}(NM,{\mathbb C})$ Gaudin models on elliptic curve

We describe the most general ${\rm GL}_{NM}$ classical elliptic finite-dimensional integrable system, which Lax matrix has $n$ simple poles on elliptic curve. For $M=1$ it reproduces the classical inhomogeneous spin chain, for $N=1$ it is the Gaudin type (multispin) extension of the spin Ruijsenaars