Differences

This shows you the differences between two versions of the page.

--- personal:blog:2017:0203_jump_for_gams_users [2018/06/18 15:11]
127.0.0.1 external edit
+++ personal:blog:2017:0203_jump_for_gams_users [2023/12/22 11:12]
antonello [Complete script]
@@ Line 236: / Line 236: @@
 <code Julia>
+# Transport example
 # Transposition in JuMP of the basic transport model used in the GAMS tutorial
 #
@@ Line 247: / Line 249: @@
 # The Scientific Press, Redwood City, California, 1988.
 # - JuMP implementation: Antonello Lobianco
-using JuMP, DataFrames
+using CSV, DataFrames, GLPK, JuMP
 # Sets
 plants  = ["seattle","san_diego"]          # canning plants
 markets = ["new_york","chicago","topeka"]  # markets
 # Parameters
 a = Dict(              # capacity of plant i in cases
@@ Line 264: / Line 266: @@
   "topeka"    => 275,
 )
 #  distance in thousands of miles
-d_table = wsv"""
+d_table = CSV.read(IOBuffer("""
 plants     new_york  chicago  topeka
 seattle    2.5       1.7      1.8
 san_diego  2.5       1.8      1.4
-"""
+"""), DataFrame, delim=" ", ignorerepeated=true,copycols=true)
 d = Dict( (r[:plants],m) => r[Symbol(m)] for r in eachrow(d_table), m in markets)
 f = 90 # freight in dollars per case per thousand miles
 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]
 # Model declaration
-trmodel = Model() # transport model
+trmodel = Model(GLPK.Optimizer) # transport model
 # Variables
 @variables trmodel begin
     x[p in plants, m in markets] >= 0 # shipment quantities in cases
 end
 # Constraints
 @constraints trmodel begin
@@ Line 293: / Line 295: @@
         sum(x[p,m] for p in plants)  >=  b[m]
 end
 # Objective
 @objective trmodel Min begin
     sum(c[p,m]*x[p,m] for p in plants, m in markets)
 end
 print(trmodel)
+optimize!(trmodel)
+status = termination_status(trmodel)
-status = solve(trmodel)
+if status == MOI.OPTIMAL
+    println("Objective value: ", objective_value(trmodel))
-if status == :Optimal
-    println("Objective value: ", getobjectivevalue(trmodel))
     println("Shipped quantities: ")
-    println(getvalue(x))
+    println(value.(x))
     println("Shadow prices of supply:")
-    [println("$p = $(getdual(supply[p]))") for p in plants]
+    [println("$p = $(dual(supply[p]))") for p in plants]
     println("Shadow prices of demand:")
-    [println("$m = $(getdual(demand[m]))") for m in markets]
+    [println("$m = $(dual(demand[m]))") for m in markets]
 else
     println("Model didn't solved")
     println(status)
 end
 # Expected result:
 # obj= 153.675