Implement CPF direction switching as in legacy CPF.

Ref: MATPOWER/matpower#23 Ref: MATPOWER/matpower#25
MATPOWER · Apr 15, 2021 · 53a8f54 · 53a8f54
1 parent b5034de
commit 53a8f54
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/lib/mp_network_acps.m b/lib/mp_network_acps.m
@@ -697,7 +697,9 @@ function cpf_add_node_balance_constraints(obj, mm, dm, mpopt)
                     end
                     if ~s.done
                         s.done = 1;
-                        s.warmstart = struct('nmt', nmt, 'dmt', dmt);
+                        dir_from_jac_eigs = isempty(find(nx.z(ad.npv+ad.npq+1:end-1) > 0, 1));
+                        s.warmstart = struct('nmt', nmt, 'dmt', dmt, ...
+                            'dir_from_jac_eigs', dir_from_jac_eigs);
                     end
                     s.evnts(i).msg = msg;
                 end