pypsa 1.2.1__tar.gz → 1.2.2__tar.gz

{pypsa-1.2.1 → pypsa-1.2.2}/.github/workflows/test-models.yml

@@ -120,7 +120,10 @@ jobs:
     - name: Install package from ref
       if: (steps.filter.outputs.src == 'true' || github.event_name == 'schedule') && env.pinned == 'false'
       run: |
+        # Allow upgrading internal packages without exclude-newer cap
+        sed -i '/^exclude-newer\s*=/d' pixi.toml
         pixi remove ${{ env.PACKAGE_NAME }}
+        pixi upgrade linopy
         pixi add --pypi --git https://github.com/${{ github.repository }}.git ${{ github.event.repository.name }} --rev ${{ github.event_name == 'pull_request' && github.event.pull_request.head.sha || github.sha }}
 
     - name: Run snakemake test workflows

{pypsa-1.2.1/pypsa.egg-info → pypsa-1.2.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pypsa
-Version: 1.2.1
+Version: 1.2.2
 Summary: Python for Power Systems Analysis
 Author-email: "PyPSA Developers, see https://docs.pypsa.org/latest/contributing/contributors.html" <t.brown@tu-berlin.de>
 License: MIT License
@@ -42,9 +42,9 @@ License-File: LICENSE
 Requires-Dist: numpy
 Requires-Dist: scipy
 Requires-Dist: pandas>=2.0
-Requires-Dist: xarray<2026.4
+Requires-Dist: xarray
 Requires-Dist: netcdf4
-Requires-Dist: linopy>=0.6.1
+Requires-Dist: linopy>=0.7.0
 Requires-Dist: matplotlib
 Requires-Dist: plotly
 Requires-Dist: pydeck>=0.6

{pypsa-1.2.1 → pypsa-1.2.2}/docs/examples/tracing-infeasibilities.ipynb

@@ -167,11 +167,7 @@
    "cell_type": "markdown",
    "id": "d3a8d32d",
    "metadata": {},
-   "source": [
-    "## Tracing the Infeasibility\n",
-    "\n",
-    "When the optimization fails due to infeasibility, PyPSA provides a convenient method to diagnose the problem. The `print_infeasibilities()` method on the optimization model will show us exactly which constraints are causing the infeasibility."
-   ]
+   "source": "## Tracing the Infeasibility\n\nWhen the optimization fails due to infeasibility, PyPSA provides a convenient method to diagnose the problem. The `format_infeasibilities()` method on the optimization model will show us exactly which constraints are causing the infeasibility."
   },
   {
    "cell_type": "code",
@@ -179,34 +175,13 @@
    "id": "4989a457",
    "metadata": {},
    "outputs": [],
-   "source": [
-    "n.model.print_infeasibilities()"
-   ]
+   "source": "print(n.model.format_infeasibilities())"
   },
   {
    "cell_type": "markdown",
    "id": "ee963ad8",
    "metadata": {},
-   "source": [
-    "### Interpreting the Infeasibility Output\n",
-    "\n",
-    "The infeasibility trace above shows us exactly what's wrong:\n",
-    "\n",
-    "1. **`Link-fix-p-lower[2015-01-01 00:00:00, DC link]: +1 Link-p[2015-01-01 00:00:00, DC link] ≥ -0`**\n",
-    "   - This constraint says the DC link power flow must be ≥ 0 (can't flow negative)\n",
-    "   - But we set the DC link capacity to 0, so Link-p value is bounded to 0\n",
-    "\n",
-    "2. **`Bus-nodal_balance[London, 2015-01-01 00:00:00]: -1 Link-p[2015-01-01 00:00:00, DC link] = 35.7962441027`**\n",
-    "   - This is London's nodal balance equation: generation ± trade ± storage = load\n",
-    "   - London has no generation, 35.8 MW load, and no AC interconnectors\n",
-    "   - The DC link (the only remaining connection) can't provide any power because its capacity is 0\n",
-    "\n",
-    "Under the hood, `print_infeasibilities()` uses Gurobi's infeasibility analysis to identify conflicting constraints.\n",
-    " **Irreducible Inconsistent Set (IIS)** identified by Gurobi shows these two constraints are fundamentally incompatible:\n",
-    "- London needs power inflow via the DC link to meet its load\n",
-    "- But the DC link is constrained to zero power\n",
-    "- This creates an impossible situation that no solution can satisfy"
-   ]
+   "source": "### Interpreting the Infeasibility Output\n\nThe infeasibility trace above shows us exactly what's wrong:\n\n1. **`Link-fix-p-lower[2015-01-01 00:00:00, DC link]: +1 Link-p[2015-01-01 00:00:00, DC link] ≥ -0`**\n   - This constraint says the DC link power flow must be ≥ 0 (can't flow negative)\n   - But we set the DC link capacity to 0, so Link-p value is bounded to 0\n\n2. **`Bus-nodal_balance[London, 2015-01-01 00:00:00]: -1 Link-p[2015-01-01 00:00:00, DC link] = 35.7962441027`**\n   - This is London's nodal balance equation: generation ± trade ± storage = load\n   - London has no generation, 35.8 MW load, and no AC interconnectors\n   - The DC link (the only remaining connection) can't provide any power because its capacity is 0\n\nUnder the hood, `format_infeasibilities()` uses Gurobi's infeasibility analysis to identify conflicting constraints.\n **Irreducible Inconsistent Set (IIS)** identified by Gurobi shows these two constraints are fundamentally incompatible:\n- London needs power inflow via the DC link to meet its load\n- But the DC link is constrained to zero power\n- This creates an impossible situation that no solution can satisfy"
   }
  ],
  "metadata": {

{pypsa-1.2.1 → pypsa-1.2.2}/docs/home/features.md

@@ -146,7 +146,7 @@ Interactive **network visualization** of SciGRID example network clustered by fe
 Interactive area plot of **electricity balance time series** in a highly-renewable sector-coupled example network, showing temporal generation (positive values) and consumption (negative values) of different technologies:
 
 <div style="width: 100%; height: 550px;">
-    <iframe src="../../../assets/interactive/carbon_management-energy_balance-area_iplot-AC-bus_carrier.html"
+    <iframe src="../../assets/interactive/carbon_management-energy_balance-area_iplot-AC-bus_carrier.html"
             width="100%" height="100%" frameborder="0" style="border: 0px solid #ccc;">
     </iframe>
 </div>

{pypsa-1.2.1 → pypsa-1.2.2}/docs/release-notes.md

@@ -6,13 +6,24 @@ SPDX-License-Identifier: CC-BY-4.0
 
 # Release Notes
 
-<!--## Upcoming Release
+## [**v1.2.2**](https://github.com/PyPSA/PyPSA/releases/tag/v1.2.2) <small>25th May 2026</small> { id="v1.2.2" }
 
+<!--
 !!! info "Upcoming Release"
 
     The features listed below have not yet been released, but will be included in the
     next update! If you would like to use these features in the meantime, you will need
-    to install the `master` branch, e.g. `pip install git+https://github.com/pypsa/pypsa`.-->
+    to install the `master` branch, e.g. `pip install git+https://github.com/pypsa/pypsa`.
+-->
+
+### Bug Fixes
+
+- Fix the sign of [Loads](./user-guide/components/loads.md) not being taken into account in the nodal balance constraint when calling [`n.optimize()`][pypsa.optimization.OptimizationAccessor.__call__]. (<!-- md:pr 1685 -->)
+
+- Fix operational constraints for non-extendable components producing `NaN` bounds when `p_nom` is infinite and `p_min_pu`/`p_max_pu` is zero. The bound now falls back to zero in this case. Relevant for linopy versions `>=0.7` where `NaN` bounds are not dropped explicitly. (<!-- md:pr 1683 -->)
+
+- Lift `xarray<2026.4` upper bound and bump `linopy>=0.7.0` floor. (<!-- md:pr 1686 -->)
+
 
 ## [**v1.2.1**](https://github.com/PyPSA/PyPSA/releases/tag/v1.2.1) <small>19th May 2026</small> { id="v1.2.1" }
 

{pypsa-1.2.1 → pypsa-1.2.2}/docs/requirements.txt

@@ -13,7 +13,7 @@ babel==2.18.0
     # via
     #   mkdocs-git-revision-date-localized-plugin
     #   mkdocs-material
-backrefs==6.2
+backrefs==7.0
     # via mkdocs-material
 beautifulsoup4==4.14.3
     # via nbconvert
@@ -27,7 +27,7 @@ cairosvg==2.9.0
     # via pypsa
 cartopy==0.25.0
     # via pypsa
-certifi==2026.2.25
+certifi==2026.4.22
     # via
     #   netcdf4
     #   pyogrio
@@ -36,13 +36,12 @@ certifi==2026.2.25
 cffi==2.0.0
     # via
     #   cairocffi
-    #   cryptography
     #   pyzmq
 cftime==1.6.5
     # via netcdf4
 charset-normalizer==3.4.7
     # via requests
-click==8.3.2
+click==8.4.0
     # via
     #   dask
     #   mkdocs
@@ -62,8 +61,6 @@ comm==0.2.3
     # via ipykernel
 contourpy==1.3.3
     # via matplotlib
-cryptography==46.0.7
-    # via google-auth
 cssselect2==0.9.0
     # via cairosvg
 cycler==0.12.1
@@ -72,7 +69,7 @@ dask==2026.3.0
     # via linopy
 debugpy==1.8.20
     # via ipykernel
-decorator==5.2.1
+decorator==5.3.1
     # via ipython
 defusedxml==0.7.1
     # via
@@ -90,9 +87,9 @@ executing==2.2.1
     # via stack-data
 fastjsonschema==2.21.2
     # via nbformat
-fonttools==4.62.1
+fonttools==4.63.0
     # via matplotlib
-fsspec==2026.3.0
+fsspec==2026.4.0
     # via dask
 geopandas==1.1.3
     # via pypsa
@@ -102,38 +99,17 @@ ghp-import==2.1.0
     #   properdocs
 gitdb==4.0.12
     # via gitpython
-gitpython==3.1.46
+gitpython==3.1.50
     # via mkdocs-git-revision-date-localized-plugin
-google-api-core==2.30.3
-    # via
-    #   google-cloud-core
-    #   google-cloud-storage
-google-auth==2.49.2
-    # via
-    #   google-api-core
-    #   google-cloud-core
-    #   google-cloud-storage
-google-cloud-core==2.5.1
-    # via google-cloud-storage
-google-cloud-storage==3.10.1
-    # via linopy
-google-crc32c==1.8.0
-    # via
-    #   google-cloud-storage
-    #   google-resumable-media
-google-resumable-media==2.8.2
-    # via google-cloud-storage
-googleapis-common-protos==1.74.0
-    # via google-api-core
 griffelib==2.0.2
     # via mkdocstrings-python
-gurobipy==13.0.1
+gurobipy==13.0.2
     # via pypsa
-highspy==1.13.1
+highspy==1.14.0
     # via
     #   pypsa
     #   tsam
-idna==3.12
+idna==3.15
     # via requests
 importlib-metadata==9.0.0 ; python_full_version < '3.12'
     # via dask
@@ -141,13 +117,11 @@ iniconfig==2.3.0
     # via pytest
 ipykernel==7.2.0
     # via mkdocs-jupyter
-ipython==9.10.1 ; python_full_version < '3.12'
-    # via ipykernel
-ipython==9.12.0 ; python_full_version >= '3.12'
+ipython==9.13.0
     # via ipykernel
 ipython-pygments-lexers==1.1.1
     # via ipython
-jedi==0.19.2
+jedi==0.20.0
     # via ipython
 jinja2==3.1.6
     # via
@@ -176,17 +150,17 @@ jupyter-core==5.9.1
     #   nbformat
 jupyterlab-pygments==0.3.0
     # via nbconvert
-jupytext==1.19.1
+jupytext==1.19.3
     # via mkdocs-jupyter
 kiwisolver==1.5.0
     # via matplotlib
 levenshtein==0.27.3
     # via pypsa
-linopy==0.6.7
+linopy==0.7.0
     # via pypsa
 locket==1.0.0
     # via partd
-lxml==6.1.0
+lxml==6.1.1
     # via mkdocs-video
 markdown==3.10.2
     # via
@@ -196,7 +170,7 @@ markdown==3.10.2
     #   mkdocstrings
     #   properdocs
     #   pymdown-extensions
-markdown-it-py==4.0.0
+markdown-it-py==4.2.0
     # via
     #   jupytext
     #   mdit-py-plugins
@@ -208,17 +182,17 @@ markupsafe==3.0.3
     #   mkdocstrings
     #   nbconvert
     #   properdocs
-matplotlib==3.10.8
+matplotlib==3.10.9
     # via
     #   cartopy
     #   pypsa
     #   salib
     #   seaborn
-matplotlib-inline==0.2.1
+matplotlib-inline==0.2.2
     # via
     #   ipykernel
     #   ipython
-mdit-py-plugins==0.5.0
+mdit-py-plugins==0.6.1
     # via jupytext
 mdurl==0.1.2
     # via markdown-it-py
@@ -226,7 +200,7 @@ mergedeep==1.3.4
     # via
     #   mkdocs
     #   mkdocs-get-deps
-mistune==3.2.0
+mistune==3.2.1
     # via nbconvert
 mkdocs==1.6.1
     # via
@@ -250,7 +224,7 @@ mkdocs-ezlinks-plugin==0.1.14
     # via pypsa
 mkdocs-get-deps==0.2.2
     # via mkdocs
-mkdocs-git-revision-date-localized-plugin==1.5.1
+mkdocs-git-revision-date-localized-plugin==1.5.2
     # via pypsa
 mkdocs-jupyter==0.26.3
     # via pypsa
@@ -272,7 +246,7 @@ mkdocstrings-python==2.0.3
     # via pypsa
 multiprocess==0.70.19
     # via salib
-narwhals==2.20.0
+narwhals==2.21.2
     # via plotly
 nbclient==0.10.4
     # via nbconvert
@@ -293,7 +267,7 @@ networkx==3.6.1
     #   tsam
 numexpr==2.14.1
     # via linopy
-numpy==2.4.4
+numpy==2.4.6
     # via
     #   bottleneck
     #   cartopy
@@ -348,11 +322,11 @@ pandas==2.3.3
     #   xarray
 pandocfilters==1.5.1
     # via nbconvert
-parso==0.8.6
+parso==0.8.7
     # via jedi
 partd==1.4.2
     # via dask
-pathspec==1.0.4
+pathspec==1.1.1
     # via
     #   mkdocs
     #   properdocs
@@ -373,31 +347,22 @@ plotly==6.7.0
     # via pypsa
 pluggy==1.6.0
     # via pytest
-polars==1.40.0
+polars==1.40.1
     # via linopy
-polars-runtime-32==1.40.0
+polars-runtime-32==1.40.1
     # via polars
 prompt-toolkit==3.0.52
     # via ipython
 properdocs==1.6.7
     # via mkdocs-redirects
-proto-plus==1.27.2
-    # via google-api-core
-protobuf==7.34.1
-    # via
-    #   google-api-core
-    #   googleapis-common-protos
-    #   proto-plus
 psutil==7.2.2
-    # via ipykernel
+    # via
+    #   ipykernel
+    #   ipython
 ptyprocess==0.7.0 ; sys_platform != 'emscripten' and sys_platform != 'win32'
     # via pexpect
 pure-eval==0.2.3
     # via stack-data
-pyasn1==0.6.3
-    # via pyasn1-modules
-pyasn1-modules==0.4.2
-    # via google-auth
 pycparser==3.0 ; implementation_name != 'PyPy'
     # via cffi
 pydeck==0.9.2
@@ -412,7 +377,7 @@ pygments==2.20.0
     #   pytest
 pygtrie==2.5.0
     # via mkdocs-ezlinks-plugin
-pymdown-extensions==10.21.2
+pymdown-extensions==10.21.3
     # via
     #   mkdocs-material
     #   mkdocstrings
@@ -426,7 +391,7 @@ pyproj==3.7.2
     # via
     #   cartopy
     #   geopandas
-pyshp==3.0.3
+pyshp==3.0.5
     # via cartopy
 pytest==9.0.3
     # via pypsa
@@ -438,7 +403,7 @@ python-dateutil==2.9.0.post0
     #   jupyter-client
     #   matplotlib
     #   pandas
-pytz==2026.1.post1
+pytz==2026.2
     # via pandas
 pyyaml==6.0.3
     # via
@@ -464,12 +429,8 @@ referencing==0.37.0
     # via
     #   jsonschema
     #   jsonschema-specifications
-requests==2.33.1
-    # via
-    #   google-api-core
-    #   google-cloud-storage
-    #   linopy
-    #   mkdocs-material
+requests==2.34.2
+    # via mkdocs-material
 rpds-py==0.30.0
     # via
     #   jsonschema
@@ -521,7 +482,7 @@ tqdm==4.67.3
     # via
     #   linopy
     #   tsam
-traitlets==5.14.3
+traitlets==5.15.0
     # via
     #   ipykernel
     #   ipython
@@ -531,18 +492,18 @@ traitlets==5.14.3
     #   nbclient
     #   nbconvert
     #   nbformat
-tsam==3.3.0
+tsam==2.3.9
     # via pypsa
 typing-extensions==4.15.0
     # via
     #   beautifulsoup4
     #   ipython
     #   referencing
-tzdata==2026.1
+tzdata==2026.2
     # via
     #   mkdocs-git-revision-date-localized-plugin
     #   pandas
-urllib3==2.6.3
+urllib3==2.7.0
     # via requests
 validators==0.35.0
     # via pypsa
@@ -550,7 +511,7 @@ watchdog==6.0.0
     # via
     #   mkdocs
     #   properdocs
-wcwidth==0.6.0
+wcwidth==0.7.0
     # via prompt-toolkit
 webencodings==0.5.1
     # via
@@ -561,5 +522,5 @@ xarray==2026.2.0
     # via
     #   linopy
     #   pypsa
-zipp==3.23.1 ; python_full_version < '3.12'
+zipp==4.1.0 ; python_full_version < '3.12'
     # via importlib-metadata

{pypsa-1.2.1 → pypsa-1.2.2}/pyproject.toml

@@ -33,9 +33,9 @@ dependencies = [
     "numpy",
     "scipy",
     "pandas>=2.0",
-    "xarray<2026.4",  # See https://github.com/PyPSA/linopy/pull/647
+    "xarray",
     "netcdf4",
-    "linopy>=0.6.1",
+    "linopy>=0.7.0",
     "matplotlib",
     "plotly",
     "pydeck>=0.6",

{pypsa-1.2.1 → pypsa-1.2.2}/pypsa/optimization/constraints.py

@@ -14,7 +14,7 @@ import linopy
 import pandas as pd
 import xarray as xr
 from linopy import merge
-from numpy import inf, isfinite, maximum, sqrt, tile
+from numpy import inf, isfinite, isinf, maximum, sqrt, tile
 from xarray import DataArray, concat, where
 
 from pypsa.common import as_index, expand_series
@@ -102,8 +102,12 @@ def define_operational_constraints_for_non_extendables(
         min_pu = min_pu.sel(snapshot=sns)
         max_pu = max_pu.sel(snapshot=sns)
 
-    lower = min_pu * nominal_fix
-    upper = max_pu * nominal_fix
+    is_inf = isinf(nominal_fix)
+    is_zero_min = min_pu == 0
+    is_zero_max = max_pu == 0
+
+    lower = (min_pu * nominal_fix).where(~(is_inf & is_zero_min), 0)
+    upper = (max_pu * nominal_fix).where(~(is_inf & is_zero_max), 0)
 
     active = c.da.active.sel(name=fix_i, snapshot=sns)
 
@@ -1144,9 +1148,8 @@ def define_nodal_balance_constraints(
             dims=["snapshot", "name"],
         )
     else:
-        loads_values = loads.da.p_set.where(
-            loads.da.active.sel(name=loads.active_assets, snapshot=sns)
-        )
+        active_mask = loads.da.active.sel(name=loads.active_assets, snapshot=sns)
+        loads_values = (-loads.da.p_set * loads.da.sign).where(active_mask)
         loads_values = loads_values.reindex(name=loads.static.index.unique("name"))
         load_buses = loads._as_xarray("bus").rename("Bus")
         if n.has_scenarios:

{pypsa-1.2.1 → pypsa-1.2.2}/pypsa/optimization/optimize.py

@@ -594,7 +594,7 @@ class OptimizationAccessor(OptimizationAbstractMixin):
             and compute_infeasibilities
             and "gurobi" in available_solvers
         ):
-            n.model.print_infeasibilities()
+            logger.info(n.model.format_infeasibilities())
 
         return status, condition
 

{pypsa-1.2.1 → pypsa-1.2.2/pypsa.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pypsa
-Version: 1.2.1
+Version: 1.2.2
 Summary: Python for Power Systems Analysis
 Author-email: "PyPSA Developers, see https://docs.pypsa.org/latest/contributing/contributors.html" <t.brown@tu-berlin.de>
 License: MIT License
@@ -42,9 +42,9 @@ License-File: LICENSE
 Requires-Dist: numpy
 Requires-Dist: scipy
 Requires-Dist: pandas>=2.0
-Requires-Dist: xarray<2026.4
+Requires-Dist: xarray
 Requires-Dist: netcdf4
-Requires-Dist: linopy>=0.6.1
+Requires-Dist: linopy>=0.7.0
 Requires-Dist: matplotlib
 Requires-Dist: plotly
 Requires-Dist: pydeck>=0.6

{pypsa-1.2.1 → pypsa-1.2.2}/pypsa.egg-info/requires.txt

@@ -1,9 +1,9 @@
 numpy
 scipy
 pandas>=2.0
-xarray<2026.4
+xarray
 netcdf4
-linopy>=0.6.1
+linopy>=0.7.0
 matplotlib
 plotly
 pydeck>=0.6

{pypsa-1.2.1 → pypsa-1.2.2}/test/test_lopf_basic_constraints.py

@@ -404,3 +404,79 @@ def test_define_fixed_operational_constraints_extendable():
     assert n.c.generators.dynamic.p["gen0"].eq(5).all()
     assert n.c.generators.dynamic.p["gen1"].eq(5).all()
     assert n.c.generators.dynamic.p["gen2"].eq(0).all()
+
+
+def test_nodal_balance_respects_load_sign():
+    """
+    The sign attribute of Load components must be considered in the
+    nodal balance constraint. With sign=-1 (default), a load consumes
+    power; with sign=+1, it represents a fixed supply injection.
+
+    Verifies both:
+      1. The nodal balance constraint RHS equals -p_set * sign.
+      2. The solved generator dispatch matches the resulting net demand.
+    """
+    n = pypsa.Network()
+    n.add("Bus", "B1")
+    n.add("Generator", "gen", bus="B1", p_nom=200, marginal_cost=10)
+    n.add("Load", "L_consume", bus="B1", p_set=100, sign=-1)
+    n.add("Load", "L_supply", bus="B1", p_set=100, sign=1)
+
+    n.optimize(solver_name="highs")
+
+    # Constraint RHS check: -p_set * sign per load, summed per bus
+    # Expected for bus B1: -(100 * -1) + -(100 * 1) = 100 - 100 = 0
+    rhs = n.model.constraints["Bus-nodal_balance"].rhs.sel(name="B1")
+    assert (rhs == 0).all().item()
+
+    # Generator stays idle
+    assert n.c.generators.dynamic.p["gen"].abs().max() < TOLERANCE
+    assert n.objective < TOLERANCE
+
+
+def test_nodal_balance_load_sign_rhs_values():
+    """
+    Verify the per-load contribution to the nodal balance RHS for both
+    sign conventions, independently of solver behavior.
+    """
+    n = pypsa.Network()
+    n.add("Bus", "B1")
+    n.add("Bus", "B2")
+    n.add("Generator", "gen1", bus="B1", p_nom=500, marginal_cost=1)
+    n.add("Generator", "gen2", bus="B2", p_nom=500, marginal_cost=1)
+    n.add("Load", "L_default", bus="B1", p_set=80)  # default sign = -1
+    n.add("Load", "L_supply", bus="B2", p_set=30, sign=1)
+
+    n.optimize.create_model()
+
+    rhs = n.model.constraints["Bus-nodal_balance"].rhs
+    # Default sign=-1: RHS contribution = -p_set * (-1) = +p_set
+    assert rhs.sel(name="B1").item() == 80.0
+    # sign=+1: RHS contribution = -p_set * (+1) = -p_set
+    assert rhs.sel(name="B2").item() == -30.0
+
+
+def test_define_fixed_operational_constraints_infinite_p_nom():
+    """Non-extendable generator with p_nom=inf and p_min_pu=0 must not produce NaN bounds.
+
+    Without the inf*0 fallback, the lower bound would be NaN and the solve would fail.
+    """
+    n = pypsa.Network()
+    n.add("Bus", "bus0")
+    n.add("Load", "load0", bus="bus0", p_set=10)
+    n.add(
+        "Generator",
+        "gen_inf",
+        bus="bus0",
+        p_nom=float("inf"),
+        p_min_pu=0,
+        p_max_pu=1,
+        marginal_cost=1,
+    )
+
+    n.optimize.create_model()
+    lower = n.model.constraints["Generator-fix-p-lower"].rhs
+    assert (lower == 0).all()
+
+    n.optimize.solve_model()
+    assert n.c.generators.dynamic.p["gen_inf"].eq(10).all()