0001 function t_loadcase(quiet)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033 if nargin < 1
0034 quiet = 0;
0035 end
0036
0037 t_begin(234, quiet);
0038
0039
0040 casefile = 't_case9_opf';
0041 matfile = 't_mat9_opf';
0042 pfcasefile = 't_case9_pf';
0043 pfmatfile = 't_mat9_pf';
0044 casefilev2 = 't_case9_opfv2';
0045 matfilev2 = 't_mat9_opfv2';
0046 pfcasefilev2 = 't_case9_pfv2';
0047 pfmatfilev2 = 't_mat9_pfv2';
0048
0049 [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
0050 MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
0051 QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
0052 [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, ...
0053 TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ...
0054 ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
0055
0056
0057 [baseMVA, bus, gen, branch, areas, gencost] = feval(casefile);
0058
0059 eval(['save ' matfile '.mat baseMVA bus gen branch areas gencost']);
0060
0061
0062 mpc = feval(casefilev2);
0063
0064 eval(['save ' matfilev2 '.mat mpc']);
0065
0066
0067
0068 tmp1 = {mpc.baseMVA, mpc.bus, mpc.gen, mpc.branch, [1 5], mpc.gencost};
0069 tmp2 = {mpc.baseMVA, mpc.bus, mpc.gen, mpc.branch, [1 5], mpc.gencost};
0070
0071 tmp1{3}(2:3, [PC1, PC2, QC1MIN, QC1MAX, QC2MIN, QC2MAX]) = zeros(2,6);
0072 tmp1{4}(1, ANGMAX) = 360;
0073 tmp1{4}(9, ANGMIN) = -360;
0074 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp1{:});
0075
0076
0077 t = 'loadcase(opf_M_file_v1) without .m extension : ';
0078 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase(casefile);
0079 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0080 t_is(bus1, bus, 12, [t 'bus']);
0081 t_is(gen1, gen, 12, [t 'gen']);
0082 t_is(branch1, branch, 12, [t 'branch']);
0083 t_is(areas1, areas, 12, [t 'areas']);
0084 t_is(gencost1, gencost, 12, [t 'gencost']);
0085
0086 t = 'loadcase(opf_M_file_v1) with .m extension : ';
0087 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase([casefile '.m']);
0088 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0089 t_is(bus1, bus, 12, [t 'bus']);
0090 t_is(gen1, gen, 12, [t 'gen']);
0091 t_is(branch1, branch, 12, [t 'branch']);
0092 t_is(areas1, areas, 12, [t 'areas']);
0093 t_is(gencost1, gencost, 12, [t 'gencost']);
0094
0095 t = 'loadcase(opf_MAT_file_v1) without .mat extension : ';
0096 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase(matfile);
0097 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0098 t_is(bus1, bus, 12, [t 'bus']);
0099 t_is(gen1, gen, 12, [t 'gen']);
0100 t_is(branch1, branch, 12, [t 'branch']);
0101 t_is(areas1, areas, 12, [t 'areas']);
0102 t_is(gencost1, gencost, 12, [t 'gencost']);
0103
0104 t = 'loadcase(opf_MAT_file_v1) with .mat extension : ';
0105 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase([matfile '.mat']);
0106 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0107 t_is(bus1, bus, 12, [t 'bus']);
0108 t_is(gen1, gen, 12, [t 'gen']);
0109 t_is(branch1, branch, 12, [t 'branch']);
0110 t_is(areas1, areas, 12, [t 'areas']);
0111 t_is(gencost1, gencost, 12, [t 'gencost']);
0112
0113
0114 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp2{:});
0115 areas = [];
0116
0117 t = 'loadcase(opf_M_file_v2) without .m extension : ';
0118 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase(casefilev2);
0119 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0120 t_is(bus1, bus, 12, [t 'bus']);
0121 t_is(gen1, gen, 12, [t 'gen']);
0122 t_is(branch1, branch, 12, [t 'branch']);
0123 t_is(areas1, areas, 12, [t 'areas']);
0124 t_is(gencost1, gencost, 12, [t 'gencost']);
0125
0126 t = 'loadcase(opf_M_file_v2) with .m extension : ';
0127 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase([casefilev2 '.m']);
0128 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0129 t_is(bus1, bus, 12, [t 'bus']);
0130 t_is(gen1, gen, 12, [t 'gen']);
0131 t_is(branch1, branch, 12, [t 'branch']);
0132 t_is(areas1, areas, 12, [t 'areas']);
0133 t_is(gencost1, gencost, 12, [t 'gencost']);
0134
0135 t = 'loadcase(opf_MAT_file_v2) without .mat extension : ';
0136 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase(matfilev2);
0137 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0138 t_is(bus1, bus, 12, [t 'bus']);
0139 t_is(gen1, gen, 12, [t 'gen']);
0140 t_is(branch1, branch, 12, [t 'branch']);
0141 t_is(areas1, areas, 12, [t 'areas']);
0142 t_is(gencost1, gencost, 12, [t 'gencost']);
0143
0144 t = 'loadcase(opf_MAT_file_v2) with .mat extension : ';
0145 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = loadcase([matfilev2 '.mat']);
0146 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0147 t_is(bus1, bus, 12, [t 'bus']);
0148 t_is(gen1, gen, 12, [t 'gen']);
0149 t_is(branch1, branch, 12, [t 'branch']);
0150 t_is(areas1, areas, 12, [t 'areas']);
0151 t_is(gencost1, gencost, 12, [t 'gencost']);
0152
0153
0154 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp1{:});
0155
0156 t = 'loadcase(opf_struct_v1) (no version): ';
0157 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = feval(casefile);
0158 clear c;
0159 c.baseMVA = baseMVA1;
0160 c.bus = bus1;
0161 c.gen = gen1;
0162 c.branch = branch1;
0163 c.areas = areas1;
0164 c.gencost = gencost1;
0165 [baseMVA2, bus2, gen2, branch2, areas2, gencost2] = loadcase(c);
0166 t_is(baseMVA2, baseMVA, 12, [t 'baseMVA']);
0167 t_is(bus2, bus, 12, [t 'bus']);
0168 t_is(gen2, gen, 12, [t 'gen']);
0169 t_is(branch2, branch, 12, [t 'branch']);
0170 t_is(areas2, areas, 12, [t 'areas']);
0171 t_is(gencost2, gencost, 12, [t 'gencost']);
0172
0173 t = 'loadcase(opf_struct_v1) (version=''1''): ';
0174 c.version = '1';
0175 [baseMVA2, bus2, gen2, branch2, areas2, gencost2] = loadcase(c);
0176 t_is(baseMVA2, baseMVA, 12, [t 'baseMVA']);
0177 t_is(bus2, bus, 12, [t 'bus']);
0178 t_is(gen2, gen, 12, [t 'gen']);
0179 t_is(branch2, branch, 12, [t 'branch']);
0180 t_is(areas2, areas, 12, [t 'areas']);
0181 t_is(gencost2, gencost, 12, [t 'gencost']);
0182
0183
0184 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp2{:});
0185
0186 t = 'loadcase(opf_struct_v2) (no version): ';
0187 clear c;
0188 c.baseMVA = baseMVA;
0189 c.bus = bus;
0190 c.gen = gen;
0191 c.branch = branch;
0192 c.areas = areas;
0193 c.gencost = gencost;
0194 [baseMVA2, bus2, gen2, branch2, areas2, gencost2] = loadcase(c);
0195 t_is(baseMVA2, baseMVA, 12, [t 'baseMVA']);
0196 t_is(bus2, bus, 12, [t 'bus']);
0197 t_is(gen2, gen, 12, [t 'gen']);
0198 t_is(branch2, branch, 12, [t 'branch']);
0199 t_is(areas2, areas, 12, [t 'areas']);
0200 t_is(gencost2, gencost, 12, [t 'gencost']);
0201
0202 t = 'loadcase(opf_struct_v2) (version=''2''): ';
0203 clear c;
0204 c.baseMVA = baseMVA;
0205 c.bus = bus;
0206 c.gen = gen;
0207 c.branch = branch;
0208 c.areas = areas;
0209 c.gencost = gencost;
0210 c.version = '2';
0211 [baseMVA2, bus2, gen2, branch2, areas2, gencost2] = loadcase(c);
0212 t_is(baseMVA2, baseMVA, 12, [t 'baseMVA']);
0213 t_is(bus2, bus, 12, [t 'bus']);
0214 t_is(gen2, gen, 12, [t 'gen']);
0215 t_is(branch2, branch, 12, [t 'branch']);
0216 t_is(areas2, areas, 12, [t 'areas']);
0217 t_is(gencost2, gencost, 12, [t 'gencost']);
0218
0219
0220
0221 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp1{:});
0222
0223 t = 'mpc = loadcase(opf_M_file_v1) without .m extension : ';
0224 mpc1 = loadcase(casefile);
0225 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0226 t_is(mpc1.bus, bus, 12, [t 'bus']);
0227 t_is(mpc1.gen, gen, 12, [t 'gen']);
0228 t_is(mpc1.branch, branch, 12, [t 'branch']);
0229 t_is(mpc1.areas, areas, 12, [t 'areas']);
0230 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0231
0232 t = 'mpc = loadcase(opf_M_file_v1) with .m extension : ';
0233 mpc1 = loadcase([casefile '.m']);
0234 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0235 t_is(mpc1.bus, bus, 12, [t 'bus']);
0236 t_is(mpc1.gen, gen, 12, [t 'gen']);
0237 t_is(mpc1.branch, branch, 12, [t 'branch']);
0238 t_is(mpc1.areas, areas, 12, [t 'areas']);
0239 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0240
0241 t = 'mpc = loadcase(opf_MAT_file_v1) without .mat extension : ';
0242 mpc1 = loadcase(matfile);
0243 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0244 t_is(mpc1.bus, bus, 12, [t 'bus']);
0245 t_is(mpc1.gen, gen, 12, [t 'gen']);
0246 t_is(mpc1.branch, branch, 12, [t 'branch']);
0247 t_is(mpc1.areas, areas, 12, [t 'areas']);
0248 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0249
0250 t = 'mpc = loadcase(opf_MAT_file_v1) with .mat extension : ';
0251 mpc1 = loadcase([matfile '.mat']);
0252 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0253 t_is(mpc1.bus, bus, 12, [t 'bus']);
0254 t_is(mpc1.gen, gen, 12, [t 'gen']);
0255 t_is(mpc1.branch, branch, 12, [t 'branch']);
0256 t_is(mpc1.areas, areas, 12, [t 'areas']);
0257 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0258
0259
0260 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp2{:});
0261
0262 t = 'mpc = loadcase(opf_M_file_v2) without .m extension : ';
0263 mpc1 = loadcase(casefilev2);
0264 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0265 t_is(mpc1.bus, bus, 12, [t 'bus']);
0266 t_is(mpc1.gen, gen, 12, [t 'gen']);
0267 t_is(mpc1.branch, branch, 12, [t 'branch']);
0268 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0269
0270 t = 'mpc = loadcase(opf_M_file_v2) with .m extension : ';
0271 mpc1 = loadcase([casefilev2 '.m']);
0272 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0273 t_is(mpc1.bus, bus, 12, [t 'bus']);
0274 t_is(mpc1.gen, gen, 12, [t 'gen']);
0275 t_is(mpc1.branch, branch, 12, [t 'branch']);
0276 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0277
0278 t = 'mpc = loadcase(opf_MAT_file_v2) without .mat extension : ';
0279 mpc1 = loadcase(matfilev2);
0280 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0281 t_is(mpc1.bus, bus, 12, [t 'bus']);
0282 t_is(mpc1.gen, gen, 12, [t 'gen']);
0283 t_is(mpc1.branch, branch, 12, [t 'branch']);
0284 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0285
0286 t = 'mpc = loadcase(opf_MAT_file_v2) with .mat extension : ';
0287 mpc1 = loadcase([matfilev2 '.mat']);
0288 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0289 t_is(mpc1.bus, bus, 12, [t 'bus']);
0290 t_is(mpc1.gen, gen, 12, [t 'gen']);
0291 t_is(mpc1.branch, branch, 12, [t 'branch']);
0292 t_is(mpc1.gencost, gencost, 12, [t 'gencost']);
0293
0294
0295 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp1{:});
0296
0297 t = 'mpc = loadcase(opf_struct_v1) (no version): ';
0298 [baseMVA1, bus1, gen1, branch1, areas1, gencost1] = feval(casefile);
0299 clear c;
0300 c.baseMVA = baseMVA1;
0301 c.bus = bus1;
0302 c.gen = gen1;
0303 c.branch = branch1;
0304 c.areas = areas1;
0305 c.gencost = gencost1;
0306 mpc2 = loadcase(c);
0307 t_is(mpc2.baseMVA, baseMVA, 12, [t 'baseMVA']);
0308 t_is(mpc2.bus, bus, 12, [t 'bus']);
0309 t_is(mpc2.gen, gen, 12, [t 'gen']);
0310 t_is(mpc2.branch, branch, 12, [t 'branch']);
0311 t_is(mpc2.areas, areas, 12, [t 'areas']);
0312 t_is(mpc2.gencost, gencost, 12, [t 'gencost']);
0313
0314 t = 'mpc = loadcase(opf_struct_v1) (version=''1''): ';
0315 c.version = '1';
0316 mpc2 = loadcase(c);
0317 t_is(mpc2.baseMVA, baseMVA, 12, [t 'baseMVA']);
0318 t_is(mpc2.bus, bus, 12, [t 'bus']);
0319 t_is(mpc2.gen, gen, 12, [t 'gen']);
0320 t_is(mpc2.branch, branch, 12, [t 'branch']);
0321 t_is(mpc2.areas, areas, 12, [t 'areas']);
0322 t_is(mpc2.gencost, gencost, 12, [t 'gencost']);
0323
0324
0325 [baseMVA, bus, gen, branch, areas, gencost] = deal(tmp2{:});
0326
0327 t = 'mpc = loadcase(opf_struct_v2) (no version): ';
0328 clear c;
0329 c.baseMVA = baseMVA;
0330 c.bus = bus;
0331 c.gen = gen;
0332 c.branch = branch;
0333 c.gencost = gencost;
0334 mpc2 = loadcase(c);
0335 t_is(mpc2.baseMVA, baseMVA, 12, [t 'baseMVA']);
0336 t_is(mpc2.bus, bus, 12, [t 'bus']);
0337 t_is(mpc2.gen, gen, 12, [t 'gen']);
0338 t_is(mpc2.branch, branch, 12, [t 'branch']);
0339 t_is(mpc2.gencost, gencost, 12, [t 'gencost']);
0340 t_ok(strcmp(mpc2.version, '2'), [t 'version']);
0341
0342 t = 'mpc = loadcase(opf_struct_v2) (version=''2''): ';
0343 clear c;
0344 c.baseMVA = baseMVA;
0345 c.bus = bus;
0346 c.gen = gen;
0347 c.branch = branch;
0348 c.gencost = gencost;
0349 c.version = '2';
0350 mpc2 = loadcase(c);
0351 t_is(mpc2.baseMVA, baseMVA, 12, [t 'baseMVA']);
0352 t_is(mpc2.bus, bus, 12, [t 'bus']);
0353 t_is(mpc2.gen, gen, 12, [t 'gen']);
0354 t_is(mpc2.branch, branch, 12, [t 'branch']);
0355 t_is(mpc2.gencost, gencost, 12, [t 'gencost']);
0356
0357
0358
0359 [baseMVA, bus, gen, branch] = feval(pfcasefile);
0360 eval(['save ' pfmatfile '.mat baseMVA bus gen branch']);
0361
0362
0363 mpc = feval(pfcasefilev2);
0364 tmp = {mpc.baseMVA, mpc.bus, mpc.gen, mpc.branch};
0365 [baseMVA, bus, gen, branch] = deal(tmp{:});
0366
0367 eval(['save ' pfmatfilev2 '.mat mpc']);
0368
0369
0370 t = 'loadcase(pf_M_file_v1) without .m extension : ';
0371 [baseMVA1, bus1, gen1, branch1] = loadcase(pfcasefile);
0372 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0373 t_is(bus1, bus, 12, [t 'bus']);
0374 t_is(gen1, gen, 12, [t 'gen']);
0375 t_is(branch1, branch, 12, [t 'branch']);
0376
0377 t = 'loadcase(pf_M_file_v1) with .m extension : ';
0378 [baseMVA1, bus1, gen1, branch1] = loadcase([pfcasefile '.m']);
0379 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0380 t_is(bus1, bus, 12, [t 'bus']);
0381 t_is(gen1, gen, 12, [t 'gen']);
0382 t_is(branch1, branch, 12, [t 'branch']);
0383
0384 t = 'loadcase(pf_MAT_file_v1) without .mat extension : ';
0385 [baseMVA1, bus1, gen1, branch1] = loadcase(pfmatfile);
0386 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0387 t_is(bus1, bus, 12, [t 'bus']);
0388 t_is(gen1, gen, 12, [t 'gen']);
0389 t_is(branch1, branch, 12, [t 'branch']);
0390
0391 t = 'loadcase(pf_MAT_file_v1) with .mat extension : ';
0392 [baseMVA1, bus1, gen1, branch1] = loadcase([pfmatfile '.mat']);
0393 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0394 t_is(bus1, bus, 12, [t 'bus']);
0395 t_is(gen1, gen, 12, [t 'gen']);
0396 t_is(branch1, branch, 12, [t 'branch']);
0397
0398 t = 'loadcase(pf_M_file_v2) without .m extension : ';
0399 [baseMVA1, bus1, gen1, branch1] = loadcase(pfcasefilev2);
0400 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0401 t_is(bus1, bus, 12, [t 'bus']);
0402 t_is(gen1, gen, 12, [t 'gen']);
0403 t_is(branch1, branch, 12, [t 'branch']);
0404
0405 t = 'loadcase(pf_M_file_v2) with .m extension : ';
0406 [baseMVA1, bus1, gen1, branch1] = loadcase([pfcasefilev2 '.m']);
0407 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0408 t_is(bus1, bus, 12, [t 'bus']);
0409 t_is(gen1, gen, 12, [t 'gen']);
0410 t_is(branch1, branch, 12, [t 'branch']);
0411
0412 t = 'loadcase(pf_MAT_file_v2) without .mat extension : ';
0413 [baseMVA1, bus1, gen1, branch1] = loadcase(pfmatfilev2);
0414 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0415 t_is(bus1, bus, 12, [t 'bus']);
0416 t_is(gen1, gen, 12, [t 'gen']);
0417 t_is(branch1, branch, 12, [t 'branch']);
0418
0419 t = 'loadcase(pf_MAT_file_v2) with .mat extension : ';
0420 [baseMVA1, bus1, gen1, branch1] = loadcase([pfmatfilev2 '.mat']);
0421 t_is(baseMVA1, baseMVA, 12, [t 'baseMVA']);
0422 t_is(bus1, bus, 12, [t 'bus']);
0423 t_is(gen1, gen, 12, [t 'gen']);
0424 t_is(branch1, branch, 12, [t 'branch']);
0425
0426 t = 'loadcase(pf_struct_v1) (no version): ';
0427 [baseMVA1, bus1, gen1, branch1] = feval(pfcasefile);
0428 clear c;
0429 c.baseMVA = baseMVA1;
0430 c.bus = bus1;
0431 c.gen = gen1;
0432 c.branch = branch1;
0433 [baseMVA2, bus2, gen2, branch2] = loadcase(c);
0434 t_is(baseMVA2, baseMVA, 12, [t 'baseMVA']);
0435 t_is(bus2, bus, 12, [t 'bus']);
0436 t_is(gen2, gen, 12, [t 'gen']);
0437 t_is(branch2, branch, 12, [t 'branch']);
0438
0439 t = 'loadcase(pf_struct_v1) (version=''1''): ';
0440 c.version = '1';
0441 [baseMVA2, bus2, gen2, branch2] = loadcase(c);
0442 t_is(baseMVA2, baseMVA, 12, [t 'baseMVA']);
0443 t_is(bus2, bus, 12, [t 'bus']);
0444 t_is(gen2, gen, 12, [t 'gen']);
0445 t_is(branch2, branch, 12, [t 'branch']);
0446
0447 t = 'loadcase(pf_struct_v2) : ';
0448 clear c;
0449 c.baseMVA = baseMVA;
0450 c.bus = bus;
0451 c.gen = gen;
0452 c.branch = branch;
0453 c.version = '2';
0454 [baseMVA2, bus2, gen2, branch2] = loadcase(c);
0455 t_is(baseMVA2, baseMVA, 12, [t 'baseMVA']);
0456 t_is(bus2, bus, 12, [t 'bus']);
0457 t_is(gen2, gen, 12, [t 'gen']);
0458 t_is(branch2, branch, 12, [t 'branch']);
0459
0460
0461
0462 t = 'mpc = loadcase(pf_M_file_v1) without .m extension : ';
0463 mpc1 = loadcase(pfcasefile);
0464 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0465 t_is(mpc1.bus, bus, 12, [t 'bus']);
0466 t_is(mpc1.gen, gen, 12, [t 'gen']);
0467 t_is(mpc1.branch, branch, 12, [t 'branch']);
0468
0469 t = 'mpc = loadcase(pf_M_file_v1) with .m extension : ';
0470 mpc1 = loadcase([pfcasefile '.m']);
0471 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0472 t_is(mpc1.bus, bus, 12, [t 'bus']);
0473 t_is(mpc1.gen, gen, 12, [t 'gen']);
0474 t_is(mpc1.branch, branch, 12, [t 'branch']);
0475
0476 t = 'mpc = loadcase(pf_MAT_file_v1) without .mat extension : ';
0477 mpc1 = loadcase(pfmatfile);
0478 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0479 t_is(mpc1.bus, bus, 12, [t 'bus']);
0480 t_is(mpc1.gen, gen, 12, [t 'gen']);
0481 t_is(mpc1.branch, branch, 12, [t 'branch']);
0482
0483 t = 'mpc = loadcase(pf_MAT_file_v1) with .mat extension : ';
0484 mpc1 = loadcase([pfmatfile '.mat']);
0485 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0486 t_is(mpc1.bus, bus, 12, [t 'bus']);
0487 t_is(mpc1.gen, gen, 12, [t 'gen']);
0488 t_is(mpc1.branch, branch, 12, [t 'branch']);
0489
0490 t = 'mpc = loadcase(pf_M_file_v2) without .m extension : ';
0491 mpc1 = loadcase(pfcasefilev2);
0492 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0493 t_is(mpc1.bus, bus, 12, [t 'bus']);
0494 t_is(mpc1.gen, gen, 12, [t 'gen']);
0495 t_is(mpc1.branch, branch, 12, [t 'branch']);
0496
0497 t = 'mpc = loadcase(pf_M_file_v2) with .m extension : ';
0498 mpc1 = loadcase([pfcasefilev2 '.m']);
0499 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0500 t_is(mpc1.bus, bus, 12, [t 'bus']);
0501 t_is(mpc1.gen, gen, 12, [t 'gen']);
0502 t_is(mpc1.branch, branch, 12, [t 'branch']);
0503
0504 t = 'mpc = loadcase(pf_MAT_file_v2) without .mat extension : ';
0505 mpc1 = loadcase(pfmatfilev2);
0506 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0507 t_is(mpc1.bus, bus, 12, [t 'bus']);
0508 t_is(mpc1.gen, gen, 12, [t 'gen']);
0509 t_is(mpc1.branch, branch, 12, [t 'branch']);
0510
0511 t = 'mpc = loadcase(pf_MAT_file_v2) with .mat extension : ';
0512 mpc1 = loadcase([pfmatfilev2 '.mat']);
0513 t_is(mpc1.baseMVA, baseMVA, 12, [t 'baseMVA']);
0514 t_is(mpc1.bus, bus, 12, [t 'bus']);
0515 t_is(mpc1.gen, gen, 12, [t 'gen']);
0516 t_is(mpc1.branch, branch, 12, [t 'branch']);
0517
0518 t = 'mpc = loadcase(pf_struct_v1) (no version): ';
0519 [baseMVA1, bus1, gen1, branch1] = feval(pfcasefile);
0520 clear c;
0521 c.baseMVA = baseMVA1;
0522 c.bus = bus1;
0523 c.gen = gen1;
0524 c.branch = branch1;
0525 mpc2 = loadcase(c);
0526 t_is(mpc2.baseMVA, baseMVA, 12, [t 'baseMVA']);
0527 t_is(mpc2.bus, bus, 12, [t 'bus']);
0528 t_is(mpc2.gen, gen, 12, [t 'gen']);
0529 t_is(mpc2.branch, branch, 12, [t 'branch']);
0530
0531 t = 'mpc = loadcase(pf_struct_v1) (version=''1''): ';
0532 c.version = '1';
0533 mpc2 = loadcase(c);
0534 t_is(mpc2.baseMVA, baseMVA, 12, [t 'baseMVA']);
0535 t_is(mpc2.bus, bus, 12, [t 'bus']);
0536 t_is(mpc2.gen, gen, 12, [t 'gen']);
0537 t_is(mpc2.branch, branch, 12, [t 'branch']);
0538
0539 t = 'mpc = loadcase(pf_struct_v2) : ';
0540 clear c;
0541 c.baseMVA = baseMVA;
0542 c.bus = bus;
0543 c.gen = gen;
0544 c.branch = branch;
0545 c.version = '2';
0546 mpc2 = loadcase(c);
0547 t_is(mpc2.baseMVA, baseMVA, 12, [t 'baseMVA']);
0548 t_is(mpc2.bus, bus, 12, [t 'bus']);
0549 t_is(mpc2.gen, gen, 12, [t 'gen']);
0550 t_is(mpc2.branch, branch, 12, [t 'branch']);
0551
0552
0553 delete([ matfile '.mat' ]);
0554 delete([ pfmatfile '.mat' ]);
0555 delete([ matfilev2 '.mat' ]);
0556 delete([ pfmatfilev2 '.mat' ]);
0557
0558 t = 'runpf(my_M_file)';
0559 mpopt = mpoption('verbose', 0, 'out.all', 0);
0560 [baseMVA3, bus3, gen3, branch3, success, et] = runpf(pfcasefile, mpopt);
0561 t_ok( success, t );
0562
0563 t = 'runpf(my_struct)';
0564 [baseMVA4, bus4, gen4, branch4, success, et] = runpf(c, mpopt);
0565 t_ok( success, t );
0566
0567 t = 'runpf result comparison : ';
0568 t_is(baseMVA3, baseMVA4, 12, [t 'baseMVA']);
0569 t_is(bus3, bus4, 12, [t 'bus']);
0570 t_is(gen3, gen4, 12, [t 'gen']);
0571 t_is(branch3, branch4, 12, [t 'branch']);
0572
0573 t = 'runpf(modified_struct)';
0574 c.gen(3,2) = c.gen(3,2) + 1;
0575 [baseMVA5, bus5, gen5, branch5, success, et] = runpf(c, mpopt);
0576 t_is(gen5(1,2), gen4(1,2) - 1, 1, t);
0577
0578 t_end;