From bd991fa3b2dde6e0107ae126b4b934e9437e3066 Mon Sep 17 00:00:00 2001
From: Brendan O'Donoghue <bodonoghue85@gmail.com>
Date: Fri, 5 Jul 2024 11:47:23 +0100
Subject: [PATCH] mpc bug test (#280)

* mpc bug test

* explicit memset to 0

* not sure what is failing

* fix n

* add some prints

* better print

* test disable res infeas

* disable other tests

* Revert "disable other tests"

This reverts commit b1a41ba75e657d593d4a5ea2605892748f2559c8.

* cleanup

* add more tests

* rename var

---------

Co-authored-by: Brendan O'Donoghue <bodonoghue@google.com>
---
 include/glbopts.h                        |   3 +-
 src/scs.c                                |  34 +++++++++++------------
 test/problems/mpc_bug.h                  |  19 +++++++++++++
 test/problems/mpc_bug1                   | Bin 0 -> 23753 bytes
 test/problems/mpc_bug2                   | Bin 0 -> 23753 bytes
 test/problems/mpc_bug3                   | Bin 0 -> 23753 bytes
 test/problems/test_prob_from_data_file.h |  22 +++++++++++++--
 test/run_tests.c                         |   3 ++
 8 files changed, 60 insertions(+), 21 deletions(-)
 create mode 100644 test/problems/mpc_bug.h
 create mode 100644 test/problems/mpc_bug1
 create mode 100644 test/problems/mpc_bug2
 create mode 100644 test/problems/mpc_bug3

diff --git a/include/glbopts.h b/include/glbopts.h
index 4d7b7fdf..ed2a15b5 100644
--- a/include/glbopts.h
+++ b/include/glbopts.h
@@ -39,7 +39,8 @@ extern "C" {
 #define WRITE_DATA_FILENAME (0)
 #define LOG_CSV_FILENAME (0)
 #define TIME_LIMIT_SECS (0.)
-
+/* Tolerance to check negativity condition for infeasibility */
+#define INFEAS_NEGATIVITY_TOL (1e-9)
 /* redefine printfs as needed */
 #if NO_PRINTING > 0     /* Disable all printing */
 #define scs_printf(...) /* No-op */
diff --git a/src/scs.c b/src/scs.c
index 9cb1ee58..7dfd092a 100644
--- a/src/scs.c
+++ b/src/scs.c
@@ -207,13 +207,13 @@ static void compute_residuals(ScsResiduals *r, scs_int m, scs_int n) {
   r->res_unbdd_a = NAN;
   r->res_unbdd_p = NAN;
   r->res_infeas = NAN;
-  if (r->ctx_tau < 0) {
+  if (r->ctx_tau < -INFEAS_NEGATIVITY_TOL) {
     nm_ax_s = NORM(r->ax_s, m);
     nm_px = NORM(r->px, n);
     r->res_unbdd_a = SAFEDIV_POS(nm_ax_s, -r->ctx_tau);
     r->res_unbdd_p = SAFEDIV_POS(nm_px, -r->ctx_tau);
   }
-  if (r->bty_tau < 0) {
+  if (r->bty_tau < -INFEAS_NEGATIVITY_TOL) {
     nm_aty = NORM(r->aty, n);
     r->res_infeas = SAFEDIV_POS(nm_aty, -r->bty_tau);
   }
@@ -587,22 +587,22 @@ static void print_summary(ScsWork *w, scs_int i, SCS(timer) * solve_timer) {
   scs_printf("\n");
 
 #if VERBOSITY > 0
-  scs_printf("Norm u = %4f, ", SCS(norm_2)(w->u, w->d->n + w->d->m + 1));
-  scs_printf("Norm u_t = %4f, ", SCS(norm_2)(w->u_t, w->d->n + w->d->m + 1));
-  scs_printf("Norm v = %4f, ", SCS(norm_2)(w->v, w->d->n + w->d->m + 1));
-  scs_printf("Norm rsk = %4f, ", SCS(norm_2)(w->rsk, w->d->n + w->d->m + 1));
-  scs_printf("Norm x = %4f, ", SCS(norm_2)(w->xys_orig->x, w->d->n));
-  scs_printf("Norm y = %4f, ", SCS(norm_2)(w->xys_orig->y, w->d->m));
-  scs_printf("Norm s = %4f, ", SCS(norm_2)(w->xys_orig->s, w->d->m));
-  scs_printf("Norm |Ax + s| = %1.2e, ", SCS(norm_2)(r->ax_s, w->d->m));
-  scs_printf("tau = %4f, ", w->u[w->d->n + w->d->m]);
-  scs_printf("kappa = %4f, ", w->rsk[w->d->n + w->d->m]);
-  scs_printf("|u - u_t| = %1.2e, ",
+  scs_printf("Norm u = %1.6e, ", SCS(norm_2)(w->u, w->d->n + w->d->m + 1));
+  scs_printf("Norm u_t = %1.6e, ", SCS(norm_2)(w->u_t, w->d->n + w->d->m + 1));
+  scs_printf("Norm v = %1.6e, ", SCS(norm_2)(w->v, w->d->n + w->d->m + 1));
+  scs_printf("Norm rsk = %1.6e, ", SCS(norm_2)(w->rsk, w->d->n + w->d->m + 1));
+  scs_printf("Norm x = %1.6e, ", SCS(norm_2)(w->xys_orig->x, w->d->n));
+  scs_printf("Norm y = %1.6e, ", SCS(norm_2)(w->xys_orig->y, w->d->m));
+  scs_printf("Norm s = %1.6e, ", SCS(norm_2)(w->xys_orig->s, w->d->m));
+  scs_printf("Norm |Ax + s| = %1.6e, ", SCS(norm_2)(r->ax_s, w->d->m));
+  scs_printf("tau = %1.6e, ", w->u[w->d->n + w->d->m]);
+  scs_printf("kappa = %1.6e, ", w->rsk[w->d->n + w->d->m]);
+  scs_printf("|u - u_t| = %1.6e, ",
              SCS(norm_diff)(w->u, w->u_t, w->d->n + w->d->m + 1));
-  scs_printf("res_infeas = %1.2e, ", r->res_infeas);
-  scs_printf("res_unbdd_a = %1.2e, ", r->res_unbdd_a);
-  scs_printf("res_unbdd_p = %1.2e, ", r->res_unbdd_p);
-  scs_printf("ctx_tau = %1.2e, ", r->ctx_tau);
+  scs_printf("res_infeas = %1.6e, ", r->res_infeas);
+  scs_printf("res_unbdd_a = %1.6e, ", r->res_unbdd_a);
+  scs_printf("res_unbdd_p = %1.6e, ", r->res_unbdd_p);
+  scs_printf("ctx_tau = %1.6e, ", r->ctx_tau);
   scs_printf("bty_tau = %1.2e\n", r->bty_tau);
 #endif
 
diff --git a/test/problems/mpc_bug.h b/test/problems/mpc_bug.h
new file mode 100644
index 00000000..e8fc2e03
--- /dev/null
+++ b/test/problems/mpc_bug.h
@@ -0,0 +1,19 @@
+#include "glbopts.h"
+#include "problems/test_prob_from_data_file.h"
+#include "scs.h"
+
+static const char *mpc_bug(void) {
+  const char *fail;
+  scs_float OPT1 = -0.473957794500; /* from scs */
+  scs_float OPT2 = -0.029336830816; /* from scs */
+  scs_float OPT3 = -0.002215217478; /* from scs */
+  fail = _test_prob_from_data("test/problems/mpc_bug1", OPT1);
+  if (fail) {
+    return fail;
+  }
+  fail = _test_prob_from_data("test/problems/mpc_bug2", OPT2);
+  if (fail) {
+    return fail;
+  }
+  return _test_prob_from_data("test/problems/mpc_bug3", OPT3);
+}
diff --git a/test/problems/mpc_bug1 b/test/problems/mpc_bug1
new file mode 100644
index 0000000000000000000000000000000000000000..989fbb3e4bb5863af6620ae8317575a1eed2a4c3
GIT binary patch
literal 23753
zcmeI32~^J68^_aFD@|#kX^|E}HOMx@)WbWJ<ki9qkzMw!QOJ;GWGh0}!M~!#Sc-<h
z41P3B7-otISxYCTP>5+}jQRg~zfaEjbvmu@TQToB`Se@vz0ZB__dfT&&rL%$HMNG~
zPd)L+!M=liN0F~Tv-y{DlFQklTtQOZx~ukhk*cX^-ink%N7VZ%*hj$y9%P^YAn&Js
zom5Wc^ppKpa(`W{2f4Vq%js+#O`QWSa+(#CemtdzoE}|v%`E67r>{5q3`lP$r%hk2
z9NNuRPQjm@T-wu0PQwzP1zk3mQ$XIvB}Yu<WT>6qGRa6z?e=6m=_1l_jlGY*$R|Us
z`&b#v$<@s`qLGQ5cHEs1e^})Cj=8-)5z7S@YV5QWnJTwSxhLMw^E2x_OyulLw-a_E
zbB&Lf>4`ixpxI$VkzO;0c()cgE&9k>N0G;M=Ec8yf7qjX!DB`Cb-MP2ugHB@fAa7V
z>9}sM!RI2+ST!+~i#(hj``y8K8FhJCkk{`!87*I6-NXL4jGDAPXsUTmMq|fE?Y90!
zM&sws-0|Zj83pcgoPF%NjH28$7PP%BqZK;i1AAx6DCe>DweHz6`qj&<ak9vqsVPw_
zL~8FyKXg4yMu!r!hCjb6BOj9u?#?1Ro!38lN-P)L;MtHHBJKC5+BXyL>xO2WjTE`e
zG?7w88f5j(cpy@?`=0afBKNM`8J#9_jd$-UCq?$mTjc)g{V3hP=4=spP-q#lS)^yH
zuKS`yo|^2L6DRWH>kl*ciyX9a@!!Gjy@a7=32QH!I}2xCOwi3}=prPmTR6r%a~95Z
zJ96ftrnAuD_Q~Y<)xCtPtAj7-Uz7<6BTRoEr_obbuiLq>?MxDOWp?~wq_d-t<t}fN
z-?h6i@Qk28VnA2HFme9ams1>sYga$pleEiHQ0VyG`&Sw_*#oZlWi6Nn8rD(#;t7U;
zr@#WxLS6Bm91H+Afd{}-;2AI(JO};^UIc#yQ^Ct%8h8~<2XBD4z&l_jm<48o_rV<S
z0r&{a2cLii;B&AL{0n>mz690YQ}{(4R9xjdRWnbZ{c2)Ea(MkQAm7&?Q*y}&)Zf?2
zlH{Xw27NTflAhZPU43e=1@*p~U*Ou*oL2t%ME$`vQ~Ewo9(dcygtnULhmV?JOl#Vw
zW<PysOh@bdb=4-`gdE2U6a6gA$RymLQ{);m`ZvQZcfct#`upIzRsF)vh-@1tB@H&C
z_RmiXznPhl)i3SWtVl5>+uiQ-b5EPn<mCf%Po|lY!aofK*eSqH0d@+oQ-GZU>=a<9
z06PWPDZowvb_%dlfSm&D6kw+SI|bM&z)k^n3b0eaI98XQs!2aQf8OZmm3*8dTNKk#
zWMTX4eKT4V(@~eSmlp#9jf?4M!JK)=PR}tdrX!C@8<(v*R$NCpxvA$T=9w1hC@C*$
zr|((TPd@zfMn|x-n2um)F&)9qVmgAI#dHKai|Gh<7U>9fR+oNMts5uI=c0Hc)eFZP
zsa`nVNcF<;MyeN%H&VTDykWg?yuth}(h0{IsZKb~NOi(-MyeBzGg6&!oRR8;;|%MB
z<BW>qi>h?tBwc@457PCA^&nk;SP#<mhxH&`e^?JCtv{>->3YLDkghkZ1L=ChI*_h6
ztOMzK!#XHwy;0HlD|y`CtP|C>-{5(#@aLQLfakla_Z2+Pl{~)68doRyv1lL3^H;i$
z<oPSzNAmoY?jw2rO81dGe|aCt@f`j`Ki`Zi&s*uZ^1PMq7kS=F_lrDlrTaymx4d8E
z_^qVztm<(r8b6+o((zN)dd%@xI&M4<rQ@cob(iC96^~a{jg!jPN1j(aE~;J!c^*~K
zcvNy6RJUGnTvPP>P5k2dAdO!<AEfb%_o34G#rx2b;#Vd8T~R+5#UI`WO5+dj11r0J
z@b>}Mwbakjc*Fa=%B~Z<&#$0=%j?%7e<~}U^L|?DN7kv-kE~OvA6ci;^_F$Y>n-n-
z%jv(e`mM-cRi!`vJ}UJO$5E+&%BwT>PZ|ADnjb3MfBrsR*?zJPO6gBY{V1|us&|eX
zQv0MjW&fAZeo5{AcmK7+IIs4$P)!ZTC~A^p1GTyefxJFAfAwp9TsHt;9jd?jG*UoK
z1Jneyz{a39_&#U`HV4hY7N7-a30i^HU`wzSXalwe+kkCBTaeG~d*Qk_=nT4meLy~+
z>WgbPupcM~`-ASF2j~fUfdjyS;Ai0H;2>}?I2s%Sjs<<eao~7x0yq(z1WpF0fK$O~
z;B;^X=m&lQ&IJ9z0B{zlY8~mIT~+JH9_{|8AB7{TN<XTKKe`w{b5Pmw=0mho)8kA{
zk1u?Tt*mvaiHh{l4(-&`^`@q-FUlHsJ&cF4_8aX{e@)w0)U^FWJ@~0U*Z^z@mUo}2
zk8)MLf9!yI{~J$*BTCv2R@FEfz~5G&HP{kt1=@h6?f)I{UQOBG*OdKxW#^kA{94)Z
z;Xlx>?(bdW(Vy!2K2_1U*A#YFbbfR~`Re+<TV4Kcgnny)nqWzJjAnRV-uTlQ<=!>_
zmesFD@K0Iec^AB|>b#7q^DCw4@dNmwHP{Ah3)+HRL1pKyl$~!XX<izme<k(*BRsFF
z{9IM#Kcsn56WA#ozmIW$AXr-du&ntzZPaH9vR=C3x;y9%egX~w9l@cX6X*jHI1KCo
z4hMUJBS0BA68sb#1@;0tKiV7DV?bwcEa(FIf_*^Vf4SnC*L{x9ydJsXInR5}|4zbv
zo=>d%{<zQiR?Z7g#eLQv>xc8ioQL(odow`J8*^Tk$DhZSb-;OA&L{h$Jm-}IaLxH;
zY5vv@&zZO^Iw#;fOnp4i0N)1<K-Mqs+qpiD*F4_rH`WpRoBhoGWk0fixSw;-UJw`z
zt^zlLk>DQiFnAig1ZILy!1wAXelY-T!0w<2I1Zc-ZUE!KvtSm;2QH@Siu${Og~P8H
zJPq1S&kKw`QH$6`zFThG*simiE^kfTmlYIEW*Gx|XSLf!FP-z`ha$fr6U}ys{pLoH
zX9$Jf)`_O+_qIA|o{FQ1U)}sBqr*3p;^40{deKJedppg%wN?b(Z5lP?R^(2~%-Zs=
znpYeR=zC!HnxtLSXW?jdTe~p&F(xd}W5)(M*QLdfSs(q6e%|N%hj)jaWZqZrdaJv!
zG+Fm#Zga0)WSC|1cyG6*bf^8)L1CL$(UArZt{5e)r}SS&g{cqTPWd-9T!oW6XyTz~
zU2<l`l1D~*gmy+8xj%kC^^*T$`uEUGt$<leD7Z~zw*$Y0P><mYZ$7#cP8*I*ZQL|q
z8`)f4s}<XS2c>R0VfW%?4Bg-3W!ZB|9DR8!d~MxJi)rbt<$?9~Eu@YU$9P*MET-(~
zE_x=~D=D;3qVK{N>!{_GIZNh;ZKa2HMuW6(MA6c@w(5&l#8BY&=k>J;W2yNDd*d@P
z!4xz)c>T17i%5S_xb=kOAkrAq>Q-R(Lh`)x>EnH}<<v^gO54Y99T`5IGWLGl*W}>1
z@zIEJQ8c)HNJ6203~BV9w?B7bEQMd0P1(*X=z|N>hY!2{6<v#1(b2(i0c8&gI1||E
zOB(w9wGA&#7Es5Q_cgWWETyFV509jMxR&mOhiH8*Y@s``Mh5MkMAD55>2p%IM^lXZ
zw5JZ5v9vm<rCt9at4UuoqSIRcQ0hA>WY*-I<#a+TEZ(o-60!&lQ9CFLq6^y+-FE2)
zke^AqUC8+vq~~LloPXM%7L3vf-R7`};w*Dds5M$mpWLaFP@u7i^knXhcgIE3&5hqa
znl@$$32DcZef|k1@968^izY3juJyM?=Qu2<u9rF=cS~DF{~QV%JvJ|xb`HFo`{T&P
zWY_s<m}}vev~0@@zfcQ5Qg0j^*d=8o-QVtKVlm5;IyVdVoUvdT)w%jxKUtHVv}$U8
z_sx6!$RQzm?@0$gGG0~KJjTP1Cbf&r80h!~<^Q1HH$~f@?oG?+{$9i^T6{C~*|BAT
zl#sk;{V}6?q;1u?(SaHB$)?R@6Wh3X^xd~c+J9;Vl8~L^a8zbVTIqv3gs%3cxmP0M
zA9_TPVpdtvz}KH=yd0@XuNCj9m3Y6b&&&H>Mc=P%xvDB()%7T=K4sUds`^#kPIa|c
z-R;(t_G_9QHPxP)ZdYxwueRA)TkWmwcE2t5zis_^Tl@32_v>Bh-@De&()!D3AbyXp
zkp|-1MR_;A&-pigL0XQV*X0+$(%!3A!S^{2$i>Ur4wvI|=gOAn^U|uSrx8~KR<vJi
z7uU!6Qm($L>(%6nm1H;9uchL48lydJkn_z<9m##Z)-C0JlQKTnE9ZOq6?~ub=hA+0
zKiFQjhuh(DhLvoW%Qvm!dicD&vi)MaxIQB#*~#_txvH}4=KA^EQdRBn`HZUD<MYMp
zYM1}tNp-j1Tt&KNUGaFc|Je_0hj|UPr$r67%c7RpXIaba<a6ViGEdFnKh`6UJNuRG
zwXUUpX1iOyCHC{*<<(aG@HlWkY~HeQ;r_IK%lp;l9qC`&cdVbb=r4bN-?I7D=ZI72
s!-EljoJnwu2v$6CiT`2ClPiIae4nos#S4r6E&<gQ{Xc5X6L9_i2SK5cEdT%j

literal 0
HcmV?d00001

diff --git a/test/problems/mpc_bug2 b/test/problems/mpc_bug2
new file mode 100644
index 0000000000000000000000000000000000000000..4fc1ca6d003a829196f0ed7c93ec2769815fb840
GIT binary patch
literal 23753
zcmeI330RG38^=#c(k5DnN-9dqln5W{o;fpRJ%ql`&RAk(t5K1iu@@m*gDFF0%Qg&V
z@LIo&F-Az!Sjsew5lWURpUHQQ^S^UlZ&z3AIW6W~SAM<k`z-f!KmX@`o^!v^kw_#u
z@}Ii$AGfyc+P0VT?Pmu6Qbuw;XVmld49~LJGl0}hW%qWZ8hW7JkHH>F&iA3;H81;2
zAJ9SVRLwZqf5neCpdjjWzOjrPvKQ%^HImUI?crvljAS(ARo?PR1~PIn8X0P@FQZYb
znmB6c$!OA(i~e=%%cxVc5Az*$WaM>qvdg%7GHUqa=%fkSGO|bra`BQgtf}z~Upbpf
z?(GxgTz~1v_xW=AHvg<=OL_Ud6YF_<$+_7mG}BW~1AWqRm2+~CcA}x2o=&pknu_P%
z<IG>n>m8gPQ=s65>r-Pi<oDf-@*Zf*d8|v9Ohx$v4TQ}`a;|QZbwfe-@ts#I?(gul
z)-;gQBjof&UBz>aaisl7N-0LOU;O!4N+$W+|NJUjO0jk81+UAM($Cq^i&LISDLUYf
zn+4CM6xpzI<g1rbO3XN=ee*9Vop4!uH7ZX^*=}2%H|9$z!)L|S<8n^_{n66ha%R4<
z?(Zq5zk5%qww!y?HXKQnmmj!3@7QrU$7e4#7vy~Q?DYSBl5^dP{EkoMJahe^*WZfg
z$D%Up$m{i8IYL9hE>0063gq{%WY=q)FK74j=gk%6`{j80Wy=|G!=#OZKWw@<O>uvi
zaPRM2InOSm&M)NrZ$j=yNB5S(jZ>lZ+iTeiA1pPUWuIj)B>JZh)=IG#LQ+S)?vve0
zkVWp&svFc&u-6LK-f+ZLaC$KJ+s4mZ2*M#B*QAlw!sUm3mp*A>C0vootaZ#Rg!{5#
zm%6(*5mG-oy7Hy3v0&Qm=W79L>Ih05zkmOV<0gBcCw^HACWB@g%3u0}Vc=QtIcTG)
zEGGkl!A;;H@GN){Oay-be*}L5e+HAl8(=b+0;Ynu!Mk7@cpuCFGr<SoL+}y!1k3?*
z!RO#hFb{kMz6J|GNgbtMG(qK6j%O9~H2SY1>X6&pk16@R{a8}>^iY!ods~u?HX77<
ztS!BC7`FQCJ{$5%$$8$hu{EvC&((bNizR&*Dhs{mX+c}8Od>}IG@~`HNtw?cH=`38
zuTmWLTad>DVPc?-6<I`@c8FPHMQ_r3XAL}SMgKdpZdIR1D<UU@bLR$Ik?YF~!c8kH
zvb*ZKX2ms2a@x~xZq@}$nzY<E>rApGDgC1(z)k^n3b0dvodWC>V5a~(1=uOTP62ib
zuv37Y0_+rErvN(z*eSqH0d@+oQ-GZU=CP{uR7LvX_47_g3i)`)v?!vZm^|0agn$-B
zbks4q;HTiwW<_)~f99N%7iL-((UFh;#$~Hc7S++itfb2mvn>mCbS``Q4!=vRpB(t-
zosM8<5gozKB07SdMRWu^i|7b;7SR#xEYuO~tSbGeTQ{Cq&xP?utQU?qV!d#@5$lEH
zjaV-nZ^U}xc*A<(c!Tv@s1uGeVx4fD5$lBGj94ceXT&<;I3v~x#~IcM#~C%p7j@~v
zQ@sDM9>n_(>p{H#upY$w59>j^|F9m4+kaRG;{AqoAl`3S2jcyPbs*kvSO?<$hILTf
zexs)OSMj{R+b62(yus^U>Cboj0k3y;pDTEst9X7@HLsrVW8pcH*RS{-$?I2qj^y<#
zK1cHU6`vz{{qi}I<2n3?alV^ZUbo_T<#j7QFY>w-pBH)EiqDI@Zuz{(@modnS=sYg
zIDfnz#q+1C{g~sgc;0v&iswyL`!2`ZN}jLEnkTjIkG!sUUevu0@;a)d`KaJLsBgdG
zxTYNUyZFWHK^(t$J&5BMpF_p*i_f9O#jgtbyS#obj6Zx16vrPv2Uc|d;O_&hYq6ij
z@rKWN72PNJoL^4=mesF?{!~>w=kv7KkE~O%A6chjKeA56`z`B~_gg+Em(hQv^;@C8
zDocO-eN^lpj-z7#lvQW!pHljxBtMjU{QP~qqT^&8lrWy+#!+a$SnnJ+#P*4G%Kk5=
z|Kj@n=l)xVabB&qP$I!=6p84yfkaCwkhd?+6~AiZx*n)_sjm3cS3;r->Vf*80cZ#|
z0Ik60pf%V6v;l2FJFq3#3bY3uz}BE6*amb0`Mv!oxb_0Q!R}xWkl#=B!gX)34=4ls
zg8e`rus=8e90>Y?{|5g94gv>*W5BWCIM5Fq4^99-11Ez1;3RM|I0c*vP6GqLK=5;L
zIv4~7gCU^0b!3Eo)vcel==UG}C>>E*`cYT>X^8o=22~wzK0rU!J<e43_`<KTRkbfQ
zQj;EB&`)*UZ>sD5qN;f}#(b#iyy1%WtJ}Gvx}6{D!cW>@Jx~WM>p9Z|^(y=P*beRf
zGoMOF6n7r1ta&trzwN-5U@OobbO1{_|GS}FbvfTxm-Bl?*P9vqTG9F8KhUn~?_Cox
zo~rsjRo=MQ7<QL;eRM$ms`|cLRsODzaqEJ5U~zejCU{=f_~VXx@0)*1>sNF5r?m0B
zBg(5gFQe}KN=bTb3O}?49l<uB6W9q<b>2$V`KIF5WiyPgxbb(!^UBK4RaX8(oF}z_
zo#OfX7w-FlCFKuGo4+$edzK*Ur3<dRf<wS=;84&590q!V!$AT^fFFV*!H>XEpcMQR
z{1_Y!egblS)C<>RL2qyz*d6o(dw_iY>WORK_c=cEe$*S!dEImV*B|$JJ+bcl;y&kF
zIWIf~_gR0eAI=kV9(DlA1%RA4=DaM=KhH1gfb+DRPYyzT&MODwn)A!z{H-&dGjUn?
zo`CZ(+IXM~HULdQ)-Rvixjl~8Jm2g$))D)g{mlMlKeB&#oHNniJa7rP3fu_BfP2AX
z;05qHcpuCK>u4x{F$EpKuAmP%9-Iqq0QZBJzzmRIxL9f`+wTJAjl5<0Y~CJv`P{sl
zBzhP5Mcux=&1etZ*m5`_W8O}(N+0Nz;j)Vgyt8FTV|J5;p3C7rv!bbg7=_<6+DX&?
z*y5>oHl8Maac6gWyWMonEy!rh!j07HUb1y-{b>5F@%Ev2V|LK}jHow~0r51j*P$6}
z&h4Te3&v<VIY-d>*obVOuQ$*S9a{_y>HHP_nBe!<kaj!Bx|i{<_P@o^q=sj*nh)4T
zW*H7o_jOrHX|7WSMQmP0$Ll@1WqxiwrCuE!p*eUP<=ocoDV+J5CLVp!@nJw5`J|^t
z8>Yuozo!k7t_LlmH%F)I2Zt=CC5|y&4&4l+4@WMz^CT^jHk_Pd&^UN2Ii#%BkL&w2
zC2czG{JJ2P9_$@p`_bfhnteBNt=9ELwDj)s(7Fi=sQtvTL+lPLqReUCjV%mUQh1NU
zehXf&qgJ<OE}k2)g&sSb4>G*HotDmW(p<D6mO{6^)Yi|7qvjjhHoF+RgyxM|vVN-0
zLNZwx+48f*d89km{%&aI0_vak@zVt9a<VtJGaPQVj?A7-9`_*rOLFtr_+-@h?KIdm
z>_DDLEa`g9IheH|jv{Z(piJ)-)bz@<kt2Tnf_{l!(caBtK4lIHz8KnJHVymkmkk9L
z^QnES2YQAxm(sa|kB=vRu$Iyy!}Pxtq9`rS+|(sEhHhU;otd<4C&l)g`piu)j#kfW
z<=l7ZYBJG_?yxo}oO+E83z_tAIi1#z*dM5~m~6tsBuAw4=*rf^y>~SXra+5S=djBG
zWIWtFG3P=M%^z(PzSV6Z#oK0`megNO-O@A;JlEYs#?pQUd*XM}osHi-nL2hc3CX7t
zhyT5VhV1-x$U^^R)Jc2m&WCQxsnd1$Q@xXy(cecS#*E8eLOXna%R2w*B64;=5z#Ym
zHZ6;K9T;vCNSX$5p&hS%N)NULTG)j2C-)|i{R8GNBaM`seWZ<c(5fjpT{rIyB)0=Q
z_nmPIq-Lw~n#cMClD|t_y06FQl=HnwuWN=u^vBfnu63eAXwjYU7blm6(t*S^>ra}`
zAwxU&`iBDMl7r(U3#a%w^zAq1hWYxTBxGK5J0Z0t{nWwj!dDNWS+}D1KlX_x<*KqH
z-?yI^yd6oz*UEB|V(*vsd0ESqw|qtGRaX7VZbwz^sk&WtwXg1es;a-L?zg)1U)}7e
zuJ%-SyK0JkHO<bNYHv-qySCV0+s08_$5Y$K^}dYneH&*<<K;9Em*Z=sf%tY|-i_~b
z{*4QW>+y3fE&!HPu5LNYa~_bZm$o0S$M2mhTA$ySR#rRpxgoH;<6^tGJ<gYM^OfDM
z9yhEaySaURHTPqH{tQ9RH#3by_xZYE3HKY7^0{#t%bAq3Jm=5F<Kl6!y=)Kn!}ZK6
z*e};_T*>Y5`|^s8i|yj}%vEG3x6AKURb@A~&+jeO)epa)QFnj*ezB_h<^Ok5)%`bD
zlWtj8Jm2hp_5<5tT}|z2QO)hLsUh~+)-XHyy>WF}r`GTv>yhW3{mS;Xtf78pyIa)~
z`}yzkYN~#C9(WuMwQOE^JgsZ_xE$Y;@wIu+#_5Fd^7r?s&5F-)&+x}bqW`*hz$1Ez
i@=5po-$&)%3iaUoe66gWSNL}cXs+=8QFESv+y6hi9d^F}

literal 0
HcmV?d00001

diff --git a/test/problems/mpc_bug3 b/test/problems/mpc_bug3
new file mode 100644
index 0000000000000000000000000000000000000000..bfed77f2de6a189d078b6ccdf7350ad672ffa20b
GIT binary patch
literal 23753
zcmeI33s{fm8^?cg=$uk%QPBahGD&QI<R+722d9xar&$cO6*ZJh4y(k%7*?heMq|uj
z>BYDV+x}r?q!rfsmy&Y~{onk)cdqO0>gxFWb+F%c<<tAV&-*<0b3fnbexBd`L`^1>
zsfj-=#UK0bJ-hc3`SP=ff2k(<JzIV+_42{<yT*`A*?if)70J!0YdWId!Qe1CgWaj*
z*dFWDu?|Y7V*1H`tGK^24Ki$&v~;1OecJ~Ww{W3%?Q-*iL{72Lb4hCMLL)Pe|1(}>
z@7(bwtwhefFlbYtiVGdi-5%*E@~G*OK58Q8Ur5S4+02C;?)i8e5xLdDIXqcp($#w^
z*F-MP-<Ompa$=!h!XuGOyb2DNi1gALSNvS0U;2C5Pegt^O)sEGWN?1Ab%DsDQ<>xQ
zMJ9V*+fYioyT#XTi(E9!!?jT4uD=XcmX=>**|(&W-lVU4SLEoMU3^Nfo68ooC}l!B
zv)Ok<Zp*j46!O%Gbls*WtN!Ce$9^?33>LZLdf|7j&z$Ic^MNfAMV>QCy!k|AlzsB-
z70;b$sL}NKCqxdp5pwpi$aVIv=KWtd(RuBN*dZc4&mHz1B+{*u_H!$dpBcEh8j6h9
zTl1c}$mN^cILkzCxxY(QMdat!31+QChJRfc+FE3=fmwPRk-GiN6AeT<eP!ic%HMQU
zt@K5%Iux6vCvtqd{6(eZ{WE@&m2zXs?K?Um_nKUZD7}7Ovi6@+PWUT(ldj0WFD;Qp
zxv2_M3c5PX8>udg-Mr0m_;3wjTmJEBgB>)4)razWEp<>A9xa|S{@?~x;k-eX^ZycC
z3hGl+4bzV`7rMo@)m~WAOxSv%-_z-9Wx~u$kApHZU!=@A{=>w}Tc4&(D%dqGVBo_P
zJFniZX(4~5$aVbs{aY?>!Ue;@^<Wy<wwe5w(O?*O5_|w!sL0E?fd1f@;C}EVcm_NR
zo&zs{zk>e-FM^lBH1G;|70duL!7T6wcniz{bHP0DHdqK2f%m`%;3KdYd<s4TUx2da
za=)m6@~ff&^-So$sGg`s_Aft%<n{7nN<*^)4fb}iB$pYw6W;N(q(|28uRZy#1r51U
z^kBGANBaE9J(b&kn9{L8m%u*<n9x=;gQyuk?P;Ce#hm*E?dfQ<r&p}^nULdb!Q0ou
zj7*{o9b(p*QAzg5{IMs^=%0gMM2w0uBeKyxb!ws+**!Wf{BCAOR=?S;Ta#)^HoM)H
z=ASmDxvR(JpGY$$`S_>_uv37Y0_+rErvN(z*eSqH0d@+oQ-GZU>=a<906PWPDZowv
zb_%dlfSm&D6kw<Dn(7Gi=ar61<>MS<UQS0b#dbM+eay@0s88C9U;P8ym(x*jz>?(C
z0jA}2<UVK9s)*$BI?Btxc;5StX_<~r-Py6z>u1&v&#zZHf}Q1b1Ut*=2zHj!5$r6d
zBiLC^N3gR@N3gTL^rP&!Q>GV=H&VTDypigK<Be1=9B-t0;dmp}3&$JQ3&$JGU-(0<
zEY3)E!f{5b6OJ=dop79y>V)HrR3{u~SSK83lpJ4_rHcX5^@sH!U4K{)()EY+AYFf0
z57PCA^-$6J!#a?zH>?Bcdc!)9t~aa$>3YLDkghkZgNoK0C5^v=$Nkkh!Si3)`wgD=
za(}*B59(`Q!Sh_f<EyA~9RNR;?IU^qO2?JwuXG>D^H;i$<oPSzNAmpTeI&<o_z(Sj
zHLg5wrTaymx6=J0&s*t!k>{;+zsU2J_lq3A6*Qi8J&tAL$MaD-eu`R;IsQt=jpw0s
z+!VF$a=fkM@v5tFQu_MH^NPns+3O(BqdFRoT8@MA)+>%{@_xUHUpybA@r&n!G=A|u
zR2si{A6ilTs-?ec>gTff!}~yK{Na6IZPyR}KES$``dJ!pc%N6>b%OW#HS}+F{aWTv
zMa6U8PfPvCI+gm7bt?5E>r}elvQBxu<$ZEB{a01LmHDf#^vB;vrT*bKD)mowb;kaw
zqCYD0Lyh~--^Xj)Pu4*t{i&!QW%f(;&T&I(pH!#p{|ef#sNH|>zjPSq)tU-rG907G
zB*zA_7IJ~SJUB1?wH2;a!O}zZ(vP|vG7V4@)B?3Z9k4ZM26g~Dg65zFXbD<@oxsju
z7tk8)3U&kE1Z_Y*w|^JcLqKP6C^!t{^QjTI9tn;DT|igR4Ri-bgJZz4;5cwR_#QX`
zoCwYYJ;7O^7x)1<8~hOT2Iqit!H>Xs;Cyfa=mYwKAA_HOexN_N5LC8~bkVM|_0t{g
z{;MCkBkD>&%8EaF7{897qT|h5Xs4manT8%;_!wJJ>rxvf>7fVOX{hT>LtS4KHSYQt
z4@K=a>`;G0+gCKS{X<LmsTHUSs)5zrXByzWy52wbM7{rxr`!=0?FZ{>91Y=bE3gyT
z8SDaDgO%<7?NP3w?C%@O{=K&Itu6do+wtK)(5~<AU9-`j`uaXq)3|2@yK6c>9Poa9
zec!Dwf2*V48lWavQ6A$BJg;v2>5cbZH~&`Eug36CRpWUdlvj3MM%nq5%JkR{e&`Bz
z1K$K~z_&p~=dBc-Z>nfswnzUe>i;`<URU|Ky2^h@^Q0!QQ#yYC!ToVyW%<LZ=I?Y+
zpDD<C>5J=rpa<9=oCG?8?}G!t$smDKz=7aYa1b~RbONV?gTWc#yCCOBhv3>1bOvXE
zLqRWa7|8pt;kf2?pW`#HM<el^=RN0t=ioliC)T|y?sLAC^TP9RpY_N3;XE<tVaK4H
z56F3A&dc)n^Z2q3I8V#@WIw#md1ZfGbADNxzqQ44CN9g)2{;ea3J)~E)}SHC`sIB)
z*T?ai$D94eI%0pbpV`0cNA?f*GXU)cfuUdoxCx8__kf4N)8HlW26zu_-c0_BA!rTu
z1Kq(7z@^|ua3A<HcnjnM7gH5^{e8jWsp*FIgLc!S2gdzn(TU`>B_m_I?ryrg^}ya+
zL9t|(J$A^g9*OkA`Hst>m|bL|+2g>dPorsc7={0-8%qoR-a0_@WCD4Ac70cN&s~&i
z@25L6WD||}Gp%D+t!TPwv}00c%uc#-YfFi2OahG^v47FJQ;9Te*-RB1+erE;F7l51
zHyi0(AM;5I-}#y@?Dcx=(Q_wt9HD=$%guP2t9K&5!<a;Bd&~Onw|!SqmfgGwk((pv
zi0bWh<5L^x>Tfe5RVHqyq704U!ijIl`_RKac|P&vo_#f1Cp&@M?zX;o$!|H89Qs7d
zf8h!W?H1E_|L<WmaO$$_g;`OwF?pW0k^eTbzOr5`-t`;0_+^UivlnrcyJw8$ppO!0
zac0!|7MGUO%FNY)E%z>?Uf!M_R!Pe#XTeZ?6P?c~eAof6WzW8#&glUwmPT%+0$bw=
zIvG1?<)=0(%h$wF;PywYw2I@Y!^ZCI&%}jN(9F;c^VLGgAS9~Ohi8LG!?R0fV9qic
zoi+IGUZ>U6Mc+zia@#Mc?fs8t<tBVZ_Kuqhr+u)4CfbE16&u8n#*ii7<u8k;sPsjY
z<GhC2{jy-{lxv^SAJJ=i**gYP&IJE6fewr5{bPS@d|?tyy*lS=>IAH$Q{NRFNqcKO
zWkrQ)eI;z6taxL?9`|A><Cm)e7q`b!oZI~S_L}juHmI|$>!h`0pc(D3-Y=X+%m`aJ
zH*YniXhrVxRa-$8;bF3aPC@j`wgV#*_58`#<f?7hc^}fBY<#xpv>ye}&<)>aA3_P1
z`6)8>wbVbWS<(ZIFG=6YO?!7jEM4F9ec^o16(pn`KRfw(D0#$Q^9Y%<ir#LuEjG`7
zHNAbQ_wkWwtLXWm$eFY5gwoD&H}ij*zMO1(AB`MdyqH#PdFC5#;Y%vo@qv9(r&I2B
zUlWUkqp9~BQKNl=S5dPozmIZivy&p`74_S^$CvDrV!u6M?@R3?iaW%)`_h~q@!8`X
zKc=D|4MwEu_|f0<v->rVUP#NYhd)eS6-Y^E*KJ5PUP3xnz18>oEG6r1b4_d#me3F1
z8|ysL3M3&X)&8iHC23uq*fV^s2Ys3zy|2JMn&h*}ipIVCJmlp_CcT!IlT~=Xs?V!i
zuBPQ{`(9nWU)S|0sy;>6tE~E!-A;YASKsY6l=d5%9Szl<hHh75v9Gb&*;wsu>~=R5
z`<q%nnreTVdcR(m{=IJftgOGB2I6vjjWiJ7F3Y>|ea^pe0qJ}Eyag8kD=XKshUGaA
z$gfwm9e$6`ooo9(pO@BEJ?dN$Skr#7U0fgMOS$^Gu2+*QR*>CXzm}5Q(MEeZAm^Kz
zx{~{RtyjtYHdTDCU(Ip`H7w8hb7{Z0A8aq%!|m{UZEM*szi(8>_3(LlZTrP`aec-L
zvXkrOb5%vz&GqxSrLx-L^BHBg$LEXn)h_?PllpGIgOYU1y5jL>|Fa+1j*bn~9`gon
zmqjD7&$5x($>+umWuA6~|5%Sa?(A2#w^JkaGuz#{iP+D7m)BVJ!{fmHux?`G!u{#m
z#QW9lHR<1*uUS8B&|m)kzGZXi=g5HYf`ieI&m=iUhsvJ}-S^{`d+C9We4nr7uNRm7
QT>`2r`+wA&C*b=34<v7y00000

literal 0
HcmV?d00001

diff --git a/test/problems/test_prob_from_data_file.h b/test/problems/test_prob_from_data_file.h
index 5e972bdb..3b993039 100644
--- a/test/problems/test_prob_from_data_file.h
+++ b/test/problems/test_prob_from_data_file.h
@@ -19,6 +19,8 @@ static const char *_test_prob_from_data(const char *file, scs_float OPT) {
   scs_float perr, derr;
   scs_int success;
   const char *fail;
+  scs_float xt_p_x;
+  scs_float *px = SCS_NULL;
 
   read_status = SCS(read_data)(file, &d, &k, &stgs);
 
@@ -28,12 +30,24 @@ static const char *_test_prob_from_data(const char *file, scs_float OPT) {
 
   stgs->eps_abs = 1e-6;
   stgs->eps_rel = 1e-6;
+  /* Force verbosity for the test */
+  stgs->verbose = 1;
 
   sol = (ScsSolution *)scs_calloc(1, sizeof(ScsSolution));
   exitflag = scs(d, k, stgs, sol, &info);
 
-  perr = SCS(dot)(d->c, sol->x, d->n) - OPT;
-  derr = -SCS(dot)(d->b, sol->y, d->m) - OPT;
+  if (d->P) {
+    /* px = Px */
+    px = (scs_float *)scs_calloc(d->n, sizeof(scs_float));
+    memset(px, 0, d->n * sizeof(scs_float));
+    SCS(accum_by_p)(d->P, sol->x, px);
+    xt_p_x = SCS(dot)(px, sol->x, d->n);
+  } else {
+    xt_p_x = 0.;
+  }
+
+  perr = 0.5 * xt_p_x + SCS(dot)(d->c, sol->x, d->n) - OPT;
+  derr = -0.5 * xt_p_x - SCS(dot)(d->b, sol->y, d->m) - OPT;
   scs_printf("primal obj error %4e\n", perr);
   scs_printf("dual obj error %4e\n", derr);
 
@@ -47,7 +61,9 @@ static const char *_test_prob_from_data(const char *file, scs_float OPT) {
   SCS(free_cone)(k);
   SCS(free_sol)(sol);
   scs_free(stgs);
-
+  if (px) {
+    scs_free(px);
+  }
   if (fail) {
     scs_printf("%s: FAILED\n", file);
   }
diff --git a/test/run_tests.c b/test/run_tests.c
index 337b5d7a..15094ce7 100644
--- a/test/run_tests.c
+++ b/test/run_tests.c
@@ -43,9 +43,11 @@ _SKIP(random_prob)
 #if NO_READ_WRITE == 0 /* reads / writes */
 #include "problems/hs21_tiny_qp_rw.h"
 #include "problems/max_ent.h"
+#include "problems/mpc_bug.h"
 #else
 _SKIP(hs21_tiny_qp_rw)
 _SKIP(max_ent)
+_SKIP(mpc_bug)
 #endif
 
 static const char *all_tests(void) {
@@ -61,6 +63,7 @@ static const char *all_tests(void) {
   mu_run_test(unbounded_tiny_qp);
   mu_run_test(random_prob);
   mu_run_test(max_ent);
+  mu_run_test(mpc_bug);
   mu_run_test(test_exp_cone);
   return 0;
 }