How to get a length of a set with condition

Problems with syntax of GAMS
Post Reply
Posts: 1
Joined: 1 year ago

How to get a length of a set with condition

Post by sadlittlebunny95 » 1 year ago

Hello everyone,

I am new to GAMS and I have a question about set in GAMS.
I defined sets: i and k. I also defined a variable z(i,k) where z(i,k) is a binary variable in which z.lo(i,k) = 0 and z.up(i,k) = 1.

set i 'samples' /1*10/;
set k 'components' /k1*k3/;
binary variable z(i,k);

I have another sets fz(i,k), fixedzero(i,k), and fixedone(i,k) which is used to fixed values of z(i,k). During the iterative solving, some z(i,k) are fixed to one, some are fixed to zero.
set fz(i,k);
set fixedzero(i,k);
set fixedone(i,k);
fz(i,k) = fixedzero(i,k);
z.fx(fz) = 0;
fz(i,k) = fixedone(i,k);
z.fx(fz) = 1;

In the end, I would like to find the size of fixed z(i,k) for a given k, i.e. given that k = 'k1', how many z(i,'k1') are fixed to one and how many are fixed to zero. For example, if index i=1,2,3 of z(i,'k1') is fixed to one, then count_fixedone('k1') should be 3. Here is how I tried, but it does not work.

set count_fixedone(k);
set count_fixedzero(k);

count_fixedone(k) = card( i$(z.l(i,k)=1) )
* if z(i,k) is fixed to one, its lower bound is 1
display count_fixedone;

This is the errors:
8 ')' expected
149 Uncontrolled set entered as constant

Thank you very much!

User avatar
Posts: 585
Joined: 3 years ago

Re: How to get a length of a set with condition

Post by Renger » 1 year ago

If I understand your question well, you just want to count the number of variables fixed to 1. The following code sums 1 every time the condition that Z is fixed to 1 is met (upperbound and lower bound equal to 1):

Code: Select all

set i 'samples' /1*10/;
set k 'components' /k1*k3/;
binary variable z(i,k);

z.l(i,"k1") = 0;
z.UP(i,"k2") = 1;
z.Lo(i,"k2") = 1;

parameter  count_fixedone(k);

count_fixedone(k) = sum(i$(Z.UP(i,k) =1 AND Z.LO(i,k) = 1), 1);

display count_fixedone;

PS. You can post properly formatted code by klicking the icon </> in the menu. That saves the person to answer your question some work when trying to reproduce your problem.
Enjoy modeling even more: Read my blog on modeling at The lazy economist

Post Reply