Re: [Help-glpk] forall constraint in MathProg

[Noli Sicad]

Hello Xypron,

It seems that there is a still a problem.

Upon closer inspection on the constraints generated by this statement,

####
s.t. habitat{p in PERIOD, n in NATURE, s in STAND} :
  sum{pp in PERIOD:pp <= p} forest[s, pp] + nh[n, p] = standarea[s];
####

It generates 300 constraints instead of 60 (i.e. 3 x 20) constraints.

Something wrong in  nh[n, p].

Again, the desire out is like this,

 S1P1 + N1P1 = 12.3
 S1P1 + S1P2 + N1P2 = 12.3
 S1P1 + S1P2 + S1P3 + N1P3 = 12.3

S2P1 + N2P1 = 35.6
S2P1 + S2P2 + N2P2 = 35.6
S2P1 + S2P2 + S2P3 + N2P3 = 34.6



It seems that it is doing like this.

S1P1 + N1P1 = 12.3
S1P1 + S1P2 + N1P2 = 12.3
S1P1 + S1P2 + S1P3 + N1P3 = 12.3

S1P1 + N2P1 = 12.3
 S1P1 + S1P2 + N2P2 = 12.3
 S1P1 + S1P2 + S1P3 + N2P3 = 12.3

This is the actual output:

Subject To
 habitat(1,1,1): + forest(1,1) + nh(1,1) = 5.64155902015045
 habitat(1,1,2): + forest(2,1) + nh(1,1) = 5.04798337584361
 habitat(1,1,3): + forest(3,1) + nh(1,1) = 5.14843434561044


This is not needed (below)

 habitat(1,2,1): + forest(1,1) + nh(2,1) = 5.64155902015045
 habitat(1,2,2): + forest(2,1) + nh(2,1) = 5.04798337584361
 habitat(1,2,3): + forest(3,1) + nh(2,1) = 5.14843434561044
 habitat(1,3,1): + forest(1,1) + nh(3,1) = 5.64155902015045
 habitat(1,3,2): + forest(2,1) + nh(3,1) = 5.04798337584361
 habitat(1,3,3): + forest(3,1) + nh(3,1) = 5.14843434561044
 habitat(1,4,1): + forest(1,1) + nh(4,1) = 5.64155902015045
 habitat(1,4,2): + forest(2,1) + nh(4,1) = 5.04798337584361
 habitat(1,4,3): + forest(3,1) + nh(4,1) = 5.14843434561044
 habitat(1,5,1): + forest(1,1) + nh(5,1) = 5.64155902015045
 habitat(1,5,2): + forest(2,1) + nh(5,1) = 5.04798337584361
 habitat(1,5,3): + forest(3,1) + nh(5,1) = 5.14843434561044


Thanks.

Noli


On 11/6/12, Noli Sicad <address@hidden> wrote:
> Hello Xypron,
>
> Yes, this is what I wanted. All the constraints are there.
>
> Nice tip as well with the param and set combo (below).
>
>> param nperiod := 3;
>> param nnature := 20;
>> param nstand := 20;
>
>> set PERIOD := {1..nperiod};
>> set NATURE := {1..nnature};
>> set STAND := {1..nperiod};
>
> Now, MathProg has easy way of implementing 'forall' constraints.
>
> Danke Schoen Xypron.
>
> It might be good to be added in the GLPK Wikibook.
>
> Robbie, would you like to add this in Wikibook?
>
> Thanks again.
>
> Regards,
>
> Noli
>
> On 11/6/12, Xypron <address@hidden> wrote:
>> Hello Noli,
>>
>> could you provide the condition you want to formulate in mathematical
>> notation?
>>
>> "forall" is used for logical expression. Did you mean "sum"?
>> Instead of "p in 1..PERIOD" you may want to write "p in {1..PERIOD}".
>> If p is a number {pp in p} does not make sense.
>>
>> Please, save the gmpl text below as test.mod, and run command.
>> glpsol -m test.mod --check --wlp test.lp
>>
>> Check if test.lp contains the constraints you meant to write.
>>
>> Best regards
>>
>> Xypron
>>
>> # file test.lp
>> param nperiod := 3;
>> param nnature := 20;
>> param nstand := 20;
>
>> set PERIOD := {1..nperiod};
>> set NATURE := {1..nnature};
>> set STAND := {1..nperiod};
>
>> param standarea {STAND} := 5 + 5 * Uniform01();
>> var forest{STAND, PERIOD};
>> var nh{NATURE, PERIOD};
>> s.t. habitat{p in PERIOD, n in NATURE, s in STAND} :
>>   sum{pp in PERIOD:pp <= p} forest[s, pp] + nh[n, p] = standarea[s];
>> end;
>>
>>
>> On 05.11.2012 19:29, Noli Sicad wrote:
>>> Hi,
>>>
>>> I am trying to translate this set of LINDO constraint (below) into
>>> MathProg. But I don't have any success.
>>>
>>> ######
>>>
>>> 37) S1P1 + N1P1 = 12.3
>>> 38) S1P1 + S1P2 + N1P2 = 12.3
>>> 39) S1P1 + S1P2 + S1P3 + N1P3 = 12.3
>>>
>>> 40) S2P1 + N2P1 = 35.6
>>> 41) S2P1 + S2P2 + N2P2 = 35.6
>>> 42) S2P1 + S2P2 + S2P3 + N2P3 = 34.6
>>>
>>> 43) S3P1 + N3P1 = 34.6
>>> 44) S3P1 + S3P2 + N3P2 = 34.6
>>> 45) S3P1 + S3P2 + S3P3 + N3P3 = 34.6
>>>
>>> .
>>> .
>>> .
>>>
>>> 94) S20P1 + N20P1 = 46.0
>>> 95) S20P1 + S20P2 + N20P2 = 46.0
>>> 96) S20P1 + S20P2 + S20P3 + N20P3 = 46.0
>>>
>>> ###########
>>>
>>> I tried to implement like this:
>>>
>>> subject to Habitat_Constraint{ p in 1..PERIOD, n in 1..NATURE, s in
>>> STAND}:
>>>    forall {pp in p} Forest[s,pp] + NH[n,p] = StandArea[s];
>>>
>>> However, I got this error below.
>>>
>>> #######
>>> Running tool: GLPK/MathProg
>>>
>>> GLPSOL: GLPK LP/MIP Solver, v4.47
>>> Parameter(s) specified in the command line:
>>>  --math Model_I_1_original_5-evenflow_even-acres_OF_alter3.mod
>>> Reading model section from
>>> Model_I_1_original_5-evenflow_even-acres_OF_alter3.mod...
>>> Model_I_1_original_5-evenflow_even-acres_OF_alter3.mod:77: domain
>>> expression has invalid type
>>> Context: ...ERIOD , n in 1 .. NATURE , s in STAND } : forall { pp in p }
>>> MathProg model processing error
>>>
>>> ########
>>>
>>> Anybody knows how to do this set of constraint properly in MathProg
>>> properly?
>>>
>>> Any suggestions?
>>>
>>> Thanks.
>>>
>>> Regards,
>>>
>>> Noli
>>>
>>> _______________________________________________
>>> Help-glpk mailing list
>>> address@hidden
>>> https://lists.gnu.org/mailman/listinfo/help-glpk
>>>
>>
>>
>