From fd2c872a2bbc462e062f50ffb7a4e749dfd527c1 Mon Sep 17 00:00:00 2001 From: Debbie Zhuang Date: Wed, 3 Nov 2021 17:30:13 -0400 Subject: [PATCH] changed test case start/end points --- mpet/mod_CCCVCPcycle.py | 8 ++++---- tests/ref_outputs/test025/params_system.cfg | 4 ++-- .../test025/sim_output/commit.diff | 15 +++++++++++++++ .../test025/sim_output/input_dict_anode.p | Bin 1327 -> 1327 bytes .../test025/sim_output/input_dict_cathode.p | Bin 1329 -> 1329 bytes .../test025/sim_output/input_dict_system.p | Bin 2722 -> 2722 bytes .../test025/sim_output/output_data.mat | Bin 22568 -> 22128 bytes .../test025/sim_output/run_info.txt | 4 ++-- tests/ref_outputs/test029/params_system.cfg | 4 ++-- .../test029/sim_output/input_dict_cathode.p | Bin 1329 -> 1329 bytes .../test029/sim_output/input_dict_system.p | Bin 2569 -> 2569 bytes .../test029/sim_output/output_data.mat | Bin 11168 -> 12160 bytes .../test029/sim_output/run_info.txt | 4 ++-- 13 files changed, 27 insertions(+), 12 deletions(-) diff --git a/mpet/mod_CCCVCPcycle.py b/mpet/mod_CCCVCPcycle.py index 24488445..00077f39 100644 --- a/mpet/mod_CCCVCPcycle.py +++ b/mpet/mod_CCCVCPcycle.py @@ -238,8 +238,8 @@ def DeclareEquations(self): self.time_counter() < dae.Constant( 0 * s)) if capfrac_cutoffs[i] is None else ( (self.ffrac_limtrode() - <= 1 - capfrac_cutoffs[i]) if limtrode == "c" else ( - self.ffrac_limtrode() >= capfrac_cutoffs[i])) + < 1 - capfrac_cutoffs[i]) if limtrode == "c" else ( + self.ffrac_limtrode() > capfrac_cutoffs[i])) # for capacity condition, cathode is capped at 1-cap_frac, anode is at cap_Frac # if end state, then we send back to state 0 and also add one to cycle_number if i == len(constraints)-1: @@ -415,9 +415,9 @@ def DeclareEquations(self): cap_cond = ( self.time_counter() < dae.Constant( 0 * s)) if capfrac_cutoffs[i] is None else ( - (self.ffrac_limtrode() >= 1 + (self.ffrac_limtrode() > 1 - capfrac_cutoffs[i]) if limtrode == "c" else ( - self.ffrac_limtrode() <= capfrac_cutoffs[i])) + self.ffrac_limtrode() < capfrac_cutoffs[i])) # voltage cutoff v_cond = (self.time_counter() < dae.Constant( 0*s)) if voltage_cutoffs[i] is None else \ diff --git a/tests/ref_outputs/test025/params_system.cfg b/tests/ref_outputs/test025/params_system.cfg index 4d659782..1730c1ae 100644 --- a/tests/ref_outputs/test025/params_system.cfg +++ b/tests/ref_outputs/test025/params_system.cfg @@ -6,8 +6,8 @@ Vmax = 5 Vmin = 2.5 Vset = 0.12 segments = [ - [1, None, 0.05, None, None, 4], - [-1, None, 0.45, None, None, 1], + [1, None, 0.06, None, None, 4], + [-1, None, 0.44, None, None, 1], ] prevDir = false totalCycle = 4 diff --git a/tests/ref_outputs/test025/sim_output/commit.diff b/tests/ref_outputs/test025/sim_output/commit.diff index 8b137891..bcc21696 100644 --- a/tests/ref_outputs/test025/sim_output/commit.diff +++ b/tests/ref_outputs/test025/sim_output/commit.diff @@ -1 +1,16 @@ +diff --git a/tests/ref_outputs/test025/params_system.cfg b/tests/ref_outputs/test025/params_system.cfg +index 4d65978..1730c1a 100644 +--- a/tests/ref_outputs/test025/params_system.cfg ++++ b/tests/ref_outputs/test025/params_system.cfg +@@ -6,8 +6,8 @@ Vmax = 5 + Vmin = 2.5 + Vset = 0.12 + segments = [ +- [1, None, 0.05, None, None, 4], +- [-1, None, 0.45, None, None, 1], ++ [1, None, 0.06, None, None, 4], ++ [-1, None, 0.44, None, None, 1], + ] + prevDir = false + totalCycle = 4 diff --git a/tests/ref_outputs/test025/sim_output/input_dict_anode.p b/tests/ref_outputs/test025/sim_output/input_dict_anode.p index 9124f55cbfae106262c176765e5f26fbf160982b..709d53a4e08e1166cbbeb57b0f56bd993f19749b 100644 GIT binary patch delta 106 zcmZ3_wVrE3En_%40|SG9ZfbgBd}5)PUjz$~>spo)UMTJtfsjfnltAVt7fSj?@B&rl mCYGcYWhUmtmsA#{7D|D2`(!$s6-sZu!nl)(QEYM_%R&J0Y#n_7 delta 106 zcmZ3_wVrE3En_$@0|P^DVo7RIW@1i!No7H5p;!b9P|PRO*{o0;%ycbF2``lJi(m%| k`{$;nC&nihN+NTU3#E{`NrlpWo3AkLWMY(@+{dyI01J!3dK0|SG9ZfbgBd}5)PUjz$~>spo)UMTJtfsjfnltAVt7fSj?@B&rl mCYGcYWhUmtmsA#{7D|D2`(!$s6-sZu#<-J-QEc)AmW2TG^c{x) delta 106 zcmdnUwUKK>J!3d80|P^DVo7RIW@1i!No7H5p;!b9P|PRO*{o0;%ycbF2``lJi(m%| k`{$;nC&nihN+NTU3#E{`NrlpWo3AnMWMY(@Jb`5)02BiqiU0rr diff --git a/tests/ref_outputs/test025/sim_output/input_dict_system.p b/tests/ref_outputs/test025/sim_output/input_dict_system.p index 001ee008855452f36dfb21c10b14f442a26da44b..c0891886379501f199e1195d0bea11b5bc5bf58a 100644 GIT binary patch delta 87 zcmZ1^x=3_`3Zv4x9f7Y~<#zb_d9#!hifFj+FP?c#U`72GKX>~(8eb3f$7oJ=W4x!# o%D})-kYAo!RH)<^!4BjkClnIJIB&(E8sq)j2 jEO0AIOfFP%xBmd-0zm{TkSNG6Pc14`_S@{v6wLtu6%8Mx diff --git a/tests/ref_outputs/test025/sim_output/output_data.mat b/tests/ref_outputs/test025/sim_output/output_data.mat index b7f16a8e1065675682996c6a3412b8bc41b40dff..01fbe23105afbef5e3cacefcc73a856761a5a2ef 100644 GIT binary patch literal 22128 zcmeI42|QKX+s98dX)cwNN-9HTibliQ(&SRoAf>K(%8<~cghG)>1C<7qp@@()P?U-; zmvVDZ6dFj5QcAP;*=O%{{8U!hDDrlETw?Dy1K@?^Q_ie zIoX(5*;}~T*_c^7xHJEP#T+sCF{eB6KQFTWCL&fx7=LygVR|_Jy}EdPVf@>5=(?Tw z&bkePb-$_`zj`jKoh@zLhaq;-NUaX71fM)L#eYpGGQ9I0HmkZe5knQNPe^2eHa5(=donJXF^6TYl?_g%( z>gwcRV@+?Ru)682;WWQ*EzSH1M{B+*Yo3@%T?2A7A{<$J_t?dOuRPVc)KsyN8Xd8FN-6_qYH3+Yg%rpZ|M# z{>pjH!@(IYBEPktF>cVm2lYqcFH1eOf0T9GH@>@f&6*O%<}C-;G^V&vBsJ8F`Eg;K zIq=DH`ORJ4a{p# z8kE#Fn`sdJ1>+L>0zD!6f_|TgJPG|aI)eO&?Z)F{VjW_A#EwVmX5sbe`b7Pk{fFLf z@IKTH{%vLEwBEyJ);u!<<}dl{>pzSVqs5UlByTYWg83USgE2Ge=8p8cAmC#HAN0-f58CH(#&QH;(Gf}56E zOXSB_Ax&d6o7S01@JIF|`Z)(vCo_xhtb85*N3bNZj*Jv?U?VZ4_O!7 z#`Te7>qFKBw{d;Qy5KghPaRtyvM#vo#ny+c3vT~peXPE>KGNvcr0tg;py9pF$N^3GoqZPQWu zF#l&p=abM`g7GR*M+at7mZuS z=#C{gdE2hL_O%5lygu_&if<|!--pprC+ICDTRJI?ixk=ZN(AeA2S(3o%txbJPKK=v zO-3i#?IUSZImRtdr4)@}`^yphx%yZM>Z7!EV0mH710;-pTz#<5Y4T&#J#dKfX>Kht zFkF1S$)g_a7_i$^U7UxU5AHG6+0H{w1oOg&jlOxK4#|({c}%TOtx(Fc%Q4T31XtUr z^UzgRFPorA{OBVQF?Gn2?Po=>(O4Lyl~Jay!G*c)QfC>W~fF-%Xy3#fJ!|`Pbm7?O@L3s2ZD@AF6Dgi5LGmc1 zZ#RjGHnfTC*N={qJQ%fWVAbk0q;M#(({#`!`L@)?BU+w|q~iyNbb&Q4DlQ*2f;!x&1mHHCWzP<%A_8l^0!N zbU}xte;8_Y9Z`y;GOKDjgAPeyzT&v(XXm0_;BVIl6EB<2D@WUgxI_*He?v2`9G(XL zjxet}DGB3}jGHtsQ#`Q*EfU+*AU&v}4O8qycyyd(q=(y|SHWM~#W%eNf-cED(@(q! zR^X!F-uH?a2K`BP5!q*E2mY$Zc}p3CzXmh>DkMOMq>}pa`HkT3y7!qw={!if#aC)A z1AqUTlmDr9XEIUM>lnzhAU!0_c)F5v=|AVPgR*I@+*61N=Sp`|SPophL22=kW{6 zW4XvFz4wTTphGfA@A9KV;P14acDw0(NE(%#tp9CfDOw>D*`Lmbq%dD` zT)Ef&fZ`q1Xz@Y=h0;0YX!lvi`(8DLsG|5>ag$v-8hS`;P6CWea=^KbPc=`MAkRlq z7MBKBv|)wZQGT4{lInLh8^PZeGd0DNpi8pYp>z4gv0PLi@~2M*=#tzT8{0G%{4I>? z96td3)xR+LQ#|O94AhTbB6EU^ES)CG(K;ku3_9-61%I>OHq<79ziy2u?wkM}lHO}R zTnHakiuAOuEuwWuVjqPuqbwjv3;d0{Ii=HM{c`jvC%;jFqxrPB(t-)WJDQZjVFC?|0_zIzB9yXTjvh-y>+A(pc_iXU@)?4rHbUq}NJXdP0lP^WHWA=#9+es4pD2%&w zc{`)FSEGXZox?ZJEJrSG(+}oX7NRU(WKAdV_vWU3XY65IlD~UCa|}LPg7T)x@r-3E z+Hk1h34WYpqlAZFqDLCqwoWX}0dz?=*20LKup}-zetwzX0?;8jP0mTbZ&N;sh)I~!5&ShBHzk|aAt|x+etiG2 zrO3qoTn(KMNnyU?c-Wp)>4>0eG-X=%_q}w>QGkNz?LGGj;d)HjKOF8)hjJ=fZ^8X3 z$?7)>oUP|eP}g^U)lWxPwBd4Y{~9_@lBaxW3Xq>h(d3VCe|oEEpCX+f$vi)kRgo!Nbm)+>aV_YOG%{P;bLaDXG~6_M zMjG6oE|(8i{tI+S4jR|RXS5QW-yVfJ*Mkm8?4vN&+kfWR!u@}ReZ0?LRE}PKG7u3f zhwF!F?jZxg-;cfC-5d$yl6)F1zVH>?|0lcIZWkF-(T2RnRDPV~4x2;j<#7MMWPfnn zXwW5je1T6*z+^6(qI|T4)+Kq@ICB;63H<%0gf{E}e|x1q*f9chNJ`JX&{y&-7p;6< zThjq_NG5mqY|rK8Bb|W~i%r1aN^hH~vY~CWt3OXc(`HJJRcQ(B1 zv8@_aKj(g0h$@dV5G>ABB;c(laGEpc;kG`mI;$gmUB*@#gi@;zE?xtn=)Vbvl|P;WH*Q zJ*f@N`*@xoajgWsiZ`jrQmbgg#JZc6beyEAR*=axw>0Fp{nm{j&?OnUXn6QVii>tG z4q2ED{Yi?3JASCk4oQ(rN$S#5 zxPR+?|KUdPH^^VnC=hf=R{uG)hs4xU6y#CbkKP`VKk|Hpe*g2E_uG~hg7-g#J!i!6 zf1Wq~`t!!`eb0J`ka}5LLBXs&V5#u$`-ME;ZtwXwUM}=GI)Bh_{QA!EgzA20|912H z@6!2^{p_pfkXHZhIpn{L)R0*4a_K#NgsI84wTL*md@{2<3zJ=*h3Qx1q)>12Ndt3L zW<_;pvrigWA9E@(&JvPiMNlSz`vh< z4*5mBzI+ZjmVM4IjJMsRt$v&I7UQ}B1|)Fz)9;4c+_)lvogwft~kU-y1IG{I5X z%M|=2x$jMpiV67p$WbEH2mB?OY|NXNFUjkF zBnN1NzrMz|MmvI!B#mbz^5Z0<$IADy0Dm>@7+sRmt9npf!QaM%j4sJ$qvx}dz~Auv zL6yhB-zF}jLvlm%yr$dW@4{C_{5(in7azRv0Q^1oZo``=;O{XfMu+6=`<{m*!C!H2 zCLfZ*e8uqu%NcL{!QX@%PiA<5zwVcgo>>k4cJe~2r+~js122ergTEx}q9^Ph4F2lQ z4;?oL{3RKs<;0JZ>^OUoswen+s|TY?GOP3c+()tSJzyB4OLDeMSUnf~b>f&W%m9D8 ztz>ja%6TOyfXFdFDs(_ik@ShosumzWvsNzd`OyJ|wY^ z!ssCF<{JzCDz$J{o&bL}e61HJfWONw^=@1S{toVP%Ig65OA>WWYpTz!L$*)Ot{4se zl9YU~iytTHtlrr#68tq4WpqjIPhFUC7W^GAmCo-^vZ+Y2xJNqpQlsT|7yQjx!03=% zrIxbE2mDRk!{kTOI^eOMHu!sOtbe2=e9v0K=0lQm?~%4X_&dCn(IF|!R~(m=`^%;X z{B>E=8h;=B4UVkLYy^LGT84QBg1??qpMK5+e@VJhgZrNbe{a>Nta}CDW0TA|6T**^ z9Qvxy^6TL5QcXseWcI9UQJ&y$bN^I+f0ECF#|6&=e>q1h{5ve-A@OyL4#|NhSIf@? ze@DeJ`H_^m($Hly_&eNwaitdcyNumVlD)M~CiVq?72h*DB(aaeIGES;2M690m}wZL zBn|HgD0=T1MZtRl#d{`DS>Ufsv|h{;@Rwxvv&VV1;O}6a;qvFeUy>0eALZyc$=#xV zRJDM=Rq>22$<2mCb=1J$>g9|s$;qD7l4IcS?*k-KRKVX=MT`!~=F+;et$uZ=vNxkc zQciMi`Cs7gW}k^i9l+nav5XE$X&q680(eh=n#Jgl6y__AzkIdl##HdPVfP3LXYkjg zk0n0j%0L5VjqR^s=Se{7x??|Q{$TJ;IDFt)B2vP z>k;y>zc^DIt_#Cj!wkV+l4;Gho0D(Up^B3i`ws_yNj@GN!;h2HsN}nnR1%khI)4F-? z1AmQ1Csc)iza;J41~^Xxe-mamY2@^-ZNq(Ug86Zhhwer3a=_m~tS-s&+%bpyfxmhx zj4sLSN=jza$5SUTn&}UWd5H9J(zBe@UwJKG>a(lk7L_)EvnPJd`_%(Ipv>J1TX5bRDwx zV01|yIuzV(-T8WCn&k6b+n9&uHns5cBk6eY@>)mmcVY{pLsIeM^!+EnUn{3$%dNrR zK1xhJBs25}Qd;2erNfL4$sc(>^3`+5FYmX%cn&G-eL5Wf=e@dLf3NPRo1=d)g5M)+gz3b#U0C8a=E|?*S!}!wIe%Q_nMU`?!0a9ob!PXTJ{edC6OT)bC(v_b z#>>LQ>tm8SxXn_agWGKWS-8J2VjXrlUSE(0n=jT!czt?5``&p+^m}&1^mF1p$TGov z^XHbmCFe{Rp)k`rY*W8FU6CqP-k4#Nr%Kg+mX>?JR-IaP)3VBJp@v|-d7$j^4X5W$ zrZAuO4D*`cGfFVu{LszXaC%-gCZ2D;Cb!=B;j*z5rjo>|exXA&1@p~2Oz-2A47!*K zPR0_B)F^8H#seYreD)4wKd93EF?SdpHO_c66j#2Y*^*p zku!>NPtfzFb=(q~?$>}0=G^?n+b6zLrz)4Nonpn)63jQx6&KB<^TCX0O^bVfc_igO zraq6>S@ucNg3br?M~;Vmzpj%RJs;MceKd^U1819n>;v>XSdy<6X&P8{n@2VHEPc7B zw+Y3x%_&z${ox;bSW~LO%sg0<4&e^h75h20ng5FC!wQ@K%IFH8|H|aCPU!qsX55XG z{D>wcPjTV%U*`|u_=`5tNuF=wjX%KWDUAXB5?S)Q3D; z*I>(mmVD@snfl4~vKs7HCfly{f77B#O;;9w;SD;^3!QFG0Ub<-k;NB=bQnd+Er`#c zbuNj9MbkQ%?rY^5){$xIdvM$Or3T;`K=q{jo2Ce8}x2`_py$ck|1? zoR>t;yQKFgD;p=LulJ*W&QJKfg!y-#pT74xq6b@_FJ~}+T_;$9{nU;-ko)^c4u$#G z-R3ZlBMMaaDIa*seU+)8$v#_##i~+TBR^TZy{|?c(AZvl%u5~q|1XhtV~Z38jM(hF zB59vAHREbWnHYs3RBLq4oWUA00!kk9zd^SXC3m-2R?7+)(V6mj(OI36%)OaH>poB2 z%h$z}cyP$;zNIvE!k}hnH0Z9`VO&P{$4qJNVCY|}N^v6Gla8sXQ$w!m9jA37#26h+ zGm+x~HshqJ5hpa>&^lUSSDw*2m_4%Wu4%terC#)?a`Q7&r?x4e8??@qEJg>@Qa13| z`ugG2fJB)=^meLunf{X2!Tg!y9bKnypOFb%zx_iAF5M3=cNvzmTQwk@)71l zZou|9a~bMIt=A*VVmWHG#FkGhdyb=Y*JxLUAr&erq+pKbeHF^lGx15#I%NSV{+rrY{)}=mkzB{c#YA)od0Qw`jXRf zlz-w_zK&7sEWQrrPw!9c1-Y=}-mXMvy*cw+5N6)ml2_q`Z+53#WCfIsH;ypusz#-kmF3a8j=793X3$w_F|$$`DkLbkfY!~a)zhc@ zW7=+z?e(}!g&O4%aJ_??8m0Z%#(xs|kzY$KzC=~TEh@1=tbHK>o~2U-^Wj+l%jzKDN{(9r!uZ>h(y=Wpft!-lqpV>k|rvZW(iGniPEGL6+$!$ z$<#?G^N=P<=v{lCy*#ztb5-{~&-1?T=bq1pueHuuzx6x6fBT$$*4e8$`t!{74aC$H zl*H!f&y!!_>gFtF>1NMd;py$JCFbSneTb*o^@U#vn*u_(PJs@xULl4ew&e(-eU;G&EdbXz94# z(aqV;(Z$|-rL&!rs}K7xcsFAR{Ba=%QMb1Azlng>;lv-Ehtqhx|3~%E-u+*m)85cIfB;xs}!Vl z1c|>=H#d79rk&$jZ*ON0=J$1@)BFqg5O{V(@i8C3i9agMX*`Gibh^KnoaW&$zK>r{ z<5B;Q>fpZq%IhUnH!l}gJ9{rLH&H9K1bbc+l^|MsI!Q~q_m>@Xu>$NYl zoAacI{%EJOb1|&k`a=^O&xtc5hMJH% z{gL~&4&@K=XX+;g>+|jFrMtbOqo=nW{pRPn)`R)qdfjaMH@zPJ@7MjA`c3-d`uQ-O zz3kZUdie48zyJBervz{RcYXeq_dTYoJ3Q(9qvMX93i*4SAP9R`;;;9!O!cbHNgeD= zaH5b(UPXw%C($3l=Or5w+nGdvTv~3uLrBS=z}a)zo$AM#bI|cnfNLPpA77-Lf8uSm zmEgqe__$}%t=ozI7#0!zygXzFfk@Wv)SsMZqR$BAsXzG{!rE?umz_O+e8D^DTa0T> zDetf=;1`)rkK@E2)#Wt4w;%b%*^lqu@qf{0zdE?DU#gQM=Q=U&INTgko|KN*mBDU* z#fa|5`T^YrJ?zH=`f-T8y7+jk-K)b$t|EhJ;=y@r%ybJ7f(BR(R-31aKyizB}<6N8AD83i^6Sidg!&#xHB<+&i z1p-rNxvF$^fptb~e>i`!PCh?cr?)p2L{!wi0V1Cq!bvjUvgy|kJATs~naHRIT}c@s z4?F7mP;^Oc*;lBi7EpbR)UEJdMd{;Qc;m{HK||_6?Awpx4NynEDmIbS!}+|B-H+3x z(UsE0$zSKP^y>I}FvZSLd-vx$(45$InAF8N@w8bX=`Uy%+aJzX!Y1Vtwd+BgxA6YF z&N`5u#P)%+zdC*UJ~fzp-D59Haa`B<$#oZeU8i?_2w#`^N!|W*==D|hd+oQs`Xva zd1B_dvOM4vU-VRCMKU->-Qqzc#@QB9dtH2rK;9$r#sc>-9@WTDIvrR5<|T#mIbTi# zyXf&D7;DYdzD6xC1ilOSTVAWOz%ojZGwtd75f9ROhcE`+O7|+5lMhsKjx8H6UJmkU zf5qt)BwI+6BI8N*8y`2o4aqBjbN|Pp~Y&WTY#Hr(xk@ZzI+_D_?qyf2H+6a zS{CxU7C2*!()&Dhb6o?d9jJaMGpHGO(mFVqi76I)&NTqdx_}hZ;dS6kEL*RO13T${ zL=Rfm9i!JHt;X&djX+`0*SR+%>VOmN&k>`~v-x{ou4)2kJkmpa!V2H+nbQcij#=6& zH>DnIqvzj>k*w2?ey=hXKGw_Ln>o+h)5+Q1*V9djY$@4Uuk+}w8+nW}7?wINSwX&^ zzgp--CRAp4OW39(hSTjRn#XA#qZ#o%0o?}l{eT{in2zhAd63RYwjMyF=k3*{AJ1NA z94mN>vJaWIvH76Ic)?poX+MaVZ$SI+^>>h-5BWY#k4HrQkRJAYz&?8X2~3Ad_0BDA zM6zxf0-V3kBtKK}!5{K{O9=k_QQhFjVXAKHppG8>y2yzWM;p{vhgN}#g0)9ntjfRu zp1b#RUgv>ujVx0qmt-(MEP+cSD!va>j%Il+axMZTW*WrgTdY3ZF}(z_iJP!H$u$G$tsew=o8v;Fj-E>7Kn zBSMOe^MUF?E511XaG8Y@&x zV9v6(+}B8dIA52|vj~Pc%UawH4g6XH1Rp=A{NY4-a$+$;0tD&&S34^~B|ZR;y40QU~X6^7^_A&qu)irPuG?j8bsAy_)ie zlT*H^Uy=V_}RGYiaR?DRb ztlzrq#*uVZAL{yVshJFO4(SFfbzsi-(_K#_LLHp5USt^++Z2M$VFvp$rC4D0jLC|J zp&rh(l(SSnPU}cU=vAnT^G5vOVFt7ELFZNR#crd@ft^*3c{tR?xn4&5Y>$J@ZC1 z=?|xnQFv#qC7eAW@+#RU6^y+eNCPkvog?~^{`d1Di0{T z%l(!Fb8Z?11Rdb~I2XHB7(8|@0x<_tr}*As_2J_1)+eKb3&3(eZuO$1G*BpbXe+6M zGg4jR+Q>zPVBp(n1z9pIAUKXEi3#;^D(mc``ff%(}7dL!}X+8*f3{Gqm zEC;DAvJ2NiU7VxZ4bO~(IlFoq2b2DAF8VxjZ>ZNJ@Or#QVF-6ASQ<2n@`006zNp{u zXk6shJyjs`QRy2W(=w3fY8HFpWgfWRER^X4bH1yvW8^?MKTf9tiI=&ni$GlZ!}R z3)zpey>w&;57fmOR3j0^Jv$$mJxetf8e0xlY|oY|gyV4vevA}xiOmJu15ZvN{o(9V zG#ivOR`gnd)WNxN zV9>iwOA5icj%$l7G;Ny6amge>geii~0{7 zz9M!sqzY&n@OAPVmVsbd?t;V5^FT_pVo(_Tyt_U2&8XpUew+a@TwPzhia_numA`$y z&+5bH2W}VZY$*T(PCI9I!q2;RfH-_^7xEh_{a2M^qNp~M223VTj>Ks}sR zlSHD)ew_EK(_%?~IE|DqUJp0V2aW^U!y6~U^SY=>Uq8X|IJNc|4LyG*7i{u<=}G#- zx#G9w&!d@-z!!}tMo0Nd!R4tzls}v(Pfl!f>0zvZ>k2HFgw4XB47~3)=Nban70wv= zgpPqZ5Be%A1eJewqQE52<&w|WtICF`p3Ay)30(ATvyotH1L5e%sGzv;l4N2 z!P%B@AXd-55PaM+vFwc+3xobH);b%nf+c*-BnkXhE@0dQSm@zRhL8Um%juO)=?fs<3ds9(Bh z!{=c;t3cHo#^~R*%7E4;o?+6(dEnjHdrd*k$v}2;^rE}r@VpYEb^fxmm3~FwY5c(@ zBMMo4c+u~|ljw~FAo4|Rw&A5T5LtOdfz-jdJ7LyPF{eVXPxzEeuqF${=Zmcx3H5O9 zH$JnP?8kZ4pK&D%>f)?hG^2peG9U2zHHdUbmjf=t<{o}H9_O}f)49&kxxg#zg&gS* zr*5)@VC1?-AgA)$EQ`^l;G8DgA5N4fClZ>&FYgJi0%0bz`LC##0k@;=9~%nuz|{}> zO|PAjf!qR3BQ7{U&O=|Wbbi=S1R{+W6y+DQ`j9dELu~8%0&rYfXt>42G_WJ#Qx`YX z!5L$`b&Ac(LhxLoYt>0T7U14eBUT9YaLRqQ5h43=9$a|x0jY~K?-_Hjh;=^LM7*i6 zlq(0v1Zqv5!SOi1q+AFcaXJ^6%|8B_)WLZrfBVfY|3^S{jS7F|*iuj$bf5Bv^Jmsy z$m?4F@N>L_eeZRxoYpl_|3B-Sf8YA%k6uqcfkWM#9AV*{nDAcC*ClXR|L=F*{O?+S z|3Mx4Iz8&CeD|v#|LA*w{@CC1?QeV^e?N2FeYfuE_-E^$|9a}$W9PoJeQpL2b?CM( z){pL=LHFyknU1K5wF&xhi1huv^B_HSdK@Csq321bAD5obz4r$9@2%R1rCjBfyDD0> zksn0lpJqBD{e02^Js;}VqkSTx@yH+Yfr$JdpIWqDb0;&?U#>zM@t4;z>o- z?-@9;z~_wP?CCWiAS=&E;aMFRbklOnXrD%qAKT{NF|q}iD$MQ@fZsFV>|9Z8Q6Sm? zjDhm`+M3!v+;OB{<3ddX_c{!agsqE{%VEL1QmsZH71kxL8dV3{E99y1I9C_@xAVF-0ioErls}xN%H4$J zy+$x~gWlsO(1D_07Aq8{x-87Bz zhqGo!({xpsv$nM4t{9l}&{DQPoG4FD9Ij^N;SO`wKfdbr5}32zq(w$^V9pDlnqHd( zb6(o1?>H6aj8l5_ojFl3=iEc1V{XBmac;6-;a>}L=GPCjc8580HJMO4IDINIszhPV z%lv%GH^ZD&-2^E;oXZx!jV9;8XqVn2e?#kppi;xm1i z^R|zsls}v{d&(8Yz?>@s9M|lIIZtqB`@_j8U(|ojVP<;}%=uhS?aTm}v(N3*No!%w zJON;>Cd^r0^D0XR=8RMIYlNyF%-OkQ!1R;w`QyB%f8@pR+6JI}XuVn}%vshxh0?*P zXc0dA^Z5qw=Ir%G2bi;7C0h^Yn}-z|j(w?BCgv+X%9Usa$u3IUlZN`@@Oy+uoF`S-x;PuWWa9N<&L_BX?%aeq zPx;-M8jo{EMV03Xn6pBb2;~na??=n1>o8~I)1{p{FlW|owm+Pl@HPUYK)|&@ovPm^04pqjYk%!JJii z#vN&cIX^H+rF3xG9D7x|9p*ereg51bFz3;_M=3p=3kNZ&ewFZ`^SHheJt-DG2!`TwEljQ?*u1c%Q9SCz4)@J*|iSpz`;g+`# z81T9VV~bM?V(_{K`M|@I33y#Y;bCy_qiS#r=8Q9Ra{HFslNvxV7tcpNm^04( zY7#V|ZP|U{@|m59f{UmmkP^a4x*S z*2O6qIgalF%z2yDfC6KfbJrwpN*8Ck&(aDdcwK{XfCl9c=kifiyb>^HA8GZ*40v5b z$x60AoSgDS{qNo%&e4K7Hy;uo>ke~XIPmo4lQ3t#qI=doFy~Lruv;*N+UiBsg<3%(=!(n9{*{y8R;eb9ntoa^Xj(V3_mB<4-9) zoF5VsseYW9v@Xu>)B(m8FlYZu0vTW6^&{ivJ)_3slsvt@A{6F)Hgu21xZIpY-Cbyq49=3LIWFt-fm>^FBTrGv99F!a(InDe^Q>niFn zXUVC4lpapc2z#m@r_4OIF3xAk0f<~)O0{#Xd+oGZolhm%vjs6V81Q_?+{vr=HWt!zmhC}~g(*#&d1HyG+44s&kw zX(tq~_vL7Vdx!T;Yyj8hRc#l6Ipd7e{=_F(-vC;4gllwzn}Owp50nnhj*kJCI};kf zFu>(kHmnXTH?5}haPH^~p!#v9US#Xy%;TEnvJ&Q8)}<(F1amIA!`8()ZNrkM%+*bx zwBtSH4`Ls5^*~@2TL&k~lM{pPcn6DU)_~;s=R?C`&US%QPLK9& z1i9nFm%WEM>sCl{dBB`;s>DZ*83}Wq5Gt&8<5g`RPP$=zDgx$wy-Ma*BFuTIsU)R? zbAJ9&Q(u_#t=Rj{J}_t9z2cM}&hGbf)ya8qo}SIt#n~C4y0c<>BY1U&kydi74#+i6 zq;zr4Pb5OSV9uikN>KiAjyBsj!yM+E#Zum=4Rdaw{o(wX^_TC~J-_`N|INB5r|Se! z|3B9Y{(aXA{?fYVH+B1wbx(JzH1+MTodBX6-S(h~+8mfp_ovh2)9G=j-K&f1V|smb zR~ za|HlPm8e~?z^LW!jNX0w^2d2k=s+DrNh7t-@`}=gbCKl@uAYDG+pmwAq)qPEM@0Mf zJKDnp$axUazWq)$kpxl~G3#MtZLX;dVHYqVrCL^k*xRHtoAhb;LiFBJI3DrE+v{TL zuho0^?JHd@qWmEqS2^x{C}Ij>l+AEnXDUgcefu3R=IN6D5P#-%a^&~*vLp8q`_N0$ z==yiM36T*&?jOeae!h;m<4{YYd9&@i!vYHlc6&MDwc1RyYnT%;a1DF^FizJn*Sqqg z+&JzVzFe7g?INEyt-~a?(EEpRvgch#oa3urTs5NazG1YF*zyW{|1i!ehuHgvaavDc zY#GA0g<$*Oyl2hNBQ4_ABD z8GoMI%ULhNypx`*5`5g!arW>y^B&bUgVY)Pc7`^ogILydSaqh7G%Od+PbNf8A-Y!EpWHWFl1OnnPWeN} z7veA1i~KU7PrcZES3I83_Rli`=J59$${mZf@TnSO2_!{U_ePe)N1}7+t4tcYuFiFZclZRwoycyG)Wn zAU-YsaDaJImf+L;(xNz0kqDWwdE2B|WkOf7+rFbhm8iaU!)tz)8qp9FVijtl(90N^ zwR^Wph!dxR+oV$hg^5#>1#^6NOZT#IX_;`AwlZ;LhlvMarb=W+#A}c`#lyBxI*5*+ zY}+S(5hs@4j5#`M(gfmAL%boWw;=y+a(s<=FWq0wN+sKCxbNqY^B`(Hc%Vq?&Xi~C zB5Ec*pP`v2LD&tR5fraDkvKBYR*}@b?R`6!9FJ)A>5;;TWh%t5y7O~M9sSh-lt0As zMbXnPM@teq#wy855~9SzR<^(KyBjEfh`;hWNY?9*-?vb;7z|V=#DD7bU#?d_?^^@* zIrv=)eFzXWii3B@EteqL>H=On6v_}%W4CpC4woYgef3@(07^s?Z)a1hq%x5f9F^bF zrr68vPxr^>mro&%ek{n{^<^^Q6+7*a-DRm>+F0(L(=uZkv4^$lS(1$+vGqyYWm4xv z_Yq15vC^+Is{EoX;lF=OarSH};=KZ(^hP~ikjPvhMIbumCS4@kj=5z=$axUk4vjoV z>iX5PbrIJKG)b0hm`Vgho3tp-l_qw+O`vp-xEfME5zlE=Uhv!`PZ$b)$Rl;mEIdg0 zLrglKu~o=m8Zqysp=^Vh3^8aX+uw*I{**t&Up~I*6M0UrhX)hgP3KsD0*LqgE@!Aj z%Mgo)j1p-bEl*^g+#L}Rr$jJAc-Pxhst_U*vv*w0RwdN6AN7duQs`xvb^V8gFiArE zu*iAgRieajtq#*W^35Q#y!7OAoLVJSG#G z67S`#2Ey~Ux(n%~-qZDKZR%PFY3j%u%~T`9oygJ#zv67i~n; z_*G#%gCvM+9c+J%(~nX95XpLR+8;^2UwwNXPrllHdp=KItM|vRi~pR@|0mAp|5Dxh id41$?9^2CnIChCnIC>q)l0N8sVI{*Lx diff --git a/tests/ref_outputs/test029/sim_output/output_data.mat b/tests/ref_outputs/test029/sim_output/output_data.mat index 554aaa5664635d85c5d797e6f3180fc550d2507d..cb6fc8fc678ce88331b1b2197004c15ef96fa491 100644 GIT binary patch literal 12160 zcmeI23piBiAIGUpq*O{2F{ni%l6A>;POGa+xm51P7%FRAb__}>x}lp=D`~rJr6kt1 z+9)YP`bSDPQ7XApR9jO1zjJ1er(R7=JiGfm+dR+A$9ZSY`JVaw=Dgl>-eW#+-Vjp24|@c+slyV>*MFKVlvH`?atu1yU~2S4QV#O1wH{Znhs4{ z&v2@Np{_n{s@7C(p|5xd&o`HV)eZ>@35ftj;YpBq41y=I3~vcd;!zH0Y;3&5*n+u| z$#SPNmowPQ-05x}e%!y{Z6PsuS`aIdkLTI_M2xo`LD|=QLG!53_j@i;r6nIEDv`#Vqf~ZPFIX13N|Nb zR@Pa|!=Z+(;P2yE%G#1@@CuoMCT(-pQir5}^Tu|ytn^>-v<=`mmn%&6u zo;bNYJloGR7#EbW%-hi!PU~@$CB;w3_MTXp_sTqFP&kf~s89G8#&y)#n;23PV{4#bsP^uO z+?KxP3R;h)OW6!2bXe5!p&2&j-u1s8V*GgR(!&^*el<2MLUzg@w`(x#g*JIFIQ3YR z)GmiF`DY}Q93M`4O$0p?XR&DF3!K&M*nVLFR zV+W_bO&7iRxGR^+?};MmgG~!x)M5*RjH|?VG<4;n-*@TaTj_3&4H!|6-d-1whR^yk z*qlXdA2)YqfDcQHSkj_ft?=%u5Bj??lD2NOFE_;~PDhr16IHu22509uxbJ;Uo2-1`D{f(1T5{#*9U`Iy{2o z$Ft-xxqXz%*=xvlly4qPnLh4096!PAA=9@vU<-7(xhT`TH4=@S8!&zzsrHIzKhJgq zpL<_g@16}ljR!JZwZP}V4=FL?l1-RYtKN$$@QHGw%jELs_v%KCiA4ZFbSu{X1> zynkPV1+NL+?E*d(qY|udf=^wo3*`1uipuRH+fg3I8?kBNGsbj;k`4HDzQ@f)X_C;P zN$Z8rhCzpOyufF^k4dQ`_?#bf$JiWvvJB}N^59d^`OMV@@QG5s=u-YI@Hs9z48r>rz|oTW#kr> zjRD}({Eh1PXW&zI@T2Sn;L}W|&GZWREVC9Nb5Wj{`Gm|x`KoBl>;>R+Q;m^sPPTgpiu*bAS9w=fGJ*()_{a>*H4b2H^ARv8r8f!6(YMrZ-k5 zg3l*!z3pVc=fofh{YdbcF=O9YG4ScMCdWV;yrGQnI6$_eOl9eqJ_n!LPl{CS!Ds#G z8ZsB<(J6Nu&B5pI6~*s!!KZLRMxieFT>r*EA{~6rcgZBTkMf>L4B3v-V9bs*HuyYE zAGtLae5y?5=Ayi;ysD_&spsQe`_Lc9VW=HTTapY3D zIE+6~wtP^sSqVN@82#@-%k&yd{jW8tzTk8CJ6!H1j6W(1688)MpD2ZMHk0irr=In^ zrJP!aJ=I|8hl9`CKKW!WN;yaNC>ijXb)!IHFZiq|3{M;lK0^oTKG+97#Uf6S+ef)0 zWFy&*a$V*F(IenNVdK2etEZn>ieK0{6@37rR@PsUr?n}E-=HUlr*1D}V(6;9Q{cn)Qx zRYwBRj&kiI&1g69Df~mN`=b5XXxN)jF2 z3!ipz0Zc>i87cp)g9ScqTN$Yb!DsujYgQk?CwoVfC>?yFEJ}{7+X6mwE9ZykfzMS} zFQiq&`9{`+o8h~`r-P!z(WBrEr7u={glI=8wxRN{1e|YdF9@ly1D~1pTghCMB1tcc zmB44|utK359dLdmXtwSzD46G|ayq`}OkDb93PQ&CaaNb_4iC zIlT3!d5Yk3nF{^%Q#gOyx+Z$=aq#)pvt_>r`1I0m_X`4VDD782CEHQ5=08$i0X{eC ztbBO}d=3?DAahZ+EmS|24(D&PrPb5}!ROg68k&aSGiJ59PDl1*Oyv|e7iC%tHy7o* zv9o482A`Tn!FRuh^EVLYEN1fM2zxw$BFFFvcg?$Gn^vtr7g`0)nSSeukm zOAPpYA0~NVR6sqpE^f1P8~8N$_O`45pD4Zdx?EHOpO>AF%PD}*;uQ8AN$`2L+H}f5 z@LA^IcUcO2q6{2=ESG3UsoOSJnUh|JHI*N8j{%?7%AsT~N>+sV!Nk4w*qIE^oI~$w zu=ut&n`6Ldrb3m>Ch&Q@Gmy+h`RkQR7G z0iS|gUm)%~`0V+sE2Hba1AfvX>)6EgAmnxfzSrl3BV7ZNw-;d)WlQ&5 zv(Zw*W53f1lsuqlus zxNagoiIa3GK1u;!C6`=TSu?!r{(`9WzuYJ8E1*c-Um#;WU#3%eBu?>2Manj@;8fX zYP-o?iq!oDhbOVfT#BX1Kl>gFQNu;YKDRLq(7>tt3sO4+$y|!@7Q6J0Hmc)-^lu|S ztqYx)ANAw?R;cUSyhXspSYlojW!qvSGv*KryfMVSb@y;voV)CT`+AxxbPRRlJ0I#1 z?UAP;$>M$97f*R~IpYs;HgCP|^P$xIsA>y$UKC}5_+y8v5kWY2yQ}%mhjRBX1hju!N5+NKFCN-`K9rgtJb>-C!4)b6hW=Jz^} z+Q06JU-kXFA#-38j=CCv*uZhV*&mtSx$YR^(BsX^c@cm@|fol=PAwQV}ZB>77W? zgGzE8^`)1SR7yDBjya;6PPfPX@3FV9mUL#SbMC$8Yk%L&&t9|l`mg!zHEaL(S{$wB zI9N?lpQO)FceI)_!N-@&QJ=+S348+i{^shuz#!k{^Xpgxl= zP@ib1Zp1J*HZz~dfIo~3r9NYmwslm5(;}p$r1}D7;YE_%WZ*^a!$(qMDUJeKTU)zX z&tk7&b2&`*5*B|khw1Gb)cF^DEF}jo26>WbYt8#gL{vSJ^8Y42c^#IEc^y>Y*blCQ z=xPK6`+ITtn14xqedXh&h|cTw-=ZuP|Bv}4{<*9m0h1lf=W_xCANr@h{TCo$_T@?I zYP`{BGLJg1i;_0wGviYINc?DDG7gFS^OpEBSv(%sm*Y)F{o2kN2`Q zTStgX(*KfA ztqK|TKn7J#^Jz*}P{z^*MEJDIwaGQ3sPW_X8y{^!)l}RXj2T08r3N=NB7dhA-#fPT zNGmGy=dCAdP(^;QrM*Hm5>*fJlc+c}tKJxo@%O4jWB$oIpCQ$V0;1EL+LD`5)~tpa zV@n|-eT(~CI}DEHvG{^H{6KFGJ0y_HAomzd=j8!i{_&nRd6-MvymEU)(*oPKr1s4y z()-_c?&r~@%JV9Rit5Xw>g)O(QOA?<87!H}G$%~l{npGtW-=X@rWO^Csvps$wv_{+=|c;}Bgh9iQGO{d*|%JikMZE%NXp#~BrPeat`HPDTC)LOps^ z|L=*@#swAGnbe}+4qkNCFRe!p5>(d=4r)U8rpbl99Mpm$udE)`$`N*>tNwJYwpC4t z7@%tIQdy6TWu~sqw5ml;t~?#JAFI*T%YjnXNmbp*IO!lN$9bF~H*}a^1M2^{`Jl3F zGg>c`%bR_q1(`0$(j?>Kth7nw*ETjHuZ!mvTG-bkyKV77DQPvx%4*#uS;cC!x^<~2 z4$d19IZm3aI+U7hpa1r41Bx-7Z>zbl8JXC&lpdWdL_K+;%b(QwP4xs{iO*ejO$A@p zy;)ui@YVZmW`ewO3sO_dvgv^K;*>s>Z&C`r;uG?oSb(q8JcH^K@O381W1%wm8hz?~ zfhqXHIcHCts2r!uUBes$@Rg$fdy)b8@^{iQWP-18Q>HqQ@o_%0KWcp)e63g?6{rQi zoZUt?Yy)39F*j}Bz<4o2m?esX^Jgo*jw$$(TU1=F`mO;5EqyzQ0AGug56sLI!&kFR zYQ8`CDhjkK@dRJCp=H*N;EQX{%pMKCv=Vo6-$HwFu95e&I}g5&pDkRY3BFFA^_#j4 zd~L|_+a&|O&a`W9)CXTUYb*PR%5h!{F=5UCU)rIi+a17H>h^1`T<|r}WnB;%ALkw` zX_b@UEALENgcA5#F57p(D)99WQ|tTBz?ZU2aY>(v;OpNLG{<;? zuMnN4YZt`ub^X4k=?3t1xozXTaPZaJrvI5>@b%(|bFdls+P;4`Pu54+ja%eoZ^VEv zznQCU^#)%fiu3mcgRhn0CMVnBew6qu)kqzD;WR2=B`U|+X6aSs2fo^it40QcFID^d z`VrvEJuf|kjE{4mf=udq@O8{0thm*#9(AyCjTeJ2gG1a~Rp9F=H(3-1=bjNMHy4Ai z2b>3c+`*TI^dZx=;Hx4`<4~a(z6|$g961cW#&4SuaTI+0EBm(F3Gh|$N=@hqzGAbk zoL7PN;7m<_&0JT6qe2EJMxwO7pmU(%lHn~TAh z!QM~lgTUN$dNTW`eJ(@~5IWI1i=E_sRud({iqA7lNjfmEdbqnEe(qKF&T- zceSrTei%ESrQr^~Dj6ZJ?T{ZXb7by}1z$H?`-tM>d^6$l?`_~quuDmD#gmW!@ho`3Ie#cmRF3mh(Fzw|@TGP-;%*A~ z5^j+$QhKl>tn^B&*@cc0^=iXLZy$J9%SXZm)8u*IaAoxxpr5T+#P$Fk3hOc8*u_@Qc{&iH( zPPeQc4K*x{-UPlJ8heE%fUk-uUhQV^g_A41=NbmSUf)hsb^~8goCRry;H$wUeW5J) zI7cxH1aUGlctp;BO zyB^M10KV*kvz1K1S5fq44JGi^ma$6|2WQEay{k*X*Ivyf+Gf&?XiUZr57&V&=gX&T zSz`ElVzB$zcV@MyU1iK~3E=C^M&;zeAx$XqaP-1<@O9`!`m>YJUYvz@BD9``H6b?B zWl!CcdX%45koE-CqO3O#E;mxEkj({(nla|3lhY~@A-^+)SPVdeZ%n@H0*9%*n$V?#zVc#qt=slAuRxc7n_uWY zzaZ}aFO2NLzc(8$fVTiX?Rnl<(s?HAc!}q=M2-2#8&)q8b6?{4V%FOgpZp#0J~2=K z{yZi*-+(RG_Bg*I9pkm^)V3Rl=E#!W#qNwTL{cBdGUa4#;$qNfm-%<~i0LM(>LIF( zuJzfz(Py@*Y%nBfmQO&FWY3Kw6!Zg3lG4T!uj7hWj&#=`nhvIKYo0r#Ykl^)cD`u6 zlqS7CyMsC0R^X^W&`gYUv)_I|hbUQFsd6`MJh7}$>2_EpqicP3`&q#{GCs{cJlF8J z8eL-e%k1S#lC_A`OAZmTPt*whyiJ*X7Y!uZo-4+Y>%{X9R7@h{(Codi_0j$7n#8;8 z4N)f&bcwX-#Tt{hP9U03rCq2WU_?kV{tWuiJCj_W3ZnWLJ-ryz|T_-!R z`_CeNjhNTAOUZ-i++RfSmfEYm9^l>MdX#UxZ{Y-GE_FPQ$gC{L&XDI*=XGC?qSvR6 zQtMJUPsmr#uNoLibk@6)(9w9NA3C__^(g0&QfuVo*HCdIi9v@qHn|R1OU3VbJxc#Y zNz@F*4aC;m_)9{?jnwy!A{?fT*0vb9N#y@i{hpq?C;NYYc#-|yE10{KT%#hl?%q4g zH+p#FBt4spK9hOWd1A9OYWMfQD!!y~$bSE+b;0+?8*-hB^v>pRxjp%&#}}I4#+mNp sjHLa%b7iozPkr