oemof.solph 0.5.2.dev0__tar.gz → 0.5.2.dev2__tar.gz

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/.bumpversion.cfg

@@ -1,5 +1,5 @@
 [bumpversion]
-current_version = 0.5.2dev0
+current_version = 0.5.2.dev2
 commit = True
 tag = True
 

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/.readthedocs.yaml

@@ -1,13 +1,15 @@
 # See https://docs.readthedocs.io/en/stable/config-file/v2.html for details
 version: 2
 build:
+  os: ubuntu-22.04
+  tools:
+    python: "3.8"
   apt_packages:
     - coinor-cbc
 sphinx:
   configuration: docs/conf.py
 formats: []
 python:
-  version: "3.8"
   install:
     - requirements: docs/requirements.txt
     - method: pip

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/MANIFEST.in

@@ -5,7 +5,6 @@ graft tests
 
 include .bumpversion.cfg
 include .coveragerc
-include .cookiecutterrc
 include .editorconfig
 
 include AUTHORS.rst

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: oemof.solph
-Version: 0.5.2.dev0
+Version: 0.5.2.dev2
 Summary: A model generator for energy system modelling and optimisation.
 Home-page: https://oemof.org
 Author: oemof developer group
@@ -9,7 +9,6 @@ License: MIT
 Project-URL: Documentation, https://oemof-solph.readthedocs.io/
 Project-URL: Changelog, https://oemof-solph.readthedocs.io/en/latest/changelog.html
 Project-URL: Issue Tracker, https://github.com/oemof/oemof-solph/issues
-Platform: UNKNOWN
 Classifier: Development Status :: 5 - Production/Stable
 Classifier: Intended Audience :: Developers
 Classifier: Intended Audience :: Science/Research
@@ -28,11 +27,29 @@ Classifier: Programming Language :: Python :: Implementation :: CPython
 Classifier: Topic :: Utilities
 Requires-Python: >=3.8
 Description-Content-Type: text/x-rst
+License-File: LICENSE
+License-File: AUTHORS.rst
+Requires-Dist: blinker
+Requires-Dist: dill
+Requires-Dist: numpy
+Requires-Dist: pandas>=1.5.3
+Requires-Dist: pyomo<7.0,>=6.6.0
+Requires-Dist: networkx
+Requires-Dist: oemof.tools>=0.4.2
+Requires-Dist: oemof.network>=0.5.0a1
 Provides-Extra: dev
+Requires-Dist: matplotlib; extra == "dev"
+Requires-Dist: nbformat; extra == "dev"
+Requires-Dist: pytest; extra == "dev"
+Requires-Dist: sphinx; extra == "dev"
+Requires-Dist: sphinx_rtd_theme; extra == "dev"
+Requires-Dist: sphinx-copybutton; extra == "dev"
+Requires-Dist: sphinx-design; extra == "dev"
+Requires-Dist: termcolor; extra == "dev"
 Provides-Extra: examples
+Requires-Dist: matplotlib; extra == "examples"
 Provides-Extra: dummy
-License-File: LICENSE
-License-File: AUTHORS.rst
+Requires-Dist: oemof; extra == "dummy"
 
 
 |tox-pytest| |tox-checks| |appveyor| |coveralls| |codecov|
@@ -327,5 +344,3 @@ AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
 SOFTWARE.
-
-

oemof.solph-0.5.2.dev2/VERSION

@@ -0,0 +1 @@
+__version__ = "0.5.2.dev2"

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/docs/conf.py

@@ -36,7 +36,7 @@ project = "oemof.solph"
 year = "2014-2023"
 author = "oemof-developer-group"
 copyright = "{0}, {1}".format(year, author)
-version = release = "0.5.2dev0"
+version = release = "0.5.2.dev2"
 
 pygments_style = "trac"
 templates_path = ["."]

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/docs/examples/index.rst

@@ -29,5 +29,4 @@ Examples
     minimal_investment.rst
     invest_nonconvex.rst
     offset_converter.rst
-    electrical.lopf.rst
     custom_constraints.rst

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/docs/usage.rst

@@ -643,8 +643,8 @@ The following example pictures a Pumped Hydroelectric Energy Storage (PHES). Bot
 
     solph.components.GenericStorage(
         label='PHES',
-        inputs={b_el: solph.flows.Flow(investment= solph.Investment(ep_costs=500))},
-        outputs={b_el: solph.flows.Flow(investment= solph.Investment(ep_costs=500)},
+        inputs={b_el: solph.flows.Flow(nominal_value=solph.Investment(ep_costs=500))},
+        outputs={b_el: solph.flows.Flow(nominal_value=solph.Investment(ep_costs=500)},
         loss_rate=0.001,
         inflow_conversion_factor=0.98, outflow_conversion_factor=0.8),
         investment = solph.Investment(ep_costs=40))
@@ -912,7 +912,7 @@ turbines.
 
     solph.components.Source(label='new_wind_pp', outputs={electricity: solph.flows.Flow(
         fix=wind_power_time_series,
-	investment=solph.Investment(ep_costs=epc, maximum=50000))})
+	nominal_value=solph.Investment(ep_costs=epc, maximum=50000))})
 
 Let's slightly alter the case and consider for already existing wind power
 capacity of 20,000 kW. We're still expecting the total wind power capacity, thus we
@@ -922,7 +922,7 @@ allow for 30,000 kW of new installations and formulate as follows.
 
     solph.components.Source(label='new_wind_pp', outputs={electricity: solph.flows.Flow(
         fix=wind_power_time_series,
-	    investment=solph.Investment(ep_costs=epc,
+	    nominal_value=solph.Investment(ep_costs=epc,
 	                                maximum=30000,
 	                                existing=20000))})
 
@@ -957,7 +957,7 @@ example of a converter:
         label='converter_nonconvex',
         inputs={bus_0: solph.flows.Flow()},
         outputs={bus_1: solph.flows.Flow(
-            investment=solph.Investment(
+            nominal_value=solph.Investment(
                 ep_costs=4,
                 maximum=100,
                 minimum=20,
@@ -1004,8 +1004,8 @@ mathematical background, like variables and constraints, which are used.
 
 .. _multi_period_mode_label:
 
-Using the multi-period (investment) mode (experimental)
--------------------------------------------------------
+Multi-period (investment) mode (experimental)
+---------------------------------------------
 Sometimes you might be interested in how energy systems could evolve in the longer-term, e.g. until 2045 or 2050 to meet some
 carbon neutrality and climate protection or RES and energy efficiency targets.
 
@@ -1020,7 +1020,7 @@ only unfolds if you look at long-term investments. Let's see how.
 First, you start by defining your energy system as you might have done before, but you
 
 * choose a longer-term time horizon (spanning multiple years, i.e. multiple periods) and
-* explicitly define the `periods` attribute of your energy system which maps time steps to the simulated period.
+* explicitly define the `periods` attribute of your energy system which lists the time steps for each period.
 
 .. code-block:: python
 
@@ -1028,10 +1028,10 @@ First, you start by defining your energy system as you might have done before, b
     import oemof.solph as solph
 
     my_index = pd.date_range('1/1/2013', periods=17520, freq='H')
-    periods = {
-        0: pd.date_range('1/1/2013', periods=8760, freq='H'),
-        1: pd.date_range('1/1/2014', periods=8760, freq='H'),
-    }
+    periods = [
+        pd.date_range('1/1/2013', periods=8760, freq='H'),
+        pd.date_range('1/1/2014', periods=8760, freq='H'),
+    ]
     my_energysystem = solph.EnergySystem(timeindex=my_index, periods=periods)
 
 If you want to use a multi-period model you have define periods of your energy system explicitly. This way,
@@ -1057,17 +1057,16 @@ and adjust to your needs:
 
         Returns
         -------
-        periods : dict
-            pd.date_ranges defining the time stamps for the respective period,
-            starting with period 0
+        periods : list
+            periods for the optimization run
         """
         years = sorted(list(set(getattr(datetimeindex, "year"))))
-        periods = {}
+        periods = []
         filter_series = datetimeindex.to_series()
         for number, year in enumerate(years):
             start = filter_series.loc[filter_series.index.year == year].min()
             end = filter_series.loc[filter_series.index.year == year].max()
-            periods[number] = pd.date_range(start, end, freq=datetimeindex.freq)
+            periods.append(pd.date_range(start, end, freq=datetimeindex.freq))
 
         return periods
 
@@ -1144,7 +1143,7 @@ Here is an example
         inputs={hydrogen_bus: solph.flows.Flow()},
         outputs={
             electricity_bus: solph.flows.Flow(
-                investment=solph.Investment(
+                nominal_value=solph.Investment(
                     maximum=1000,
                     ep_costs=1e6,
                     lifetime=30,
@@ -1178,7 +1177,7 @@ This would mean that for investments in the particular period, these values woul
         inputs={hydrogen_bus: solph.flows.Flow()},
         outputs={
             electricity_bus: solph.flows.Flow(
-                investment=solph.Investment(
+                nominal_value=solph.Investment(
                     maximum=1000,
                     ep_costs=[1e6, 1.1e6],
                     lifetime=30,
@@ -1269,6 +1268,13 @@ Besides the `invest` variable, new variables are introduced as well. These are:
     * You can specify periods of different lengths, but the frequency of your timeindex needs to be consistent. Also,
       you could use the `timeincrement` attribute of the energy system to model different weightings. Be aware that this
       has not yet been tested.
+    * For now, both, the `timeindex` as well as the `timeincrement` of an energy system have to be defined since they
+      have to be of the same length for a multi-period model.
+    * You can choose whether to re-evaluate assets at the end of the optimization horizon. If you set attribute
+      `use_remaining_value` of the energy system to True (defaults to False), this leads to the model evaluating the
+      difference in the asset value at the end of the optimization horizon vs. at the time the investment was made.
+      The difference in value is added to or subtracted from the respective investment costs increment,
+      assuming assets are to be liquidated / re-evaluated at the end of the optimization horizon.
     * Also please be aware, that periods correspond to years by default. You could also choose
       monthly periods, but you would need to be very careful in parameterizing your energy system and your model and also,
       this would mean monthly discounting (if applicable) as well as specifying your plants lifetimes in months.
@@ -1466,12 +1472,11 @@ This nonlinearity is linearised in the
         inputs={b_diesel: solph.flows.Flow()},
         outputs={
             b_el: solph.flows.Flow(
-                nominal_value=None,
                 variable_costs=0.04,
                 min=0.2,
                 max=1,
                 nonconvex=solph.NonConvex(),
-                investment=solph.Investment(
+                nominal_value=solph.Investment(
                     ep_costs=90,
                     maximum=150, # required for the linearization
                 ),

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/docs/whatsnew/v0-5-1.rst

@@ -23,8 +23,8 @@ New features
 
 * Add option to run multi-period (dynamic) investment models with oemof.solph as an experimental feature:
     * You can change from standard model to multi-period model by defining the newly introduced `periods`
-      attribute of your energy system. Be aware that it is experimental as of now. `periods` is a dictionary
-      mapping the periods you want to model (usually years) to pandas.date_range objects.
+      attribute of your energy system. Be aware that it is experimental as of now. `periods` is a list
+      of the periods you want to model (usually years) given as pandas.date_range objects.
     * Add attributes `periods` to :class:`oemof.solph._energy_system.EnergySystem`.
     * Introduce new Pyomo Sets `PERIODS` and `TIMEINDEX` in :class:`oemof.solph.models.Model`.
     * Index all investment-related variables with `PERIODS` and flow variable with `TIMEINDEX`, which

oemof.solph-0.5.2.dev2/docs/whatsnew/v0-5-2.rst

@@ -0,0 +1,52 @@
+v0.5.2 (????)
+-------------
+
+API changes
+###########
+
+* New bool attribute `use_remaining_value` of `oemof.solph.EnergySystem`
+* Use list for period definition in multi-period investment optimization.
+
+New features
+############
+
+* Allow for evaluating differences in the remaining vs. the original value
+  for multi-period investments.
+* Allow to define minimum up- and down-time per time step
+
+Documentation
+#############
+
+Bug fixes
+#########
+
+* Fix handling of investment annuities and fixed costs for multi-period models:
+  Limit to costs that occur within the optimization horizon to prevent a
+  bias towards investments happening earlier in the optimization horizon.
+* Fix bugs in multi-period documentation.
+* Fix y intersect of OffsetConverter
+* Fix minimum uptime being relevant for initial downtime (and vice versa).
+* Fix duplicated discounting of fixed costs for multi-period investment
+
+Testing
+#######
+
+Other changes
+#############
+
+* Improved compatibility with upcoming oemof.network
+
+Known issues
+############
+
+* Documentation, example and API of OffsetConverter are inconsistent
+  and might not work as expected.
+
+Contributors
+############
+
+* Patrik Schönfeldt
+* Johannes Kochems
+* Julian Endres
+* Hendrik Huyskens
+* Raul Ciria Aylagas

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/pyproject.toml

@@ -19,5 +19,5 @@ exclude = '''
 '''
 
 [build-system]
-requires = ["setuptools<60.0", "wheel"]
+requires = ["setuptools", "wheel"]
 build-backend = "setuptools.build_meta"

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/setup.cfg

@@ -28,7 +28,7 @@ python_files =
 	tests.py
 addopts = 
 	-ra
-	--strict
+	--strict-markers
 	--ignore=docs/conf.py
 	--ignore=setup.py
 	--ignore=ci

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/setup.py

@@ -29,7 +29,7 @@ long_description = "%s" % (
 
 setup(
     name="oemof.solph",
-    version="0.5.2dev0",
+    version="0.5.2.dev2",
     license="MIT",
     description=(
         "A model generator for energy system modelling and optimisation."

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/src/oemof/solph/__init__.py

@@ -1,4 +1,4 @@
-__version__ = "0.5.2dev0"
+__version__ = "0.5.2.dev2"
 
 from . import buses
 from . import components

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/src/oemof/solph/_energy_system.py

@@ -62,6 +62,10 @@ class EnergySystem(es.EnergySystem):
         For a standard model, periods are not (to be) declared, i.e. None.
         A list with one entry is derived, i.e. [0].
 
+    use_remaining_value : bool
+        If True, compare the remaining value of an investment to the
+        original value (only applicable for multi-period models)
+
     kwargs
     """
 
@@ -71,14 +75,17 @@ class EnergySystem(es.EnergySystem):
         timeincrement=None,
         infer_last_interval=None,
         periods=None,
-        **kwargs,
+        use_remaining_value=False,
+        groupings=None,
     ):
         # Doing imports at runtime is generally frowned upon, but should work
         # for now. See the TODO in :func:`constraint_grouping
         # <oemof.solph.groupings.constraint_grouping>` for more information.
         from oemof.solph import GROUPINGS
 
-        kwargs["groupings"] = GROUPINGS + kwargs.get("groupings", [])
+        if groupings is None:
+            groupings = []
+        groupings = GROUPINGS + groupings
 
         if not (
             isinstance(timeindex, pd.DatetimeIndex)
@@ -133,10 +140,9 @@ class EnergySystem(es.EnergySystem):
             else:
                 msg = (
                     "Ensure that your timeindex and timeincrement are "
-                    "consistent.\nIf you are not considering non-equidistant "
-                    "timeindices, consider only specifying a timeindex."
+                    "consistent."
                 )
-                warnings.warn(msg, debugging.SuspiciousUsageWarning)
+                warnings.warn(msg, debugging.ExperimentalFeatureWarning)
 
         elif timeindex is not None and timeincrement is None:
             df = pd.DataFrame(timeindex)
@@ -151,16 +157,19 @@ class EnergySystem(es.EnergySystem):
         if timeincrement is not None and (pd.Series(timeincrement) <= 0).any():
             msg = (
                 "The time increment is inconsistent. Negative values and zero "
-                "is not allowed.\nThis is caused by a inconsistent "
+                "are not allowed.\nThis is caused by a inconsistent "
                 "timeincrement parameter or an incorrect timeindex."
             )
             raise TypeError(msg)
 
         super().__init__(
-            timeindex=timeindex, timeincrement=timeincrement, **kwargs
+            groupings=groupings,
+            timeindex=timeindex,
+            timeincrement=timeincrement,
         )
 
-        if periods is not None:
+        self.periods = periods
+        if self.periods is not None:
             msg = (
                 "CAUTION! You specified the 'periods' attribute for your "
                 "energy system.\n This will lead to creating "
@@ -170,48 +179,71 @@ class EnergySystem(es.EnergySystem):
                 "now. If you find anything suspicious or any bugs, "
                 "please report them."
             )
-            warnings.warn(msg, debugging.SuspiciousUsageWarning)
-        self.periods = periods
-        self._extract_periods_years()
-        self._extract_periods_matrix()
+            warnings.warn(msg, debugging.ExperimentalFeatureWarning)
+            self._extract_periods_years()
+            self._extract_periods_matrix()
+            self._extract_end_year_of_optimization()
+            self.use_remaining_value = use_remaining_value
 
     def _extract_periods_years(self):
-        """Map simulation years to the respective period based on time indices
+        """Map years in optimization to respective period based on time indices
 
-        Returns
-        -------
-        periods_years: dict
-            the simulation year of the start of each a period,
-            relative to the start of the optimization run and starting with 0
+        Attribute `periods_years` of type list is set. It contains
+        the year of the start of each period, relative to the
+        start of the optimization run and starting with 0.
         """
         periods_years = [0]
-        if self.periods is not None:
-            start_year = self.periods[0].min().year
-            for k, v in enumerate(self.periods):
-                if k >= 1:
-                    periods_years.append(v.min().year - start_year)
+        start_year = self.periods[0].min().year
+        for k, v in enumerate(self.periods):
+            if k >= 1:
+                periods_years.append(v.min().year - start_year)
 
-            self.periods_years = periods_years
+        self.periods_years = periods_years
 
     def _extract_periods_matrix(self):
         """Determines a matrix describing the temporal distance to each period.
+
+        Attribute `periods_matrix` of type list np.array is set.
         Rows represent investment/commissioning periods, columns represent
         decommissioning periods. The values describe the temporal distance
         between each investment period to each decommissioning period.
+        """
+        periods_matrix = []
+        period_years = np.array(self.periods_years)
+        for v in period_years:
+            row = period_years - v
+            row = np.where(row < 0, 0, row)
+            periods_matrix.append(row)
+        self.periods_matrix = np.array(periods_matrix)
+
+    def _extract_end_year_of_optimization(self):
+        """Extract the end of the optimization in years
+
+        Attribute `end_year_of_optimization` of int is set.
+        """
+        duration_last_period = self.get_period_duration(-1)
+        self.end_year_of_optimization = (
+            self.periods_years[-1] + duration_last_period
+        )
+
+    def get_period_duration(self, period):
+        """Get duration of a period in full years
+
+        Parameters
+        ----------
+        period : int
+            Period for which the duration in years shall be obtained
 
         Returns
         -------
-        period_distance_matrix: np.array
-
+        int
+            Duration of the period
         """
-        periods_matrix = []
-        if self.periods is not None:
-            period_years = np.array(self.periods_years)
-            for v in period_years:
-                row = period_years - v
-                row = np.where(row < 0, 0, row)
-                periods_matrix.append(row)
-            self.periods_matrix = np.array(periods_matrix)
+        return (
+            self.periods[period].max().year
+            - self.periods[period].min().year
+            + 1
+        )
 
 
 def create_time_index(

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/src/oemof/solph/_models.py

@@ -132,12 +132,21 @@ class BaseModel(po.ConcreteModel):
             "constraint_groups", []
         )
 
-        self._constraint_groups += [
-            i
-            for i in self.es.groups
-            if hasattr(i, "CONSTRAINT_GROUP")
-            and i not in self._constraint_groups
-        ]
+        if self.is_cellular:
+            for es in energysystem:
+                self._constraint_groups += [
+                    i
+                    for i in es.groups
+                    if hasattr(i, "CONSTRAINT_GROUP")
+                    and i not in self._constraint_groups
+                ]
+        else:
+            self._constraint_groups += [
+                i
+                for i in self.es.groups
+                if hasattr(i, "CONSTRAINT_GROUP")
+                and i not in self._constraint_groups
+            ]
 
         self.flows = self.es.flows()
         if self.is_cellular:
@@ -393,12 +402,11 @@ class Model(BaseModel):
         """Add all basic sets to the model, i.e. NODES, TIMESTEPS and FLOWS.
         Also create sets PERIODS and TIMEINDEX used for multi-period models.
         """
-        self.nodes = self.es.nodes
+        self.nodes = list(self.es.nodes)
         if self.is_cellular:
             # collect all nodes from the child cells
             for cell in self.ec:
-                for node in cell.nodes:
-                    self.nodes.append(node)
+                self.nodes.extend(cell.nodes)
         # create set with all nodes
         self.NODES = po.Set(initialize=[n for n in self.nodes])
 

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/src/oemof/solph/_options.py

@@ -13,7 +13,6 @@ SPDX-FileCopyrightText: Johannes Kochems
 SPDX-License-Identifier: MIT
 
 """
-
 from oemof.solph._plumbing import sequence
 
 
@@ -62,7 +61,7 @@ class Investment:
         Units lifetime, given in years; only applicable for multi-period
         models
     age : int, :math:`a`
-        Units start age, given in years at the beginning of the simulation;
+        Units start age, given in years at the beginning of the optimization;
         only applicable for multi-period models
     interest_rate : float, :math:`ir`
         Interest rate for calculating annuities when investing in a particular
@@ -184,14 +183,14 @@ class NonConvex:
         from the actual output.
     inactivity_costs : numeric (iterable or scalar)
         Costs associated with not operating the flow.
-    minimum_uptime : numeric (1 or positive integer)
-        Minimum time that a flow must be greater then its minimum flow after
-        startup. Be aware that minimum up and downtimes can contradict each
-        other and may lead to infeasible problems.
-    minimum_downtime : numeric (1 or positive integer)
-        Minimum time a flow is forced to zero after shutting down.
-        Be aware that minimum up and downtimes can contradict each
-        other and may to infeasible problems.
+    minimum_uptime : numeric or list of numeric (1 or positive integer)
+        Minimum number of time steps that a flow must be greater then its
+        minimum flow after startup. Be aware that minimum up and downtimes
+        can contradict each other and may lead to infeasible problems.
+    minimum_downtime : numeric or list of numeric (1 or positive integer)
+        Minimum number of time steps a flow is forced to zero after
+        shutting down. Be aware that minimum up and downtimes can
+        contradict each other and may to infeasible problems.
     maximum_startups : numeric (0 or positive integer)
         Maximum number of start-ups in the optimization timeframe.
     maximum_shutdowns : numeric (0 or positive integer)
@@ -199,12 +198,12 @@ class NonConvex:
     initial_status : numeric (0 or 1)
         Integer value indicating the status of the flow in the first time step
         (0 = off, 1 = on). For minimum up and downtimes, the initial status
-        is set for the respective values in the edge regions e.g. if a
-        minimum uptime of four timesteps is defined, the initial status is
-        fixed for the four first and last timesteps of the optimization period.
-        If both, up and downtimes are defined, the initial status is set for
-        the maximum of both e.g. for six timesteps if a minimum downtime of
-        six timesteps is defined besides a four timestep minimum uptime.
+        is set for the respective values in the beginning e.g. if a
+        minimum uptime of four timesteps is defined and the initial status is
+        set to one, the initial status is fixed for the four first timesteps
+        of the optimization period. Otherwise if the initial status is set to
+        zero and the first timesteps are fixed for the number of minimum
+        downtime steps.
     negative_gradient_limit : numeric (iterable, scalar or None)
         the normed *upper bound* on the positive difference
         (`flow[t-1] < flow[t]`) of two consecutive flow values.
@@ -232,8 +231,8 @@ class NonConvex:
             custom_attributes = {}
 
         self.initial_status = initial_status
-        self.minimum_uptime = minimum_uptime
-        self.minimum_downtime = minimum_downtime
+        self.minimum_uptime = sequence(minimum_uptime)
+        self.minimum_downtime = sequence(minimum_downtime)
         self.maximum_startups = maximum_startups
         self.maximum_shutdowns = maximum_shutdowns
 
@@ -247,11 +246,7 @@ class NonConvex:
         for attribute, value in custom_attributes.items():
             setattr(self, attribute, value)
 
-    @property
-    def max_up_down(self):
-        """Return maximum of minimum_uptime and minimum_downtime.
-
-        The maximum of both is used to set the initial status for this
-        number of time steps within the edge regions.
-        """
-        return max(self.minimum_uptime, self.minimum_downtime)
+        if initial_status == 0:
+            self.first_flexible_timestep = self.minimum_downtime[0]
+        else:
+            self.first_flexible_timestep = self.minimum_uptime[0]

{oemof.solph-0.5.2.dev0 → oemof.solph-0.5.2.dev2}/src/oemof/solph/_plumbing.py

@@ -61,7 +61,7 @@ class _Sequence(UserList):
     --------
     >>> s = _Sequence(default=42)
     >>> len(s)
-    0
+    1
     >>> s[1]
     42
     >>> s[2]
@@ -79,7 +79,7 @@ class _Sequence(UserList):
     def __init__(self, *args, **kwargs):
         self.default = kwargs["default"]
         self.default_changed = False
-        self.highest_index = -1
+        self.highest_index = 0
         super().__init__(*args)
 
     def __getitem__(self, key):