refactor:add more file

rspower/examples/test_newtonpf.txt

@@ -6,6 +6,8 @@
 #include ../lib/dsbus_dv.txt
 #include ../lib/make_jac.txt
 #include ../lib/make_sdzip.txt
+#include ../lib/make_sbus.txt
 
 Ybus = make_y_bus(baseMVA, bus, branch);
-return make_sdzip(baseMVA, bus, 1);
\ No newline at end of file
+r = make_sbus(baseMVA, bus, gen, 1, 1, 1);
+return r;
\ No newline at end of file

rspower/lib/make_sbus.txt

@@ -29,5 +29,11 @@
 // define named indices into bus, gen matrices
 
 fn make_sbus(baseMVA, bus, gen, mpopt, Vm, Sg) {
-    nb = size(bus, 1);
+    nb = size(bus, 0);
+    Sd = make_sdzip(baseMVA, bus, mpopt);
+    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?
+    ngon = size(on, 0);
+    Cg = sparse(gbus, range(0,ngon - 1), 1, nb, ngon);
+    return Cg;
 }
\ No newline at end of file