refactor: little change

rspower/lib/runpf.txt

@@ -22,18 +22,22 @@ fn runpf() {
     // initial state
     //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
+    vcb = set(ones(size(V0)), pq, zeros(size(pq)));   // create mask of voltage-controlled buses
     k = find(get_multi(vcb, gbus));            // in-service gens at v-c buses
     gbus_index = get_multi(gbus, k);            // indices of in-service generators at v-c buses
     on_k = get_multi(on, k);
     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);
+    V_init = set(V0, gbus_k, gen_vg ./ abs(voltage_g) .* voltage_g);
 
     Ybus = make_y_bus(baseMVA, bus, branch);
     repeat = 1;
+    while repeat {
+        // compute bus power injections
+        Sbus = make_sbus(baseMVA, bus, gen, 1, 1, 1);
+        repeat = 0;
+    }
 
-    Sbus = make_sbus(baseMVA, bus, gen, 1, 1, 1);
-    return pq;
+    return V_init;
 }
\ No newline at end of file