Re: [Help-glpk] Set pair e.g. ( (i, j) in T) not properly implemented in matrix generation (i.e. CPLEX LP formulation, etc.) in particular data set (e.g. repeating i's in (i, j) in T

[Noli Sicad]

Hi Andrew,

> I see nothing wrong. For example, there is
> Adjacent_Constraint(1,34,19,6) generated for (34,19) in ADJACENT which
> is listed below (33,21). For (33,22), (33,23), etc., corresponding
> Adjacent_Constraints are not generated because the condition
> isHarvestable[t,i1,j] and isHarvestable[t,i2,j] specified for
> Adjacent_Constraint is not met.

I think the condition isHarvestable[t,i1,j] and isHarvestable[t,i2,j]
specified for Adjacent_Constraint is met for stands (i1,i2) (e.g.
(22,33), (23, 33))

Period 1
   (1,21,5)
   (1,22,7)
   (1,23,7)
   (1,27,5)
   (1,31,5)
   (1,32,6)
   (1,33,5)

Period 2
   (2,21,6)
   (2,22,8)
   (2,23,8)
   (2,26,5)
   (2,27,6)
   (2,31,6)
   (2,32,7)
   (2,33,6)

This can be verify by placing this script in the MathProg model

display HARVEST;

The result is:

#####################################################################################
Display statement at line 138
HARVEST:
   (1,1,6)
   (1,2,5)
   (1,3,5)
   (1,10,6)
   (1,11,6)
   (1,17,5)
   (1,18,6)
   (1,19,6)
   (1,20,6)
   (1,21,5)
   (1,22,7)
   (1,23,7)
   (1,27,5)
   (1,31,5)
   (1,32,6)
   (1,33,5)
   (1,34,6)
.
.
.
  (2,20,7)
   (2,21,6)
   (2,22,8)
   (2,23,8)
   (2,26,5)
   (2,27,6)
   (2,31,6)
   (2,32,7)
   (2,33,6)
   (2,34,7)
   (2,35,7)
.
##########

This is the CPLEX LP conversion.

 Adjacent_Constraint(1,20,34,6): + Y(1,20,6) + Y(1,34,6) <= 1
 Adjacent_Constraint(1,20,77,6): + Y(1,20,6) + Y(1,77,6) <= 1
 Adjacent_Constraint(1,21,33,5): + Y(1,21,5) + Y(1,33,5) <= 1
 Adjacent_Constraint(1,32,72,6): + Y(1,32,6) + Y(1,72,6) <= 1
 Adjacent_Constraint(1,32,73,6): + Y(1,32,6) + Y(1,73,6) <= 1
 Adjacent_Constraint(1,33,21,5): + Y(1,21,5) + Y(1,33,5) <= 1
 Adjacent_Constraint(1,34,19,6): + Y(1,19,6) + Y(1,34,6) <= 1
 Adjacent_Constraint(1,34,20,6): + Y(1,20,6) + Y(1,34,6) <= 1
 Adjacent_Constraint(1,34,65,6): + Y(1,34,6) + Y(1,65,6) <= 1
 Adjacent_Constraint(1,34,77,6): + Y(1,34,6) + Y(1,77,6) <= 1
 Adjacent_Constraint(1,35,38,6): + Y(1,35,6) + Y(1,38,6) <= 1
 Adjacent_Constraint(1,35,39,6): + Y(1,35,6) + Y(1,39,6) <= 1

There is no Adjacent_Constraint for

   (1,21,5)
   (1,22,7)
   (1,23,7)
   (1,27,5)
   (1,31,5)
   (1,32,6)


Noli