Linear Programming with Python (Pulp) - Scheduling

Question

I'm trying to optimize a schedule by formulating it as an (integer) linear problem. I've run into some problems with the coding part (mostly due to Pulp), which is fairly separate from the formulation.

The Problem:

Schedule n ships over a period of 5 years (timesteps in months), subject to the following constraints:

one and only one ship must be on "Local Patrol" (state = 'C') each month
one and only one ship must go on "Extended Patrol" (state = 'A') each year, which is a four month block from July to October inclusive
each ship must have a 4 month long inspection once during the 5 year period (state = '4')
each ship must have an 8 month long break once during the 5 year period (state = '8')
each ship must have regular maintenance approx every 5 months (state = 'M')
a ship can be in only one state during a given month, but it can also have no state (open/vacant, '_')

Objective Function formulation:

Define slack variables that count the number of times (and the amount) the schedule deviates from the 5 month ('M') maintenance plan.

I can explain the setup more if needed, but for now I have some basic questions that might be the source of my problems.

Here is one of my constraints (blockList stores the start of the 4 and 8 month block for each ship):

# Each ship must have one 8block and one m4block in a 5 year schedule, which must be between 2 and 3 years apart
for n in ship_list:
    prob += blockList[n][0] - blockList[n][1] >= 24
    prob += blockList[n][0] - blockList[n][1] <= 36

I want to take the absolute value of the LHS of the above expressions, but I can't because the objects are initialized as pulp.LpVariable. Any recommendations for a workaround?

Anyways, that may or may not fix my problem. If it doesn't, my next question is for anyone familiar with pulp. How can I know if the solution true, or if there was some sort of elasticity added by the solver in the background? (Running the code as is with that constraint gives incorrect results, and using a different solver like GLPK gives an error).

pchtsp · Accepted Answer · 2020-06-08 02:18:38Z

1

Here are a few suggestions:

Taking Absolute Values

In PuLP, you can do that by relying on the underlying Python commands to do that.

Something along the lines of:

if blockList[n][0] > blockList[n][1]:
    prob += blockList[n][0] - blockList[n][1] >= 24
else:
    prob += blockList[n][1] - blockList[n][0] >= 24

To see if the Solver that to "add elasticity" you'd examine the value of pulp.constants.LpStatusOptimal. If this value is 1, then the problem was solved to optimality. Note that the typical practice is to add a dummy slack or a surplus variable, give it a small penalty in the objective function. If in the solution that dummy variable is non-zero, that means that the problem needed the extra 'elasticity.'
Finally, my main suggestion is that you start small, with a minimal example that works in PuLP. You can even start with one of the case studies here. At each step, write down your LP to a file and examine it. This will right away tell you which constraints or variables are off, and you can alter that.
```
class pulp.LpProblem(name='outfile.lp', sense=1)
```

Hope this gets you moving forward.

edited Jun 8, 2020 at 2:18

pchtsp

8741 gold badge7 silver badges20 bronze badges

answered Sep 13, 2013 at 21:46

Ram Narasimhan

22.6k5 gold badges54 silver badges56 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

user2777139 Over a year ago

Thanks Ram, I ended up using cases to get it to work. I wasn't able to use a generic ">" operation though.

KHAN irfan Over a year ago

Hi @Ram. Kindly answer this question datascience.stackexchange.com/questions/49792/…

Collectives™ on Stack Overflow

Linear Programming with Python (Pulp) - Scheduling

1 Answer 1

2 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

2 Comments

Your Answer

Sign up or log in

Post as a guest

Related