How do I model a band matrix in GAMS?

A band matrix can be modeled efficiently in GAMS using circular leads (++) and circular lags (–), see chapter lag and lead operator of the GAMS User’s Guide for more information. Using circular leads and lags, the first and last members of the set are assumed to be adjacent, so as to form a circular sequence of members. These operators only work with ordered sets.

Find below an example, where the generation of a model with the usage of circular leads (equation *eq*) was ~50 times faster than the original formulation using *sum* - statements and lots of comparisons (equation *eq1* and *eq2*). Both models are identical.

```
set k / k1*k300/
t / t1*t700/;
alias (t,tt,ll,t2);
parameter l(k); l(k) = uniformInt(1,6);
variable y(k,t), z(k,t), obj;
equation defobj;
defobj.. obj =e= 0;
set ktl(k,t),
kl(k,ll);
kl(k,ll) = ord(ll) <= l(k);
equation eq;
eq(k,t).. sum(kl(k,ll), y(k,t++(ord(ll)-1))) =l= z(k,t++(l(k)-1));
model m1 /eq, defobj/;
solve m1 min obj using lp;
*******************************************************************
*original formulation
equation eq1, eq2;
eq1(k,t)$(ord(t)<=card(t)-l(k)+1)..
sum(t2$((ord(t2)>=ord(t)) and(ord(t2)<ord(t)+l(k))),y(k,t2))
=l=
sum(t2$(ord(t2)=l(k)+ord(t)-1),Z(k,t2));
eq2(k,t)$(ord(t)> card(t)-l(k)+1) ..
sum(t2$(ord(t2)>=ord(t)),Y(k,t2)) +
sum(t2$(ord(t2)<l(k)-card(t)+ord(t)),Y(k,t2))
=l=
sum(t2$(ord(t2)=l(k)-card(t)+ord(t)-1),Z(k,t2));
model m2 /defobj, eq1, eq2/;
solve m2 min obj using lp;
```