Merge branch 'bd-dev' into smarter_solver

Author: BDonnot

.gitignore

@@ -275,3 +275,6 @@ test.txt
 test_output.txt
 test_pypower_fdpf.py
 lightsim2grid/tests/_grid2op_for_test/
+bug_sparselu
+bug_sparselu_eigen.cpp
+test_segfault.sh

CHANGELOG.rst

@@ -20,6 +20,9 @@ Change Log
 
 [0.7.6] 2023-xx-yy
 --------------------
+- [BREAKING] now able to retrieve `dcSbus` with a dedicated method (and not with the old `get_Sbus`).
+  If you previously used `gridmodel.get_Subus()` to retrieve the Sbus used for DC powerflow, please use
+  `gridmodel.get_dcSbus()` instead.
 - [FIXED] now voltage is properly set to 0. when shunts are disconnected
 - [FIXED] now voltage is properly set to 0. when storage units are disconnected
 - [FIXED] a bug where non connected grid were not spotted in DC
@@ -31,6 +34,7 @@ Change Log
   one.
 - [ADDED] possibility to set / retrieve the names of each elements of the grid.
 - [ADDED] embed in the generator models the "non pv" behaviour. (TODO need to be able to change Q from python side)
+- [ADDED] computation of PTPF (Power Transfer Distribution Factor) is now possible
 - [IMPROVED] now performing the new grid2op `create_test_suite` 
 - [IMPROVED] now lightsim2grid properly throw `BackendError`
 

lightsim2grid/tests/test_SameResPP.py

@@ -219,8 +219,8 @@ def _aux_test(self, pn_net):
         Vdc = backend._grid.dc_pf(Vinit, max_iter, tol_this)
         backend._grid.reactivate_result_computation()
         backend._grid.tell_topo_changed()
-        Ydc_me = copy.deepcopy(backend._grid.get_Ybus())
-        Sdc_me = copy.deepcopy(backend._grid.get_Sbus())
+        Ydc_me = copy.deepcopy(backend._grid.get_dcYbus())
+        Sdc_me = copy.deepcopy(backend._grid.get_dcSbus())
         assert np.max(np.abs(V_init_ref[pp_vect_converter] - Vdc[:nb_sub])) <= 100.*self.tol,\
             f"\t Error for the DC approximation: resulting voltages are different " \
             f"{np.max(np.abs(V_init_ref[pp_vect_converter] - Vdc[:nb_sub])):.5f}pu"

lightsim2grid/tests/test_init_from_pypowsybl.py

@@ -125,12 +125,12 @@ def test_compare_pp(self):
         mat_pp = self.ref_samecase.get_dcYbus().todense()
         assert max_ <= self.tol_eq, f"error for dcYbus: {max_:.2e}"
         # check Sbus without slack
-        Sbus_ordered = self.gridmodel.get_Sbus()[reorder].reshape(-1)
+        Sbus_ordered = self.gridmodel.get_dcSbus()[reorder].reshape(-1)
         if slack_id > 0:
-            max_ = np.abs(Sbus_ordered[:slack_id] - self.ref_samecase.get_Sbus()[:slack_id]).max()
+            max_ = np.abs(Sbus_ordered[:slack_id] - self.ref_samecase.get_dcSbus()[:slack_id]).max()
             assert max_ <= self.tol_eq, f"error for dc Sbus: {max_:.2e}"
-        if slack_id != self.gridmodel.get_Sbus().shape[0] - 1:
-            max_ = np.abs(Sbus_ordered[(slack_id+1):] - self.ref_samecase.get_Sbus()[(slack_id+1):]).max()
+        if slack_id != self.gridmodel.get_dcSbus().shape[0] - 1:
+            max_ = np.abs(Sbus_ordered[(slack_id+1):] - self.ref_samecase.get_dcSbus()[(slack_id+1):]).max()
             assert max_ <= self.tol_eq, f"error for dc Sbus: {max_:.2e}"
 
         # same in AC

lightsim2grid/tests/test_ptdf.py

@@ -32,7 +32,10 @@ def setUp(self) -> None:
             warnings.filterwarnings("ignore")
             self.gridmodel = init(self.case)
         self.V_init = 1. * self.gridmodel.get_bus_vn_kv()
-        self.gridmodel.change_solver(self.get_solver_type())
+        solver_type = self.get_solver_type()
+        if solver_type not in self.gridmodel.available_solvers():
+            self.skipTest("Solver type not supported on this platform")
+        self.gridmodel.change_solver(solver_type)
         self.gridmodel.dc_pf(self.V_init, 1, 1e-8)
         self.dcYbus = 1.0 * self.gridmodel.get_dcYbus()
         self.dcSbus = 1.0 * self.gridmodel.get_dcSbus().real