Differences
This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
personal:blog:2017:0203_jump_for_gams_users [2017/02/07 10:26] antonello [Declaration of the model] |
personal:blog:2017:0203_jump_for_gams_users [2023/12/22 11:12] antonello [Complete script] |
||
---|---|---|---|
Line 236: | Line 236: | ||
<code Julia> | <code Julia> | ||
+ | # Transport example | ||
+ | |||
# Transposition in JuMP of the basic transport model used in the GAMS tutorial | # Transposition in JuMP of the basic transport model used in the GAMS tutorial | ||
# | # | ||
Line 247: | Line 249: | ||
# The Scientific Press, Redwood City, California, 1988. | # The Scientific Press, Redwood City, California, 1988. | ||
# - JuMP implementation: | # - JuMP implementation: | ||
- | + | ||
- | using JuMP, DataFrames | + | using CSV, DataFrames, GLPK, JuMP |
+ | |||
# Sets | # Sets | ||
plants | plants | ||
markets = [" | markets = [" | ||
+ | |||
# Parameters | # Parameters | ||
a = Dict( # capacity of plant i in cases | a = Dict( # capacity of plant i in cases | ||
Line 264: | Line 266: | ||
" | " | ||
) | ) | ||
+ | |||
# distance in thousands of miles | # distance in thousands of miles | ||
- | d_table = wsv""" | + | d_table = CSV.read(IOBuffer(""" |
plants | plants | ||
seattle | seattle | ||
san_diego | san_diego | ||
- | """ | + | """ |
d = Dict( (r[: | d = Dict( (r[: | ||
+ | |||
f = 90 # freight in dollars per case per thousand miles | f = 90 # freight in dollars per case per thousand miles | ||
+ | |||
c = Dict() # transport cost in thousands of dollars per case ; | c = Dict() # transport cost in thousands of dollars per case ; | ||
[ c[p,m] = f * d[p,m] / 1000 for p in plants, m in markets] | [ c[p,m] = f * d[p,m] / 1000 for p in plants, m in markets] | ||
+ | |||
# Model declaration | # Model declaration | ||
- | trmodel = Model() # transport model | + | trmodel = Model(GLPK.Optimizer) # transport model |
+ | |||
# Variables | # Variables | ||
@variables trmodel begin | @variables trmodel begin | ||
x[p in plants, m in markets] >= 0 # shipment quantities in cases | x[p in plants, m in markets] >= 0 # shipment quantities in cases | ||
end | end | ||
+ | |||
# Constraints | # Constraints | ||
@constraints trmodel begin | @constraints trmodel begin | ||
Line 293: | Line 295: | ||
sum(x[p,m] for p in plants) | sum(x[p,m] for p in plants) | ||
end | end | ||
+ | |||
# Objective | # Objective | ||
@objective trmodel Min begin | @objective trmodel Min begin | ||
sum(c[p, | sum(c[p, | ||
end | end | ||
+ | |||
print(trmodel) | print(trmodel) | ||
+ | |||
+ | optimize!(trmodel) | ||
+ | status = termination_status(trmodel) | ||
- | status = solve(trmodel) | + | if status == MOI.OPTIMAL |
- | + | println(" | |
- | if status == :Optimal | + | |
- | println(" | + | |
println(" | println(" | ||
- | println(getvalue(x)) | + | println(value.(x)) |
println(" | println(" | ||
- | [println(" | + | [println(" |
println(" | println(" | ||
- | [println(" | + | [println(" |
+ | |||
else | else | ||
println(" | println(" | ||
println(status) | println(status) | ||
end | end | ||
- | + | ||
# Expected result: | # Expected result: | ||
# obj= 153.675 | # obj= 153.675 |