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 _______________________________________________ Help-glpk mailing list [email protected] https://lists.gnu.org/mailman/listinfo/help-glpk
