refactor: little change

rspower/examples/test_newtonpf.txt

@@ -12,9 +12,9 @@
 // get generator status
 // bus_gen_status = zeros(size(bus, 1), 1);
 // bus_gen_status(gen(:, GEN_BUS)) = gen(:, GEN_STATUS) > 0;
-nb = size(bus, 1);
-ng = size(gen, 1);
-Cg = sparse(slice(gen, [0], [GEN_BUS-1,GEN_BUS]), (1:ng)', slice(gen, [0], [GEN_STATUS-1,GEN_STATUS]) > 0, nb, ng);  // gen connection matrix
+nb = size(bus, 0);
+ng = size(gen, 0);
+Cg = sparse(slice(gen, [0], [GEN_BUS-1,GEN_BUS]), range(0,ng - 1), slice(gen, [0], [GEN_STATUS-1,GEN_STATUS]) > 0, nb, ng);  // gen connection matrix
                                         // element i, j is 1 if, generator j at bus i is ON
 bus_gen_status = Cg * ones(ng, 1);      // number of generators at each bus that are ON
 
@@ -29,9 +29,9 @@ pq  = find(bus_type == PQ || ~~bus_gen_status);   // PQ bus indices
 // V0    = ones(size(bus, 1), 1);            // flat start
 V0  = slice(bus, [0], [VM-1,VM]) .* exp(c(0,1) * pi/180 * slice(bus, [0], [VA-1,VA]));
 vcb = ones(size(V0));           // create mask of voltage-controlled buses
-vcb(pq) = 0;                    // exclude PQ buses
-k = find(vcb(gbus));            // in-service gens at v-c buses
-V0(gbus(k)) = gen(on(k), VG) ./ abs(V0(gbus(k))).* V0(gbus(k));
+//vcb(pq) = 0;                    // exclude PQ buses
+//k = find(vcb(gbus));            // in-service gens at v-c buses
+//V0(gbus(k)) = gen(on(k), VG) ./ abs(V0(gbus(k))).* V0(gbus(k));
 
 Ybus = make_y_bus(baseMVA, bus, branch);
 r = make_sbus(baseMVA, bus, gen, 1, 1, 1);

rspower/lib/newtonpf.txt

@@ -54,5 +54,4 @@ fn newtonpf(Ybus, Sbus, V0, ref, pv, pq, mpopt) {
 
     // evaluate F(x0)
     mis = V .* conj(Ybus * V) - Sbus(Vm);
-   ·
 }
\ No newline at end of file