First of all, MIPS is definitely not a very robust LP solver, so you should use
another solver if one is available. It appears (from the results of rundcopf)
that you do have the Optimization Toolbox, i.e. linprog available. I suggest
trying the linprog dual-simplex solver.
mpopt = mpoption('linprog.Algorithm', 'dual-simplex', 'most.solver', 'OT');
Hope this helps,
Ray
> On Jul 20, 2017, at 10:23 PM, Alain Golderman <[email protected]> wrote:
>
> Dear all,
>
> I'm starting to use the MOST software and I was trying to launch a simple
> test case that was already working with Matpower (DC OPF)
>
> The test case is constitue of 42 buses with 37 loads, 1 generator and 61
> lines.
>
> The only difference between the test case that I run in matpower and in MOST
> is the insertion of a simple load profile in the 37 buses.
>
> """"""""""
> function loadprofile = load_profile
>
> [CT_LABEL, CT_PROB, CT_TABLE, CT_TBUS, CT_TGEN, CT_TBRCH, CT_TAREABUS, ...
> CT_TAREAGEN, CT_TAREABRCH, CT_ROW, CT_COL, CT_CHGTYPE, CT_REP, ...
> CT_REL, CT_ADD, CT_NEWVAL, CT_TLOAD, CT_TAREALOAD, CT_LOAD_ALL_PQ, ...
> CT_LOAD_FIX_PQ, CT_LOAD_DIS_PQ, CT_LOAD_ALL_P, CT_LOAD_FIX_P, ...
> CT_LOAD_DIS_P, CT_TGENCOST, CT_TAREAGENCOST, CT_MODCOST_F, ...
> CT_MODCOST_X] = idx_ct;
>
>
> loadprofile = struct( ...
> 'type', 'mpcData', ...
> 'table', CT_TLOAD, ...
> 'rows', [2 3 4 5 6 7 8 9 10 11 12 14 15 16 17 18 20 21 22 24 25 26 27 28
> 29 30 32 33 34 35 36 37 38 39 40 41 42], ...
> 'col', CT_LOAD_ALL_PQ, ...
> 'chgtype', CT_REP, ...
> 'values', [] );
>
>
> loadprofile.values(:, 1, :) = [
> 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62
> 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62
> 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62
> 1.62 1.62 1.62 1.62 1.62 1.62 1.62 1.62;
> 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16
> 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16
> 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16
> 2.16 2.16 2.16 2.16 2.16 2.16 2.16 2.16;
>
> """"""""""
> I run the DC OPF with the MIPS solver, which gives a correct and reasonable
> solution when using it with matpower. However, in the MOST application, the
> solver fails to find a solution to the problem.
>
> "" The error message """
>
> MATPOWER Interior Point Solver -- MIPS, Version 1.2.2, 16-Dec-2016
> (using built-in linear solver)
> Numerically Failed
>
> Did not converge in 3 iterations.
>
> ============================================================================
> - MOST: LP solver 'MIPS' failed with exit flag = -1
> You can query the workspace to debug.
> When finished, type the word "return" to continue.
>
> - Post-processing results.
> - MOST: Done.
>
> """""""""""
>
> I turned the problem up and down and I still can't get rid of this problem,
> and I was hoping that I will find a solution or a part of it by emailing you.
>
> Please accept, the assurance of my highest consideration
>
> Sincerely yours,
>
> Alain Golderman
>
>
>
>
> PS :
>
> Results of the rundcopf, its converge !
>
> ""
> MATPOWER Version 6.0, 16-Dec-2016 -- DC Optimal Power Flow
> The interior-point algorithm uses a built-in starting point;
> ignoring user-supplied X0.
> Optimization terminated.
>
> Converged in 0.58 seconds
> Objective Function Value = 7.99 $/hr
