LightSim2Grid 0.13.0__tar.gz → 0.13.1__tar.gz

{lightsim2grid-0.13.0 → lightsim2grid-0.13.1}/CHANGELOG.rst

@@ -28,7 +28,30 @@ TODO: in `main.cpp` check the returned policy of pybind11 and also the `py::call
 TODO: a cpp class that is able to compute (DC powerflow) ContingencyAnalysis and TimeSeries using PTDF and LODF
 TODO: integration test with pandapower (see `pandapower/contingency/contingency.py` and import `lightsim2grid_installed` and check it's True)
 
-[0.13.0] 2026-04-xx
+[0.13.1]  2026-04-21
+--------------------
+- [BREAKING] when loading a powergrid from pypowsybl with "use_buses_for_sub" tagged
+  and disconnected element on the grid will now raise a RuntimeError. Before there were
+  some "automatic" bahaviour to try to find a possible bus which could lead to 
+  error afterwards.
+- [BREAKING] adding a more precise information about linear solvers. The "refactor" timings
+  are now also available in solver.get_timers_jacobian() which now returns a tuple of size 10
+- [FIXED] an issue where disconnected powerlines could be tagged as "fakely connected"
+- [FIXED] some issues when loading a grid from pypowsybl in case of disconnected elements.
+- [FIXED] remove the undefined behaviour while maintaining compile time check to prevent 
+  wrong conversion from different bus labelling.
+- [IMPROVED] the CI to allow automatic push on pypi on new version tag (introduced in version 0.12.0)
+- [IMPROVED] reduce code duplication between ContingencyAnalysis and TimeSerie (cpp side)
+- [IMPROVED] handling of branches disconnected at only one side: less code duplication and
+  it should be working with TimeSeries and ContingencyAnalysis
+- [IMPROVED] speed (DC mode): avoid the systematic call to "refactor" when Ybus is not changed
+  when using DC approximation.
+- [IMPROVED] simplify the future integration of other linear solvers and the logic when linear_solvers
+  are called by decoupling "refactor" steps from "solve" steps (they used to be all under the same
+  "solve" method).
+- [IMPROVED] added "final" and "override" key-words on some methods in src/element_container
+  
+[0.13.0] 2026-04-15
 --------------------
 - [PENDING DEPRECATION] the cpp module (lightsim2grid_cpp) will not be usable directly anymore.
   This means that calls like "from lightsim2grid_cpp import XXX" will not work. To replace them 

{lightsim2grid-0.13.0 → lightsim2grid-0.13.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.2
 Name: lightsim2grid
-Version: 0.13.0
+Version: 0.13.1
 Summary: LightSim2grid implements a c++ backend targeting the Grid2Op platform.
 Author-Email: Benjamin Donnot <benjamin.donnot@rte-france.com>
 Classifier: Development Status :: 4 - Beta

{lightsim2grid-0.13.0 → lightsim2grid-0.13.1}/lightsim2grid/__init__.py

@@ -6,7 +6,7 @@
 # SPDX-License-Identifier: MPL-2.0
 # This file is part of LightSim2grid, LightSim2grid implements a c++ backend targeting the Grid2Op platform.
 
-__version__ = "0.13.0"
+__version__ = "0.13.1"
 
 __all__ = [
     "newtonpf",

{lightsim2grid-0.13.0 → lightsim2grid-0.13.1}/lightsim2grid/gridmodel/from_pypowsybl/_from_pypowsybl.py

@@ -22,25 +22,35 @@ PP_BUG_RATIO_TAP_CHANGER = version.parse("1.9")
 PYPOWSYBL_VER = version.parse(pypo.__version__)
 
 
-def _aux_get_bus(vl_df, bus_df, df, conn_key="connected", bus_key="bus_id", vl_key="voltage_level_id"):
+def _aux_get_bus(vl_df, bus_df, first_bus_per_vl, el_type, df, conn_key="connected", bus_key="bus_id", vl_key="voltage_level_id"):
     if df.shape[0] == 0:
         # no element of this type so no problem
         return np.zeros(0, dtype=int), np.ones(0, dtype=bool), np.zeros(0, dtype=int)
     # retrieve which elements are disconnected 
     mask_disco = ~df[conn_key]
-    
+    if mask_disco.any() and first_bus_per_vl is None:
+        raise RuntimeError(f"Some elements of type {el_type} are disconnected (no bus attached to them) and you "
+                           "are in 'use pypowsybl bus as grid2op substation' mode. The converter cannot assign "
+                           f"substation for {el_type}:\n {mask_disco[mask_disco].index.to_numpy()}\n, this is not supported. "
+                           "Please upgrade the pypowsybl grid (*eg* grid.xiidm of the grid2op env) "
+                           "and manually reconnect any non connected elements there.")
     # retrieve the bus where the element are
     tmp_bus_id = df[bus_key].copy()
-    # element disconnected are, by default assigned to first bus of their substation
-    tmp_disco = vl_df.loc[df.loc[mask_disco, vl_key], "vl_id"].values
-    bus_el_disco = bus_df.iloc[tmp_disco]
-    tmp_bus_id[mask_disco] = bus_el_disco.index
+    
+    if mask_disco.any():
+        # element disconnected are, by default assigned to first bus of their substation
+        el_disco = vl_df.loc[df.loc[mask_disco, vl_key]].index
+        tmp_bus_id.loc[mask_disco] = first_bus_per_vl.loc[el_disco, "first_bus_name"].values
+        
+    # assign bus id
     bus_id = bus_df.loc[tmp_bus_id.values]["bus_global_id"].values.copy()
-    # deactivate the element not on the main component
-    # wrong_component = bus_df.loc[tmp_bus_id.values]["connected_component"].values != 0
-    # mask_disco[wrong_component] = True
-    # assign bus -1 to disconnected elements
-    bus_id[mask_disco] = -1
+        
+    if mask_disco.any():
+        # deactivate the element not on the main component
+        # wrong_component = bus_df.loc[tmp_bus_id.values]["connected_component"].values != 0
+        # mask_disco[wrong_component] = True
+        # assign bus -1 to disconnected elements
+        bus_id[mask_disco] = -1
     
     sub_id = bus_df.loc[tmp_bus_id.values]["glop_sub_id"].values
     return bus_id, mask_disco.values, sub_id
@@ -166,8 +176,6 @@ def init(net : pypo.network.Network,
         voltage_levels = net.get_voltage_levels()
     
     all_buses_vn_kv = voltage_levels.loc[bus_df["voltage_level_id"].values]["nominal_v"].values
-    sub_unique = None
-    sub_unique_id = None
     nb_bus_per_vl = bus_df[["voltage_level_id", "name"]].groupby("voltage_level_id").count()
     
     if buses_for_sub is not None and buses_for_sub:
@@ -192,6 +200,7 @@ def init(net : pypo.network.Network,
         bus_df["glop_sub_id"] = np.arange(n_sub_ls)  # np.concatenate([np.arange(n_sub_ls) for _ in range(n_busbar_per_sub)])
         sub_names = bus_df.index.values.astype(str)
         voltage_levels["vl_id"] = bus_df[["voltage_level_id", "bus_global_id"]].groupby("voltage_level_id").min()
+        first_bus_per_vl = None
     else:        
         # the "substation" in lightsim2grid
         voltage_levels["nb_bus_per_vl"] = nb_bus_per_vl["name"]        
@@ -230,6 +239,36 @@ def init(net : pypo.network.Network,
         sub_names = voltage_levels.index.values.astype(str)
         bus_df["glop_sub_id"] = voltage_levels.loc[bus_df["voltage_level_id"].values, "glop_sub_id"].values
         
+        # retrieve the "first bus of every substation"
+        # this is used to connected reconnected element
+        first_bus_per_vl = bus_df[["glop_sub_id", "voltage_level_id", "name"]].groupby("glop_sub_id").first().reset_index().set_index("voltage_level_id")
+        first_bus_per_vl["first_bus_name"] = [el[0] for el in first_bus_per_vl["name"].values]
+        first_bus_per_vl.drop(["glop_sub_id", "name"], axis=1, inplace=True)
+        
+        # make sure every voltage level as a "first bus"
+        # I would raise an error in the case a voltage is fully disconnected but...
+        check_grid_ok = voltage_levels.index.isin(first_bus_per_vl.index)
+        if np.any(~check_grid_ok):
+            warnings.warn("There are some voltage levels without any connected buses. "
+                          f"Check voltage levels {voltage_levels[~check_grid_ok].index}")
+            # name_added_bus = voltage_levels[~check_grid_ok].index
+            nm_vl_without_bus = voltage_levels.loc[~check_grid_ok, ["vl_id", "glop_sub_id"]]
+            to_add = voltage_levels.loc[~check_grid_ok, ["vl_id", "glop_sub_id"]].copy()
+            to_add.index = to_add.index.astype(str) + "added_bus"
+            to_add["name"] = [[el] for el in to_add.index.astype(str)]
+            to_add["v_mag"] = np.nan
+            to_add["v_angle"] = np.nan
+            to_add["connected_component"] = 99999
+            to_add["synchronous_component"] = 99999
+            to_add["voltage_level_id"] = nm_vl_without_bus.index
+            to_add["orig_id"] = 9999
+            to_add["local_id"] = 1
+            to_add["bus_global_id"] = to_add["glop_sub_id"]
+            for added_bus_nm in to_add.index:
+                bus_df.loc[added_bus_nm] = to_add.loc[added_bus_nm]
+            for vl_bus_added in nm_vl_without_bus.index:
+                first_bus_per_vl.loc[vl_bus_added, "first_bus_name"] = vl_bus_added + "added_bus"
+            
     # all_buses_vn_kv = np.concatenate([all_buses_vn_kv for _ in range(n_busbar_per_sub)])
     model.init_bus(n_sub_ls,
                    n_busbar_per_sub_ls,
@@ -260,7 +299,7 @@ def init(net : pypo.network.Network,
     max_q = max_q_aux.astype(np.float32)
     min_q[~np.isfinite(min_q)] = min_float_value
     max_q[~np.isfinite(max_q)] = max_float_value
-    gen_bus, gen_disco, gen_sub = _aux_get_bus(voltage_levels, bus_df, df_gen)
+    gen_bus, gen_disco, gen_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "gen", df_gen)
 
     # dirty fix for when regulating elements are not the same
     bus_reg = copy.deepcopy(df_gen["regulated_element_id"].values)
@@ -294,7 +333,7 @@ def init(net : pypo.network.Network,
         df_load = net.get_loads().sort_index()
     else:
         df_load = net.get_loads()
-    load_bus, load_disco, load_sub = _aux_get_bus(voltage_levels, bus_df, df_load)
+    load_bus, load_disco, load_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "load", df_load)
     model.init_loads(df_load["p0"].values,
                      df_load["q0"].values,
                      load_bus
@@ -326,8 +365,8 @@ def init(net : pypo.network.Network,
     line_x = df_line_pu["x"].values
     line_h_or = (df_line_pu["g1"].values + 1j * df_line_pu["b1"].values)
     line_h_ex = (df_line_pu["g2"].values + 1j * df_line_pu["b2"].values)
-    lor_bus, lor_disco, lor_sub = _aux_get_bus(voltage_levels, bus_df, df_line, conn_key="connected1", bus_key="bus1_id", vl_key="voltage_level1_id")
-    lex_bus, lex_disco, lex_sub = _aux_get_bus(voltage_levels, bus_df, df_line, conn_key="connected2", bus_key="bus2_id", vl_key="voltage_level2_id")
+    lor_bus, lor_disco, lor_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "line (side 1)", df_line, conn_key="connected1", bus_key="bus1_id", vl_key="voltage_level1_id")
+    lex_bus, lex_disco, lex_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "line (side 2)", df_line, conn_key="connected2", bus_key="bus2_id", vl_key="voltage_level2_id")
     model.init_powerlines_full(line_r,
                                line_x,
                                line_h_or,
@@ -392,8 +431,8 @@ def init(net : pypo.network.Network,
             # bug in per unit view in both python and java
             ratio[has_r_tap_changer] = 1. * ratio_tap_changer.loc[df_trafo_pu.loc[has_r_tap_changer].index, "rho"].values
 
-    tor_bus, tor_disco, tor_sub = _aux_get_bus(voltage_levels, bus_df, df_trafo, conn_key="connected1", bus_key="bus1_id", vl_key="voltage_level1_id")
-    tex_bus, tex_disco, tex_sub = _aux_get_bus(voltage_levels, bus_df, df_trafo, conn_key="connected2", bus_key="bus2_id", vl_key="voltage_level2_id")
+    tor_bus, tor_disco, tor_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "trafo (side 1)", df_trafo, conn_key="connected1", bus_key="bus1_id", vl_key="voltage_level1_id")
+    tex_bus, tex_disco, tex_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "trafo (side 2)", df_trafo, conn_key="connected2", bus_key="bus2_id", vl_key="voltage_level2_id")
     model.init_trafo(trafo_r,
                      trafo_x,
                      trafo_h,
@@ -415,7 +454,7 @@ def init(net : pypo.network.Network,
     else:
         df_shunt = net.get_shunt_compensators()
         
-    sh_bus, sh_disco, sh_sub = _aux_get_bus(voltage_levels, bus_df, df_shunt)    
+    sh_bus, sh_disco, sh_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "shunts", df_shunt)    
     shunt_kv = voltage_levels.loc[df_shunt["voltage_level_id"].values]["nominal_v"].values
     model.init_shunt(df_shunt["g"].values * shunt_kv**2,
                      -df_shunt["b"].values * shunt_kv**2,
@@ -433,8 +472,8 @@ def init(net : pypo.network.Network,
     else:
         df_dc = net.get_hvdc_lines()
         df_sations = net.get_vsc_converter_stations()
-    hvdc_bus_from_id, hvdc_from_disco, hvdc_sub_from_id  = _aux_get_bus(voltage_levels, bus_df, df_sations.loc[df_dc["converter_station1_id"].values]) 
-    hvdc_bus_to_id, hvdc_to_disco, hvdc_sub_to_id = _aux_get_bus(voltage_levels, bus_df, df_sations.loc[df_dc["converter_station2_id"].values]) 
+    hvdc_bus_from_id, hvdc_from_disco, hvdc_sub_from_id  = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "hvdc (side 1)", df_sations.loc[df_dc["converter_station1_id"].values]) 
+    hvdc_bus_to_id, hvdc_to_disco, hvdc_sub_to_id = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "hvdc (side 2)", df_sations.loc[df_dc["converter_station2_id"].values]) 
     loss_percent = np.zeros(df_dc.shape[0])  # TODO 
     loss_mw = np.zeros(df_dc.shape[0])  # TODO
     model.init_dclines(hvdc_bus_from_id,
@@ -460,7 +499,7 @@ def init(net : pypo.network.Network,
         df_batt = net.get_batteries().sort_index()
     else:
         df_batt = net.get_batteries()
-    batt_bus, batt_disco, batt_sub = _aux_get_bus(voltage_levels, bus_df, df_batt)
+    batt_bus, batt_disco, batt_sub = _aux_get_bus(voltage_levels, bus_df, first_bus_per_vl, "storage", df_batt)
     model.init_storages(df_batt["target_p"].values,
                         df_batt["target_q"].values,
                         batt_bus

lightsim2grid-0.13.1/lightsim2grid/tests/case_14_iidm/grid copy.xiidm

@@ -0,0 +1,147 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<iidm:network xmlns:iidm="http://www.powsybl.org/schema/iidm/1_10" xmlns:slt="http://www.powsybl.org/schema/iidm/ext/slack_terminal/1_5" id="ieee14cdf" caseDate="1993-08-19T00:00:00.000Z" forecastDistance="0" sourceFormat="IEEE-CDF" minimumValidationLevel="STEADY_STATE_HYPOTHESIS">
+    <iidm:substation id="S1">
+        <iidm:voltageLevel id="VL1" nominalV="135.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B1" name="Bus 1     HV" v="143.1" angle="0.0"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="B1-G" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="232.4" targetV="143.1" targetQ="-16.9" bus="B1" connectableBus="B1">
+                <iidm:minMaxReactiveLimits minQ="-1.7976931348623157E308" maxQ="1.7976931348623157E308"/>
+            </iidm:generator>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S2">
+        <iidm:voltageLevel id="VL2" nominalV="135.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B2" name="Bus 2     HV" v="141.075" angle="-4.98"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="B2-G" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="40.0" targetV="141.075" targetQ="42.4" bus="B2" connectableBus="B2">
+                <iidm:minMaxReactiveLimits minQ="-40.0" maxQ="50.0"/>
+            </iidm:generator>
+            <iidm:load id="B2-L" loadType="UNDEFINED" p0="21.7" q0="12.7" bus="B2" connectableBus="B2"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S3">
+        <iidm:voltageLevel id="VL3" nominalV="135.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B3" name="Bus 3     HV" v="136.35" angle="-12.72"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="B3-G" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="0.0" targetV="136.35" targetQ="23.4" bus="B3" connectableBus="B3">
+                <iidm:minMaxReactiveLimits minQ="0.0" maxQ="40.0"/>
+            </iidm:generator>
+            <iidm:load id="B3-L" loadType="UNDEFINED" p0="94.2" q0="19.0" bus="B3" connectableBus="B3"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S4">
+        <iidm:voltageLevel id="VL4" nominalV="135.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B4" name="Bus 4     HV" v="137.565" angle="-10.33"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B4-L" loadType="UNDEFINED" p0="47.8" q0="-3.9" bus="B4" connectableBus="B4"/>
+        </iidm:voltageLevel>
+        <iidm:voltageLevel id="VL7" nominalV="14.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B7" name="Bus 7     ZV" v="14.868" angle="-13.37"/>
+            </iidm:busBreakerTopology>
+        </iidm:voltageLevel>
+        <iidm:voltageLevel id="VL9" nominalV="12.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B9" name="Bus 9     LV" v="12.672" angle="-14.94"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B9-L" loadType="UNDEFINED" p0="29.5" q0="16.6" bus="B9" connectableBus="B9"/>
+            <iidm:shunt id="B9-SH" sectionCount="1" voltageRegulatorOn="false" bus="B9" connectableBus="B9">
+                <iidm:shuntLinearModel bPerSection="0.13194444444444445" maximumSectionCount="1"/>
+            </iidm:shunt>
+        </iidm:voltageLevel>
+        <iidm:twoWindingsTransformer id="T4-7-1" r="0.0" x="0.4098752" g="0.0" b="0.0" ratedU1="132.03" ratedU2="14.0" bus1="B4" connectableBus1="B4" voltageLevelId1="VL4" bus2="B7" connectableBus2="B7" voltageLevelId2="VL7"/>
+        <iidm:twoWindingsTransformer id="T4-9-1" r="0.0" x="0.8008992" g="0.0" b="0.0" ratedU1="130.815" ratedU2="12.0" bus1="B4" connectableBus1="B4" voltageLevelId1="VL4" bus2="B9" connectableBus2="B9" voltageLevelId2="VL9"/>
+    </iidm:substation>
+    <iidm:substation id="S5">
+        <iidm:voltageLevel id="VL5" nominalV="135.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B5" name="Bus 5     HV" v="137.7" angle="-8.78"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B5-L" loadType="UNDEFINED" p0="7.6" q0="1.6" bus="B5" connectableBus="B5"/>
+        </iidm:voltageLevel>
+        <iidm:voltageLevel id="VL6" nominalV="12.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B6" name="Bus 6     LV" v="12.84" angle="-14.22"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="B6-G" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="0.0" targetV="12.84" targetQ="12.2" bus="B6" connectableBus="B6">
+                <iidm:minMaxReactiveLimits minQ="-6.0" maxQ="24.0"/>
+            </iidm:generator>
+            <iidm:load id="B6-L" loadType="UNDEFINED" p0="11.2" q0="7.5" bus="B6" connectableBus="B6"/>
+        </iidm:voltageLevel>
+        <iidm:twoWindingsTransformer id="T5-6-1" r="0.0" x="0.36290880000000003" g="0.0" b="0.0" ratedU1="125.82000000000001" ratedU2="12.0" bus1="B5" connectableBus1="B5" voltageLevelId1="VL5" bus2="B6" connectableBus2="B6" voltageLevelId2="VL6"/>
+    </iidm:substation>
+    <iidm:substation id="S8">
+        <iidm:voltageLevel id="VL8" nominalV="20.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B8" name="Bus 8     TV" v="21.8" angle="-13.36"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="B8-G" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="0.0" targetV="21.8" targetQ="17.4" bus="B8" connectableBus="B8">
+                <iidm:minMaxReactiveLimits minQ="-6.0" maxQ="24.0"/>
+            </iidm:generator>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S10">
+        <iidm:voltageLevel id="VL10" nominalV="12.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B10" name="Bus 10    LV" v="12.611999999999998" angle="-15.1"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B10-L" loadType="UNDEFINED" p0="9.0" q0="5.8" bus="B10" connectableBus="B10"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S11">
+        <iidm:voltageLevel id="VL11" nominalV="12.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B11" name="Bus 11    LV" v="12.684" angle="-14.79"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B11-L" loadType="UNDEFINED" p0="3.5" q0="1.8" bus="B11" connectableBus="B11"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S12">
+        <iidm:voltageLevel id="VL12" nominalV="12.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B12" name="Bus 12    LV" v="12.66" angle="-15.07"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B12-L" loadType="UNDEFINED" p0="6.1" q0="1.6" bus="B12" connectableBus="B12"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S13">
+        <iidm:voltageLevel id="VL13" nominalV="12.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B13" name="Bus 13    LV" v="12.600000000000001" angle="-15.16"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B13-L" loadType="UNDEFINED" p0="13.5" q0="5.8" bus="B13" connectableBus="B13"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="S14">
+        <iidm:voltageLevel id="VL14" nominalV="12.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="B14" name="Bus 14    LV" v="12.432" angle="-16.04"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="B14-L" loadType="UNDEFINED" p0="14.9" q0="5.0" bus="B14" connectableBus="B14"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:line id="L1-2-1" r="3.5320050000000003" x="10.7837325" g1="0.0" b1="1.448559670781893E-4" g2="0.0" b2="1.448559670781893E-4" bus1="B1" connectableBus1="B1" voltageLevelId1="VL1" bus2="B2" connectableBus2="B2" voltageLevelId2="VL2"/>
+    <iidm:line id="L1-5-1" r="9.8469675" x="40.64904" g1="0.0" b1="1.3497942386831276E-4" g2="0.0" b2="1.3497942386831276E-4" bus1="B1" connectableBus1="B1" voltageLevelId1="VL1" bus2="B5" connectableBus2="B5" voltageLevelId2="VL5"/>
+    <iidm:line id="L2-3-1" r="8.563927499999998" x="36.0800325" g1="0.0" b1="1.2016460905349794E-4" g2="0.0" b2="1.2016460905349794E-4" bus1="B2" connectableBus1="B2" voltageLevelId1="VL2" bus2="B3" connectableBus2="B3" voltageLevelId2="VL3"/>
+    <iidm:line id="L2-4-1" r="10.5905475" x="32.13432" g1="0.0" b1="9.327846364883403E-5" g2="0.0" b2="9.327846364883403E-5" bus1="B2" connectableBus1="B2" voltageLevelId1="VL2" bus2="B4" connectableBus2="B4" voltageLevelId2="VL4"/>
+    <iidm:line id="L2-5-1" r="10.379137499999999" x="31.68963" g1="0.0" b1="9.492455418381344E-5" g2="0.0" b2="9.492455418381344E-5" bus1="B2" connectableBus1="B2" voltageLevelId1="VL2" bus2="B5" connectableBus2="B5" voltageLevelId2="VL5"/>
+    <iidm:line id="L3-4-1" r="12.2125725" x="31.170217499999993" g1="0.0" b1="3.511659807956104E-5" g2="0.0" b2="3.511659807956104E-5" bus1="B3" connectableBus1="B3" voltageLevelId1="VL3" bus2="B4" connectableBus2="B4" voltageLevelId2="VL4"/>
+    <iidm:line id="L4-5-1" r="2.4330375" x="7.6745475" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B4" connectableBus1="B4" voltageLevelId1="VL4" bus2="B5" connectableBus2="B5" voltageLevelId2="VL5"/>
+    <iidm:line id="L6-11-1" r="0.13677119999999998" x="0.286416" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B6" connectableBus1="B6" voltageLevelId1="VL6" bus2="B11" connectableBus2="B11" voltageLevelId2="VL11"/>
+    <iidm:line id="L6-12-1" r="0.1769904" x="0.3683664" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B6" connectableBus1="B6" voltageLevelId1="VL6" bus2="B12" connectableBus2="B12" voltageLevelId2="VL12"/>
+    <iidm:line id="L6-13-1" r="0.09525600000000001" x="0.18758879999999997" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B6" connectableBus1="B6" voltageLevelId1="VL6" bus2="B13" connectableBus2="B13" voltageLevelId2="VL13"/>
+    <iidm:line id="L7-8-1" r="0.0" x="0.49322000000000005" g1="0.0" b1="-0.8689254867430936" g2="0.0" b2="0.6082478407201655" bus1="B7" connectableBus1="B7" voltageLevelId1="VL7" bus2="B8" connectableBus2="B8" voltageLevelId2="VL8"/>
+    <iidm:line id="L7-9-1" r="0.0" x="0.1848168" g1="0.0" b1="0.7729662176660504" g2="0.0" b2="-0.9017939206103922" bus1="B7" connectableBus1="B7" voltageLevelId1="VL7" bus2="B9" connectableBus2="B9" voltageLevelId2="VL9"/>
+    <iidm:line id="L9-10-1" r="0.04580639999999999" x="0.12168" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B9" connectableBus1="B9" voltageLevelId1="VL9" bus2="B10" connectableBus2="B10" voltageLevelId2="VL10"/>
+    <iidm:line id="L9-14-1" r="0.18303840000000002" x="0.3893472" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B9" connectableBus1="B9" voltageLevelId1="VL9" bus2="B14" connectableBus2="B14" voltageLevelId2="VL14"/>
+    <iidm:line id="L10-11-1" r="0.118152" x="0.27658079999999996" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B10" connectableBus1="B10" voltageLevelId1="VL10" bus2="B11" connectableBus2="B11" voltageLevelId2="VL11"/>
+    <iidm:line id="L12-13-1" r="0.3181248" x="0.28782719999999995" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B12" connectableBus1="B12" voltageLevelId1="VL12" bus2="B13" connectableBus2="B13" voltageLevelId2="VL13"/>
+    <iidm:line id="L13-14-1" r="0.2461392" x="0.5011488" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="B13" connectableBus1="B13" voltageLevelId1="VL13" bus2="B14" connectableBus2="B14" voltageLevelId2="VL14"/>
+    <iidm:extension id="VL1">
+        <slt:slackTerminal id="B1-G"/>
+    </iidm:extension>
+</iidm:network>

lightsim2grid-0.13.1/lightsim2grid/tests/dont_run_test_reco_deco_shunt.py

@@ -0,0 +1,53 @@
+import sys
+import grid2op
+from grid2op.tests.BaseBackendTest import (BaseTestNames,  # noqa: E402
+                                           BaseTestLoadingCase,
+                                           BaseTestLoadingBackendFunc, 
+                                           BaseTestTopoAction,
+                                           BaseTestEnvPerformsCorrectCascadingFailures,
+                                           BaseTestChangeBusAffectRightBus,
+                                        #    BaseTestShuntAction,
+                                           BaseTestResetEqualsLoadGrid,
+                                           BaseTestVoltageOWhenDisco,
+                                           BaseTestChangeBusSlack,
+                                           BaseIssuesTest,
+                                           BaseStatusActions,
+                                           MakeBackend, AlwaysLegal, CompleteAction, AmbiguousAction
+                                           )
+
+# gdb python --eval-command "run test_reco_deco_shunt.py"
+from lightsim2grid import LightSimBackend
+
+backend2 = LightSimBackend()
+print("MAKE HERE")
+env_change_q = grid2op.make(
+    "rte_case14_realistic",
+    test=True,
+    gamerules_class=AlwaysLegal,
+    action_class=CompleteAction,
+    backend=backend2,
+    _add_to_name="_3"
+)
+assert not env_change_q.done
+print("MAKE END")
+
+param = env_change_q.parameters
+param.NO_OVERFLOW_DISCONNECTION = True
+env_change_q.change_parameters(param)
+print("RESET HERE")
+env_change_q.reset(seed=0, options={"time serie id": 0})
+print("RESET END")
+
+# obs_disco_sh, *_ = env_change_q.step(
+#             env_change_q.action_space({"shunt": {"set_bus": [(0, -1)]}})
+#         )
+# obs_disco_sh, *_ = env_change_q.step(
+#             env_change_q.action_space({"shunt": {"set_bus": [(0, -1)]}})
+#         )
+
+
+obs_co_bus2_sh_alone, reward, done, info = env_change_q.step(
+    env_change_q.action_space({"shunt": {"set_bus": [(0, 2)]}})
+)
+print(f"{done = }")
+env_change_q.reset(seed=0, options={"time serie id": 0})