Pyomo sum inside sum with various index
I would like to do this:
but I don't know exacly how, this is what I've tried:
def R2(model, da, gr, cu):
return sum(sum(model.Dap[asi, pe] for asi in model.A[asi, gr, cu]) for pe in model.Pd[pe, da]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
pyomo
asked Nov 18 at 23:06
Moltkerrimo
11
add a comment |
I would like to do this:
but I don't know exacly how, this is what I've tried:
def R2(model, da, gr, cu):
return sum(sum(model.Dap[asi, pe] for asi in model.A[asi, gr, cu]) for pe in model.Pd[pe, da]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
pyomo
asked Nov 18 at 23:06
Moltkerrimo
11
add a comment |
I would like to do this:
but I don't know exacly how, this is what I've tried:
def R2(model, da, gr, cu):
return sum(sum(model.Dap[asi, pe] for asi in model.A[asi, gr, cu]) for pe in model.Pd[pe, da]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
pyomo
asked Nov 18 at 23:06
Moltkerrimo
11
I would like to do this:
but I don't know exacly how, this is what I've tried:
def R2(model, da, gr, cu):
return sum(sum(model.Dap[asi, pe] for asi in model.A[asi, gr, cu]) for pe in model.Pd[pe, da]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
pyomo
pyomo
asked Nov 18 at 23:06
Moltkerrimo
11
asked Nov 18 at 23:06
Moltkerrimo
11
edited Nov 18 at 23:45
asked Nov 18 at 23:06
Moltkerrimo
11
asked Nov 18 at 23:06
Moltkerrimo
11
asked Nov 18 at 23:06
Moltkerrimo
11
11
add a comment |
add a comment |
1 Answer
1
active
oldest
votes
I am not 100% I understand your mathematical formulation.
Based on the way you write the constraint, it seems you have multiple Indexed Sets (i.e. sets which are themselves indexed over other sets), which is often a symptom that the formulation is a bit clunky, although it is sometimes necessary.
If what you are trying to express is a system of constraints, i.e. one constraint for each element of DA, GR, CU, and the sets Pd and A are indexed over DA, GR, CU:
Then you should write
def R2(model, d, g, c):
return sum(sum(model.Dap[p] for a in model.A[g, c])
for p in model.Pd[d]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
Otherwise, you should try to make the formulation more clear.
answered Nov 21 at 1:01
Giorgio Balestrieri
645
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53366334%2fpyomo-sum-inside-sum-with-various-index%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
I am not 100% I understand your mathematical formulation.
Based on the way you write the constraint, it seems you have multiple Indexed Sets (i.e. sets which are themselves indexed over other sets), which is often a symptom that the formulation is a bit clunky, although it is sometimes necessary.
If what you are trying to express is a system of constraints, i.e. one constraint for each element of DA, GR, CU, and the sets Pd and A are indexed over DA, GR, CU:
Then you should write
def R2(model, d, g, c):
return sum(sum(model.Dap[p] for a in model.A[g, c])
for p in model.Pd[d]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
Otherwise, you should try to make the formulation more clear.
answered Nov 21 at 1:01
Giorgio Balestrieri
645
add a comment |
I am not 100% I understand your mathematical formulation.
Based on the way you write the constraint, it seems you have multiple Indexed Sets (i.e. sets which are themselves indexed over other sets), which is often a symptom that the formulation is a bit clunky, although it is sometimes necessary.
If what you are trying to express is a system of constraints, i.e. one constraint for each element of DA, GR, CU, and the sets Pd and A are indexed over DA, GR, CU:
Then you should write
def R2(model, d, g, c):
return sum(sum(model.Dap[p] for a in model.A[g, c])
for p in model.Pd[d]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
Otherwise, you should try to make the formulation more clear.
answered Nov 21 at 1:01
Giorgio Balestrieri
645
add a comment |
I am not 100% I understand your mathematical formulation.
Based on the way you write the constraint, it seems you have multiple Indexed Sets (i.e. sets which are themselves indexed over other sets), which is often a symptom that the formulation is a bit clunky, although it is sometimes necessary.
If what you are trying to express is a system of constraints, i.e. one constraint for each element of DA, GR, CU, and the sets Pd and A are indexed over DA, GR, CU:
Then you should write
def R2(model, d, g, c):
return sum(sum(model.Dap[p] for a in model.A[g, c])
for p in model.Pd[d]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
Otherwise, you should try to make the formulation more clear.
answered Nov 21 at 1:01
Giorgio Balestrieri
645
I am not 100% I understand your mathematical formulation.
Based on the way you write the constraint, it seems you have multiple Indexed Sets (i.e. sets which are themselves indexed over other sets), which is often a symptom that the formulation is a bit clunky, although it is sometimes necessary.
If what you are trying to express is a system of constraints, i.e. one constraint for each element of DA, GR, CU, and the sets Pd and A are indexed over DA, GR, CU:
Then you should write
def R2(model, d, g, c):
return sum(sum(model.Dap[p] for a in model.A[g, c])
for p in model.Pd[d]) <= 1
model.R2 = Constraint(model.DA, model.GR, model.CU, rule=R2)
Otherwise, you should try to make the formulation more clear.
answered Nov 21 at 1:01
Giorgio Balestrieri
645
answered Nov 21 at 1:01
Giorgio Balestrieri
645
answered Nov 21 at 1:01
Giorgio Balestrieri
645
answered Nov 21 at 1:01
Giorgio Balestrieri
645
645
add a comment |
add a comment |
Thanks for contributing an answer to Stack Overflow!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Some of your past answers have not been well-received, and you're in danger of being blocked from answering.
Please pay close attention to the following guidance:
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53366334%2fpyomo-sum-inside-sum-with-various-index%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
