I’ve just made some changes to the dev version to handle this situation gracefully. It now prints out a messages telling you that your base and target cases are identical and returns after 0 iterations, rather than getting stuck in an infinite loop.
Ray > On Jan 15, 2015, at 5:25 PM, Abhyankar, Shrirang G. <abhy...@mcs.anl.gov> > wrote: > > Arun, > You have the same generation/load for both the base (mpcb) and target > (mpct) cases. What are you seeing is that runcpf is solving the base case > over and over again since there is zero power transfer. Use different > generation/load for mpcb and mpct and run the continuation power flow again. > See the MATPOWER User's Manual > <http://www.pserc.cornell.edu//matpower/manual.pdf> for a CPF example on a > test 9-bus case. > > Shri > > From: arun s nair <arunsnair...@gmail.com <mailto:arunsnair...@gmail.com>> > Reply-To: MATPOWER discussion forum <matpowe...@list.cornell.edu > <mailto:matpowe...@list.cornell.edu>> > Date: Thu, 15 Jan 2015 13:02:57 -0700 > To: <MATPOWER-L@cornell.edu <mailto:MATPOWER-L@cornell.edu>> > Subject: Matpower- Continuation Power Flow > > Hi All > > I am trying to use matpower to evaluate IEEE 39 bus system and I tried to run > an example code for runcpf (Continuation Power Flow) > > define_constants; > mpopt = mpoption('out.all',0,'verbose',2); > mpopt = mpoption('cpf.stop_at','nose','cpf.step',0.2); > mpopt = mpoption(mpopt,'cpf.plot.level',2); > mpcb = loadcase(t_case9_pfv2); % load base case > mpct = mpcb; % set up target case with > mpct.gen(:, [PG QG]) = mpcb.gen(:, [PG QG]) > mpct.bus(2,[PD]) = mpcb.bus(2, [PD]) > results = runcpf(mpcb, mpct, mpopt); > results.cpf > > The program is keep on running and is not converging. Can anybody please help > on this. > > If anyone can send me a running example for runcpf of matpower , it will be > of great help and if anybody has worked on IEEE 39 bus using runcpf, please > let me know > > > Thanking you all in advance. > > Regards > > Arun Nair > >