On 01/25/2018 03:47 PM, Ray Zimmerman wrote:
Hi Baljinnyam,
Thanks for your offer to contribute. This would be a nice addition to
MATPOWER.
I suggest starting by adding two new MATPOWER options, each with
default values of 0 …
opf.current_balance = 0 (complex power balance) or 1 (complex current
balance)
opf.v_cartesian = 0 (polar) or 1 (Cartesian)
Then find all the places in the existing OPF code where the current
(default) formulation is assumed and add conditionals to implement the
other variations. You might find it useful to use the MATLAB Profiler
with a simple runopf() call to map out the functions that are called.
I think much of the work will be in opf_setup(), where we’ll need to
use different variables and constraints.
Another thing to do early in the process is to duplicate the
t/t_opf_mips.m file for each of the 3 new variations. The initial goal
will be to get all of those tests to pass for each of the versions.
Of course, you will want to begin all of this on your own fresh fork
of the current MATPOWER development branch using the guidelines in the
Contributor’s Guide
<https://github.com/MATPOWER/matpower/blob/master/CONTRIBUTING.md#submitting-additions-or-modifications-to-the-code>.
Don’t hesitate to ask if you have questions or need any help along the
way.
Thanks again,
Ray
On Jan 24, 2018, at 4:40 AM, Baljinnyam Sereeter
<b.seree...@tudelft.nl <mailto:b.seree...@tudelft.nl>> wrote:
Dear Matpower developers,
I use Matpower for my work to solve OPF problems. As we know that the
power balance (mismatch) equations in Polar coordinates are used as
equality constraints in the OPF formulation in Matpower. However,
depending on the formulation of equality constraints (power or
current mismatch) and specification of the coordinates (Polar or
Cartesian), OPF problem can be formulated mathematically in four
different ways for a single physical formulation. I am working on
other three versions since some advantages can be obtained for the
derivatives or properties of matrix equations of Newton update step.
I have finished other versions theoretically and implemented on
Matlab. So I was wondering if you would like to include the other
versions in Matpower.
If you are interested, I would like to discuss how we can design the
implementation. Thanks in advance.
Regards,
Baljinnyam S
Dear Ray,
Thank you for your suggestions. I will begin with the current MATPOWER
development branch and try to implemented as you suggested. I am not
sure when I finish it since I have to work on my PhD. But I will try to
finish as soon as possible. I will contact you when I have a question.
Thank you.
Regards,
Baljinnyam