How do I model a band matrix?

[aileen]

How do I model a band matrix in GAMS?

[aileen]

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.

Code: Select all

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;