0001 function t_opf_dc_bpmpd(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 num_tests = 23;
0038
0039 t_begin(num_tests, quiet);
0040
0041 [PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
0042 VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
0043 [GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
0044 MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
0045 QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
0046 [F_BUS, T_BUS, BR_R, BR_X, BR_B, RATE_A, RATE_B, RATE_C, ...
0047 TAP, SHIFT, BR_STATUS, PF, QF, PT, QT, MU_SF, MU_ST, ...
0048 ANGMIN, ANGMAX, MU_ANGMIN, MU_ANGMAX] = idx_brch;
0049
0050 casefile = 't_case9_opf';
0051 if quiet
0052 verbose = 0;
0053 else
0054 verbose = 0;
0055 end
0056
0057 t0 = 'DC OPF (BPMPD): ';
0058 mpopt = mpoption('out.all', 0, 'verbose', verbose);
0059 mpopt = mpoption(mpopt, 'opf.dc.solver', 'BPMPD');
0060
0061
0062 if have_fcn('bpmpd')
0063
0064 ib_data = [1:BUS_AREA BASE_KV:VMIN];
0065 ib_voltage = [VM VA];
0066 ib_lam = [LAM_P LAM_Q];
0067 ib_mu = [MU_VMAX MU_VMIN];
0068 ig_data = [GEN_BUS QMAX QMIN MBASE:APF];
0069 ig_disp = [PG QG VG];
0070 ig_mu = (MU_PMAX:MU_QMIN);
0071 ibr_data = (1:ANGMAX);
0072 ibr_flow = (PF:QT);
0073 ibr_mu = [MU_SF MU_ST];
0074 ibr_angmu = [MU_ANGMIN MU_ANGMAX];
0075
0076
0077 load soln9_dcopf;
0078
0079
0080 t = t0;
0081 [baseMVA, bus, gen, gencost, branch, f, success, et] = rundcopf(casefile, mpopt);
0082 t_ok(success, [t 'success']);
0083 t_is(f, f_soln, 3, [t 'f']);
0084 t_is( bus(:,ib_data ), bus_soln(:,ib_data ), 10, [t 'bus data']);
0085 t_is( bus(:,ib_voltage), bus_soln(:,ib_voltage), 3, [t 'bus voltage']);
0086 t_is( bus(:,ib_lam ), bus_soln(:,ib_lam ), 3, [t 'bus lambda']);
0087 t_is( bus(:,ib_mu ), bus_soln(:,ib_mu ), 2, [t 'bus mu']);
0088 t_is( gen(:,ig_data ), gen_soln(:,ig_data ), 10, [t 'gen data']);
0089 t_is( gen(:,ig_disp ), gen_soln(:,ig_disp ), 3, [t 'gen dispatch']);
0090 t_is( gen(:,ig_mu ), gen_soln(:,ig_mu ), 3, [t 'gen mu']);
0091 t_is(branch(:,ibr_data ), branch_soln(:,ibr_data ), 10, [t 'branch data']);
0092 t_is(branch(:,ibr_flow ), branch_soln(:,ibr_flow ), 3, [t 'branch flow']);
0093 t_is(branch(:,ibr_mu ), branch_soln(:,ibr_mu ), 2, [t 'branch mu']);
0094
0095
0096
0097
0098
0099
0100 mpc = loadcase(casefile);
0101 mpc.A = sparse([1;1;1;2;2;2],[10;11;13;11;12;14],[-1;1;-1;1;-1;-1],2,14);
0102 mpc.u = [0; 0];
0103 mpc.l = [-Inf; -Inf];
0104 mpc.zl = [0; 0];
0105
0106 mpc.N = sparse([1;2], [13;14], [1;1], 2, 14);
0107 mpc.fparm = ones(2,1) * [1 0 0 1];
0108 mpc.H = sparse(2,2);
0109 mpc.Cw = [1000;1];
0110
0111 t = [t0 'w/extra constraints & costs 1 : '];
0112 [r, success] = rundcopf(mpc, mpopt);
0113 t_ok(success, [t 'success']);
0114 t_is(r.gen(1, PG), 116.15974, 5, [t 'Pg1 = 116.15974']);
0115 t_is(r.gen(2, PG), 116.15974, 5, [t 'Pg2 = 116.15974']);
0116 t_is(r.var.val.z, [0; 0.3348], 4, [t 'user vars']);
0117 t_is(r.cost.usr, 0.3348, 4, [t 'user costs']);
0118
0119
0120 mpc = loadcase(casefile);
0121 mpc.A = sparse([1;1;1;2;2;2],[19;20;25;20;21;26],[-1;1;-1;1;-1;-1],2,26);
0122 mpc.u = [0; 0];
0123 mpc.l = [-Inf; -Inf];
0124 mpc.zl = [0; 0];
0125
0126 mpc.N = sparse([1;2], [25;26], [1;1], 2, 26);
0127 mpc.fparm = ones(2,1) * [1 0 0 1];
0128 mpc.H = sparse(2,2);
0129 mpc.Cw = [1000;1];
0130
0131 t = [t0 'w/extra constraints & costs 2 : '];
0132 [r, success] = rundcopf(mpc, mpopt);
0133 t_ok(success, [t 'success']);
0134 t_is(r.gen(1, PG), 116.15974, 5, [t 'Pg1 = 116.15974']);
0135 t_is(r.gen(2, PG), 116.15974, 5, [t 'Pg2 = 116.15974']);
0136 t_is(r.var.val.z, [0; 0.3348], 4, [t 'user vars']);
0137 t_is(r.cost.usr, 0.3348, 4, [t 'user costs']);
0138
0139 t = [t0 'infeasible : '];
0140
0141 mpc = loadcase(casefile);
0142 mpc.A = sparse([1;1], [10;11], [1;1], 1, 14);
0143 mpc.u = Inf;
0144 mpc.l = 600;
0145 [r, success] = rundcopf(mpc, mpopt);
0146 t_ok(~success, [t 'no success']);
0147 else
0148 t_skip(num_tests, 'BPMPD_MEX not available');
0149 end
0150
0151 t_end;