Ray, I understand there are two ways of calculating generation shift factor,
1. From makeB function, one could calculate the B prime matrix( Bp), then one can remove the row/column associated with the reference bus. THE TRANSFORMER TAP IS IGNORED while Bp is formed. Therefore, we can compute the flow on the lines by injecting 1 MW at a bus i by Bp\dP where dP(i)=1 only at bus i, zeros everywhere else, this flow is called generation shift factor. 2. We repeat the procedure above by setting 1MW injection at bus i, the injection at all other buses set to 0. Then call rundcpf, we could end up with the line flow which we is essentially generation shift factor. But the RUNDCPF DOES CONSIDER TAP. Which one is correct in calculating generation shift factor? Thanks Cuong
