From 06c29fed2cb4e817310359ae1494a02413cdbb07 Mon Sep 17 00:00:00 2001 From: Jan Benda Date: Tue, 25 Mar 2025 11:30:42 +0100 Subject: [PATCH] removed absolute values from motivation --- motivation.pdf | Bin 198014 -> 198710 bytes motivation.py | 4 ++-- susceptibility1.tex | 42 ++++++++++++++++++++---------------------- 3 files changed, 22 insertions(+), 24 deletions(-) diff --git a/motivation.pdf b/motivation.pdf index f2b2e4a098502ee1180a0c6d9dca7bf1066197da..c54abd0de70a8b6adacba25bfb7237624595e8c0 100644 GIT binary patch delta 20364 zcmX7vbzD_V6URkRx|B|@gmiazBi-HIomW~3=`LxwbX>ZTED+rJ+>zCQ|bvK*@%J1fi@k{zlm%$9)v&9yH_uYqm7@4d!~pRCY~JIO_&J@}BPqV)jn*dymP0>crhu zBU6~5h!kh@`O1GT7~uY_O2Rk-7%#tH;))IebCQJFByntNC2@7~K4OaasYsuJeE34S zQrlUh8qe^$Gl*)VnCd-&U+fymUI)!l8xFzD262*z8s5wORpxy)d@k68n7I`<8ZgMe z7)K=&1XymFm=Z3W`Pwd{tRHjzRMI~T`8`yYxvS+&azO9YKlGm`ik2#P`qwz#tCIp~ zGHr;W_(@d8BkV%*30WnN_fGNRR9g&C1Ae}aG-FOu5pDnVJKvj&^X@Q2?D<0f{i`8c z4>kFg6d=TmgbDQd&u8sFUlnF~L4;r3!VIF)o)#2EzVJj57A@LsKgM4_#S`&KEMx6m z^tdGxGXpQ$WGs5}1-HIj^rw^xUzn)`@&>*M?$7*b%K7L&-O~Rja~J0=pYMCaZq4+H zpocf;Vc$wfKG4|_!JrAiYVI}1(bV0)LC=OMj>RqK7kJA@^p+`5zj?|Fj_Vdhk@=U^ zM;F}HHC96r1UzlOJy_sFbafV?AFzDCg^9jt#+?HfXLWQ;TA2G0?|oL8O}j4D0v=g^uYvMWT)bpJzc)Z|i9cniy%Ekl z5tMMgpNEWN0FzLX*gG2uLPG~x-a(!C^&|R42mah2+WhtbS_!jB-Nn7W7+jy0$ma$Y zh2U<*Jb}M<%1Ulf7od}qk}@VO^m z@dU5o5MkDAKz~fdkfu`x?5}TXAD0#hcg@!AXQy@_WCi$Zi7f&&bzB0TQe>m;LJ1L8^L<5ndovb%rd4Ox_5d zVmK!|62H#aK#OMEAHa_o>%}v`%O(TzjsmYAwXFPJ5xniZ3d$#Q0#mady-pXtwXI>T zHkfq%j0-FjZMW9@krD$~HU>V7-1>RlqTJ?A6h)Jno=zNqU$sNNJ0@P-kbnhXRlQ0}(M!-zPw`_^1&d3nf}(*a%;lkgg!uqI z4mG6y{$Pay%$zVt$lj?J1DHAbh&>HoKX?++KRe|^mH)X-(`{d#Uo{$tJ`5U;BqXtn zS4YL*P1J|uZl7Uv2Ob=1SVofkgh39~+dr!TrU&D;Tx*a{5)8O_K#bXus`nId3W{)5 z(R9BS@HmD}G>{_H^ATmxQm$NS0~VO}xH$I#`b#IETS_w>8hek|;8P@KZ1aCS1R|}N zpm^q-2ontrX6mp!%lLuk6fd{NXTwR01RW$l-IDXGU>j-~2$x%rDH(jFd5qnRM z(fPK!leJy?#;$-Q%-ziCcLly7B`&CKHWZBj9|nH2G-jAdD+VN!t8e$k#Ws+~r{$e+ zi~`2MLUQL<;7ix;15A0EjndE!=GrRjr;$cJKNh2j;JF>)EMH0VeN$}+C1;(`0P#BV zsb3_c;unLCMh3Z66R%IEzqL{v#~Ib&{ddK;Dt;YAqP8I|#ibYHUN(raFaEo(P_gNy z**Xf!E#O2E4^r#k4r9*I?UlD-E+``Uko>F#YAhN~e}-Pi}U( z|7mIJ$!28hiihg5Ts=1*sVEbyPf^n2v~uOw^Up7=Jbn&@1jltNXV}JaSZg`FX3cXh>%fsZQ$H z-%YCIxg5Iss?R^6dCuENYb-!QzPGK(cdA%7?UYOY=y;cl)X75euO=QJEc>%*1*C4K zJ{XYiV@&!()ZlOAKkDLQEc!h&^19iwNLXexnxn)MDubWXZl(SN@7|P%N*pt;A@V6* zrRSx&7SOAJ9EOiEcucf)ere3UWqe`@J>EWsm~k8`{&;6EZ&!|OkNplGjXrgTp^=g> zKa{g%*kzf;K7n0+l5_h90+>cKvyR5EI!LTk*mRUl$g!9~)pHkD$Sd_Z?JXkv^qZvo z9-a^gws?tL5cRY5`-|A1M?nzwe*N?b`qMT z*n4+OOgwi^I<9m`csD{A-Ksc#GSucVEI(fQy$O7|TfgNyU9RKO5PxE!8E;e3aNT-CY+$i2ZEGKrk=ao3t+^iBC8JR^h3nKz4MO5R?35Hy4a zPU_Umla4Lxf_(MXFmN92IGp?XKgS{uqWd@0J^F-6UXn8;ItBgmTQXM^xOC!v$WxD# zM56yX$v^kW}sud2|BhS;86k$U|%2;vQdV`})5@e#@k z?_FA<8@YHv0SZ`O>Xce0x+0Iw4jJA;mTSta0WVUwxLw!#{lZ(Okf){&o*5e(GdCR7 z>X1KS;aPRpmc)ssY!!#4T_YYn=YzVV?Li}(b$F5Bb?s{}HvEaptL}6W3`ql2iq-2r zDQam6ICX{f!aI~`3zSytCgOyFt7;J|jK>it=`L%jyc#Z5;|0~f z5kHuj%FBUXMYg)SNBOonQ@Q!^=mV7=;Z~G9U7+ zE(Zr!b9@A^Rl}ws5Iz{(VpJa)hh~mQ#A^~xil)Z8J5&h$$1F^53@@1l*U*K?1JBMk z)JO5{$`ar_SNUN@NYa=dVc}^e{AUgMoJ@T-HHfC7jRWY%qkplepH*G?zV@#@GE|wG zm|i8%I_A1B1d&3>)Xs*NM#InS62Au}=0OHx%N^!XzFK{>c2|>>EQ-JPk{l-L(xnM% zK8iiZdj!Qj*5Hrut=0&@)4w4Oq#S|AFFL%V_)QH?>(W%1yo*Is-YDree^Ynx-n=2~ z1{vQFG;gY>x}3)IF;(g=s;fBG+Ub+M8mt$gPrUX0+sNRFc4Dly(8M2*FW5Y$9L8?+ zM;e~SgSjUX%UM5EjhFM0{H32-W87KU@)k9)lR~p^9QLBQ;VkF-{>3-%{nSsI^!dZP z(r+mX7O(t(&Q@1HtDz^Qp&G zG1}YPnrM7K1aogJ!@fmJ7fA&-;uFcFAu5cEAGNiI-{$ULHX{nT&%^hUQD1Ju!;_c~ z60Xd%ip!l5D%Lsx%JH1ME|GU%wC=h7G&9?QH8pSIzU@=TYT0f3LqY{I42CFuvvD*wKWc1*(E5%wYDx3TEzDELj$|(A%V}^{O0KQ!;WC~ zgIlGUX$9MfDWvYrf2i0G?tZ(hkycIAC2iLzT*H6%Ylwpa%S;`=yi_NYp84&To9HlE zFC^3Z6$nt_6B{w>`IM!a*zkSJD_=QM0#7&fB6pcJwqIiry@MSc{U|`4tX^mXE206h zZ{3KzY^K!(&ZDYbB4QE|NmW<)Zi#y`HYudv}T#f|ibW(Kjoy z?P!8RPSPxpCw~xs;b5#=rJKRCvRJZ8I)WKCm@OekHJ)ea;szsG)jVI*JSDYR_x@-% zsX8Bbu=y?9?%{MzHV5+?zwvOlIPxcqCW)0M_CCmhPQWmHDRul;4<`J)V}rc22|7@b zqlLj2@9GkLFV0aj68zP;_+$?dxK&DG?e6DJnk{BDmVIgZx7pu8wH%3s<(}WjQp;vV zIa{{9pY)q4hcBn-@6)cin9)8N{6nZ++!#iQ@N;wMUfCuN`pJ?mHt!};RU@G< zcxl8h>`O29_Fvi#&XVshU^iJYzFv$!Daj+nvMXbg+JD74OtX^1*&lls9nP;78%`|O z<3`mpzHf>?k1Bv478^QzAJ$etMyDZU2G08@aTz*MxQ_kppE7N7)!aU_DME~nb20?(-j^li(dwPS-DP_sN2?xk1%0g{2E3LL`cfUp?Kyj9hxhu8MxR_|x-fSw4TC()oFBir7Y^;JB8kyQwc0h02qL74N|?Pq@fP%P|s za;%|1{jUT++vE(?`Y3x4iWy@7KfO`=no(`S4``(tnU)>T4?y9LtwhS|uL}qhNSfgd zU>@@~sG+)~NpOz|0%7zL$NhgEn=To@3+Dus_-Czn(@>2FaOK1XDE~s-a2w-QmCqU3 z`~Ch-_!EXMY*uL_6Hpmd2kh35-vHTFl|Ma={EPzV*{**CIX#}VfD@&JAv+R))<>}_ zR@}x1v^R4wHR?zpFb}M7UloU2u17)uBgY4+rg6H%;RMw&=kCfx>V{_%L+$%Kc(~b3hBd0nj8U67wZh zQzs;b5FA|GTEn98A0M>sbj9xD{cL}W`H2+S_~dMxm_#&XgKF#3E2isSAfom}X#ym1 zw%MUM{xibaEm?jkspuZ_W*;et#YmX2qYwBT5Q#S-YFV4Fh%%VIedt;aq3H}yxjMH^ zwe|M($NZ8+g+gl5b3*sm2HAw`6el z#-(|E;8%}I@eF?0-)IJ25JF7+N$y_y8U#+7F5f`_CONDsRWSO`JNyz$FDS331uQV_ z=lOB-1zS4x5bziB0L{6WEfS^vn~IouohAvsFaYE>7@gn)cE#wZ9ilkVn1X`$>!xj@ zEP!xzLQ?}$niSL=0M<8_0MN0NT0l@eFfu_ACXr002>{pjV<53^X>#kQ0L{3-J7=DE zfGr=J|KLU6ypqYUCn}Gfe9s70LjyjoGGSpc22a_7lkF(2nI+=wcsx77w2um3o}b?w z%O}y6kPibe5+m33IycmdA@fv&9Zy;Wo#mqrc zO8jF$XWex9qeO%wKu(04TK_=Z?i59+#Yj`Oe0&8L$eD_60O~gRJ|YILn6UOW26zRY zI2%17woa>yM|<6kR_8kbS9P zW**X(Q#AuMMqrWQ-kW*!UtbHR$&|pSHn$w=MFo7Nu zVn#FV!4Iz(GcZ{%di1%|1fhaeYgB0bH3#us?>*}h0V~vY4#(kN$m;pCY{pEb`w0AchyH zd#N7Q)|6t;a`0`7Wf8=MoJhZtZv|_P3AvFbFrkl@-irr-@E8!GD ze}^EwL#JPDF?6f^(iHmhdp4d^V%p|JnR7r<&`H%g^6WtM;c4lxl>yhZcFW9H@P}g~ zxvB65O_?PImC#}Nw*F$fL?}m>dY>|V!P}GgNm`1xvN9}Dow8l$wzb;tzO&8RV|B}9 zei|%NQyE;1JG@a1IO^}TwQorZR9USq6UulT@EqJfoS}WUn1H_3PJO}kP71;;}X|w;R9 zZG@n@k2D44M7jh=Y|YT5LQ`4=Q_><^7XGQrLKNNIp>vFxc6!Zx)lb^7-&}{(sd6gD z@ZbsC@Vz$boc7d;lKky2_GVA0Edu2#-vuz@@Tqa-6`7Iy?5;d`|3YujF1M8^WSx1Z zNHh0mV*ZrL!pKpzqn22LWY5el!x~KHCSA{9E)8!LZKp)=ZgHYH%e-&`k~XupYd7i3 zc_Zi8oNj7&ZS5wQ6>ur!|AU2Xsfsy;zQz>l#qCy_qOla=WG0{|z?}-2`%47D{E^q} zpIe$5Kv{_pafl1$A(Cx0VOguuKV8T&QztoF)yK9h@fzBZ!_Kb7l;II~ zlx}jg#N+H(|0`~dMOC{5=1%bVT3wN>+(G2qnsW zAN8WGdrPW=QZ#yeb`Un0WeZI*u1D_9^8R4wqYBf(iMPtm0;!vkZLW`PvM6ffkTbdU zQ$;RR4VGYyXow6n9-}H|Q*3!6>xb#Ds_Z=Ovdh(V%&nra(*s~-JyHIX{pu-_tsJgh?_)J*4hzv_HqmW|T{L!71Po{S#9kJiG z__Q%3AeRSWwAuq7wpo5d3*Gb^>smqZ*~=LBi^n}paopC3GN=~%>4a#7#cN0_!$DUC{yP9@&jRQZ-) z4#{50lzseS!|3=dg;z{Wd#x|MU4XB07F0LC;01HrsM-q_0W*+x@K^gJ7cjK{Z|Kw$^%k!6+jHF4Lbe^_aB z4TO}j4F-W{p&b7xQt8#PX9iMHk~xYBZ8=rR6+BUFfgg@sF1=e{Sh?ds3w{lJs%hV;iilhQYz-KbnCC@OpI(W) zVt1xYemwy7jo3(Lo1AX{{t%?aBKK-D9-Ip34}^X6aHK$ZJsUu#h*{9)ivTtNenMpGuwD6a1)E zz_fWMd3^-zwR_=H51)W)0f*R>Tu0ed9OysbSx2&Lo99%s=ufrNK=g;V(fxtn0+mHc zqIuo{dn?O|QC22U0Er^}#UOt%4mtxYCLrV10Jf=qnRl0)`_^ard?d(kFws`Y?D>ycxbZ5Zc4XCHFBqa+}VfhavOv`!d4@IqyiL_O&5u*n|Y7sVRB+7 zeV$c69Ui-Q%#m6b{{IHopxyS_FHtZr-#L}EfNJO&3ZN>)%C#ouUTbzk;=F*MDO12_ z@Z`qaK0Yjg4uB5917E>w%)VY?$PbfY%y>Pgch$k{zr!IR~yH{Iq__`@wyS$)x(RNZ_;mqcF(aN}r_?eEXgF878E0;6;4$ z45;`O6-8iJCH~zun8Uu1$;hfs`0`1`lb~Pn1y`Co#1^=&GK+214MAQ@wqR|irE4@m z_&R?~9*ZDjcgl-3|6WEIz7=b+3lvZch-2lkNqa!ZCAmJG)$x&N1s&!^VHbZQ@n`v{ z#Orr@$4Bxb6fAD>Ay&_C6(sg7XLlC9$p=iGx<}uII^XjqrHn?=^lohW?N7m`ABVUS zg{eVgqcQx$%0ofsyQg?Hbz_zNl80v)>q&|~|2sneOd4J5$v25IQFz2%6ORV`HP1G1 z)U^$iyQ#L-g;*tCB+(;c#VAlo`-K3>w}#Io2H*sem56>=yZu~}=8~85uJm(H6lnx( z{<&SCyv|?IoAp+_;3KgLJwz?QX=#Oxh7(HeC4%dWu1UX1fSZ)rEFkejpq4@un15~G ztbQDY4{!~uP?qxEAdalw4TUr#ZE3yy&!eZr(V*s1gk&0^WDWks&Oj=0NRP(%;SKoW zMz8rslN`{uaQ<;kgTNUq&_~^mgYZm12FHMUzN$;k18QLB3*yJ}K|wl2rPf2Cl%J5IHAwLtvV zrn>nbtCMzlEnPG16$SGuhUroggz=cfc36w&7MMnhIK`(`eAEy8Hl=z?A?rJwtdvsoakr9r_BuCmv7hBdlSdsy75KUnhrybg5q zh^Kd&oVe8(`obEg8h;zA(N$S6X7?;#NPb9NfK;zYv%U~6{gSt%(FmrPuwE$T5mz+x zvDBHghcKfyYlG;iemWxcoADu%&87)>hm1n6r`hm6^tbl^m~{W{|^N8;Q{fvX}qe%2Y|pm*Uezwo{vTYABEWfaGgj zlnz?QkULo0H!Bm%d^VXRiyBKfO2py$i=2FD2-&P9L=X9# zQGJ%vO5;O4FmV7r!~9Ep?OB+1NqJCM_o-mqsDdDQKZNS*%0zf|M#0CCYtN|#BLkm>A)v*SMJ|c`8{cj?aB)dij(zeM@gy(?GDQy6aFs$HNUQ4 zIcS%F>jsbT!LVj>FD_K+c0%jLF|p>V&yt2@!Tu>eiEAq6NVC@{adMu5cfINtv>uqg z=!CL9?a9uq3 zaz^LM(OhJYJtaG|szmswlM&2#F0uJ#4!U+shhWeqE~YZ_%w~Xt&CfbX)eH)9t2rF| z*2bjZF%02-)|HdSL){Xeb@fn-;>x7wlfhXbK9ha zvVueA9n{g9)119y+4p(dv@LVZ)dDp1LTyFsBMr~3lOc5~x|7Lg$6XZQj_n20`N=vX z<)5^O#w*9n5p9umPJ@aqU4N-O-=km2`Ic0q<05mJo+OtNvPrgNbcZUB=^W`5jw~>B zCC&~hE`pRwiC67h%0}+u4kw{s#iDTM)4>Tur87|`NWQ5rcs~Wqy%096_M#)_QW(sF z0*4**20wI=Ca_X1r>0yFJMPt591uE8I)~oDmHcYYOjds9bt%8X6@J!J@iDq2cESBi zUdWsm?$C9UCQ?d@&}P)dkCRw+XSNHan1ksL742R_CZv?ibV>X|L(`Y%_g67K|6plq z*~Uvqzf2F#P1t2~N%xL$zt=% z(RLrN?p^iY;?&T!{jY3EeqmT9n{`g|d84^uxKl^o`4nB+oDVoEjhuL?j`DP;q?sHy zaLG5n{wBeWv)vk*kQ?fN{C*OEV%kEDUWn(e5~p9QuUALi5CVE%zK%ssJ{m6? zM+YMx?-&OHk6Oxf>1|k7`iq6q9xNcIt;tmVQ@B&U!9T5dsyAC5URKND?e8Dc&v8Y* z{)1FPZcfEfofz5P=Dz~n(=;Yuf^fGbMXm%jG2q~%uhlycfg5fAmU7q2GMI6kUzhhM zJ9PO<=kQ9+E++0Or1kw)aGLZl_k!KQFRFKPGpktv4_j#7hn5uo6-SVfVHmbl1qt8R z85J;2?vL6b3qs9OFGqFCHQ*{7xd8!M3vmAR0dwlAIBj~mst zw0v;a@{V)A^iDarVxksrt7y*bt(rX`Odk7u@>TGlsO?NxWtlrFGB8vlEs&smVP~Gs z<$eDE*yKieAivFj*$i*WGzF`Tu@F0w-n{w>YPaXPcnQsXEp6vaxxjeC;iaO3ER~qa$|HS>vGFBh1bRbQz7G6P4{tk@^z-G*4Z`?zbv6y2{H}3F01}rdRG4;gN85 zKut)$qJy655Dot2{PX5Ki1+aHi-+d@7=FYXv`1>KpUVy}LN$-X(PD$;pttE1W2L$bUO*p z4xokNBC(;-C_ot|W~e?oRl|-EK>UnF-%+T|9*PkrsqUcEG8~rH45*+!0(o)YvDbU9 z72Nwl+{b4{mfX*+wC%Xjf_r`7N$R@YY@n(mul+JV6l8yR`cmq;Vc_R;6g^=#Je!o- zB3Kf~_}EopI?+)nMRyGx84O-RbGMs-d&bVC@{`&9AxpRLL}S88gzS#Q&xRT`*)IHM zF2v6WM$~3dc$KL&bP>u3hjIy!O|+m0`IDi?S=U6i-=gW+$g*`)q(m z#d!(7j1>^h2~p}k2$jH%+0MB=QmI{Tz7NxmF@Bb=z?ZOO+$;MtJs+i2;_EMCcKLuK z8TNPl=Ps`oRhJxn_!t4qEnQ%6C*^WN96eGm>rpI`|Ih33g`@Up_NXShmsB#&MHbBW zb9q}yA@rjE_QNxbGy@9mN+|hGrLcB{Zx5AfswypE=hA@7A>({~v~!U;t&ah0TU*^( z;^++j0Zq%%{jPaK-L>%pf+OO2pNmeV9loccO^^I`$%Ppj|Lnu~7aOJen+551*Io=I zb}>tNmSPbA4Jr-BLegc9e5^7k$H5cxRuV7cjn;_mBYkGBt(wwPGNrNw4k;VD&SZ4| zYTYt9ePD^C>2XlzTSCWTIDU3g`zS_|K0)d`{v;m9f5KrHDs31-z0J>N@52ALZHK7m zalP~V-1~fTdTicN4T<$0$Jqp;9(Z6-GFqpe4~ks2HTt+>J-}2?I-g@pg#t!T4OqCg zhX1TcXkpN1`_!d7)SeMy{3jwIZaRC}u&Ed7aiNwhoV{QOqjOGRgLDd04n-bo_p2pt zhSuAoy0ltHek@Ru2{U7~sglo(Dq)^N`s3~~cU`)-@WJ6b%Kj+CV^4E*3a7x!px$+y zLsPq+y%)hx%_em+vz~7l$_-|`wX2L8-eg8yWCx8m6!>tFBV^1~NbfyRB}Q%5Z?Vqe zkki7lr}yQ)8tEvtPOagpCz2`Hsf5sxjN9a-iQvt zty%Q|$MJ6E2@y0*hX{6SkarsN2|qGaoc1k;Hke%`79APVJu@ZSqt-mA7ebx7!X9Y} zZyZ#esd}%rM$KgGC(Aq*cR&jysmB(B7KS=MBAz7f@fI24jk@dKBo~_Lw3PPT91w0( z_pGn(urOyb`Z6jNyDnRAdns7cMk}TwnI2lD87FDEZZAZUhJ#TOjy4!LZtHZ__FP!n z=P3EzrzSV_>fx!VqM~aV^tz zjB)1FC5YyG8m@5p1+kiiMNK~Lsvq9km{e19+%&$S%;TFbP=NQ1V&f=Cg*mU;_$$K` zOnk7&a4>K%jW4*5)a=Qw#5L#EV8Gh|*%Vyuooi}>DFz>+^&;!kUr>TWRW=-)k_hF&$vWngL^E3lIn$dXR&&shz2)+n>DnXR$sIrjJc z$tRtUtzmrE2ujtuiWoo^IqoaCXuN+G+M1fScP}2Ue9|S^hmxY&VyLabRn9$3MOjFt zm>q1HErP}1PnZUqW4oM*A5zOt?lkwn6iUhV{pZAv&6U-A}2qF!`425WIeP!4K@9ZvnI8^ zG;n@qiN8JV{0Vbs)@I~&TdcxCOL?R)q{(EVccB%oE>y(lZ*8?4j? zf#lll1P4dfVH7gUB9oa?F8h6~*nCO4jNp8PzkLM}Lp*!RUOF2>K?-)MmK@IBEdx2m zm~6soR*6_VcDl-_#Geb570v}{d4eYU39hhN^!?YgtmIkhrsm7rAdxSVKe2R|SFXCv zo~Kk5Q=6m| zj#|`zdyL~UFO0n8Q!)p~7V+b|TemlQ#8s(Ji&xjPJhjOcvB00%fww8fvrc4Gi?q-! zkL=0Cm;1{WQh}$z4tWPETgHJ67Iv?PwMv0Cet`kn;&&P~DJWfWc|*SYVc?g?J2o2qCYgRzQdMyudy=;mSg-F02~9^Qm1JtLEA z&(O4Z&806i?O^gl-w4l#z*AK>XQT4ctt8)Y)ibDh{4Aw~NsFsv)4OR&e;YclU}+`m zKy9hn!RFzZcEz#KS+K7(4{;a`RlEm&(ZHfp!qkepm_rQGHpp5914W0 z`4Uj;bTV$Zx>Tmex(oE%wp+NI;7@L{V{2v9f9sp z&tjR2RqpsjIu*R|#LUk;dXC%?5P3#n`&+A_$}jg19C|V5kT!+nWaVdOaqIImT6t(^ zJX_6B1$W@m0u=fx7Pv;-+SO~6k`<@ck&@>3{S}%jF{`zN7?+>ITv8jb~?s=G4 zcM0V~e(5dpgEY$D$A$HWZ)vC1pDYavmhm=kzPx)lP-q5YmGgh-%zoK>;%vX87t5t! zk2?-|>HFE6S9Uevt8r{rJn~VkGGWZD&;3cqd+X{twQq5*yyDSMx!lNl@4YzaQG=9D z+B?$tn;3J{`;}roOM+j#pI{Mc&mTtO+1c=k1i;Oe7$&e@B$Khr9)akv?g(2)?K_0P zQ=wy-ZCvniGKNz@Ki%{9W)3m%-WUpUCQrh{TWBz5Qd>(jtWSqYSi2yhN9QbA>-a(T zm0!q>PSeTnwl?(T3k#O?6)Vw2`o3$rK9bg7mkR2ME1?x+anvuHii1R-aJ7C^)@g=z z45ux?G$+F$=JUYD?N{7yI6v7~2+2+bpZYtVo_r3;LUIb-{e66QfD&?+(fpTm4wt+ud5zf>RXfvBC5n zEV;&}5?q?Ls2Q)a@RrX_*0y4fIpvFPsTJjfRv7bmnOFdl?1zxv*(1SV0x$eSyxOi{ z{Z=o*>MJK^X00F4IvXD*RL|qG+hUTS&2V>&Axe1nkh0Q9>-Qh(j=JQ+3%WaeMx?NYw6>4Su0*3av5q09S_#$7MsTa9R)9!p%FA^M73576%A73fC%o z;a_c4vyoHiq{`pDyXtEVteX4n36HyBJY3!^!S`OxokcpXgQFKwkFqwGovqI_e8Rw= z|9ne%e$-Ud{L^1UxwQ8>LuCi5I`33-yNq>;s;u~I4}{BI)#YVu;%4+J?AXY^B!GPt zwSj+~O`r<@@5l7vz1nzZDm_EkM@rV0~+1$+hULK=(*JO)wwlAYw(DDp+=o9Fp;Vs~5h z)dIBv(mJ{l`xQ(*PzTe44HKU8yg8T##Zb+!*lTdlZ0+r=gt&FheOBsV+NAz|uYE_N zAC8LMYETS35=x>cOkiI?Uwp-spRTNwSIk^ovrofnFpg- zk+tvt1`LIo&zj%wiSFEk4U4Y!39)EsR0B=)ekF|$=HLn4KLlgaG zoa^0If3ZOOB;H~z8n`i&Bi4iX<4xS%g9U6})M6Z9EXYD~ZJ(f*+13gJ8q;ben?lKm zqMLpwVJ1Gmlqqjwk$uz~H;86IRDm#)oPN)nmCa0fVlxk;D><#M?d%L2{zNCM7Jf2v zTs;m>5!Q-!(LIm?jHJf3&luo;RHcy{akS-C2hWKuZ6Di6L;IKVTA~k0FZ_siWSc4% zI=Ua;KOH<5J2jSc3Pw>LhcaHxwh$1R)1NRz!l9|lfAQL8Vf+;+M0p1J{3*z&zvM%e zL8g4qc=iZV zJo8o586Njt7j_g8@*S`AoQ&@aKztxK*hCJi?b%O^u^Gw5lZE;2S=xn|!<Dal;(wfYLhWIgMm8kcqZ0Eb_U4;-V4X_{mI8KPjdrHFLzm){!CxYzlj|L z+8*DhS}{!tsZwgj%vv#t4+ia(jmoxj!BeI;Bh{;1wlg0a95`%Chjw!RZ2yAq1Y2QeTVvGe-eab?{69pjS{90xA?i3e>Li58Yu-)MJGr!7FA+0AuOWIew zbfJNUUGE^iLi4w}*7WIG@rJxo%4Q$EgP6ao4(K~DJ>z~zvFmF8!@r?>I8dm_olMVJ zR&9~W1ZB7c&5OmSdsCCMIOrY1VB=iRRt_2sJi{D*`wx`;DPHJjCfFCs;#&n;F zF(j#hQrEhufN_aF`t!K{g0v`?q|Xpt?jYGcwbZ@P)f()>9;W8KO~K`rK{d(sO)aHR z#XwMUKtPR3tVnJrQBIlUOjEgI*XQwQq%?RMMRE%t{6}A4mAj+~THUfip4d8ZQ~Ojq z3%0m#!SL8jf}8Z+j$4$M*z@5HTGkbp=2CrB z9{!eFlbIl?iFtO`))bOIlB!jqz`8j9iHH+N=nePCpZ&C!R3c(QdEqR4CQ_2Lh1(ZP z@7Aj?8&;r|LfJ*B%yr3`RV1oSYWRjR(XHeaoboK6MYC1MsoQ7TgbHz;~r19A0JB`t)| zq_paOFjpCnX}FC_$qH|@z6x@X=X03pygY(sdq@!YERq^rmKNUMF*7F^zAn9yz7y8JntU7>pEVaur$$6=>nFn3iT>e>^mdx+)itC3y-ZaY{fr0S8G~2 zE4~nFN_d%i2xC7_effyl>I8K$eAiA}2lC&6%=S#p9j0ICsGbXT_9ZHuvzYeo8d#<* ziCZa>4b}V%PU}AupVR8WY1>Zf4Q0BEFe=*u9&ei}t5O(_cZ!$3c+2IswGBH!N(rkq zQ%r~FbC=#To)sTYO)$Zf(-eEENqEynUU!defC-@A$dR4Sf^P{ZtW5yO&CaL9vai1pe)or_0Wy3=!W3W zhhC>gBqm+)X+5Eo(-{xCWNRZ4)mt*{{W)Lp7!jW+9nH>c%$&>$kUTe3fO;bvRN?-{ zfG}Zek@mJOUrg5oUnQ(&=-Krga1d#&i9Z};Yx=_Xur>(zaI=vA;qk8Trt?m&2PIpa z@FtIH&|Z-9o7vZ}80)csu>-Mc z>^Ntfl~s=Mo_NaSqJlP^`V9Z-hsw?6m&a)@_(wwoSkt7O6wzgIhTSz$K%*e_b*@vD z%{b+#CyFb(#dX0f(EnL?5kIAwB8E%uSm^1tjj^gfT84Z4lRl7WxNu6TEP3sVd>G7M zk!<(zC8R0025BTl%{Pc@?_l_)*LLBreL=^e`NK-N!-uC>l0 zjcq!8kMOzPtoRv$cJXmh(@D(sJi__Wp-ACNl;wB=RpXPYV!Opq&fCG!=58+&%Lk*l zx0LjuvBjLxJhGs*GT^c%5gvH{?l>6`<;;jZl22j_BZo5}kE#5T?NMCk1F}&S94VMLA|S$uYZYzs>F#$LxfM(33+ee9%$b8$IT>31C>es;aJsR|0p+S18)v#A|9YR8{4;_{agvY@kT^u=I8R+wFsq_28r|9v%5u zNKkNT)Vo z=GkW^1ZSm7s}Jtz4Ir!O6U>g<={WIE#q{8=GyZQQjD8vteclc&9pAo)XYq9P&d{_> zSV_^cfP5X=RQM`J^i8^ZJeLGf`#dcM(GpR1;+FZ~TF?fql)|S#-m5vDF9ygKT9fsf zdZ$|@bnfb{C)EjgHx=MoL2Ui3(++75>bWiZ;^%ZLZ|J)7K8ZpH*c#6W1_?3`3~&Z- zn5}KS%1(cgG?^Q{$U=gTAVVN%(ina0b}Y-vjJO z9fQqYCfaV-62I=-Zh;)}D$2VAgOrip>@<3{zjxP^l~D%0Y#T|p&X=*iNSG4r3D-F~ zLAmJPW$F9m)A3*j#wYtB8M4g95GS62MU0g z=0e_%g+{pZ*8|GS%yA98-!wx5<|%Tw?DM)=K)=P0_jj75tlQ&#%tkg_9}1C-kU4p& ztK@~3A$!oX(!CoO2CZs4fSUUomV^COw_0oILX3vuZzD9OMoR5$2z{NA>?0t&N411liMOeS6U;$090Jn%t|DF+xZuIVCJ7piss~|wEX#BIrE-@kc4dWCRjOw%hi;d|E>g}SvI@Vs_9f;$7U3oi4$=>9c~OAtl8X zlYSHH;ooAtZV$^eW^M^}S4al=J^7%BUYc3|6KV{=qpw(gBh{Aa3dC!Gl0c2VHgt3jq zxUDZAT|Z-go$9Jk+O|;03XfQVSE&>a7L5vvEscyV22KfhuLc;#Q&A_IrP5Pt~9(E;;carocLz!4}kYFEc_ tEK)^zH-jOB+wbj|em{{@V>e}e!3 delta 19604 zcmW(+bzBuq6cr_;6$z9Sf=g!@` z=bSsQHWKMM6zR{TF>t<|ry4`rHyl%YM{`_>D>UT~G%?ic`E4TAe`GCK(!*4$l=$bFC)JMbtd z{{CA?5Ov?RGdP&r=xzy?XR_Nspl2ElF!9jNg3j=hdZKKYOg^*9l@pl#%a$`2yf31C;i(HcSuCKD+B! z(vU7Ca_Y;0&QLh=W`pLUM2U#~Ju5bfH|h844(-V?&-bsRx$6$~7HvsQvENgP`mUYA zZ7n16zq(0=_d?061;EeT3&jTy2ce9uDg1ENXKPC;f7U@RYWRgm zo_Q*|_#dptYE)(KrLym1%=)PSRi{>(LWV1ky!>84Y!+N-&pVMS0ao4H99JHNJQYk` z76m*ArRiOb9k0&Mo_9v3Ue|2vEnxQ#$m5=U0%xunqbgP8a`;GaVsC%yC18F^1(1T$ zV?Ea*vh~SKUbfG~Goi~YhF>J?Tr)_ok}4Xi^1;ECNpLHg(^ydI%FBqCQ@7KiJoJo1#jT5kU^ok6pe2pae(QBY%L_Z*HQ7H(WK z{yYb+dpCPUz#+hYoK&d+4FYlq;Yp9{vg_4;wsW_-mYRgFr% z4CA|Im|9g%r=+5T1G^lGXKt|xz$W;9WqrVNHENu-UY1E7+}jtdG(Rv5)k$+>acovL zd^MrY=X3(R6K$376V6f}7>bvgmj1t-- zEU+8P9UeTo5hYaQuRNQM9t0;@iM~M7oeGJ_q}yWIV)6kOUnJ9l>--z#)BJM^F%>BY zE>x<9j}mkc_uo4(&g4zYgUkH4?N@!rqw3pOS&m*$ED#5Gwq#n72eN~@TS3?C+5z3M z7WD(5{rhVrAxN%6*D*Y=7S>{SBI@H5+PXLA>djB@p|QeUrRxE1Jo0L` zsZR)DW1UfddjClA>cscFgYjKSW=q<7^#XaQudpH#Di|*?f74V5OaO2NLiIpssT|kG z97EQ;>4Gy~@Vs+!@a+-37JbF~b1bF0j{aX@-*xY5!vt9M%`(r>7;s?3!+Sq@SemB) zz;!DAZLa04o;CtbQt^s$$dCK=C~}XhDk6AJ1qAINNAa$hVCWRBl$uvg3xX?->g0A( zpJ99OC;(~%4H|d^C_M6*#|JddJrI=avT7A2Wk_dP&pWd3jcGN>8=&0^j{JN*%eU?C zV9K2z>pIjNSEGK2=zlAH{{N-iFRayv^Ir8;kcdpLHQ1nnprRxE39eojz6QTD)agoz zCty&(${XH?GjDvl05{yifXp7=1qY@pz+wtgcXB3>BqPPRVPR)WNca#65isqBgZvPD zJT=Cq`Mudr60rZnY+x9!Yt1*S);#ahK_2`>nBpn+eWXG`MQ1gPlt(&wZa1LyEMzk| z)I#M^z$8QE4eSJS3dWOZ92Ng(6cp&CVZBruMji+ye){QSvm5A*R*b~eE-~QnOgxES zrN3vr34A|CwSMM#5+H%`#<+357X$y3{2yP$W}nNYc&fD2)Mg;&`vzf8sMd(HfXTiY zhBCxYw`ve^=5<=>`kjQLcl?)G_Nt>!-GMp*)_pRQi1(UV$yU}LGlezl=huQ`Sh7cXU^U2Rx3kQ4Yb=Ga;gsyBww1_Ka($C2q- zYLIs*w;btCsMh!C?S&Fc>&M?N@rQ$-AWMgw7;$hd%1P?5q}oclk6DIY(shN;I3mEW zj&WGH`8(~gUo)kUtf)0v`!-ovcaBblsXb8z(em(|ltv(|p;-qmpmj=LY;q-;;Pmo_W5K03on;(`EBQ;SHsW5XY@ z(QF9(@{{x!a*1W=*>6lgP3`XbP z2x)0(NqN9>HWLkTF?n8tW|10!8Bc?_RNGdDRUFwyAz`Q4Xpzz)+23aH(?_L0TxQDm z`fwpCe!anh6B`Vgc@Y$?V2~4nGEJyL60G9m6w)>cfy8&qg_3=V zfrpsZx5zGo1u^KTi@&S*raY67J(R8J5oTNq)I|3eM84|HF8x`xETTpr!K?JaGSZqZ zJ0pmeCCC1SUpg(HEasJOvjjHD;lUJ+o(E~Nj!t9IyXb`!)}+*vwk`AA9F^vvY?=}O!+#BH zvL|0Ip%emJ*l91qffC+Y$%vHxxP2-5O7tr_`3orr`EQ#0$yc-Jfd$QiikVqOov2r% z+ipT`9(n&{arA+J|gqPZhe3Jxh8arcJy@ zh-}SflurMqynjDXUmIYo?agCkYB3GzVQ&as@l@4%sf&dp1DUKMurg>uh2CqNP9<5_ zL$-Ab#(5UnD&tsvWrxfM4xvX2xlc(MSC##yR@x1;eoS#QBN{C?n0|k}3C^#~uN1`Q zGgO`u8J{Wh-rlzoOWJVPRKh69)T`EH>cD=qFu7gnn_(>_p`xZUIu-gM_s!iAKi5q2 z%NNb|O!}rSJz!JlYHOTj6w`<5aO7#KLPwLPxGE!*&Yw@|GN(rDpnx>n&PW&8xD_{E zZrl7C-;g}NOYia1Oyx*+UvOS|LNgV#-Q|@0TWgzNk?d^Y!vnwzo|ZeHLe4KZlsvog(~|R3!2%S$eE+gas9rL zOsxwn7QvpgUz1N=!hD9$_Tr{=N2d9tEJy9wYx!?qW4`?KU$n?dSTpaZ!hg*Y7?08S zZbV1MrlNJ_)U?R_CoAIw4(~@R!x|`@Ch%KmWljW7JPEStFj^#eR&#_%3fKvYm)NqX z#D{l*RQbkih7cQ7$lauKX_?sW*#WW&Y=QCVzmqu~X)-A}hKO!9d6ZEs#Ud)*@87=? z=`9=hH{c*N2Ak5#j_lT47}P)%rU=4%6)Vq=yV@E1AxQ4{GJfbautGx`3@K zN#UyJg65WgBCX?VtFB{>HCAf_!?g|Me3Wj+^stZ1jq&;W8}CEtk~IVqe>Y|7uau#B z1v1gGlU2{8L78A{XxAFZN5&x?|KrEBXKO`mak4}dm3gSlEuVxV%O+8mo`G0w&*Gr2 znjW}#C^=c53E5(omGi=g{S`2{@R#cca^v>qeMaMoa8JIFv{*}Am?Dc6R;mB#_^n8> zV@agHAQZHyj63p(+Ac*!qJES3(f56?dj2d-rH9TBlym+uN-2EzSz;xcFjyn7P>b}9pOPw&Uj912P z9@~Bv2g#siLdTV^C@FgK4cfQlFav2?U)1`+;GTyQG93zj0Xul+q9<=6-n_q2PUTLz z&k6~ODesLi0M^({B|onIi<=J#%u2>NYKDnPjav6dyYUtNSbdFeS$7YoD%07S-grfX zIz^K{eW(y!sbK3an%DFlfGZ}7F>qnTt~%DuJ{zMY%(pesrEsq()N$h&GJ;^QqQ#WB z5Fyozi>}>;?<8A(jQ9{U=Yk9Tyqu)FPrm-bghHsZ2xZO@#%h~{jGnqrOjpyg_cu48 z**Q`3G-&|$MCsi~i%@=H0>1_9wQ^v9x<5^J+bWr$dYvtLh@wY7cFMl;z%yD>(RSEH zpz}e)sa7g(`grv9j|9h9w~a@p1D%}6hXs6y5rZ4hsP1y;*bhr@Mdtz>bj`iYZ=K(g z=PS_BZ$0gv9?kp+k?sk1Ku;J`?;$9$sFhtwIF$=>_M_Vhr@L~dC)E9gzBOb%q25Oa zS#D~ecJAfuUANwRI50lu=AkNMOxL+iY02 zNxQxXsuoYI#I9fMwd8vK*_Hc#E!6+J)WUHLQ-`SS_iZk$b|Y|_R$) zX^|4QzI8i*P`y&4!;EMioFm6^efrW>2f|PFSQeMlzm{eXqV~nn{WDq+d7I)Q(>U>f z>=Z!v)D8__qX1QoAEUEus;`s%c?`dC;R)&i6&9O1jW$qv5GSf_(y+V*F|~~7-$I<5 z*Y45PZ`o_t-=U%#zb`GsJOqhFjuIO%Rvqx#9~;SlrWO-UtV1EC;r>eyNRT}7`@1%N zUfxPjQ6nT8=is^K6!;xzaRMOQqSCJUrvWZZT42Ba3;)0uPO_i+C~ZRm#Ou^qe8!LU zC_D;9$$B@v^=smKbcs`Ay|%77f6ivDvR^<&NM0$QZ--{N859Gu$U~6IF?F+h-sN71 z=xu*^^&c3%1W|UKQSJ7w_-xA|JtP6z<=-~oJXd`?qE%EGfZLe$Gl9R0IGLTUG+_4} z!YK!mSN#YXqn~Hv2;@kKWq$2)})VULXNkPDl7#F8)kmkk|b^VJnsz{r}8sJru*4xJKzkZ8=3$8c?ID$#AX>BizlXt0#1U8 zVBz4u8f6cv3MsHMg7u?d6g=YmE!(T!Re6k9%+m)-UVa^zx<(%4 z@?aaaN_9dg*8eh{%+?LbS^%FnpM2%tAo>sB!(;%#)ga4O0kjVNN~K}Flp7?Qa#}}k zLC#AM+pUnfU!C@9W+Tqx*8j?tErqT|zsu`dh8-R8l8C&CJp|}Ag9;bZF|XG72UrE6 zCFi`h*|t+q%7fMka-<*CwYvt95JlC*NRabjgZkGXw{^UD!64!_ zhyB;RDWD6+3Z1oDc8`1yc&48mR6?r#>^f0}jBU1X&hS9HtUH&%t)cdpQadaBMo`a? zPpx$?J)=d5K;MyfdOg-%Soe^R0X+n2hPvKy8 zPIw>;E1*OwfrYbirt>=PlCKLD7~DXWocPt;Q|#!P;duc(3q3Gz9!b-kca3sVXsrR9 zXE$)RHkDcuWa_Ungvi9)PwHd*SM7s4Y2Zio|G6QvN;A?Kd}k;dFHG5`;K1C#Jc;dN z^!rsE7^4zNB};ov7s>-r_&Z>eI8h8)*hhlO9=XtAp?KhK-?iTt%zQc<5_q88iJlZd zD021XCzvtR#Zx%h6W>LkhZ>rPAaA#WE$?NhSf9a`Lm`SnE?;o2sH9#tjWD^v4c{N% zzD`s|fiJUo=5)=kSv^qG5GU0(ZgYK=3owK$bk;%bAtw@}5J#P2Bap9DdcaD147%y`%O#nNh$4XKqzGIQ;+MRv&j?on10&?)ejf0&PsBEd$5}_Ph%~ zt8D1Bk^|Eodh?;X}6id&Jx_0Bs`^kCfWeRKbRh6yRn zSuvw~O?!!Fjx|6z7uY{GoWwY62ozY1^wMfXQSiHb;X;j3E!O@@fO1VFhP05Tr2&|n z%=uX>C|{FPa^lLXZXzzYlu@@6!1pEr5O*0#`V93GZOEYwOaGnCev_zy+X$5FxJk^! zOR^|aTH4BQV0MVOPj&=;F{}}w7YEIa-4v?a%#1p5$1{Jr5Rd>8DUfA`@M}38% zv#d|$)Fh2f@3PM;YE@iP^U2bPc1nTzl#i-q{Hq?kL>}3x!XKW65)~-op&Ctu@42#( zs}n_qhPD#d9!2fIwSgXcown3ps#7%;LRn8kx4!+uX&RV?Sk$dn#cex;NO`OI;xvcj z5{EEdhN0r=tV>HB!&Q5xi^_}b(4v3cR2}nOOwC>2+scdNE_Kl2Tz;^0`%(=Z#9`5U z`Lv4ijE;{g2%+Vw)Td1!5#3>lFTpkT)jK29@O8Sz8N8@ zdG4&{p?%|Bb1hS~&<7qYEjo;zK}Czj_Jw{P+HhU+hL{px1%P6~okM6`wjj3F-qx&o zV!X75l}yu(_~DBY#Y9Z5o@F0v(eL@S-B#gQ;7sk%<1W?D$6?=osXDv zbl#1MT^6eFy7b8>YQskpl5a9a{>-1Moz_Bxk8MoN(7r&*%==mwL*3MtzR+nE>rvXA zOt~wWouWIL%XK|qZO44}qXr7b;n--{5?v8Y>k00_p(rYIdpG&H)W4u8g( z&_lssp<@2)5?-MtojqK}nn^9sZQ6#OcuO{%G$`qO>~dy=B2RC6bmj`Ka_aI@e}VxT93H<-!IIk=tdmiQb``)x}8$jExP5SSd$jzd$3tT;zmVzC#^ z?fNy6KjYBq(zGa@E3}{;{0o^!sUU%q%Y-U1reSfTgtfML0K^(dYnTCFI+P;>7jg*lCYXP zIPbk0)nUq0=CE;^@Ci_%6hC$EiY|N?0Rwk?vMLDzY1G&NG;u4sv@QTF4Wa&49X#hb zH3{nwpgs9!w+iBY$F3_7r_yo?Se9O5J3R10D9?>d7S)JfdBttqrBH?B-;ELK0HP6o z;}wW&j<{;hG8_nA)f^+EEG3v|?m>c~Tl?M|VJr?KOb&?7M zP{E!-Zq5+&AOH?w;sfc6z90DEWekf-qye$BLGD|3d#4o=9RPa*#UBMw(<1U~a z@l>Jy!b7ScC7>4R9~?3X)r_+?vp+tro ztTQ%xpAM4TfFuLjAxtgBd9~v74b%s^DbxTs6abi}cElz(&|dfDmwHQYhFL)XjyQ|X zPXm34ExwvGvIXMbs?&@{W8G*tWffNdQ0M)=+v3D_0Qz(83-h) zUzqe+z%mZ&WBP<=&Ij(tC{%zpMn_3+0|XsWsDrGPdcxkTT=750J9 z7*67!{q^CJSNW?%;_uSs=4lV1)kAtAJ-z-iUWMol7rS;Yx zaG{;)4Ixs%P8SoKT{hFu2Xzkg0fV$m5UFaT{nj+l$OEgubbI$3W^O#j$dwW7D&jC_ zrz|WVYbM0_(t8ZDzKu=Uo>akckH(t$wAbh#qww$OO_@-VvuU1-`2^ zi{S)eofvz14~bD_SF(ci(1q)uG-+Jsbs`wROi+w3!3>fSCtVvvR{L3y#=@8LD2}Z& zaf$aByp>}h5zSfAo^j9FcH_yzuH=Z-I6Oz;A%D|Ixh?+Yk+RVZRYF`N#H|n%000=O zzb5SWAQPSGZ(#zx!!#tHSK~U`UFG_iWSpp?fK`}id>bwndQ2THg#juI=+)g@oxv?# zzCFS9lpY`#Os#@#Bmc$8UZd$ox@JNPT}lKaP=OMYk-MCq2(rZ7Js+Qh_? zE}&T%*f?_>J#;Th#>bb!&}++Su-IrTHW12SK_%ERw_N`I*jj>10GPE-9xhCo*KK4( z%u6SsSqO4Q^TF|uj9Mfmb4Jc#;Nzy}a5@k4n;q`&zK=WaHnQhl0BSAkHRk2g@iJp&u(a}UW&+2u%0095X89N_NJ-(NmFs25`WMx4 zf1tAZPRoWaN@Rb**EugLD}9^RKMZSnD2b_PM5lqlq;nhl4MrEAv^iI^vPXb?K>2{( zRuEZc>|bhK>(taYuQc2Q#V!524+%I4pG!hvt`D~yW<8$8i2+4Twrg&)Z6F)Bx=>~a z!< z=eRR9d}FaN_g&N4&>X2ro3tvq=wK%h%N=omvi@G{ARxxE9v=lKeR=U#(l8=fg{6@3 zJ}(hUSHeM72_S1p)|FB5BXv`B$fudglrWTb&lyXqY{?54zU(05QVYIM!BwA^+wmgM zurtfywegfoyA;lRZ(`|bsuS0OzPWpK8@z|&!J{bS9^<5_5jnoYs2V19C4`=TV%k;t zJNMJzwd>?OLd~av>C~Z;g3&jzHcm$+D``#Sf8W-Tlc=1rSL@sH*m%+@n=>Q&jVoIx zAk?KDTC>5sU4K_q1Y#jou}7^xeC+4zJ+TeuZUPZ1CDr{R*qmZV>ol9le2}Jg4^Qk{ zx3180zOCJV?^{0SRO?6Oj9ClA0UF*IS+Z0L;ufk!nn@bVY#YC!=a`z}&LRmmZADi< zO5UdBrs1ezJ))k%m2f)M?WVli+$i471DaOt$D;^r)z2a}Qp_#{CAnctD~_LS71S zkbQq7EA-e9{XDkexzcWs?piS(GST|A7Fo1W(v4S!h2mYgwwM?OB;Dm2MVgS#B64izVD&~keYNrLL#nL`^EUDwwaG{E zO7B}wnkf`yJ>x ztwsn#{BL4(E!o=ivHsrtsi@xZpz4t#bZXHk3$e`?5Z5B_^2&#+QG zx9_=@b9Kc#&F=Y~o!$YMW<@qZs9Gt(b7ULZ6PYeJIW^A=Nmho@al)=do@4&eSFBw$ zE^8;D8}EWjfA+vOBO6aAgG;_uk4AWX1gV-G5XVh;UxUj-2Q#vonHIt@<#eGBMKSlV zwJ(~Stu6aqm};l=nqJY^DseS=reA(RB1?t4ePLUTWk;$~YK9B=$fjdxqo7}+!4Xb% zYV>UIhF>eHv-W;?z!-f_ua#(W6&EgBtDE@|hp zK@vzmjZh=)*;<+zy?HkwH|?}T4*NY4PR!1Z3ly?KOb+2t333xf`IB0C2cxijw&V-e zXtdAEvdqihRlr-&Ca=)ZbCtyMFKhjx7oSxUDB5Y zv{(+_lCN+5+9y}eLSdGji3|`}IF|e32)KVGn|&@j{%SDiBF-c#i#x+OPwVa~xKt*l zfu??7yrPeA155I$%opZ4dty6n^&H|jZ@-!Z%WwH>6S3726UaUIOPewCWb{3}!Ol%E znNcbPsdX`s0>iI(T9OpA<{sxRdg_efMSE@Wcyn3;1gdY!zL=!ywA>x zXgXCo+SOEb0dPd6>*%DUVKI8qv{1+~VbMfDn5pDn?G5vC@2`P02lI%@Ysq>I>F2q} zEUVfNRzKSlT8f2S!Xr{!g!V8GGcklD2IVdFX+C;+JhfU(RJ4mXDm%E-VYk6KXwx!7 zq8D229@i>m%C~3f8)D~r7wRCotuPm4WJ`jL%ViW@WG_23n3BuF5@e>GEQrcv7UgrIkt~Y1K*WKsOaN&NMDlj$zYp zL9bj9uJPSU6|?7(c`IcN?Mih-WUCc=DY@}*9aM)F2zhS7jB;Jo=WBPbjdAtOr9G3U zY7Rcre(t=eZ=gYGw?iF-r0SWGCL})`q)%BlMwp$-h#~fMMfJH` zf&7}z!P5JI7;hxFS)_!Ayi++(ItY^vj-XU+Ta^Rf+FBCn_C&pZ+dy4v5&bc@)tjMn z*p%H?>O)1)mzrt(x&FLO*lj+#6+;^{)qE#xpPlZgp_GNU3f(BkhnqG1NF|8|04YbB zW#tXoV~uPmkE>N7zDJ8q-j~6Wcm7}?oIudfi2`h+NUd-hd~EmS zCmR@hI3f^_Z`)w{I9nn`Adw6KB0|~?k7r^H$Vb(TKHtYii+S--HZGcyrCS~84^7WNA(-52dLdNxIb!ubZ7N05t7vh z0Z0w^={3-iolpJRmpcz5W_H2byfuTsS)}we$Y#>R=vG{{K?Ps!pkU*y8Si<7=;%Us zFfxPvpOxquXGD(@Tk2W7SXNq(0UCBwtbypaKcKWdbW0HN|Fzw*G3kKUtGr1?fmYwq zeR;tvh*|l3ti0d>&VILi%2e7tOLPg$*9mt^g@i+2wpqC6M7r@P$TI-?2{a&x9<(B{ zAkO;>l5M59E4SS;kd>sG*)PGVn}GUsY?zA2`*N^*@9Z9d_%0_8dC?l-!|3j%Eh)OB zb9qlH2mx#J~mM%W@CI7svM|%+>!2$D#?*Fmo!rM0WTsW$ho}QCoH~y!=K? z;U|pnb(w^TIIh@gEH7nugtzj7?r!pNj5w`Utg1~MfO{&zZ}EererN%$@%V5gWb0}( zw>Y(Gtq9$CPqR0yPy@aU$Yr0TJoGS8qH@{Ymcw*PJVmmGd&MuveG(-vZSF|7s#no z^XrogkR6LMpfjzP!ljwlEl(q`u=IVQJkLU^nTWZ7qj)*jfK5>|cH3)ioOeObY;GEz zbu8g?1w~O>Qx$0bFh>c(#Kg;HrjIo8BWJX8Vg0@i+b?t5=JD+h7G)ZHr&WQL@FGu*b%;SmxXdA~^2l=p(09BJ8|P7bSG3O8nvxb9 zsm@%b0aY&`^kHd03Y8LWI;@b1GgiV7CdCnDmykUjuk631J0R*tOy8@FRZ=3{QaN-7 z;`e>6lso#a!}&NaA?2~5xwXV^O$oUA%|dwXj~MRN;)xOY@|zQ^)diUs{?NX6YE(MF zw7q+If+rfFg>tHsu?z7CJ<^k(@+_^XGrEY&KhmRpPL!&Ho%5oZn6FUmpgR!mTcxGzI(&XGY0Gwz??z9N{>z;NF6 z^O=#?J4w52E6%Uu;nRqsUe*-GtA@s>#KI*O3iuvx1{4tg&ecq7Gr?k1Vt&46xE`f2 zco$+kzL?VbHU7F>YQ4t%RxP0w-S*Db@jbq2V@(#9<7zN9%0iQFZ+vq&OD@)Ug?*G? zY}4Af>*k)WbLu-~nI_1%?H#H^^sY5zuk=)R3U~ z-g93&z;*8K)E&|IgI-4l9mRR+XlA(7{G8xa|8TJZs^Ac%3sJN7f{et%Tqm%uL5Vin zZmW^D2I_aqv>EO|fX#q+#Mn}?#eqZ87m`1V>E_~+348@p6%%N_Mlw^R3;<#|b^tef zAo){oc&ryyJ96J(k?vO=(oJhEu886jR3%yMBjcdP$N7=1Or3|ZNaQT)->VoN>5XU7 zXvTByZgcEof4P$km!kjL<3A*o65Ofn0c7$C);;U^wvFW#6PhWP1EDb}#NW+oEya>X zT!v8LW+DvKftkv!m>=1lB|v47y{Do7L-(9%i!+*hwuvIo-t^umrj%E(w^}>1@unGW z$LFGjDMoaugUTTfH+_T4FsbZmPr1Ed-0*pW4K%1C+ds~uoz7$Gl9Sb)kF~jsbar^Z z=~|6DJk%X9*q66rrjDhIluEo}loi;>qJsvxv7M7*nRX$(Rjg(S)4;I32`$Z3_jIl3 z@`zr+{mC$Yxqa9}S?`CTtIBISID&yP1uHF`z3}s_P%o7vD6{BbJ-09~`X;=N@*=me zo*flzi)iP(`-W1(ww^a0Xe{$LqCBN1wal=E+u>$d78lITxkU*7iiDGBd0{-qWz;cE zj3i5DLYf`ridLnlTL2`kY3G6fIhm%R%*sW~A!f6-&yRs3Rp(pkE~;Zm(KMyxb%}j- z%eT_W3W6W#bi=NkYhuVG7llSs-#hGgvtaOGYBA)zpXkp0+0V77;HEj?3!6UAJZ8Ch z=^Q28MQr+`Xq=kJalYOkhwH1f3~%2<$u5#4fkk0C>BbgT8XHh55g5{}zIbb_@n$NZv&?J+c%^Q=z80 zz^AxTYy#W_!fNDV6v5*9X4?=W<@=}=13ZhK{8dVAFI8ft92pLJla{hpWsa@l6kMB@ zy~CgV!J&s0>uNxfcuuv>dCO9(9CYv2lTDIN!;j(?Zw4%;gubU}a&4Wq8H*}Ogdy>M zTwQJD(Y|Gg@yMGq2}5}dwW)#}^9kFA`VR8C+{|I4kR+}MU9;o2AJ5CD6V_%Vp_bmo z1mB4c7{n18MW;7ZAwL}w!BYq4(xse}15n=h0;+J9#xH<#%acj!oYYM3;PX+rAW#0R z)8hO7>*V?KAi2GxHt+jRx%B7#m0S+=TCBbeW;VBnJ$|na56=&(@`z=9`pnQ$zo;(? z)1No;SG3!Jhdr?OkzD#^b2``lK|omh%$@b-!EpWrLr8rtTuj=gTZdzSB000w$9n6N z9mmm`R6vQCht!9m@}5-REsbr(pc<`Jb;48U`1T4mS13eyy_QokI?~sOQXMgdUGdbxP z)_t^Lh#@F%3>Q!!#kFrb4>D3XKmWjB6f<2^FuNOlk}&1woNYOu3-cpM8OKd~a(?fd zrm8Sab>Ih+kxp-^%zuX12XOaq=*0R=Zx7h2!kKbxqB>bN21QtKBwAwn{#-s3?YOIB z;sV8Ql$Q6-7M#>6r`am?k-Ng>tl4Xn++E`-w+{ef66J?T zACzxz$;;euimr|}vt%9p17Vk##~F_Y#?(+c+`kH1k*Y})NhFmr4HR6GNbVkXMFSqW z%)8QG3?-K@KbD!fL@JSMQo!A>e~RK_qi!N{MUiyi-$jIyB@)Yj1txot%%JA^aIJd{7HEcSyf!45f1 zN+Xg>oW_LM{GAsIN~6vFo0|!{9uFY&c<|V*$-)z0Io3Jmxc%GoJShGGX9ZX`P0ivw z&!%!x9mgn7*U27aXFXuS<$cAQD!qsl9Wb?3w|F2FIQ5mR_#*rz!M;RLA)jJSy z+x8L!+-`8+Fl8vNagrNZ{ptpMX*$o^^e0Q_jf@8iup<@Kf9!Yo4r2a{(|$a*c|L=L zz*p*IPGPOTfdpQn%%t>}I>X;jdHUc;H>B`6x@_g^@}12dDvV$KRe#waIIeo0YLlH- zJe{PgeJBPldR}(gJfAmKp05a&`CdA|Ctrqgq~1)>Xk9wX^ZmJum(m5=+xXTmb5wse zW|kxM^iJA@px>T@Svey>XRIHGDcpuyYn_d2VVZRCO%9)3ym#{9j z76rMW>P&Lz3Z{#`5-TZeJ43E#P&-Sax0hcp%SY;g!`gV>SrzRGy#!9}$$)Qh5*G%r zp<4;iIJdVo2_KPvl-TnDdyYx75}sy(GKV{{H)^s|jFYm4=`d#p%(EFvE|~*)Gv;OQf*srR$)4L43MlIh+JDh$eP^x=(xM-QC;{4;fK}@ki zA*_S+k~yET|L+`g3bY`vjOLE}921-dIl-nybdlYd-p+li zP9R-IY4?y$T+fvB{nk&Vl;DI-OCat*i1)01_6IMR|Gz%hTA#W?L%ad4n{P{gE#<4O z@HbP9Sl8Z~H$T1DeRt4%i4lKN@xABrdc@6RYFiF}X=X2p(r`0?_ia>PIwDjk4w=R2 zv{@uF9PJ{{!NNy|zkaUHXSov0%~62WyDRaibn85M-g(+UCD+Pf&Y)=*>bjzdp&kEI zuxoOr2vDe_`=F3ZVP09UUs%u585zgFPatUK_p&oKCevx`Xk$3m^*Y0=4`P#%NKkRZ zhx1;(h8x*s@LqVCz$y>_16P4>sOx>Tzdo2?gQ2}=thfz4v^}RV=n%BpN_CHjXgn|Y z-LIVLv@%N2Z>etg>L%xBTnGG3sJ`P;-;PYN0KUR-E-3B%yz&36bTXdtjj;Y1?pdkz zRMPz=PWAi%@x3hMtj1?m38GsCFg&^*N^{zEF|58IYpR4_?hn}PT5yz}l4u#qi2r^m zMtuLC%88wAN$EMT^*q)!KHRZyFuG-|TunfJ59Yk={6QW6(#b?Q>Q}%Vp=?>385EQm zCx9JyKZsoJXVd23)bHW#2Ln9`G5Ya1pT8!DEp`SKSD=r{lBZ%N>d$@|f66$2?^i@d zI1$C6prOHhM?uPpAZhMF9-=Yc9&3{RW7wkNds?VuS|U|;2*PkK%k8)Qnk&pXWufti zC+v_K874^k-(k31ys97=LePzWIXC$c>ZmN31D!_OaAM%+fE8p2#=OG}x7niEWq8J_nlW5@b%^S{(ChIpz*@w?46kTr3 z_;4x#HRBI%*qiy}mrVQs7kwj9W$5^vou6q3ZOf+{R-IK7FLCC;%9KMO1@NAfr^bt4 z)-Oaja{GA5NI!LTkx;S@C4O<T_`iFXY7&Hcxt_n0|cb)9~`MZg6(HUcJD7Yz?a ztkUet9BU^pR@Ae@+E9;N3erZ{3?96CBJ$c-zd}vS*^<*C0tchseW8#9KFk^71fII) z?O;=(ov+)q)p1w$V{biXe<+Q7%Ge{+X$%$Bl|jY_~ZMo zlfLJl_qp!texB#PpMT!#df!=wskQnOij+17;x!Ny9(u8Q-efN8H`rpHZC9Azk|DOq zOoFF-;cE4Rx8YvoV22S+v*8Tt-F>^GKi_Zg!?zrG_g4rrA=2OBXm-c@r zKh@#OT$Q@jE;%?idr~cKLRfU&@wJ-HxJiJCtsG^fwP+l+Bt)y<#GBxWwS(URAHOI@ z<@)3ZL!O$}>KH2GQ?rv8R$BN#&qpBb-Y}izrhrKT?8%H>G0~K^CHCBdc)Zy(SoM23 zRG?1IA|>UjCfy8GaT%omBhvC@3lA;fSlWe8Z*wG559CgY*(mKeCV{lk-X}7k6s^=9 zN;dwy65`PEFy)ngWqctvLTJ5+voejDo@bpj$}Qmbi`Fl?MkgQdcXed#i_Ym&3r~u& z4`ZDUKflPEGbJT(Tu52^hu|wQeV=|v8m$Of6HnN6QG`nXXdb9#)5802M>AC5n-#-E zk&H=wR@Yec5Hm>fbeWO$?JK8b`@rpeJopz@yfx9-;0ABTayu4|O3-_-xV_;;z{LI! zY?ar3clYKHC$P^p4Rf~~UwY9ox0!!-WRuSx{m!u8ZdcB!DUOt0t=+FNK0U1oa1|E0 zB`x31^93go9w`DD-4oc}PrkZ$XuDX6{SK1%xEGA2@6o;8cj$k$IQjCf1O2PI<9m9n zO5Hf7&G;IKf8)J!$Kf+n*}Ctg`dysz`Lb!j;$`HI_xt8bPRoX<^z+1wNSqYm{1x}- zmQ|(i66@2-qu=B@1hm;HD}#N)LqALyx^@H(GrL}tLZk+OwzL}Ro1E*z-10^hxmU>z zr6YDtdg!PYUy17+v2z-&)HiYzH$@jUp^)&gWtB~+r1}X{1u^}ROY^Zj!f5o3QF6oO z2_gC8;ezF-6PbK`u#sq))K|tPH&(RExV$Ox&mEE9zGP-`-&;J4<&bHUB|70U0%|=q zs;2tRbRVF7kK-5ktV35u0)!;3ebpTwC2rMv8IyMB1N;H)T+tosU^-g%4j31;ZT7sL z%gUvw*+E;CM-ysIwwx+NX6b_2Eh`v3d@5g)W`9m~=~i#<;L%*q@4`F-=%l;nG0`S| zRnw9ilX=BRSmJKMbj743EnG8i#JFv)O;lta6aplKZ5sbw-|afSKqTIm#+Dg*I1;lo zwAfD@<+sN;PB3Dd#e}t6w*?tkARQlecNI*|R7=-a`K6uwbVU9P*~K?6P95`mggNhz z?&EAXo6jh@{NK&wx;80Uy`9aE{=z8p{urb%JT2oV%D=kD1=oW9Jdh^E$85Z9Uv9+zZa z$7I_fX~)|P>dX1ZNx z$h(`nv{$juvp&8*#`VrEbh^}egU11WsB6rT^vD0y3lH*{ap*KI--?U${@^t8T_eg{ zD&(S^N6WYCqBz+n*s;8&%WC(n?CFNltB}UG9eIZH;ZU~*t<3(Foo69f?e=?|Me5|M zBWXJM?cG0I5l*k!hFZ>^U!#@3`tJ6!b{r`3z!yMa^)}&Dn*)>*pnBPU@p2{bkZvqm zo$1S{N|*Rb(e1UwDx;qKIo+O_s3~^eJ$Hazxca@t#i$X35&CAAZ!y5QO+0 zdKYFu$8MXy4$SVr4Km?weQm}s5mWo^48z-P_O1QXC@|d90Z%d7Dv?vpo$5K^%1Yn0 zQ%5OSE7_ZB<|6itFu@}50i}iHe4jsj==VZB(!(B`;G8pqHCIC|d&oxGu0eb84RoGn z2eifR=nFv&YS}cfk8BX=G;l5MPw%F5>?`L)ZhZWq36`i4Rq3i=JN1@wGR}g;tb3ZE zkaeHqYseg`i_dqa4_o8HOVs2upQTF9EB8P^e%1LQ!>5nX4|ZArK63V7r~!4ju|c=j z%?;V&6^twPAOq)hZ57MO`N}STb;??exqQ$A8MD1shiJ=;yHGjPG&Hf;;KwZa#nXOp{TB z#hc_@^*Cni9|~}0>+C7Rrblkx-NR7OHDjsCIW`OCLv|sWOHOPv_WjCoW@aQHy64Bf zKwjR=JKYW}$Y0ZN|6A$0Y-g9H$lU)wqzFq&noo6O>}${Zh|j&lH68s#&}zhkgFS|& z5!ueEvxc5dZp>8Ma=}NNq8Y3A`|l3V%8#Js=C~E!y4J6{ZJOVwXp?R!>b!4TP<8&I1yL8dF^RR zZ2!saa0mp7fuPn{Xj;@(&~aY~1gQ#M57i)5sr3_rL88%^UoZ>;z7a!#5wKq|)c@gO z5HwwJ5PAy&1wqq@;-G^70*To$7Ys+BVCbLwfstq^a?RU%{8t|sj)bBhzhHq3fD~e>X&;F)-A|;CQJp3Aya z2pqX)@aMr0%$mK87#sygA~x{gC>UgeV+aNXK_O|w5+KX%C>RVTrLJyaXDRg$Z#)g| diff --git a/motivation.py b/motivation.py index bac6148..bd8f48c 100644 --- a/motivation.py +++ b/motivation.py @@ -106,7 +106,7 @@ def motivation_all_small(dev_desired = '1', ylim=[-1.25, 1.25], c1=10, devs=['2' arrays_original[3]], spikes_pure, cell, grid0, mean_type, group_mean, xlim=xlim, row=1, array_chosen=array_chosen, - color0_burst=color0_burst, color01=color01, color02=color02,ylim_log=(-15, 3), + color0_burst=color0_burst, color01=color01, color02=color02,ylim_log=(-22, 3), color012=color012,color012_minus = color01_2,color0=color0) ########################################################################## @@ -285,4 +285,4 @@ def plot_pictograms(ax_s, grid00): if __name__ == '__main__':#2.5 motivation_all_small(dev_desired = '1', c1=10, devs=['05'], save=True, end='all', - sorted_on='LocalReconst0.2NormAm')# \ No newline at end of file + sorted_on='LocalReconst0.2NormAm')# diff --git a/susceptibility1.tex b/susceptibility1.tex index 29f296e..821122f 100644 --- a/susceptibility1.tex +++ b/susceptibility1.tex @@ -295,20 +295,19 @@ \newcommand{\bdiff}{\ensuremath{|\boneabs{} - \btwoabs{}|}}%diff of both beat frequencies -\newcommand{\signalnoise}{$s_\xi(t)$}%su\right m - -\newcommand{\bsumb}{$\bsum{}=\fbase{}$}%su\right m -\newcommand{\btwob}{$\Delta f_{2}=\fbase{}$}%sum -\newcommand{\boneb}{$\Delta f_{1}=\fbase{}$}%sum -\newcommand{\bsumbtwo}{$\bsum{}=2 \fbase{}$}%sum -\newcommand{\bsumbc}{$\bsum{}=\fbasecorr{}$}%sum -\newcommand{\bsume}{$\bsum{}=\feod{}$}%sum -\newcommand{\bsumehalf}{$\bsum{}=\feod{}/2$}%sum +\newcommand{\signalnoise}{$s_\xi(t)$} + +\newcommand{\bsumb}{$\bsum{}=\fbase{}$} +\newcommand{\btwob}{$\Delta f_{2}=\fbase{}$} +\newcommand{\boneb}{$\Delta f_{1}=\fbase{}$} +\newcommand{\bsumbtwo}{$\bsum{}=2 \fbase{}$} +\newcommand{\bsumbc}{$\bsum{}=\fbasecorr{}$} +\newcommand{\bsume}{$\bsum{}=\feod{}$} +\newcommand{\bsumehalf}{$\bsum{}=\feod{}/2$} \newcommand{\bdiffb}{$\bdiff{}=\fbase{}$}%diff of both beat frequencies \newcommand{\bdiffbc}{$\bdiff{}=\fbasecorr{}$}%diff of both beat frequencies \newcommand{\bdiffe}{$\bdiff{}=\feod{}$}%diff of both -\newcommand{\bdiffehalf}{$\bdiff{}=\feod{}/2$}%diff of both -%beat frequencies +\newcommand{\bdiffehalf}{$\bdiff{}=\feod{}/2$}%diff of both beat frequencies %%%%%%%%%%%%%%%%%%%%%%% Frequency combinations \newcommand{\fsum}{\ensuremath{f_{1} + f_{2}}}%sum @@ -425,7 +424,7 @@ We like to think about signal encoding in terms of linear relations with unique The transfer function used to describe linear properties of a system is the first-order term of a Volterra series. Higher-order terms successively approximate nonlinear features of a system \citep{Rieke1999}. Second-order kernels have been used in the time domain to predict visual responses in catfish \citep{Marmarelis1972}. In the frequency domain, second-order kernels are known as second-order response functions or susceptibilities. Nonlinear interactions of two stimulus frequencies generate peaks in the response spectrum at the sum and the difference of the two. Including higher-order terms of the Volterra series, the nonlinear nature of spider mechanoreceptors \citep{French2001}, mammalian visual systems \citep{Victor1977, Schanze1997}, locking in chinchilla auditory nerve fibers \citep{Temchin2005}, and bursting responses in paddlefish \citep{Neiman2011} have been demonstrated. -Noise in nonlinear systems, however, linearizes the system's response properties \citep{Yu1989, Chialvo1997}. Vice versa, in the limit to small stimuli, nonlinear systems can be well described by linear response theory \citep{Roddey2000, Doiron2004, Rocha2007, Sharafi2013}. With increasing stimulus amplitude, the contribution of the second-order kernel of the Volterra series becomes more relevant. For these weakly nonlinear responses analytical expressions for the second-order susceptibility have been derived for leaky-integrate-and-fire (LIF) \citep{Voronenko2017} and theta model neurons \citep{Franzen2023}. In the suprathreshold regime, where the LIF generates a baseline firing rate in the absence of an external stimulus, the linear response function has a peak at the baseline firing rate and its harmonics (\subfigrefb{fig:lifresponse}{A}) and the second-order susceptibility shows very distinct ridges of elevated nonlinear responses, exactly where two stimulus frequencies equal or add up to the neuron's baseline firing rate (\subfigrefb{fig:lifresponse}{B}). In experimental data, such structures in the second-order susceptibility have not been reported yet. +Noise in nonlinear systems, however, linearizes the system's response properties \citep{Yu1989, Chialvo1997}. Vice versa, in the limit to small stimuli, nonlinear systems can be well described by linear response theory \citep{Roddey2000, Doiron2004, Rocha2007, Sharafi2013}. With increasing stimulus amplitude, the contribution of the second-order kernel of the Volterra series becomes more relevant. For these weakly nonlinear responses analytical expressions for the second-order susceptibility have been derived for leaky-integrate-and-fire (LIF) \citep{Voronenko2017} and theta model neurons \citep{Franzen2023}. In the suprathreshold regime, where the LIF generates a baseline firing rate in the absence of an external stimulus, the linear response function has a peak at the baseline firing rate and its harmonics (\subfigrefb{fig:lifresponse}{A}) and the second-order susceptibility shows very distinct ridges of elevated nonlinear responses, exactly where one of two stimulus frequencies equals or both add up to the neuron's baseline firing rate (\subfigrefb{fig:lifresponse}{B}). In experimental data, such structures in the second-order susceptibility have not been reported yet. Here we search for such weakly nonlinear responses in electroreceptors of the two electrosensory systems of the wave-type electric fish \textit{Apteronotus leptorhynchus}, i.e. the tuberous (active) and the ampullary (passive) electrosensory system. The p-type electroreceptors of the active system (P-units) are driven by the fish's high-frequency, quasi-sinusoidal electric organ discharges (EOD) and encode disturbances of it \citep{Bastian1981a}. The electroreceptors of the passive system are tuned to lower-frequency exogeneous electric fields such as caused by muscle activity of prey \citep{Kalmijn1974}. As different animals have different EOD-frequencies, being exposed to stimuli of multiple distinct frequencies is part of the animal's everyday life \citep{Benda2020,Henninger2020} and weakly nonlinear interactions may occur in the electrosensory periphery. In communication contexts \citep{Walz2014, Henninger2018} the EODs of interacting fish superimpose and lead to periodic amplitude modulations (AMs or beats) of the receiver's EOD. Nonlinear mechanisms in P-units, enable encoding of AMs in their time-dependent firing rates \citep{Bastian1981a, Walz2014, Middleton2006, Barayeu2023}. When multiple animals interact, the EOD interferences induce second-order amplitude modulations referred to as envelopes \citep{Yu2005, Fotowat2013, Stamper2012Envelope} and saturation nonlinearities allow also for the encoding of these in the electrosensory periphery \citep{Savard2011}. Field observations have shown that courting males were able to react to the extremely weak signals of distant intruding males despite the strong foreground EOD of the nearby female \citep{Henninger2018}. Weakly nonlinear interactions at particular combinations of signals can be of immediate relevance in such settings as they could boost detectability of the faint signals \citep{Schlungbaum2023}. @@ -437,29 +436,28 @@ Here we search for such weakly nonlinear responses in electroreceptors of the tw } \end{figure*} -Theoretical work on leaky integrate-and-fire and conductance-based models suggests a distinct structure of the second-order response function for neurons with low levels of intrinsic noise driven in the supra-threshold regime with low stimulus amplitudes (\figrefb{fig:lifresponse}, \citealp{Voronenko2017}). Here, we explored a large set of recordings of P-units and ampullary cells of the active and passive electrosensory systems of the brown ghost knifefish \textit{Apteronotus leptorhynchus} to search for such weakly nonlinear responses in real neurons. Our work is supported by simulations of LIF-based models of P-unit spiking. We start with demonstrating the basic concepts using example P-units and models. +Theoretical work on leaky integrate-and-fire and conductance-based models suggests a distinct structure of the second-order response function for neurons with low levels of intrinsic noise driven in the supra-threshold regime with low stimulus amplitudes (\figrefb{fig:lifresponse}, \citealp{Voronenko2017}). Here, we explored a large set of recordings of P-units and ampullary cells of the active and passive electrosensory systems of the brown ghost knifefish \textit{Apteronotus leptorhynchus} to search for such weakly nonlinear responses in real neurons. Additional simulations of LIF-based models of P-unit spiking put the experimental findings into context. We start with demonstrating the basic concepts using example P-units and models. \subsection{Nonlinear responses in P-units stimulated with two frequencies} Without external stimulation, a P-unit is driven by the fish's own EOD alone (with EOD frequency \feod{}) and spontaneously fires action potentials at the baseline rate \fbase{}. Accordingly, the power spectrum of the baseline activity has a peak at \fbase{} (\subfigrefb{fig:motivation}{A}). Superposition of the receiver's EOD with an EOD of another fish with frequency $f_1$ results in a beat, a periodic amplitude modulation of the receiver's EOD. The frequency of the beat is given by the difference frequency $\Delta f_1 = f_1 - \feod$ between the two fish. P-units encode this beat in their firing rate \citep{Bastian1981a, Barayeu2023} and consequently the power spectrum of the response has a peak at the beat frequency (\subfigrefb{fig:motivation}{B}). A second peak at the first harmonic of the beat frequency is indicative of a nonlinear process that here is easily identified by the clipping of the P-unit's firing rate at zero. Pairing the fish with another fish at a higher beat frequency $\Delta f_2 = f_2 - \feod > \Delta f_1$ results in a weaker response with a single peak in the response power spectrum, suggesting a linear response (\subfigrefb{fig:motivation}{C}). The weaker response to this beat can be explained by the beat tuning of the cell \citep{Walz2014}. Note that $\Delta f_2$ has been deliberately chosen to match the recorded P-unit's baseline firing rate. -When stimulating with both foreign signals simultaneously, additional peaks appear in the response power spectrum at the sum \bsum{} and the difference frequency \bdiff{} (\subfigrefb{fig:motivation}{D}). Thus, the cellular response is not equal to the sum of the responses to the two beats presented separately. These peaks at the sum and the difference of the two stimulus frequencies are a hallmark of nonlinear interactions that by definition are absent in linear systems. +When stimulating with both foreign signals simultaneously, additional peaks appear in the response power spectrum at the sum $\Delta f_1 + \Delta f_2$ and the difference frequency $\Delta f_2 - \Delta f_1$ (\subfigrefb{fig:motivation}{D}). Thus, the cellular response is not equal to the sum of the responses to the two beats presented separately. These peaks at the sum and the difference of the two stimulus frequencies are a hallmark of nonlinear interactions that by definition are absent in linear systems. \subsection{Linear and weakly nonlinear regimes} -\begin{figure*}[tp] + +\begin{figure*}[p] \includegraphics[width=\columnwidth]{nonlin_regime.pdf} - %\includegraphics[width=\columnwidth]{regimes/regimes.pdf} - \notejb{Update caption to new figure.}\\ - \caption{\label{fig:nonlin_regime} Linear and nonlinear responses of a model P-unit in a three-fish setting in dependence on increasing stimulus amplitudes. The model P-unit was stimulated with two sinewaves of equal amplitude (contrast) at difference frequencies $\bone=30$\,Hz and $\btwo=130$\,Hz relative the receiver's EOD frequency. \btwo{} was set to match the baseline firing rate \fbase{} of the P-unit. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D} Top: the stimulus, an amplitude modulation of the receiver's EOD resulting from the stimulation with the two sine waves. The contrasts of both beats increase from \panel{A} to \panel{D} as indicated. Middle: Spike raster of the model P-unit response. Bottom: power spectrum of the firing rate estimated from the spike raster with a Gaussian kernel ($\sigma=1$\,ms). \figitem{A} At very low stimulus contrasts the response is linear. The only peaks in the response spectrum are at the two stimulating beat frequencies (green and purple marker). \figitem{B} At higher stimulus contrasts, nonlinear responses appear at the sum and the difference of the stimulus frequencies (orange and red marker). \figitem{C} At even stronger stimulation additional peaks appear in the power spectrum. \figitem{D} At a contrast of 10\,\% the response at the sum of the stimulus frequencies almost disappears. \figitem{E} Amplitude of the linear (\bone{}, \btwo{}) and nonlinear (\bdiff{}, \bsum{}) responses of the model P-unit as a function of beat contrast (thick lines). Thin lines indicate the initial linear and quadratic dependence on stimulus amplitude for the linear and nonlinear responses, respectively.} + \caption{\label{fig:nonlin_regime} Linear and nonlinear responses of a model P-unit in a three-fish setting in dependence on increasing stimulus amplitudes. The model P-unit was stimulated with two sinewaves of equal amplitude (contrast) at difference frequencies $\bone=30$\,Hz and $\btwo=130$\,Hz relative the receiver's EOD frequency. \btwo{} was set to match the baseline firing rate \fbase{} of the P-unit. The model used has the cell identifier 2013-01-08-aa (table~\ref{modelparams}). \figitem{A--D} Top: the stimulus, an amplitude modulation of the receiver's EOD resulting from the stimulation with the two sine waves. The contrasts of both beats increase from \panel{A} to \panel{D} as indicated. Middle: Spike raster of the model P-unit response. Bottom: power spectrum of the firing rate estimated from the spike raster with a Gaussian kernel ($\sigma=1$\,ms). \figitem{A} At low stimulus contrasts the response is linear. The only peaks in the response spectrum are at the two stimulating beat frequencies (green and purple marker) and their amplitudes increase linearly with stimulus contrast (thin lines). \figitem{B} At intermediate stimulus contrasts, nonlinear responses appear at the sum and the difference of the stimulus frequencies (orange and red marker). Their amplitudes grow quadraticlaly with stimulus constrast (thin lines). \figitem{C} At stronger stimulation the amplitudes of these nonlinear repsonses deviate from the quadratic dependency on stimulus contrast. \figitem{D} At higher stimulus contrasts additional peaks appear in the power spectrum. \figitem{E} Amplitude of the linear (at $\Delta f_1$ and $\Delta f_2$) and nonlinear (at $\Delta f_2 - \Delta f_1$ and $\Delta f_1 + \Delta f_2$) responses of the model P-unit as a function of beat contrast (thick lines). Thin lines indicate the initial linear and quadratic dependence on stimulus amplitude for the linear and nonlinear responses, respectively.} \end{figure*} -The stimuli used in \figref{fig:motivation} had the same not-small amplitude. Whether this stimulus conditions falls into the weakly nonlinear regime as in \citet{Voronenko2017} is not clear. In order to illustrate how the responses to two beat frequencies develop over a range of amplitudes we use a stochastic leaky-integrate-and-fire (LIF) based P-unit model fitted to a specific electrophysiologically measured cell \citep{Barayeu2023}. +The stimuli used in \figref{fig:motivation} had the same not-small amplitude. Whether this stimulus condition falls into the weakly nonlinear regime as in \citet{Voronenko2017} is not clear. In order to illustrate how the responses to two beat frequencies develop over a range of amplitudes we use a stochastic leaky-integrate-and-fire (LIF) based P-unit model fitted to a specific electrophysiologically measured cell \citep{Barayeu2023}. -At very low stimulus contrasts (in the example cell less than approximately 0.5\,\% relative to the receiver's EOD amplitude) the spectrum has small peaks only at the beat frequencies (note, $\Delta f_2 = f_{base}$, \subfigref{fig:nonlin_regime}{A}, green and purple). The amplitudes of these peaks initially increase linearly with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines), an indication of the linear response at lowest stimulus amplitudes. +At very low stimulus contrasts (in the example cell less than approximately 1.5\,\% relative to the receiver's EOD amplitude) the spectrum has small peaks only at the beat frequencies (note, $\Delta f_2 = f_{base}$, \subfigref{fig:nonlin_regime}{A}, green and purple). The amplitudes of these peaks initially increase linearly with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines), an indication of the linear response at lowest stimulus amplitudes. -This linear regime is followed by the weakly nonlinear regime. In addition to peaks at the stimulus frequencies, peaks at the sum and the difference of the stimulus frequencies appear in the response spectrum (\subfigref{fig:nonlin_regime}{B}, orange and red). The amplitudes of these two peaks initially increase quadratically with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines). +This linear regime is followed by the weakly nonlinear regime (in the example cell approximately between 1.5\,\% and 4\,\% stimulus contrasts). In addition to peaks at the stimulus frequencies, peaks at the sum and the difference of the stimulus frequencies appear in the response spectrum (\subfigref{fig:nonlin_regime}{B}, orange and red). The amplitudes of these two peaks initially increase quadratically with stimulus amplitude (\subfigref{fig:nonlin_regime}{E}, thin lines). -At higher stimulus amplitudes (\subfigref{fig:nonlin_regime}{C \& D}) additional peaks appear in the response spectrum. The linear response and the weakly-nonlinear response start to deviate from their linear and quadratic dependence on amplitude (\subfigref{fig:nonlin_regime}{E}). The responses may even decrease for intermediate stimulus contrasts (\subfigref{fig:nonlin_regime}{D}). At high stimulus contrasts, additional nonlinearities in the system, in particular clipping of the firing rate, shape the responses. +At higher stimulus amplitudes the linear response and the weakly-nonlinear response start to deviate from their linear and quadratic dependence on amplitude (\subfigrefb{fig:nonlin_regime}{C \& E}) and additional peaks appear in the response spectrum (\subfigrefb{fig:nonlin_regime}{D}). At high stimulus contrasts, additional nonlinearities in the system, in particular clipping of the firing rate, shape the responses. For this example, we have chosen two specific stimulus (beat) frequencies. One of these matching the P-unit's baseline firing rate. In the following, however, we are interested in how the nonlinear responses depend on different combinations of stimulus frequencies in the weakly nonlinear regime. For the sake of simplicity we will drop the $\Delta$ notation event though P-unit stimuli are beats.