Here it is:
>> t_most_w_ds(0, 'DEFAULT', 2)
1..1
=============================================================================
MATPOWER Optimal Scheduling Tool -- MOST Version 1.0.1
A multiperiod stochastic secure OPF with unit commitment
----- Built on MATPOWER -----
by Carlos E. Murillo-Sanchez, Universidad Nacional de Colombia--Manizales
and Ray D. Zimmerman, Cornell University
(c) 2012-2018 Power Systems Engineering Research Center (PSERC)
=============================================================================
- Building indexing structures.
- Building expected storage-tracking mechanism.
- Building constraint submatrices.
- Building DC flow constraints.
- Splitting storage injections into charge/discharge.
- Building CCV constraints for piecewise-linear costs.
- Building contingency reserve constraints.
- Building ramping transitions and reserve constraints.
- Building storage constraints.
- Building dynamical system constraints.
- Building cost structures.
- Assembling full set of costs.
- Assembling full set of constraints.
- Assembling full set of variable bounds.
- Calling QP solver.
============================================================================
Gurobi Version 8.1.0 -- interior point QP solver
Academic license - for non-commercial use only
Optimize a model with 190492 rows, 69548 columns and 502789 nonzeros
Model has 14508 quadratic objective terms
Coefficient statistics:
Matrix range [3e-04, 4e+02]
Objective range [1e-02, 2e+03]
QObjective range [7e+01, 3e+02]
Bounds range [1e-01, 1e+02]
RHS range [1e-02, 8e+00]
Presolve removed 96168 rows and 16036 columns
Presolve time: 0.37s
Presolved: 94324 rows, 77872 columns, 309792 nonzeros
Presolved model has 14508 quadratic objective terms
Ordering time: 0.00s
Barrier statistics:
Free vars : 7844
AA' NZ : 1.199e+06
Factor NZ : 2.491e+07 (roughly 300 MBytes of memory)
Factor Ops : 5.025e+10 (roughly 1 second per iteration)
Threads : 2
Objective Residual
Iter Primal Dual Primal Dual Compl Time
0 3.21127763e+09 -5.60608500e+08 5.05e+04 3.44e+04 2.09e+06 7s
1 1.23510591e+09 -5.18031778e+08 1.67e+04 1.14e+04 8.06e+05 9s
2 3.09251550e+08 -4.70853904e+08 3.98e+03 2.71e+03 2.19e+05 10s
3 4.86452629e+07 -3.87386002e+08 5.57e+02 3.79e+02 3.44e+04 12s
4 5.03195119e+06 -2.84694039e+08 2.50e+01 1.70e+01 3.34e+03 14s
5 3.01371313e+06 -1.14411122e+08 9.80e-01 6.67e-01 7.82e+02 15s
6 2.46579202e+06 -1.25694347e+07 8.03e-02 5.47e-02 9.54e+01 17s
7 1.67961319e+06 -3.34425330e+06 2.73e-03 1.86e-03 3.17e+01 19s
8 1.42578260e+06 1.83460040e+05 6.39e-05 4.36e-05 7.84e+00 21s
9 1.24063177e+06 9.38673372e+05 1.01e-06 3.83e-06 1.91e+00 23s
10 1.17154487e+06 1.00539424e+06 4.91e-08 2.98e-07 1.05e+00 25s
11 1.13851434e+06 1.11144905e+06 1.30e-08 5.70e-08 1.71e-01 27s
12 1.13408494e+06 1.11747899e+06 6.49e-09 2.85e-08 1.05e-01 28s
13 1.13116265e+06 1.12115682e+06 1.31e-09 5.75e-09 6.31e-02 30s
14 1.12943362e+06 1.12393229e+06 3.08e-10 1.36e-09 3.47e-02 32s
15 1.12786553e+06 1.12669317e+06 2.50e-09 4.93e-09 7.40e-03 34s
16 1.12769444e+06 1.12702149e+06 3.85e-09 2.35e-09 4.25e-03 36s
17 1.12761652e+06 1.12716655e+06 2.31e-09 1.41e-09 2.84e-03 38s
18 1.12748736e+06 1.12738832e+06 1.06e-08 4.33e-11 6.25e-04 41s
19 1.12746902e+06 1.12741631e+06 1.89e-07 2.08e-11 3.33e-04 44s
20 1.12746151e+06 1.12742771e+06 4.26e-05 1.87e-08 2.13e-04 46s
21 1.12744909e+06 1.12744595e+06 1.40e-04 4.47e-08 2.04e-05 48s
22 1.12744908e+06 1.12744606e+06 1.39e-04 4.24e-08 1.92e-05 50s
23 1.12744899e+06 1.12744667e+06 1.00e-04 2.69e-08 1.26e-05 51s
24 1.12744850e+06 1.12744708e+06 1.13e-04 1.84e-08 8.28e-06 53s
25 1.12744936e+06 1.12744714e+06 1.03e-04 1.66e-08 7.56e-06 54s
26 1.12744878e+06 1.12744739e+06 8.27e-05 1.34e-08 4.97e-06 56s
Barrier solved model in 26 iterations and 56.38 seconds
Optimal objective 1.12744878e+06
============================================================================
- MOST: QP solved successfully.
- Post-processing results.
- MOST: Done.
not ok 1 - dynamical system state (Z)
index got expected abs(got - exp)
--------------- ---------------- ---------------- ----------------
(13,24) 0.0424164 0.0425214 0.000104919 *
max diff @ (13,24) = 0.000104919 > allowed tol of 0.0001
Ran 1 of 1 tests: 0 passed, 1 failed
Elapsed time 68.55 seconds.
> On 20 Nov 2018, at 16:28, Ray Zimmerman <[email protected]> wrote:
>
> I’ve been unable to reproduce this failure so far, even on the same version
> of Matlab. The only difference is my OS is still High Sierra 10.13.6.
>
> Would you mind downloading yet another version of t_most_w_ds.m from the same
> URL
> <https://raw.githubusercontent.com/MATPOWER/most/master/lib/t/t_most_w_ds.m>,
> then run ...
>
> t_most_w_ds(0, 'DEFAULT', 2)
>
> … and send the result? This should run t_most_w_ds exactly as is run when you
> type test_most, but with verbose output so I can see what’s going on.
>
> Thanks for helping to debug these test failures.
>
> Ray
>
>
>
>> On Nov 20, 2018, at 9:50 AM, Jack Heron <[email protected]
>> <mailto:[email protected]>> wrote:
>>
>>
>> Hi Ray,
>>
>> Progress, but still something failed I think. This is what it gave me:
>>
>> >> test_most
>> t_most_3b_1_1_0........ok
>> t_most_3b_3_1_0........ok
>> t_most_30b_1_1_0.......ok
>> t_most_30b_3_1_0.......ok
>> t_most_fixed_res.......ok
>> t_most_30b_1_1_0_uc....ok
>> t_most_uc..............ok (156 of 260 skipped)
>> t_most_suc.............ok (111 of 185 skipped)
>> t_most_w_ds............not ok
>> ##### Ran 1 of 1 tests: 0 passed, 1 failed
>> Ran 615 of 615 tests: 347 passed, 1 failed, 267 skipped
>> Elapsed time 95.96 seconds.
>>
>>
>> Anything else I should try?
>>
>> Thanks,
>>
>> Jack
>>
>>
>>
>>> On 20 Nov 2018, at 15:43, Ray Zimmerman <[email protected]
>>> <mailto:[email protected]>> wrote:
>>>
>>> Oops … sorry. That new test file introduced a new bug, which should be
>>> fixed in this version
>>> <https://raw.githubusercontent.com/MATPOWER/most/master/lib/t/t_most_w_ds.m>.
>>>
>>> But, no worries with your simulations, these bugs were in the tests, not in
>>> MOST itself.
>>>
>>> Ray
>>>
>>>
>>>> On Nov 20, 2018, at 7:17 AM, Jack Heron <[email protected]
>>>> <mailto:[email protected]>> wrote:
>>>>
>>>>
>>>> Hi Ray,
>>>>
>>>> Thanks for your help. A fresh Matlab session did indeed reveal the Gurobi
>>>> version number. I’ve also got my hands on CPLEX and the latest version of
>>>> Matlab in effort to fix this. My new ‘mpver' is shown below.
>>>>
>>>> I updated the most/lib/t/t_most_w_ds.m file as you described. This time
>>>> ’test_most' ran quickly, although it produced a new set of errors. The
>>>> result is also below.
>>>>
>>>> I’m proceeding with my simulation despite this test, but if there is a fix
>>>> it would give me confidence going forward. Please do let me know if
>>>> there’s anything I should try.
>>>>
>>>> Thanks,
>>>>
>>>> Jack
>>>>
>>>>
>>>> >> test_most
>>>> t_most_3b_1_1_0........ok
>>>> t_most_3b_3_1_0........ok
>>>> t_most_30b_1_1_0.......ok
>>>> t_most_30b_3_1_0.......ok
>>>> t_most_fixed_res.......ok
>>>> t_most_30b_1_1_0_uc....ok
>>>> t_most_uc..............ok (104 of 260 skipped)
>>>> t_most_suc.............ok (74 of 185 skipped)
>>>> t_most_w_ds............Not enough input arguments.
>>>>
>>>> Error in t_most_w_ds (line 52)
>>>> if isempty(solver)
>>>>
>>>> Error in t_run_tests (line 61)
>>>> feval( test_names{k}, ~verbose );
>>>>
>>>> Error in test_most (line 57)
>>>> all_ok = t_run_tests( tests, verbose );
>>>>
>>>>
>>>>
>>>> >> mpver
>>>>
>>>> MATPOWER Version 7.0b1 31-Oct-2018
>>>> MATLAB Version 9.5 24-May-2018 Release: (R2018b)
>>>> Optimization Toolbox Version 8.2 24-May-2018 Release: (R2018b)
>>>> MIPS Version 1.3 30-Oct-2018
>>>> MOST Version 1.0.1 30-Oct-2018
>>>> SDP_PF -- not installed --
>>>> YALMIP -- not installed --
>>>> BPMPD_MEX -- not installed --
>>>> CLP -- not installed --
>>>> CPLEX Version 12.8.0.0
>>>> GLPK -- not installed --
>>>> Gurobi Version 8.1.0
>>>> IPOPT -- not installed --
>>>> KNITRO -- not installed --
>>>> MINOPF -- not installed --
>>>> MOSEK -- not installed --
>>>> PARDISO -- not installed --
>>>> PDIPMOPF -- not installed --
>>>> SCPDIPMOPF -- not installed --
>>>> SDPT3 -- not installed --
>>>> SeDuMi -- not installed --
>>>> TRALMOPF -- not installed --
>>>> Architecture: MACI64
>>>>
>>>> MATPOWER 7.0b1 is distributed under the 3-clause BSD License.
>>>> Please see the LICENSE file for details.
>>>>
>>>>
>>>>> On 19 Nov 2018, at 17:30, Ray Zimmerman <[email protected]
>>>>> <mailto:[email protected]>> wrote:
>>>>>
>>>>> Hi Jack,
>>>>>
>>>>> A few comments.
>>>>>
>>>>> - I notice that the code in that test used to select the best solver
>>>>> available had some errors, causing it to attempt to use MIPS (known to
>>>>> fail on this one) rather than Gurobi. I think I’ve fixed that problem in
>>>>> the latest version on the MOST GitHub repository. Try replacing your
>>>>> most/lib/t/t_most_w_ds.m file with this one
>>>>> <https://raw.githubusercontent.com/MATPOWER/most/master/lib/t/t_most_w_ds.m>
>>>>> and see if that works.
>>>>>
>>>>> - You may also have an issue with the installation of Gurobi on your
>>>>> machine. It seems odd to me that mpver does not print the version number,
>>>>> as it does for me. If typing mpver in a fresh Matlab session does not
>>>>> print the Gurobi version number, I’d be interested in following up with
>>>>> you to uncover whether there’s bug in MATPOWER or its an issue with your
>>>>> installation or what.
>>>>>
>>>>> - Finally, it’s always a challenge with numerical optimization to find
>>>>> the right options and tolerances to make tests pass across a bunch of
>>>>> different architectures, Matlab versions, solvers, etc. and this test is
>>>>> particularly troublesome, which is why I now moved it last. In the end,
>>>>> if everything else passes except this test, I wouldn’t worry about it.
>>>>>
>>>>> Ray
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>> On Nov 19, 2018, at 5:31 AM, Jack Heron <[email protected]
>>>>>> <mailto:[email protected]>> wrote:
>>>>>>
>>>>>>
>>>>>> Hi,
>>>>>>
>>>>>> I have an update. This time I let ’test_most’ run for almost 45 minutes,
>>>>>> and it finally came up with the following warnings and error message.
>>>>>> Does anyone know what this could be?
>>>>>>
>>>>>> Any help would be much appreciated.
>>>>>>
>>>>>> Thanks,
>>>>>>
>>>>>> Jack
>>>>>>
>>>>>>>> test_most
>>>>>> t_most_3b_1_1_0........ok
>>>>>> t_most_3b_3_1_0........ok
>>>>>> t_most_30b_1_1_0.......ok
>>>>>> t_most_30b_3_1_0.......ok
>>>>>> t_most_fixed_res.......ok
>>>>>> t_most_30b_1_1_0_uc....ok
>>>>>> t_most_uc..............ok (156 of 260 skipped)
>>>>>> t_most_suc.............ok (111 of 185 skipped)
>>>>>> t_most_w_ds............> In mplinsolve (line 70)
>>>>>> In mips (line 467)
>>>>>> In qps_mips (line 180)
>>>>>> In qps_matpower (line 252)
>>>>>> In most (line 2029)
>>>>>> In t_most_w_ds (line 247)
>>>>>> In t_run_tests (line 61)
>>>>>> In test_most (line 57)
>>>>>> Warning: Matrix is close to singular or badly scaled. Results may be
>>>>>> inaccurate. RCOND = 1.450050e-16.
>>>>>>> In mplinsolve (line 70)
>>>>>> In mips (line 467)
>>>>>> In qps_mips (line 180)
>>>>>> In qps_matpower (line 252)
>>>>>> In most (line 2029)
>>>>>> In t_most_w_ds (line 247)
>>>>>> In t_run_tests (line 61)
>>>>>> In test_most (line 57)
>>>>>> Warning: Matrix is close to singular or badly scaled. Results may be
>>>>>> inaccurate. RCOND = 1.437052e-16.
>>>>>>> In mplinsolve (line 70)
>>>>>> In mips (line 467)
>>>>>> In qps_mips (line 180)
>>>>>> In qps_matpower (line 252)
>>>>>> In most (line 2029)
>>>>>> In t_most_w_ds (line 247)
>>>>>> In t_run_tests (line 61)
>>>>>> In test_most (line 57)
>>>>>> Warning: Matrix is close to singular or badly scaled. Results may be
>>>>>> inaccurate. RCOND = 1.169742e-19.
>>>>>>> In mplinsolve (line 70)
>>>>>> In mips (line 467)
>>>>>> In qps_mips (line 180)
>>>>>> In qps_matpower (line 252)
>>>>>> In most (line 2029)
>>>>>> In t_most_w_ds (line 247)
>>>>>> In t_run_tests (line 61)
>>>>>> In test_most (line 57)
>>>>>> Warning: Matrix is close to singular or badly scaled. Results may be
>>>>>> inaccurate. RCOND = 6.418324e-28.
>>>>>> Reference to non-existent field 'Z'.
>>>>>>
>>>>>> Error in t_most_w_ds (line 252)
>>>>>> t_is(mdo.results.Z, s.Z, 4, t);
>>>>>>
>>>>>> Error in t_run_tests (line 61)
>>>>>> feval( test_names{k}, ~verbose );
>>>>>>
>>>>>> Error in test_most (line 57)
>>>>>> all_ok = t_run_tests( tests, verbose );
>>>>>>
>>>>>>
>>>>>>
>>>>>>> On 16 Nov 2018, at 16:43, Jack Heron <[email protected]
>>>>>>> <mailto:[email protected]>> wrote:
>>>>>>>
>>>>>>> Hi,
>>>>>>>
>>>>>>> I have just now installed matpower. When I run ‘test_matpower'
>>>>>>> everything is OK. However when I then run ‘test_most’, the operation
>>>>>>> gets stuck and seems to go on indefinitely - until I’m forced to
>>>>>>> terminate execution.
>>>>>>>
>>>>>>> I’ve tried several times, fresh install, etc. But it always gets stuck
>>>>>>> at the same place (’t_most_w_ds’). I have the latest version of Gurobi
>>>>>>> installed, and an older Matlab Optimisation toolbox.
>>>>>>>
>>>>>>> Can anyone help?
>>>>>>>
>>>>>>> This is the command window text I’m faced with:
>>>>>>>
>>>>>>>>> test_most
>>>>>>> t_most_3b_1_1_0........ok
>>>>>>> t_most_3b_3_1_0........ok
>>>>>>> t_most_30b_1_1_0.......ok
>>>>>>> t_most_30b_3_1_0.......ok
>>>>>>> t_most_fixed_res.......ok
>>>>>>> t_most_30b_1_1_0_uc....ok
>>>>>>> t_most_uc..............ok (156 of 260 skipped)
>>>>>>> t_most_suc.............ok (111 of 185 skipped)
>>>>>>> t_most_w_ds............Operation terminated by user during mplinsolve
>>>>>>> (line 70)
>>>>>>>
>>>>>>>
>>>>>>> In mips (line 467)
>>>>>>> dxdlam = mplinsolve([M dg; dg' sparse(neq, neq)], [-N; -g],
>>>>>>> opt.linsolver, mplinsolve_opt);
>>>>>>>
>>>>>>> In qps_mips (line 180)
>>>>>>> [x, f, eflag, output, lambda] = mips(p);
>>>>>>>
>>>>>>> In qps_matpower (line 252)
>>>>>>> qps_mips(H, c, A, l, u, xmin, xmax, x0, mips_opt);
>>>>>>>
>>>>>>> In most (line 2029)
>>>>>>> mdo.QP.lambda ] = qps_matpower( mdi.QP.H, mdi.QP.C, ...
>>>>>>>
>>>>>>> In t_most_w_ds (line 247)
>>>>>>> mdo = most(mdi, mpopt);
>>>>>>>
>>>>>>> In t_run_tests (line 61)
>>>>>>> feval( test_names{k}, ~verbose );
>>>>>>>
>>>>>>> In test_most (line 57)
>>>>>>> all_ok = t_run_tests( tests, verbose );
>>>>>>>
>>>>>>>
>>>>>>> ————
>>>>>>>
>>>>>>> This is my output to ‘mpver’ :
>>>>>>>
>>>>>>>>> mpver
>>>>>>>
>>>>>>> MATPOWER Version 7.0b1 31-Oct-2018
>>>>>>> MATLAB Version 9.1 25-Aug-2016 Release:
>>>>>>> (R2016b)
>>>>>>> Optimization Toolbox Version 7.5 25-Aug-2016 Release:
>>>>>>> (R2016b)
>>>>>>> MIPS Version 1.3 30-Oct-2018
>>>>>>> MOST Version 1.0.1 30-Oct-2018
>>>>>>> SDP_PF -- not installed --
>>>>>>> YALMIP -- not installed --
>>>>>>> BPMPD_MEX -- not installed --
>>>>>>> CLP -- not installed --
>>>>>>> CPLEX -- not installed --
>>>>>>> GLPK -- not installed --
>>>>>>> Gurobi Version
>>>>>>> IPOPT -- not installed --
>>>>>>> KNITRO -- not installed --
>>>>>>> MINOPF -- not installed --
>>>>>>> MOSEK -- not installed --
>>>>>>> PARDISO -- not installed --
>>>>>>> PDIPMOPF -- not installed --
>>>>>>> SCPDIPMOPF -- not installed --
>>>>>>> SDPT3 -- not installed --
>>>>>>> SeDuMi -- not installed --
>>>>>>> TRALMOPF -- not installed --
>>>>>>> Architecture: MACI64
>>>>>>>
>>>>>>> MATPOWER 7.0b1 is distributed under the 3-clause BSD License.
>>>>>>> Please see the LICENSE file for details.
>>>>>>>
>>>>>>> ————
>>>>>>>
>>>>>>> I’m running MacOS Mojave version 10.14
>>>>>>>
>>>>>>> Thanks,
>>>>>>>
>>>>>>> Jack
>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>
>>>
>>
>
