LightSim2Grid 0.12.1__tar.gz → 0.12.2__tar.gz

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2/LightSim2Grid.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: LightSim2Grid
-Version: 0.12.1
+Version: 0.12.2
 Summary: LightSim2Grid implements a c++ backend targeting the Grid2Op platform.
 Home-page: https://github.com/grid2op/lightsim2grid/
 Author: Benjamin Donnot

{lightsim2grid-0.12.1/LightSim2Grid.egg-info → lightsim2grid-0.12.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: LightSim2Grid
-Version: 0.12.1
+Version: 0.12.2
 Summary: LightSim2Grid implements a c++ backend targeting the Grid2Op platform.
 Home-page: https://github.com/grid2op/lightsim2grid/
 Author: Benjamin Donnot

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/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.12.1"
+__version__ = "0.12.2"
 
 __all__ = [
     "newtonpf",

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/lightsim2grid/gridmodel/from_pandapower/_aux_add_dc_line.py

@@ -38,7 +38,7 @@ def _aux_add_dc_line(model, pp_net, pp_to_ls):
     branch_from_id = pp_bus_to_ls(pp_net.dcline["from_bus"].values, pp_to_ls)
     branch_to_id = pp_bus_to_ls(pp_net.dcline["to_bus"].values, pp_to_ls)
     p_mw = -pp_net.dcline["p_mw"].values
-    if np.any(~np.isfinite(p_mw)):
+    if (~np.isfinite(p_mw)).any():
         warnings.warn("Some non finite values are found for p_mw, they have been replaced by 0.")
         p_mw[~np.isfinite(p_mw)] = 0.
         
@@ -47,20 +47,20 @@ def _aux_add_dc_line(model, pp_net, pp_to_ls):
     vm_or_pu = pp_net.dcline["vm_from_pu"].values
     vm_ex_pu = pp_net.dcline["vm_to_pu"].values
     
-    min_q_or = pp_net.dcline["min_q_from_mvar"].values
-    if np.any(~np.isfinite(min_q_or)):
+    min_q_or = pp_net.dcline["min_q_from_mvar"].values.copy()
+    if (~np.isfinite(min_q_or)).any():
         min_q_or[~np.isfinite(min_q_or)] = -SOME_KIND_OF_INF_FOR_PMIN_PMAX
         
-    max_q_or = pp_net.dcline["max_q_from_mvar"].values
-    if np.any(~np.isfinite(max_q_or)):
+    max_q_or = pp_net.dcline["max_q_from_mvar"].values.copy()
+    if (~np.isfinite(max_q_or)).any():
         max_q_or[~np.isfinite(max_q_or)] = SOME_KIND_OF_INF_FOR_PMIN_PMAX
         
-    min_q_ex = pp_net.dcline["min_q_to_mvar"].values
-    if np.any(~np.isfinite(min_q_ex)):
+    min_q_ex = pp_net.dcline["min_q_to_mvar"].values.copy()
+    if (~np.isfinite(min_q_ex)).any():
         min_q_ex[~np.isfinite(min_q_ex)] = -SOME_KIND_OF_INF_FOR_PMIN_PMAX
         
-    max_q_ex = pp_net.dcline["max_q_to_mvar"].values
-    if np.any(~np.isfinite(max_q_ex)):
+    max_q_ex = pp_net.dcline["max_q_to_mvar"].values.copy()
+    if (~np.isfinite(max_q_ex)).any():
         max_q_ex[~np.isfinite(max_q_ex)] = +SOME_KIND_OF_INF_FOR_PMIN_PMAX
                           
     model.init_dclines(branch_from_id.astype(np.int32),

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/lightsim2grid/gridmodel/from_pandapower/_aux_add_load.py

@@ -29,7 +29,7 @@ def _aux_add_load(model, pp_net, pp_to_ls):
 
     ratio = 1.0
     if "scaling" in pp_net.load:
-        ratio = pp_net.load["scaling"].values
+        ratio = pp_net.load["scaling"].values.copy()
         ratio[~np.isfinite(ratio)] = 1.0
         
     model.init_loads(pp_net.load["p_mw"].values * ratio,

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/lightsim2grid/gridmodel/from_pandapower/_aux_add_sgen.py

@@ -35,7 +35,7 @@ def _aux_add_sgen(model, pp_net, pp_to_ls):
         return
 
     if "min_p_mw" in pp_net.sgen:
-        min_p_mw = pp_net.sgen["min_p_mw"].values
+        min_p_mw = pp_net.sgen["min_p_mw"].values.copy()
     else:
         min_p_mw = np.zeros(pp_net.sgen.shape[0]) - SOME_KIND_OF_INF_FOR_PMIN_PMAX
     if np.any(~np.isfinite(min_p_mw)):
@@ -43,7 +43,7 @@ def _aux_add_sgen(model, pp_net, pp_to_ls):
     min_p_mw[~np.isfinite(min_p_mw)] = 0.
 
     if "max_p_mw" in pp_net.sgen:
-        max_p_mw = pp_net.sgen["max_p_mw"].values
+        max_p_mw = pp_net.sgen["max_p_mw"].values.copy()
     else:
         max_p_mw = np.zeros(pp_net.sgen.shape[0]) + SOME_KIND_OF_INF_FOR_PMIN_PMAX
     if np.any(~np.isfinite(max_p_mw)):
@@ -51,7 +51,7 @@ def _aux_add_sgen(model, pp_net, pp_to_ls):
     max_p_mw[~np.isfinite(max_p_mw)] = 0.
 
     if "min_q_mvar" in pp_net.sgen:
-        min_q_mvar = pp_net.sgen["min_q_mvar"].values
+        min_q_mvar = pp_net.sgen["min_q_mvar"].values.copy()
     else:
         min_q_mvar = np.zeros(pp_net.sgen.shape[0]) - SOME_KIND_OF_INF_FOR_PMIN_PMAX
     if np.any(~np.isfinite(min_q_mvar)):
@@ -59,7 +59,7 @@ def _aux_add_sgen(model, pp_net, pp_to_ls):
     min_q_mvar[~np.isfinite(min_q_mvar)] = 0.
 
     if "max_q_mvar" in pp_net.sgen:
-        max_q_mvar = pp_net.sgen["max_q_mvar"].values
+        max_q_mvar = pp_net.sgen["max_q_mvar"].values.copy()
     else:
         max_q_mvar = np.zeros(pp_net.sgen.shape[0]) + SOME_KIND_OF_INF_FOR_PMIN_PMAX
     if np.any(~np.isfinite(max_q_mvar)):

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/lightsim2grid/gridmodel/from_pandapower/_aux_add_slack.py

@@ -96,7 +96,7 @@ def _aux_add_slack(
             # so I assume i need to add as many generators as number of slack bus
             nb_slack = len(slack_bus_ids)
             if "slack_weight" in pp_net.ext_grid:
-                slack_coeff = 1.0 * pp_net.ext_grid["slack_weight"].values
+                slack_coeff = 1.0 * pp_net.ext_grid["slack_weight"].values.copy()
             else:
                 slack_coeff = np.ones(nb_slack)
             slack_coeff_norm = slack_coeff / slack_coeff.sum()

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/lightsim2grid/gridmodel/from_pandapower/_aux_add_trafo.py

@@ -40,7 +40,7 @@ def _aux_add_trafo(
                            "Some pp_net.trafo[\"parallel\"] != 1 it is not handled by lightsim yet.")
 
     # fix the missing values
-    tap_neutral = 1.0 * pp_net.trafo["tap_neutral"].values
+    tap_neutral = pp_net.trafo["tap_neutral"].values.copy()
     if np.any(~np.isfinite(tap_neutral)):
         warnings.warn("There were some Nan in the pp_net.trafo[\"tap_neutral\"], they have been replaced by 0")
     tap_neutral[~np.isfinite(tap_neutral)] = 0.
@@ -48,17 +48,17 @@ def _aux_add_trafo(
     if np.any(tap_neutral != 0.):
         raise RuntimeError("lightsim converter supposes that tap_neutral is 0 for the transformers")
 
-    tap_step_pct = 1.0 * pp_net.trafo["tap_step_percent"].values
+    tap_step_pct = pp_net.trafo["tap_step_percent"].values.copy()
     if np.any(~np.isfinite(tap_step_pct)):
         warnings.warn("There were some Nan in the pp_net.trafo[\"tap_step_percent\"], they have been replaced by 0")
     tap_step_pct[~np.isfinite(tap_step_pct)] = 0.
 
-    tap_pos = 1.0 * pp_net.trafo["tap_pos"].values
+    tap_pos = pp_net.trafo["tap_pos"].values.copy()
     if np.any(~np.isfinite(tap_pos)):
         warnings.warn("There were some Nan in the pp_net.trafo[\"tap_pos\"], they have been replaced by 0")
     tap_pos[~np.isfinite(tap_pos)] = 0.
 
-    shift_ = 1.0 * pp_net.trafo["shift_degree"].values
+    shift_ = pp_net.trafo["shift_degree"].values.copy()
     if np.any(~np.isfinite(tap_pos)):
         warnings.warn("There were some Nan in the pp_net.trafo[\"shift_degree\"], they have been replaced by 0")
     shift_[~np.isfinite(shift_)] = 0.
@@ -81,7 +81,7 @@ def _aux_add_trafo(
             raise RuntimeError("Ideal phase shifters are not modeled. Please remove all 3-winding trafos "
                                "with \"tap_changer_type\" set to \"Ideal\".")
             
-    tap_angles_ = 1.0 * pp_net.trafo["tap_step_degree"].values
+    tap_angles_ = pp_net.trafo["tap_step_degree"].values
     if np.any(~np.isfinite(tap_angles_)):
         warnings.warn("There were some Nan in the pp_net.trafo[\"tap_step_degree\"], they have been replaced by 0")
     tap_angles_[~np.isfinite(tap_angles_)] = 0.

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/lightsim2grid/lightSimBackend.py

@@ -1165,10 +1165,10 @@ class LightSimBackend(Backend):
         self.name_sub = pp_cls.name_sub
         self._grid.set_substation_names(self.name_sub)
 
-        self.prod_pu_to_kv = self._init_pp_backend.prod_pu_to_kv
-        self.load_pu_to_kv = self._init_pp_backend.load_pu_to_kv
-        self.lines_or_pu_to_kv = self._init_pp_backend.lines_or_pu_to_kv
-        self.lines_ex_pu_to_kv = self._init_pp_backend.lines_ex_pu_to_kv
+        self.prod_pu_to_kv = self._init_pp_backend.prod_pu_to_kv.copy()
+        self.load_pu_to_kv = self._init_pp_backend.load_pu_to_kv.copy()
+        self.lines_or_pu_to_kv = self._init_pp_backend.lines_or_pu_to_kv.copy()
+        self.lines_ex_pu_to_kv = self._init_pp_backend.lines_ex_pu_to_kv.copy()
 
         # TODO storage check grid2op version and see if storage is available !
         if self.__has_storage:
@@ -1364,13 +1364,14 @@ class LightSimBackend(Backend):
         # TODO speed optimization here to read it "better"
         cls = type(self)
         # handle object in topo vect
-        res[cls.load_pos_topo_vect] = cls.global_bus_to_local(np.array([el.bus_id for el in self._grid.get_loads()]),
-                                                                        cls.load_to_subid)
-        res[cls.gen_pos_topo_vect] = cls.global_bus_to_local(np.array([el.bus_id for el in self._grid.get_generators()]),
-                                                                        cls.gen_to_subid)
+        res[cls.load_pos_topo_vect] = cls.global_bus_to_local(
+            np.array([el.bus_id for el in self._grid.get_loads()]), cls.load_to_subid)
+        res[cls.gen_pos_topo_vect] = cls.global_bus_to_local(
+            np.array([el.bus_id for el in self._grid.get_generators()]), cls.gen_to_subid)
         if self.__has_storage:
-            res[cls.storage_pos_topo_vect] = cls.global_bus_to_local(np.array([el.bus_id for el in self._grid.get_storages()]),
-                                                                                cls.storage_to_subid)
+            res[cls.storage_pos_topo_vect] = cls.global_bus_to_local(
+                np.array([el.bus_id for el in self._grid.get_storages()]), cls.storage_to_subid)
+            
         lor_glob_bus = np.concatenate((np.array([el.bus1_id for el in self._grid.get_lines()]),
                                         np.array([el.bus1_id for el in self._grid.get_trafos()])))
         res[cls.line_or_pos_topo_vect] = cls.global_bus_to_local(lor_glob_bus, cls.line_or_to_subid)
@@ -1382,7 +1383,8 @@ class LightSimBackend(Backend):
         # handle shunts (not in topo vect)
         if cls.shunts_data_available:
             self.sh_bus.flags.writeable = True
-            self.sh_bus[:] = cls.global_bus_to_local(self._grid.get_shunts().get_bus_id(),
+            sh_bus_global = np.asarray([el.bus_id for el in self._grid.get_shunts()])
+            self.sh_bus[:] = cls.global_bus_to_local(sh_bus_global,
                                                      cls.shunt_to_subid).copy()
             self.sh_bus.flags.writeable = False
         
@@ -1402,14 +1404,6 @@ class LightSimBackend(Backend):
             self._init_pp_backend.__class__ = type(self._init_pp_backend).init_grid(type(self))
         self._backend_action_class = _BackendAction.init_grid(type(self))
         
-        # self._init_action_to_set = self._backend_action_class()
-        # try:
-        #     _init_action_to_set = self.get_action_to_set()
-        # except TypeError:
-        #     # I am in legacy grid2op version...
-        #     _init_action_to_set = _dont_use_get_action_to_set_legacy(self)
-            
-        # self._init_action_to_set += _init_action_to_set
         if self.prod_pu_to_kv is not None:
             assert np.isfinite(self.prod_pu_to_kv).all()
         if self.load_pu_to_kv is not None:
@@ -1421,6 +1415,13 @@ class LightSimBackend(Backend):
         if self.__has_storage and self.n_storage > 0 and self.storage_pu_to_kv is not None:
             assert np.isfinite(self.storage_pu_to_kv).all()
 
+        if self._dist_slack_non_renew and not type(self).redispatching_unit_commitment_availble:
+            raise Grid2OpException("You asked to distribute the slack on renewable generators, "
+                                   "yet the environment does not seem to know which generators are "
+                                   "renewables or not. Make sure env.redispatching_unit_commitment_availble is True "
+                                   "(this include having a 'prods_charac.csv' properly formatted and "
+                                   "at the right location)")
+            
     def _count_object_per_bus(self):
         # should be called only when self.topo_vect and self.shunt_topo_vect are set
         # todo factor that more properly to update it when it's modified, and not each time
@@ -1491,11 +1492,14 @@ class LightSimBackend(Backend):
                 self._need_islanding_detection = True
         
         # update the injections
+        gen_v = backendAction.prod_v.values.copy()
+        gen_v[~backendAction.prod_v.changed] = 0.  # sometimes there were nan here, which lead to a warnings
+        gen_v /= self.prod_pu_to_kv
         try:
             self._grid.update_gens_p(backendAction.prod_p.changed,
                                      backendAction.prod_p.values)
             self._grid.update_gens_v(backendAction.prod_v.changed,
-                                     backendAction.prod_v.values / self.prod_pu_to_kv)
+                                     gen_v)
             self._grid.update_loads_p(backendAction.load_p.changed,
                                       backendAction.load_p.values)
             self._grid.update_loads_q(backendAction.load_q.changed,

{lightsim2grid-0.12.1 → lightsim2grid-0.12.2}/setup.py

@@ -14,7 +14,7 @@ import warnings
 from pybind11.setup_helpers import Pybind11Extension, build_ext
 
 
-__version__ = "0.12.1"
+__version__ = "0.12.2"
 KLU_SOLVER_AVAILABLE = False
 
 # Try to link against SuiteSparse (if available)