ltbams 1.0.8__py3-none-any.whl → 1.0.10__py3-none-any.whl

tests/test_io.py

@@ -1,7 +1,7 @@
 """
 Test file input/output.
 """
-
+import os
 import unittest
 import numpy as np
 
@@ -9,6 +9,9 @@ import ams
 
 
 class TestMATPOWER(unittest.TestCase):
+    """
+    Test IO functions for MATPOWER and PYPOWER.
+    """
 
     def setUp(self):
         self.mpc5 = ams.io.matpower.m2mpc(ams.get_case('matpower/case5.m'))
@@ -34,6 +37,12 @@ class TestMATPOWER(unittest.TestCase):
         np.testing.assert_array_equal(system5.GCost.c2.v,
                                       np.zeros(system5.StaticGen.n))
 
+        # Check if Area is added
+        self.assertGreater(system5.Area.n, 0)
+
+        # Check if Zone is added
+        self.assertGreater(system5.Zone.n, 0)
+
         system14 = ams.system.System()
         # Test gentype length check
         mpc14 = self.mpc14.copy()
@@ -41,11 +50,41 @@ class TestMATPOWER(unittest.TestCase):
         with self.assertRaises(ValueError, msg='gentype length check failed!'):
             ams.io.matpower.mpc2system(mpc14, system14)
 
+        system14 = ams.system.System()
         ams.io.matpower.mpc2system(self.mpc14, system14)
         # In case14.m, the gencost has type 2 cost model, with 3 parameters.
         np.testing.assert_array_less(np.zeros(system14.StaticGen.n),
                                      system14.GCost.c2.v,)
 
+    def test_system2mpc(self):
+        """Test conversion from AMS System to MPC."""
+        system5 = ams.system.System()
+        ams.io.matpower.mpc2system(self.mpc5, system5)
+        mpc5 = ams.io.matpower.system2mpc(system5)
+
+        self.assertEqual(mpc5['baseMVA'], self.mpc5['baseMVA'])
+
+        # Bus
+        # type, PD, QD, GS,BS, VM, VA. BASE_KV, VMAX, VMIN
+        bus_cols = [1, 2, 3, 4, 5, 7, 8, 9, 11, 12]
+        np.testing.assert_array_equal(mpc5['bus'][:, bus_cols],
+                                      self.mpc5['bus'][:, bus_cols])
+
+        # Branch, Gen, Gencost, can have minor differences but is okay
+
+        # String type data
+        np.testing.assert_array_equal(mpc5['gentype'], self.mpc5['gentype'])
+        np.testing.assert_array_equal(mpc5['genfuel'], self.mpc5['genfuel'])
+        np.testing.assert_array_equal(mpc5['bus_name'], self.mpc5['bus_name'])
+
+        # Area quantity
+        self.assertEqual(np.unique(mpc5['bus'][:, 6]).shape[0],
+                         np.unique(self.mpc5['bus'][:, 6]).shape[0])
+
+        # Zone quantity
+        self.assertEqual(np.unique(mpc5['bus'][:, 10]).shape[0],
+                         np.unique(self.mpc5['bus'][:, 10]).shape[0])
+
     def test_gencost1(self):
         """Test when gencost is type 1."""
         mpcgc1 = self.mpc14.copy()
@@ -54,3 +93,88 @@ class TestMATPOWER(unittest.TestCase):
         system = ams.system.System()
         ams.io.matpower.mpc2system(mpcgc1, system)
         self.assertEqual(system.GCost.n, 5)
+
+    def test_mpc2m(self):
+        """Test conversion from MPC to M file."""
+        mpc5 = ams.io.matpower.m2mpc(ams.get_case('matpower/case5.m'))
+        mpc14 = ams.io.matpower.m2mpc(ams.get_case('matpower/case14.m'))
+
+        # Test conversion to M file
+        mfile5 = ams.io.matpower.mpc2m(mpc5, './case5out.m')
+        mfile14 = ams.io.matpower.mpc2m(mpc14, './case14out.m')
+
+        # Check if the files exist
+        self.assertTrue(os.path.exists(mfile5))
+        self.assertTrue(os.path.exists(mfile14))
+
+        mpc5read = ams.io.matpower.m2mpc(mfile5)
+        mpc14read = ams.io.matpower.m2mpc(mfile14)
+
+        # Check if the numerical values are the same
+        for key in mpc5:
+            if key in ['bus_name', 'gentype', 'genfuel']:
+                continue
+            np.testing.assert_array_almost_equal(
+                mpc5[key], mpc5read[key], decimal=5,
+                err_msg=f"Mismatch in {key} when converting case5.m"
+            )
+        for key in mpc14:
+            if key in ['bus_name', 'gentype', 'genfuel']:
+                continue
+            np.testing.assert_array_almost_equal(
+                mpc14[key], mpc14read[key], decimal=5,
+                err_msg=f"Mismatch in {key} when converting case14.m"
+            )
+
+        # Clean up the generated files
+        os.remove(mfile5)
+        os.remove(mfile14)
+
+    def test_system2m(self):
+        """Test conversion from AMS System to M file."""
+        s5 = ams.load(ams.get_case('matpower/case5.m'),
+                      no_output=True)
+        mfile5 = './case5out.m'
+        ams.io.matpower.write(s5, mfile5)
+
+        # Clean up the generated file
+        os.remove(mfile5)
+
+
+class TestSystemExport(unittest.TestCase):
+    """
+    Test system export functions.
+    """
+
+    def test_system_to(self):
+        """Test conversion from AMS System to several formats."""
+        ss = ams.load(ams.get_case('matpower/case5.m'),
+                      no_output=True)
+
+        mpc = ss.to_mpc()
+
+        self.assertIsInstance(mpc, dict)
+        self.assertIn('bus', mpc)
+        self.assertIn('branch', mpc)
+        self.assertIn('gen', mpc)
+        self.assertIn('gencost', mpc)
+        self.assertIn('baseMVA', mpc)
+        self.assertIn('bus_name', mpc)
+        self.assertIn('gentype', mpc)
+        self.assertIn('genfuel', mpc)
+
+        ss.to_m("./case5_out.m")
+        self.assertTrue(os.path.exists("./case5_out.m"))
+        os.remove("./case5_out.m")
+
+        ss.to_json("./case5_out.json")
+        self.assertTrue(os.path.exists("./case5_out.json"))
+        os.remove("./case5_out.json")
+
+        ss.to_xlsx("./case5_out.xlsx")
+        self.assertTrue(os.path.exists("./case5_out.xlsx"))
+        os.remove("./case5_out.xlsx")
+
+        ss.to_raw("./case5_out.raw", overwrite=True)
+        self.assertTrue(os.path.exists("./case5_out.raw"))
+        os.remove("./case5_out.raw")

tests/test_omodel.py

@@ -70,7 +70,7 @@ class TestOModelrepr(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True,
                            default_config=True,
                            no_output=True,

tests/test_report.py

@@ -20,7 +20,7 @@ class TestReport(unittest.TestCase):
 
     def setUp(self) -> None:
         self.ss = ams.main.load(
-            ams.get_case("5bus/pjm5bus_demo.xlsx"),
+            ams.get_case("5bus/pjm5bus_demo.json"),
             default_config=True,
             no_output=True,
         )
@@ -125,7 +125,7 @@ class TestReport(unittest.TestCase):
         Test report with RTEDES solved.
         """
         self.ss.files.no_output = False
-        self.ss.RTEDES.run()
+        self.ss.RTEDES.run(solver='SCIP')
         self.assertTrue(self.ss.report())
         self.assertTrue(os.path.exists(self.expected_report))
 
@@ -180,7 +180,7 @@ class TestReport(unittest.TestCase):
         Test report with EDES solved.
         """
         self.ss.files.no_output = False
-        self.ss.EDES.run()
+        self.ss.EDES.run(solver='SCIP')
         self.assertTrue(self.ss.report())
         self.assertTrue(os.path.exists(self.expected_report))
 
@@ -199,7 +199,7 @@ class TestReport(unittest.TestCase):
         Test report with UC solved.
         """
         self.ss.files.no_output = False
-        self.ss.UC.run()
+        self.ss.UC.run(solver='SCIP')
         self.assertTrue(self.ss.report())
         self.assertTrue(os.path.exists(self.expected_report))
 
@@ -218,7 +218,7 @@ class TestReport(unittest.TestCase):
         Test report with UCDG solved.
         """
         self.ss.files.no_output = False
-        self.ss.UCDG.run()
+        self.ss.UCDG.run(solver='SCIP')
         self.assertTrue(self.ss.report())
         self.assertTrue(os.path.exists(self.expected_report))
 
@@ -237,7 +237,7 @@ class TestReport(unittest.TestCase):
         Test report with UCES solved.
         """
         self.ss.files.no_output = False
-        self.ss.UCES.run()
+        self.ss.UCES.run(solver='SCIP')
         self.assertTrue(self.ss.report())
         self.assertTrue(os.path.exists(self.expected_report))
 

tests/test_routine.py

@@ -104,7 +104,7 @@ class TestSetOptzValueACOPF(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True,
                            default_config=True,
                            no_output=True,
@@ -144,7 +144,7 @@ class TestSetOptzValueDCOPF(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True,
                            default_config=True,
                            no_output=True,

tests/test_rtn_acopf.py

@@ -0,0 +1,75 @@
+import unittest
+
+import ams
+
+
+class TestACOPF(unittest.TestCase):
+    """
+    Test routine `ACOPF`.
+    """
+
+    def setUp(self) -> None:
+        self.ss = ams.load(ams.get_case('matpower/case14.m'),
+                           setup=True, no_output=True, default_config=True)
+
+    def test_init(self):
+        """
+        Test initialization.
+        """
+        self.ss.ACOPF.init()
+        self.assertTrue(self.ss.ACOPF.initialized, "ACOPF initialization failed!")
+
+    def test_trip_gen(self):
+        """
+        Test generator tripping.
+        """
+        stg = 2
+        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
+
+        self.ss.ACOPF.update()
+        self.ss.ACOPF.run()
+        self.assertTrue(self.ss.ACOPF.converged, "ACOPF did not converge under generator trip!")
+        self.assertAlmostEqual(self.ss.ACOPF.get(src='pg', attr='v', idx=stg),
+                               0, places=6,
+                               msg="Generator trip does not take effect!")
+
+        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1)
+
+    def test_trip_line(self):
+        """
+        Test line tripping.
+        """
+        self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
+
+        self.ss.ACOPF.update()
+        self.ss.ACOPF.run()
+        self.assertTrue(self.ss.ACOPF.converged, "ACOPF did not converge under line trip!")
+        self.assertAlmostEqual(self.ss.ACOPF.get(src='plf', attr='v', idx='Line_3'),
+                               0, places=6,
+                               msg="Line trip does not take effect!")
+
+        self.ss.Line.alter(src='u', idx='Line_3', value=1)
+
+    def test_set_load(self):
+        """
+        Test setting and tripping load.
+        """
+        # --- run ACOPF ---
+        self.ss.ACOPF.run()
+        pgs = self.ss.ACOPF.pg.v.sum()
+
+        # --- set load ---
+        self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
+        self.ss.ACOPF.update()
+
+        self.ss.ACOPF.run()
+        pgs_pqt = self.ss.ACOPF.pg.v.sum()
+        self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
+
+        # --- trip load ---
+        self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0)
+        self.ss.ACOPF.update()
+
+        self.ss.ACOPF.run()
+        pgs_pqt2 = self.ss.ACOPF.pg.v.sum()
+        self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")

tests/test_rtn_dcopf.py

@@ -10,7 +10,7 @@ class TestDCOPF(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True, default_config=True, no_output=True)
         # decrease load first
         self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])

tests/test_rtn_dcopf2.py

@@ -10,7 +10,7 @@ class TestDCOPF2(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True, default_config=True, no_output=True)
         # decrease load first
         self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])

tests/test_rtn_ed.py

@@ -11,7 +11,7 @@ class TestED(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True, default_config=True, no_output=True)
         # decrease load first
         self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
@@ -102,7 +102,7 @@ class TestEDDG(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True, default_config=True, no_output=True)
         # decrease load first
         self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
@@ -193,7 +193,7 @@ class TestEDES(unittest.TestCase):
     """
 
     def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.xlsx"),
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
                            setup=True, default_config=True, no_output=True)
         # decrease load first
         self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
@@ -217,7 +217,7 @@ class TestEDES(unittest.TestCase):
         loc_offtime = np.array([0, 2])
         self.ss.EDTSlot.ug.v[loc_offtime, stg_uid] = 0
 
-        self.ss.EDES.run()
+        self.ss.EDES.run(solver='SCIP')
         self.assertTrue(self.ss.EDES.converged, "ED did not converge under generator trip!")
         pg_pv1 = self.ss.EDES.get(src='pg', attr='v', idx=stg)
         np.testing.assert_almost_equal(np.zeros_like(loc_offtime),
@@ -232,7 +232,7 @@ class TestEDES(unittest.TestCase):
         self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
         self.ss.EDES.update()
 
-        self.ss.EDES.run()
+        self.ss.EDES.run(solver='SCIP')
         self.assertTrue(self.ss.EDES.converged, "ED did not converge under generator trip!")
         pg_pv1 = self.ss.EDES.get(src='pg', attr='v', idx=stg)
         np.testing.assert_array_less(np.zeros_like(pg_pv1), pg_pv1,
@@ -246,7 +246,7 @@ class TestEDES(unittest.TestCase):
         Test line tripping.
         """
         self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
-        self.ss.EDES.run()
+        self.ss.EDES.run(solver='SCIP')
         self.assertTrue(self.ss.EDES.converged, "EDES did not converge!")
         plf_l3 = self.ss.EDES.get(src='plf', attr='v', idx='Line_3')
         np.testing.assert_almost_equal(np.zeros_like(plf_l3),
@@ -259,14 +259,14 @@ class TestEDES(unittest.TestCase):
         """
         Test setting and tripping load.
         """
-        self.ss.EDES.run()
+        self.ss.EDES.run(solver='SCIP')
         pgs = self.ss.EDES.pg.v.sum()
 
         # --- set load ---
         self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
         self.ss.EDES.update()
 
-        self.ss.EDES.run()
+        self.ss.EDES.run(solver='SCIP')
         pgs_pqt = self.ss.EDES.pg.v.sum()
         self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
 
@@ -274,6 +274,6 @@ class TestEDES(unittest.TestCase):
         self.ss.PQ.alter(src='u', idx='PQ_2', value=0)
         self.ss.EDES.update()
 
-        self.ss.EDES.run()
+        self.ss.EDES.run(solver='SCIP')
         pgs_pqt2 = self.ss.EDES.pg.v.sum()
         self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")

tests/test_rtn_opf.py

@@ -0,0 +1,142 @@
+"""
+Test wrapper routines of gurobi_optimodes.
+"""
+import unittest
+
+import ams
+from ams.shared import skip_unittest_without_gurobi_optimods
+
+
+class TestOPF(unittest.TestCase):
+    """
+    Test routine `OPF` in DC.
+    """
+
+    def setUp(self) -> None:
+        self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
+                           setup=True, default_config=True, no_output=True)
+        # decrease load first
+        self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
+
+    def test_init(self):
+        """
+        Test initialization.
+        """
+        self.ss.OPF.init()
+        self.assertTrue(self.ss.OPF.initialized, "OPF initialization failed!")
+
+    @skip_unittest_without_gurobi_optimods
+    def test_trip_gen(self):
+        """
+        Test generator tripping.
+        """
+        stg = 'PV_1'
+        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
+
+        self.ss.OPF.update()
+        self.ss.OPF.run(opftype='DC')
+        self.assertTrue(self.ss.OPF.converged, "OPF in DC did not converge under generator trip!")
+        self.assertAlmostEqual(self.ss.OPF.get(src='pg', attr='v', idx=stg),
+                               0, places=6,
+                               msg="Generator trip does not take effect!")
+
+        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1)  # reset
+
+    @skip_unittest_without_gurobi_optimods
+    def test_trip_gen_ac(self):
+        """
+        Test generator tripping.
+        """
+        stg = 'PV_1'
+        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
+
+        self.ss.OPF.update()
+        self.ss.OPF.run(opftype='AC')
+        self.assertTrue(self.ss.OPF.converged, "OPF in AC did not converge under generator trip!")
+        self.assertAlmostEqual(self.ss.OPF.get(src='pg', attr='v', idx=stg),
+                               0, places=6,
+                               msg="Generator trip does not take effect!")
+
+        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1)  # reset
+
+    @skip_unittest_without_gurobi_optimods
+    def test_trip_line(self):
+        """
+        Test line tripping.
+        """
+        self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
+
+        self.ss.OPF.update()
+        self.ss.OPF.run(opftype='DC')
+        self.assertTrue(self.ss.OPF.converged, "OPF in DC did not converge under line trip!")
+        self.assertAlmostEqual(self.ss.OPF.get(src='plf', attr='v', idx='Line_3'),
+                               0, places=6,
+                               msg="Line trip does not take effect!")
+
+        self.ss.Line.alter(src='u', idx='Line_3', value=1)  # reset
+
+    @skip_unittest_without_gurobi_optimods
+    def test_trip_line_ac(self):
+        """
+        Test line tripping.
+        """
+        self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
+
+        self.ss.OPF.update()
+        self.ss.OPF.run(opftype='AC')
+        self.assertTrue(self.ss.OPF.converged, "OPF in AC did not converge under line trip!")
+        self.assertAlmostEqual(self.ss.OPF.get(src='plf', attr='v', idx='Line_3'),
+                               0, places=6,
+                               msg="Line trip does not take effect!")
+
+        self.ss.Line.alter(src='u', idx='Line_3', value=1)  # reset
+
+    @skip_unittest_without_gurobi_optimods
+    def test_set_load(self):
+        """
+        Test setting and tripping load.
+        """
+        # --- run OPF ---
+        self.ss.OPF.run(opftype='DC')
+        pgs = self.ss.OPF.pg.v.sum()
+
+        # --- set load ---
+        self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
+        self.ss.OPF.update()
+
+        self.ss.OPF.run(opftype='DC')
+        pgs_pqt = self.ss.OPF.pg.v.sum()
+        self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
+
+        # --- trip load ---
+        self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
+        self.ss.OPF.update()
+
+        self.ss.OPF.run(opftype='DC')
+        pgs_pqt2 = self.ss.OPF.pg.v.sum()
+        self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
+
+    @skip_unittest_without_gurobi_optimods
+    def test_set_load_ac(self):
+        """
+        Test setting and tripping load.
+        """
+        # --- run OPF ---
+        self.ss.OPF.run(opftype='AC')
+        pgs = self.ss.OPF.pg.v.sum()
+
+        # --- set load ---
+        self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
+        self.ss.OPF.update()
+
+        self.ss.OPF.run(opftype='AC')
+        pgs_pqt = self.ss.OPF.pg.v.sum()
+        self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
+
+        # --- trip load ---
+        self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
+        self.ss.OPF.update()
+
+        self.ss.OPF.run(opftype='AC')
+        pgs_pqt2 = self.ss.OPF.pg.v.sum()
+        self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")

tests/test_rtn_pflow.py

@@ -145,75 +145,3 @@ class TestPFlow(unittest.TestCase):
         self.ss.PFlow.run()
         pgs_pqt2 = self.ss.PFlow.pg.v.sum()
         self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
-
-
-class TestACOPF(unittest.TestCase):
-    """
-    Test routine `ACOPF`.
-    """
-
-    def setUp(self) -> None:
-        self.ss = ams.load(ams.get_case('matpower/case14.m'),
-                           setup=True, no_output=True, default_config=True)
-
-    def test_init(self):
-        """
-        Test initialization.
-        """
-        self.ss.ACOPF.init()
-        self.assertTrue(self.ss.ACOPF.initialized, "ACOPF initialization failed!")
-
-    def test_trip_gen(self):
-        """
-        Test generator tripping.
-        """
-        stg = 2
-        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
-
-        self.ss.ACOPF.update()
-        self.ss.ACOPF.run()
-        self.assertTrue(self.ss.ACOPF.converged, "ACOPF did not converge under generator trip!")
-        self.assertAlmostEqual(self.ss.ACOPF.get(src='pg', attr='v', idx=stg),
-                               0, places=6,
-                               msg="Generator trip does not take effect!")
-
-        self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1)
-
-    def test_trip_line(self):
-        """
-        Test line tripping.
-        """
-        self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
-
-        self.ss.ACOPF.update()
-        self.ss.ACOPF.run()
-        self.assertTrue(self.ss.ACOPF.converged, "ACOPF did not converge under line trip!")
-        self.assertAlmostEqual(self.ss.ACOPF.get(src='plf', attr='v', idx='Line_3'),
-                               0, places=6,
-                               msg="Line trip does not take effect!")
-
-        self.ss.Line.alter(src='u', idx='Line_3', value=1)
-
-    def test_set_load(self):
-        """
-        Test setting and tripping load.
-        """
-        # --- run ACOPF ---
-        self.ss.ACOPF.run()
-        pgs = self.ss.ACOPF.pg.v.sum()
-
-        # --- set load ---
-        self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
-        self.ss.ACOPF.update()
-
-        self.ss.ACOPF.run()
-        pgs_pqt = self.ss.ACOPF.pg.v.sum()
-        self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
-
-        # --- trip load ---
-        self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0)
-        self.ss.ACOPF.update()
-
-        self.ss.ACOPF.run()
-        pgs_pqt2 = self.ss.ACOPF.pg.v.sum()
-        self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")