I should have also mentioned that you can already accomplish what you are after
by specifying a TXFR matrix for the input rather than a SLACK matrix.
Specifically, if D is the normalized slack matrix you want to use, you can
simply do …
txfr = eye(nb, nb) - D ./ sum(D);
PTDF = makePTDF(mpc, [],
First of all, thank you for your carefully written and detailed report. The
bottom line is … you’re right.
This looks like a classic example of the aphorism “if you didn’t test it, it
doesn’t work”. This undocumented* feature was never covered in the test suite,
so this bug has been there since