ltbams 1.0.9__py3-none-any.whl → 1.0.10__py3-none-any.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- ams/__init__.py +0 -1
- ams/_version.py +3 -3
- ams/cases/5bus/pjm5bus_demo.json +1324 -0
- ams/core/__init__.py +1 -0
- ams/core/common.py +30 -0
- ams/core/model.py +1 -1
- ams/core/symprocessor.py +1 -1
- ams/extension/eva.py +1 -1
- ams/interface.py +40 -24
- ams/io/matpower.py +31 -17
- ams/io/psse.py +278 -1
- ams/main.py +2 -2
- ams/models/group.py +2 -1
- ams/models/static/pq.py +7 -3
- ams/opt/param.py +1 -2
- ams/routines/__init__.py +2 -3
- ams/routines/acopf.py +5 -108
- ams/routines/dcopf.py +8 -0
- ams/routines/dcpf.py +1 -1
- ams/routines/ed.py +4 -2
- ams/routines/grbopt.py +150 -0
- ams/routines/pflow.py +2 -2
- ams/routines/pypower.py +631 -0
- ams/routines/routine.py +4 -10
- ams/routines/uc.py +2 -2
- ams/shared.py +26 -43
- ams/system.py +118 -2
- docs/source/api.rst +2 -0
- docs/source/getting_started/install.rst +9 -6
- docs/source/images/dcopf_time.png +0 -0
- docs/source/images/educ_pie.png +0 -0
- docs/source/release-notes.rst +21 -47
- {ltbams-1.0.9.dist-info → ltbams-1.0.10.dist-info}/METADATA +87 -47
- {ltbams-1.0.9.dist-info → ltbams-1.0.10.dist-info}/RECORD +54 -71
- {ltbams-1.0.9.dist-info → ltbams-1.0.10.dist-info}/WHEEL +1 -1
- tests/test_1st_system.py +1 -1
- tests/test_case.py +14 -14
- tests/test_export_csv.py +1 -1
- tests/test_interface.py +24 -2
- tests/test_io.py +50 -0
- tests/test_omodel.py +1 -1
- tests/test_report.py +6 -6
- tests/test_routine.py +2 -2
- tests/test_rtn_acopf.py +75 -0
- tests/test_rtn_dcopf.py +1 -1
- tests/test_rtn_dcopf2.py +1 -1
- tests/test_rtn_ed.py +9 -9
- tests/test_rtn_opf.py +142 -0
- tests/test_rtn_pflow.py +0 -72
- tests/test_rtn_pypower.py +315 -0
- tests/test_rtn_rted.py +8 -8
- tests/test_rtn_uc.py +18 -18
- ams/pypower/__init__.py +0 -8
- ams/pypower/_compat.py +0 -9
- ams/pypower/core/__init__.py +0 -8
- ams/pypower/core/pips.py +0 -894
- ams/pypower/core/ppoption.py +0 -244
- ams/pypower/core/ppver.py +0 -18
- ams/pypower/core/solver.py +0 -2451
- ams/pypower/eps.py +0 -6
- ams/pypower/idx.py +0 -174
- ams/pypower/io.py +0 -604
- ams/pypower/make/__init__.py +0 -11
- ams/pypower/make/matrices.py +0 -665
- ams/pypower/make/pdv.py +0 -506
- ams/pypower/routines/__init__.py +0 -7
- ams/pypower/routines/cpf.py +0 -513
- ams/pypower/routines/cpf_callbacks.py +0 -114
- ams/pypower/routines/opf.py +0 -1803
- ams/pypower/routines/opffcns.py +0 -1946
- ams/pypower/routines/pflow.py +0 -852
- ams/pypower/toggle.py +0 -1098
- ams/pypower/utils.py +0 -293
- ams/routines/cpf.py +0 -65
- ams/routines/dcpf0.py +0 -196
- ams/routines/pflow0.py +0 -113
- tests/test_rtn_dcpf.py +0 -77
- {ltbams-1.0.9.dist-info → ltbams-1.0.10.dist-info}/entry_points.txt +0 -0
- {ltbams-1.0.9.dist-info → ltbams-1.0.10.dist-info}/top_level.txt +0 -0
tests/test_rtn_opf.py
ADDED
@@ -0,0 +1,142 @@
|
|
1
|
+
"""
|
2
|
+
Test wrapper routines of gurobi_optimodes.
|
3
|
+
"""
|
4
|
+
import unittest
|
5
|
+
|
6
|
+
import ams
|
7
|
+
from ams.shared import skip_unittest_without_gurobi_optimods
|
8
|
+
|
9
|
+
|
10
|
+
class TestOPF(unittest.TestCase):
|
11
|
+
"""
|
12
|
+
Test routine `OPF` in DC.
|
13
|
+
"""
|
14
|
+
|
15
|
+
def setUp(self) -> None:
|
16
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
17
|
+
setup=True, default_config=True, no_output=True)
|
18
|
+
# decrease load first
|
19
|
+
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
20
|
+
|
21
|
+
def test_init(self):
|
22
|
+
"""
|
23
|
+
Test initialization.
|
24
|
+
"""
|
25
|
+
self.ss.OPF.init()
|
26
|
+
self.assertTrue(self.ss.OPF.initialized, "OPF initialization failed!")
|
27
|
+
|
28
|
+
@skip_unittest_without_gurobi_optimods
|
29
|
+
def test_trip_gen(self):
|
30
|
+
"""
|
31
|
+
Test generator tripping.
|
32
|
+
"""
|
33
|
+
stg = 'PV_1'
|
34
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
35
|
+
|
36
|
+
self.ss.OPF.update()
|
37
|
+
self.ss.OPF.run(opftype='DC')
|
38
|
+
self.assertTrue(self.ss.OPF.converged, "OPF in DC did not converge under generator trip!")
|
39
|
+
self.assertAlmostEqual(self.ss.OPF.get(src='pg', attr='v', idx=stg),
|
40
|
+
0, places=6,
|
41
|
+
msg="Generator trip does not take effect!")
|
42
|
+
|
43
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1) # reset
|
44
|
+
|
45
|
+
@skip_unittest_without_gurobi_optimods
|
46
|
+
def test_trip_gen_ac(self):
|
47
|
+
"""
|
48
|
+
Test generator tripping.
|
49
|
+
"""
|
50
|
+
stg = 'PV_1'
|
51
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
52
|
+
|
53
|
+
self.ss.OPF.update()
|
54
|
+
self.ss.OPF.run(opftype='AC')
|
55
|
+
self.assertTrue(self.ss.OPF.converged, "OPF in AC did not converge under generator trip!")
|
56
|
+
self.assertAlmostEqual(self.ss.OPF.get(src='pg', attr='v', idx=stg),
|
57
|
+
0, places=6,
|
58
|
+
msg="Generator trip does not take effect!")
|
59
|
+
|
60
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1) # reset
|
61
|
+
|
62
|
+
@skip_unittest_without_gurobi_optimods
|
63
|
+
def test_trip_line(self):
|
64
|
+
"""
|
65
|
+
Test line tripping.
|
66
|
+
"""
|
67
|
+
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
68
|
+
|
69
|
+
self.ss.OPF.update()
|
70
|
+
self.ss.OPF.run(opftype='DC')
|
71
|
+
self.assertTrue(self.ss.OPF.converged, "OPF in DC did not converge under line trip!")
|
72
|
+
self.assertAlmostEqual(self.ss.OPF.get(src='plf', attr='v', idx='Line_3'),
|
73
|
+
0, places=6,
|
74
|
+
msg="Line trip does not take effect!")
|
75
|
+
|
76
|
+
self.ss.Line.alter(src='u', idx='Line_3', value=1) # reset
|
77
|
+
|
78
|
+
@skip_unittest_without_gurobi_optimods
|
79
|
+
def test_trip_line_ac(self):
|
80
|
+
"""
|
81
|
+
Test line tripping.
|
82
|
+
"""
|
83
|
+
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
84
|
+
|
85
|
+
self.ss.OPF.update()
|
86
|
+
self.ss.OPF.run(opftype='AC')
|
87
|
+
self.assertTrue(self.ss.OPF.converged, "OPF in AC did not converge under line trip!")
|
88
|
+
self.assertAlmostEqual(self.ss.OPF.get(src='plf', attr='v', idx='Line_3'),
|
89
|
+
0, places=6,
|
90
|
+
msg="Line trip does not take effect!")
|
91
|
+
|
92
|
+
self.ss.Line.alter(src='u', idx='Line_3', value=1) # reset
|
93
|
+
|
94
|
+
@skip_unittest_without_gurobi_optimods
|
95
|
+
def test_set_load(self):
|
96
|
+
"""
|
97
|
+
Test setting and tripping load.
|
98
|
+
"""
|
99
|
+
# --- run OPF ---
|
100
|
+
self.ss.OPF.run(opftype='DC')
|
101
|
+
pgs = self.ss.OPF.pg.v.sum()
|
102
|
+
|
103
|
+
# --- set load ---
|
104
|
+
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
|
105
|
+
self.ss.OPF.update()
|
106
|
+
|
107
|
+
self.ss.OPF.run(opftype='DC')
|
108
|
+
pgs_pqt = self.ss.OPF.pg.v.sum()
|
109
|
+
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
110
|
+
|
111
|
+
# --- trip load ---
|
112
|
+
self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
|
113
|
+
self.ss.OPF.update()
|
114
|
+
|
115
|
+
self.ss.OPF.run(opftype='DC')
|
116
|
+
pgs_pqt2 = self.ss.OPF.pg.v.sum()
|
117
|
+
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
118
|
+
|
119
|
+
@skip_unittest_without_gurobi_optimods
|
120
|
+
def test_set_load_ac(self):
|
121
|
+
"""
|
122
|
+
Test setting and tripping load.
|
123
|
+
"""
|
124
|
+
# --- run OPF ---
|
125
|
+
self.ss.OPF.run(opftype='AC')
|
126
|
+
pgs = self.ss.OPF.pg.v.sum()
|
127
|
+
|
128
|
+
# --- set load ---
|
129
|
+
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
|
130
|
+
self.ss.OPF.update()
|
131
|
+
|
132
|
+
self.ss.OPF.run(opftype='AC')
|
133
|
+
pgs_pqt = self.ss.OPF.pg.v.sum()
|
134
|
+
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
135
|
+
|
136
|
+
# --- trip load ---
|
137
|
+
self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
|
138
|
+
self.ss.OPF.update()
|
139
|
+
|
140
|
+
self.ss.OPF.run(opftype='AC')
|
141
|
+
pgs_pqt2 = self.ss.OPF.pg.v.sum()
|
142
|
+
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
tests/test_rtn_pflow.py
CHANGED
@@ -145,75 +145,3 @@ class TestPFlow(unittest.TestCase):
|
|
145
145
|
self.ss.PFlow.run()
|
146
146
|
pgs_pqt2 = self.ss.PFlow.pg.v.sum()
|
147
147
|
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
148
|
-
|
149
|
-
|
150
|
-
class TestACOPF(unittest.TestCase):
|
151
|
-
"""
|
152
|
-
Test routine `ACOPF`.
|
153
|
-
"""
|
154
|
-
|
155
|
-
def setUp(self) -> None:
|
156
|
-
self.ss = ams.load(ams.get_case('matpower/case14.m'),
|
157
|
-
setup=True, no_output=True, default_config=True)
|
158
|
-
|
159
|
-
def test_init(self):
|
160
|
-
"""
|
161
|
-
Test initialization.
|
162
|
-
"""
|
163
|
-
self.ss.ACOPF.init()
|
164
|
-
self.assertTrue(self.ss.ACOPF.initialized, "ACOPF initialization failed!")
|
165
|
-
|
166
|
-
def test_trip_gen(self):
|
167
|
-
"""
|
168
|
-
Test generator tripping.
|
169
|
-
"""
|
170
|
-
stg = 2
|
171
|
-
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
172
|
-
|
173
|
-
self.ss.ACOPF.update()
|
174
|
-
self.ss.ACOPF.run()
|
175
|
-
self.assertTrue(self.ss.ACOPF.converged, "ACOPF did not converge under generator trip!")
|
176
|
-
self.assertAlmostEqual(self.ss.ACOPF.get(src='pg', attr='v', idx=stg),
|
177
|
-
0, places=6,
|
178
|
-
msg="Generator trip does not take effect!")
|
179
|
-
|
180
|
-
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1)
|
181
|
-
|
182
|
-
def test_trip_line(self):
|
183
|
-
"""
|
184
|
-
Test line tripping.
|
185
|
-
"""
|
186
|
-
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
187
|
-
|
188
|
-
self.ss.ACOPF.update()
|
189
|
-
self.ss.ACOPF.run()
|
190
|
-
self.assertTrue(self.ss.ACOPF.converged, "ACOPF did not converge under line trip!")
|
191
|
-
self.assertAlmostEqual(self.ss.ACOPF.get(src='plf', attr='v', idx='Line_3'),
|
192
|
-
0, places=6,
|
193
|
-
msg="Line trip does not take effect!")
|
194
|
-
|
195
|
-
self.ss.Line.alter(src='u', idx='Line_3', value=1)
|
196
|
-
|
197
|
-
def test_set_load(self):
|
198
|
-
"""
|
199
|
-
Test setting and tripping load.
|
200
|
-
"""
|
201
|
-
# --- run ACOPF ---
|
202
|
-
self.ss.ACOPF.run()
|
203
|
-
pgs = self.ss.ACOPF.pg.v.sum()
|
204
|
-
|
205
|
-
# --- set load ---
|
206
|
-
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
|
207
|
-
self.ss.ACOPF.update()
|
208
|
-
|
209
|
-
self.ss.ACOPF.run()
|
210
|
-
pgs_pqt = self.ss.ACOPF.pg.v.sum()
|
211
|
-
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
212
|
-
|
213
|
-
# --- trip load ---
|
214
|
-
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0)
|
215
|
-
self.ss.ACOPF.update()
|
216
|
-
|
217
|
-
self.ss.ACOPF.run()
|
218
|
-
pgs_pqt2 = self.ss.ACOPF.pg.v.sum()
|
219
|
-
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
@@ -0,0 +1,315 @@
|
|
1
|
+
"""
|
2
|
+
Test wrapper routines of PYPOWER.
|
3
|
+
"""
|
4
|
+
import unittest
|
5
|
+
|
6
|
+
import ams
|
7
|
+
from ams.shared import skip_unittest_without_PYPOWER
|
8
|
+
|
9
|
+
|
10
|
+
class TestDCPF1(unittest.TestCase):
|
11
|
+
"""
|
12
|
+
Test routine `DCPF1`.
|
13
|
+
"""
|
14
|
+
|
15
|
+
def setUp(self) -> None:
|
16
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
17
|
+
setup=True, default_config=True, no_output=True)
|
18
|
+
# decrease load first
|
19
|
+
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
20
|
+
|
21
|
+
def test_init(self):
|
22
|
+
"""
|
23
|
+
Test initialization.
|
24
|
+
"""
|
25
|
+
self.ss.DCPF1.init()
|
26
|
+
self.assertTrue(self.ss.DCPF1.initialized, "DCPF1 initialization failed!")
|
27
|
+
|
28
|
+
@skip_unittest_without_PYPOWER
|
29
|
+
def test_trip_gen(self):
|
30
|
+
"""
|
31
|
+
Test generator tripping.
|
32
|
+
"""
|
33
|
+
stg = 'PV_1'
|
34
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
35
|
+
|
36
|
+
self.ss.DCPF1.update()
|
37
|
+
self.ss.DCPF1.run()
|
38
|
+
self.assertTrue(self.ss.DCPF1.converged, "DCPF1 did not converge under generator trip!")
|
39
|
+
self.assertAlmostEqual(self.ss.DCPF1.get(src='pg', attr='v', idx=stg),
|
40
|
+
0, places=6,
|
41
|
+
msg="Generator trip does not take effect!")
|
42
|
+
|
43
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1) # reset
|
44
|
+
|
45
|
+
@skip_unittest_without_PYPOWER
|
46
|
+
def test_trip_line(self):
|
47
|
+
"""
|
48
|
+
Test line tripping.
|
49
|
+
"""
|
50
|
+
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
51
|
+
|
52
|
+
self.ss.DCPF1.update()
|
53
|
+
self.ss.DCPF1.run()
|
54
|
+
self.assertTrue(self.ss.DCPF1.converged, "DCPF1 did not converge under line trip!")
|
55
|
+
self.assertAlmostEqual(self.ss.DCPF1.get(src='plf', attr='v', idx='Line_3'),
|
56
|
+
0, places=6,
|
57
|
+
msg="Line trip does not take effect!")
|
58
|
+
|
59
|
+
self.ss.Line.alter(src='u', idx='Line_3', value=1) # reset
|
60
|
+
|
61
|
+
@skip_unittest_without_PYPOWER
|
62
|
+
def test_set_load(self):
|
63
|
+
"""
|
64
|
+
Test setting and tripping load.
|
65
|
+
"""
|
66
|
+
# --- run DCPF1 ---
|
67
|
+
self.ss.DCPF1.run()
|
68
|
+
pgs = self.ss.DCPF1.pg.v.sum()
|
69
|
+
|
70
|
+
# --- set load ---
|
71
|
+
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
|
72
|
+
self.ss.DCPF1.update()
|
73
|
+
|
74
|
+
self.ss.DCPF1.run()
|
75
|
+
pgs_pqt = self.ss.DCPF1.pg.v.sum()
|
76
|
+
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
77
|
+
|
78
|
+
# --- trip load ---
|
79
|
+
self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
|
80
|
+
self.ss.DCPF1.update()
|
81
|
+
|
82
|
+
self.ss.DCPF1.run()
|
83
|
+
pgs_pqt2 = self.ss.DCPF1.pg.v.sum()
|
84
|
+
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
85
|
+
|
86
|
+
|
87
|
+
class TestPFlow1(unittest.TestCase):
|
88
|
+
"""
|
89
|
+
Test routine `PFlow1`.
|
90
|
+
"""
|
91
|
+
|
92
|
+
def setUp(self) -> None:
|
93
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
94
|
+
setup=True, default_config=True, no_output=True)
|
95
|
+
# decrease load first
|
96
|
+
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
97
|
+
|
98
|
+
def test_init(self):
|
99
|
+
"""
|
100
|
+
Test initialization.
|
101
|
+
"""
|
102
|
+
self.ss.PFlow1.init()
|
103
|
+
self.assertTrue(self.ss.PFlow1.initialized, "PFlow1 initialization failed!")
|
104
|
+
|
105
|
+
@skip_unittest_without_PYPOWER
|
106
|
+
def test_trip_gen(self):
|
107
|
+
"""
|
108
|
+
Test generator tripping.
|
109
|
+
"""
|
110
|
+
stg = 'PV_1'
|
111
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
112
|
+
|
113
|
+
self.ss.PFlow1.update()
|
114
|
+
self.ss.PFlow1.run()
|
115
|
+
self.assertTrue(self.ss.PFlow1.converged, "PFlow1 did not converge under generator trip!")
|
116
|
+
self.assertAlmostEqual(self.ss.PFlow1.get(src='pg', attr='v', idx=stg),
|
117
|
+
0, places=6,
|
118
|
+
msg="Generator trip does not take effect!")
|
119
|
+
|
120
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1) # reset
|
121
|
+
|
122
|
+
@skip_unittest_without_PYPOWER
|
123
|
+
def test_trip_line(self):
|
124
|
+
"""
|
125
|
+
Test line tripping.
|
126
|
+
"""
|
127
|
+
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
128
|
+
|
129
|
+
self.ss.PFlow1.update()
|
130
|
+
self.ss.PFlow1.run()
|
131
|
+
self.assertTrue(self.ss.PFlow1.converged, "PFlow1 did not converge under line trip!")
|
132
|
+
self.assertAlmostEqual(self.ss.PFlow1.get(src='plf', attr='v', idx='Line_3'),
|
133
|
+
0, places=6,
|
134
|
+
msg="Line trip does not take effect!")
|
135
|
+
|
136
|
+
self.ss.Line.alter(src='u', idx='Line_3', value=1) # reset
|
137
|
+
|
138
|
+
@skip_unittest_without_PYPOWER
|
139
|
+
def test_set_load(self):
|
140
|
+
"""
|
141
|
+
Test setting and tripping load.
|
142
|
+
"""
|
143
|
+
# --- run PFlow1 ---
|
144
|
+
self.ss.PFlow1.run()
|
145
|
+
pgs = self.ss.PFlow1.pg.v.sum()
|
146
|
+
|
147
|
+
# --- set load ---
|
148
|
+
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
|
149
|
+
self.ss.PFlow1.update()
|
150
|
+
|
151
|
+
self.ss.PFlow1.run()
|
152
|
+
pgs_pqt = self.ss.PFlow1.pg.v.sum()
|
153
|
+
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
154
|
+
|
155
|
+
# --- trip load ---
|
156
|
+
self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
|
157
|
+
self.ss.PFlow1.update()
|
158
|
+
|
159
|
+
self.ss.PFlow1.run()
|
160
|
+
pgs_pqt2 = self.ss.PFlow1.pg.v.sum()
|
161
|
+
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
162
|
+
|
163
|
+
|
164
|
+
class TestDCOPF1(unittest.TestCase):
|
165
|
+
"""
|
166
|
+
Test routine `DCOPF1`.
|
167
|
+
"""
|
168
|
+
|
169
|
+
def setUp(self) -> None:
|
170
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
171
|
+
setup=True, default_config=True, no_output=True)
|
172
|
+
# decrease load first
|
173
|
+
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
174
|
+
|
175
|
+
def test_init(self):
|
176
|
+
"""
|
177
|
+
Test initialization.
|
178
|
+
"""
|
179
|
+
self.ss.DCOPF1.init()
|
180
|
+
self.assertTrue(self.ss.DCOPF1.initialized, "DCOPF initialization failed!")
|
181
|
+
|
182
|
+
@skip_unittest_without_PYPOWER
|
183
|
+
def test_trip_gen(self):
|
184
|
+
"""
|
185
|
+
Test generator tripping.
|
186
|
+
"""
|
187
|
+
stg = 'PV_1'
|
188
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
189
|
+
|
190
|
+
self.ss.DCOPF1.update()
|
191
|
+
self.ss.DCOPF1.run()
|
192
|
+
self.assertTrue(self.ss.DCOPF1.converged, "DCOPF did not converge under generator trip!")
|
193
|
+
self.assertAlmostEqual(self.ss.DCOPF1.get(src='pg', attr='v', idx=stg),
|
194
|
+
0, places=6,
|
195
|
+
msg="Generator trip does not take effect!")
|
196
|
+
|
197
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1) # reset
|
198
|
+
|
199
|
+
@skip_unittest_without_PYPOWER
|
200
|
+
def test_trip_line(self):
|
201
|
+
"""
|
202
|
+
Test line tripping.
|
203
|
+
"""
|
204
|
+
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
205
|
+
|
206
|
+
self.ss.DCOPF1.update()
|
207
|
+
self.ss.DCOPF1.run()
|
208
|
+
self.assertTrue(self.ss.DCOPF1.converged, "DCOPF1 did not converge under line trip!")
|
209
|
+
self.assertAlmostEqual(self.ss.DCOPF1.get(src='plf', attr='v', idx='Line_3'),
|
210
|
+
0, places=6,
|
211
|
+
msg="Line trip does not take effect!")
|
212
|
+
|
213
|
+
self.ss.Line.alter(src='u', idx='Line_3', value=1) # reset
|
214
|
+
|
215
|
+
@skip_unittest_without_PYPOWER
|
216
|
+
def test_set_load(self):
|
217
|
+
"""
|
218
|
+
Test setting and tripping load.
|
219
|
+
"""
|
220
|
+
# --- run DCOPF1 ---
|
221
|
+
self.ss.DCOPF1.run()
|
222
|
+
pgs = self.ss.DCOPF1.pg.v.sum()
|
223
|
+
|
224
|
+
# --- set load ---
|
225
|
+
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
|
226
|
+
self.ss.DCOPF1.update()
|
227
|
+
|
228
|
+
self.ss.DCOPF1.run()
|
229
|
+
pgs_pqt = self.ss.DCOPF1.pg.v.sum()
|
230
|
+
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
231
|
+
|
232
|
+
# --- trip load ---
|
233
|
+
self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
|
234
|
+
self.ss.DCOPF1.update()
|
235
|
+
|
236
|
+
self.ss.DCOPF1.run()
|
237
|
+
pgs_pqt2 = self.ss.DCOPF1.pg.v.sum()
|
238
|
+
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
239
|
+
|
240
|
+
|
241
|
+
class TestACOPF1(unittest.TestCase):
|
242
|
+
"""
|
243
|
+
Test routine `ACOPF1`.
|
244
|
+
"""
|
245
|
+
|
246
|
+
def setUp(self) -> None:
|
247
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
248
|
+
setup=True, default_config=True, no_output=True)
|
249
|
+
# decrease load first
|
250
|
+
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
251
|
+
|
252
|
+
def test_init(self):
|
253
|
+
"""
|
254
|
+
Test initialization.
|
255
|
+
"""
|
256
|
+
self.ss.ACOPF1.init()
|
257
|
+
self.assertTrue(self.ss.ACOPF1.initialized, "ACOPF1 initialization failed!")
|
258
|
+
|
259
|
+
@skip_unittest_without_PYPOWER
|
260
|
+
def test_trip_gen(self):
|
261
|
+
"""
|
262
|
+
Test generator tripping.
|
263
|
+
"""
|
264
|
+
stg = 'PV_1'
|
265
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
266
|
+
|
267
|
+
self.ss.ACOPF1.update()
|
268
|
+
self.ss.ACOPF1.run()
|
269
|
+
self.assertTrue(self.ss.ACOPF1.converged, "ACOPF1 did not converge under generator trip!")
|
270
|
+
self.assertAlmostEqual(self.ss.ACOPF1.get(src='pg', attr='v', idx=stg),
|
271
|
+
0, places=6,
|
272
|
+
msg="Generator trip does not take effect!")
|
273
|
+
|
274
|
+
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=1) # reset
|
275
|
+
|
276
|
+
@skip_unittest_without_PYPOWER
|
277
|
+
def test_trip_line(self):
|
278
|
+
"""
|
279
|
+
Test line tripping.
|
280
|
+
"""
|
281
|
+
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
282
|
+
|
283
|
+
self.ss.ACOPF1.update()
|
284
|
+
self.ss.ACOPF1.run()
|
285
|
+
self.assertTrue(self.ss.ACOPF1.converged, "ACOPF1 did not converge under line trip!")
|
286
|
+
self.assertAlmostEqual(self.ss.ACOPF1.get(src='plf', attr='v', idx='Line_3'),
|
287
|
+
0, places=6,
|
288
|
+
msg="Line trip does not take effect!")
|
289
|
+
|
290
|
+
self.ss.Line.alter(src='u', idx='Line_3', value=1) # reset
|
291
|
+
|
292
|
+
@skip_unittest_without_PYPOWER
|
293
|
+
def test_set_load(self):
|
294
|
+
"""
|
295
|
+
Test setting and tripping load.
|
296
|
+
"""
|
297
|
+
# --- run ACOPF1 ---
|
298
|
+
self.ss.ACOPF1.run()
|
299
|
+
pgs = self.ss.ACOPF1.pg.v.sum()
|
300
|
+
|
301
|
+
# --- set load ---
|
302
|
+
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.1)
|
303
|
+
self.ss.ACOPF1.update()
|
304
|
+
|
305
|
+
self.ss.ACOPF1.run()
|
306
|
+
pgs_pqt = self.ss.ACOPF1.pg.v.sum()
|
307
|
+
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
308
|
+
|
309
|
+
# --- trip load ---
|
310
|
+
self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
|
311
|
+
self.ss.ACOPF1.update()
|
312
|
+
|
313
|
+
self.ss.ACOPF1.run()
|
314
|
+
pgs_pqt2 = self.ss.ACOPF1.pg.v.sum()
|
315
|
+
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|
tests/test_rtn_rted.py
CHANGED
@@ -11,7 +11,7 @@ class TestRTED(unittest.TestCase):
|
|
11
11
|
"""
|
12
12
|
|
13
13
|
def setUp(self) -> None:
|
14
|
-
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.
|
14
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
15
15
|
setup=True, default_config=True, no_output=True)
|
16
16
|
# decrease load first
|
17
17
|
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
@@ -107,7 +107,7 @@ class TestRTEDDG(unittest.TestCase):
|
|
107
107
|
"""
|
108
108
|
|
109
109
|
def setUp(self) -> None:
|
110
|
-
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.
|
110
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
111
111
|
setup=True, default_config=True, no_output=True)
|
112
112
|
# decrease load first
|
113
113
|
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
@@ -203,7 +203,7 @@ class TestRTEDES(unittest.TestCase):
|
|
203
203
|
"""
|
204
204
|
|
205
205
|
def setUp(self) -> None:
|
206
|
-
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.
|
206
|
+
self.ss = ams.load(ams.get_case("5bus/pjm5bus_demo.json"),
|
207
207
|
setup=True, default_config=True, no_output=True)
|
208
208
|
# decrease load first
|
209
209
|
self.ss.PQ.set(src='p0', attr='v', idx=['PQ_1', 'PQ_2'], value=[0.3, 0.3])
|
@@ -224,7 +224,7 @@ class TestRTEDES(unittest.TestCase):
|
|
224
224
|
self.ss.StaticGen.set(src='u', idx=stg, attr='v', value=0)
|
225
225
|
|
226
226
|
self.ss.RTEDES.update()
|
227
|
-
self.ss.RTEDES.run()
|
227
|
+
self.ss.RTEDES.run(solver='SCIP')
|
228
228
|
self.assertTrue(self.ss.RTEDES.converged, "RTEDES did not converge under generator trip!")
|
229
229
|
self.assertAlmostEqual(self.ss.RTEDES.get(src='pg', attr='v', idx=stg),
|
230
230
|
0, places=6,
|
@@ -240,7 +240,7 @@ class TestRTEDES(unittest.TestCase):
|
|
240
240
|
self.ss.Line.set(src='u', attr='v', idx='Line_3', value=0)
|
241
241
|
|
242
242
|
self.ss.RTEDES.update()
|
243
|
-
self.ss.RTEDES.run()
|
243
|
+
self.ss.RTEDES.run(solver='SCIP')
|
244
244
|
self.assertTrue(self.ss.RTEDES.converged, "RTEDES did not converge under line trip!")
|
245
245
|
self.assertAlmostEqual(self.ss.RTEDES.get(src='plf', attr='v', idx='Line_3'),
|
246
246
|
0, places=6,
|
@@ -253,14 +253,14 @@ class TestRTEDES(unittest.TestCase):
|
|
253
253
|
"""
|
254
254
|
Test setting and tripping load.
|
255
255
|
"""
|
256
|
-
self.ss.RTEDES.run()
|
256
|
+
self.ss.RTEDES.run(solver='SCIP')
|
257
257
|
pgs = self.ss.RTEDES.pg.v.sum()
|
258
258
|
|
259
259
|
# --- set load ---
|
260
260
|
self.ss.PQ.set(src='p0', attr='v', idx='PQ_1', value=0.05)
|
261
261
|
self.ss.RTEDES.update()
|
262
262
|
|
263
|
-
self.ss.RTEDES.run()
|
263
|
+
self.ss.RTEDES.run(solver='SCIP')
|
264
264
|
pgs_pqt = self.ss.RTEDES.pg.v.sum()
|
265
265
|
self.assertLess(pgs_pqt, pgs, "Load set does not take effect!")
|
266
266
|
|
@@ -268,6 +268,6 @@ class TestRTEDES(unittest.TestCase):
|
|
268
268
|
self.ss.PQ.set(src='u', attr='v', idx='PQ_2', value=0)
|
269
269
|
self.ss.RTEDES.update()
|
270
270
|
|
271
|
-
self.ss.RTEDES.run()
|
271
|
+
self.ss.RTEDES.run(solver='SCIP')
|
272
272
|
pgs_pqt2 = self.ss.RTEDES.pg.v.sum()
|
273
273
|
self.assertLess(pgs_pqt2, pgs_pqt, "Load trip does not take effect!")
|