@e-zolotarevskaya I cannot assign you. No idea why.

Are you planning to build a stochastic programming API? I am not too much familiar with SP, however, I think it would require a whole new type of model (similar to the ExpansionModel)?

With regard to the looping over blocks we do this here already for the "normal" decision variables to create our results dataframe, maybe this helps?

That is extremely helpful indeed! @e-zolotarevskaya, have a look! Right now we hope that we can do this without the need for a full additional model. This is because the mpi-sppy package is built to work on top of (variations of) vanilla pyomo models. So we hope that we can work on top of vanilla oemof models. We'll see.

Sorry for the late response to this one:
I'm not so very familiar with stochastical programming either, but I think, it could be quite similar to the multi-period modeling I have just started to revisit, see #793 resp. #750.

I've discussed parts of it with @p-snft quite a while ago:

• Basically, you would have to index everything with a dedicated scenario index, i.e. to create another set in the models class. As @simnh suggested, you could create a dedicated model class for that. In addition, you usually would define a weighting, i.e. a propability distribution for your uncertain scenario parameters.
• An alternative could be to integrate everything into the existing (models and class) infrastructure and make the single-scenario consideration a special case which would require less duplicates, but might slightly impact performance. I think, the latter might not be much of an issue, but I haven't yet had the possibility to check that.
• Anyways, I think, this would be a very great feature for the community, but you propably have to touch quite a lot and it's probably not possible or at least not straightforward to use third-party packages for the purpose of stochastic modeling.

I agree that we definitely should discuss it in the upcoming meeting oemof/oemof#96

After discussing it at the dev meeting oemof/oemof#96, this should be easier to add than I originally assumed. Here is a try to summarize. @simnh and @gnn: Please correct me if I'm wrong and please add to it.

• For your special use case, trying to add a (rather generalized) component that encapsulates your constraints would be the first step.
• We basically have a flow variable in the model (oemof.solph.models.Model) that it indexed by inflow, outflow and a timestep model.flow[i, o, t].
• You have to filter the flows for the ones that shall be fixed in the first stage. You could do so, by introducing an attribute, such as stochastic and then use a constraint grouping to group the respective flows.
• You have to interchange the "standard" pyomo by mpi-sppy and provide it with the list of flows that shall be fixed.
• Until our pyomo version issue isn't resolve, you'd have to use an earlier version of mpi-sppy that is compatible with pyomo version >= 5.7.0 and < 5.7.3