I suspect you have isolated buses or islands that do not have a swing bus. Have you tried I-iv in MatPower faq on power flow divergence.
http://www.pserc.cornell.edu/matpower/#pfconvergence Shri On Sep 2, 2016, at 6:51 AM, Gamze Dogan <[email protected]<mailto:[email protected]>> wrote: Dear all, I have a problem when using runpf which does not converge. Therefore, I tried to check if the system load exceeded the steady-state loading limit (as mentionned on the website) define_constants; mpcbase = loadcase('casefile'); mpcbase.bus(:, PD) = 0; mpcbase.bus(:, QD) = 0; mpcbase.gen(:, PG) = 0; mpctarget = loadcase('casefile'); results = runcpf(mpcbase, mpctarget); results.cpf.max_lam But I have a problem using the function runcpf, it gives me the following error message: Warning: Matrix is singular to working precision. > In cpf_corrector at 124 In runcpf at 257 step 1 : lambda = NaN, corrector did not converge in 10 iterations ??? NaN's cannot be converted to logicals. Error in ==> printpf at 175 nzld = find((bus(:, PD) | bus(:, QD)) & bus(:, BUS_TYPE) ~= NONE); Error in ==> runcpf at 392 printpf(results, 1, mpopt); I do not understand what the problem is. Any suggestion is wellcome. Thank you, G. Dogan
