How to use the summation in Julia?

Specific Domains Optimization (Mathematical)
,
1

I need to put this objective function in Julia, but I don’t know how to build a sum that has 3 indexes. Can someone help me?
FunObjJulia

I started doing like this

using JuMP, Cbc
Model1 = Model(with_optimizer(Cbc.Optimizer))
CR = [11; 11; 11; 11; 11]
CO = [50; 50; 50; 50; 50]

#sets
a = 80 #atividades
t = 18 #peírodo de tempo
w = 5 #centro de trabalho
#variable
@variable(Model1,r[1:a,1:t], lower_bound=0)
@variable(Model1,o[1:a,1:t], lower_bound=0)
#objective function
@objective(Model1,Min,sum(CR[w]*(r[a,t]+o[a,t])+CO[w]*o[a,t]))
1 Like
2

The JuMP repo has a number of examples that should be helpful. For example: JuMP.jl/cannery.jl at v0.21.2 · jump-dev/JuMP.jl · GitHub

For your example, you want something like this:

using JuMP
model = Model()

CR = [11; 11; 11; 11; 11]
CO = [50; 50; 50; 50; 50]

A = 1:80 #atividades
T = 1:18 #peírodo de tempo
W = 1:5 #centro de trabalho

@variable(model, r[A, T] >= 0)
@variable(model, o[A, T] >= 0)

@objective(
    model, 
    Min, 
    sum(CR[w] * (r[a, t] + o[a, t]) + CO[w] * o[a, t] for w in W, a in A, t in T)
)
2 Likes
3

Hi @RaquelSantos, I noticed from your code that you may be Brazilian. In that case (or if you can understand brazilian portuguese), I have created some tutorials, including a few for JuMP: https://www.youtube.com/playlist?list=PLOOY0eChA1uyjKRQDeR4-LfsObNqslIUh

5 Likes
4

Oi @abelsiqueira, vou olhar sim.
Obrigada pela ajuda!!

5

@odow, @abelsiqueira Thanks for the help!

I tried to combine all the equations in a single function as in the image, only the first two sums were compiled.

FOJulia4
FOJulia4791×128 14.6 KB

I started doing like this

using JuMP, Cbc
Model1 = Model(with_optimizer(Cbc.Optimizer))
CR = [11; 11; 11; 11; 11]
CO = [50; 50; 50; 50; 50]
PS = [4951400; 4951400; 14154430]

A = 1:80 #atividades
T = 1:18 #peírodo de tempo
W = 1:5 #centro de trabalho
I = 1:5 #centro de trabalho 

@variable(Model1, r[A,T] >= 0)
@variable(Model1, o[A,T] >= 0)
@variable(Model1,wr[W,T] >= 0)
@variable(Model1,wo[W,T] >= 0) 
@variable(Model1,d[I] >= 0) 

@objective(Model1,Min,sum(CR[w]*(r[a,t]+o[a,t])+CO[w]*o[a,t] for w in W, a in A, t in T)+sun(CR[w]*(wr[w,t]+wo[w,t])+CO[w]*wo[w,t] for w in W, t in T)+sun(PS[i]*(1-d[i]) for i in I))
6

it looks like you accidently wrote sun instead of sum for the third one

4 Likes