refactor: little change

rspower/lib/runpf.txt

@@ -15,8 +15,12 @@ fn runpf() {
     pv  = find(bus_type == PV  && bus_gen_status);   // PV bus indices
     pq  = find(bus_type == PQ || ~~bus_gen_status);   // PQ bus indices
 
+    // generator info
+    on = find(slice(gen, [0], [GEN_STATUS-1, GEN_STATUS]) > 0);  // which generators are on?
+    gbus = get_multi(slice(gen, [0], [GEN_BUS-1, GEN_BUS]), on) - 1; // what buses are they at?
+
     // initial state
-    // V0    = ones(size(bus, 1), 1);            // flat start
+    //V0    = ones(nb, 1);            // flat start
     V0  = slice(bus, [0], [VM-1,VM]) .* exp(c(0,1) * pi/180 * slice(bus, [0], [VA-1,VA]));
     vcb = set(ones(size(V0)), qp, zeros(size(pq)));   // create mask of voltage-controlled buses
     k = find(get_multi(vcb, gbus));            // in-service gens at v-c buses
@@ -25,11 +29,11 @@ fn runpf() {
     gen_vg = get_multi(slice(gen, [0], [VG-1,VG]), on_k);
     gbus_k = get_multi(gbus, k);
     voltage_g = get_multi(V0, gbus_k);
-    V0 = set(V0, gbus_k, gen_vg ./ abs(voltage_g) .* voltage_g);
+    //V0 = set(V0, gbus_k, gen_vg ./ abs(voltage_g) .* voltage_g);
 
     Ybus = make_y_bus(baseMVA, bus, branch);
     repeat = 1;
 
     Sbus = make_sbus(baseMVA, bus, gen, 1, 1, 1);
-    return V0;
+    return pq;
 }
\ No newline at end of file