From 83def215a8e1529029febc6a936da7a2b05205ac Mon Sep 17 00:00:00 2001 From: Yuning Date: Tue, 5 Mar 2024 11:17:05 +0100 Subject: [PATCH] fix test --- test/test.png | Bin 23080 -> 0 bytes test/testfidelity.jl | 137 +++++++++++++++++++++---------------------- 2 files changed, 68 insertions(+), 69 deletions(-) delete mode 100644 test/test.png diff --git a/test/test.png b/test/test.png deleted file mode 100644 index 598033d3294e9c078c9215177e1eb5fb23be4d0c..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 23080 zcmXtA1z1&0w?1@9NQp?7bV_$ggLHRycOwXbbR!|%jev9tNJ&d~cQ@SO`=5JkOQ0g-BSRn%R4GX@B?ts62K=B$ga&_OH06{Cet|cVkr0DCzkK`AS{M(3 zkU*rwgjGD!4>HV!u%zDgnRA8&qe(%8g?57il=5U=e-b*t@$USxSW9QbZ)1Mknh7@= zB>7rcNKA-^EQcCV(moqjh9?J-znfWB5F#uzBv6QC5vr~h zH9&|YrXCs&1v=2>-46p{k{CM}Yye}x^m(nmP5Ki|kzB=oJH z-acG?`>P8_M@qWd8}VjX`^HD$`|uIHHxDLm#MFp=YxCd9$;nwN7Yusb(ptG4L~s$_ zV_I;1FsfyYQ)D$w&DpVB3G>lRe)Ey^<4YQ9>XZ+bald}u4yCdU4i1JKCNt{RG&D5S z);4>b=4eC)gTcaq_qn;b85xm%6>IW;e!9NrsAjd%$Ly6h$(q62@D+(t}X~Kuh5&D8|?-=-&#>0A0I!zr%ZnD6OV_^&j^+u z=9=A)U%h%2MaWz4va9s|{d-8q+2&w0v7ms<&Zp`$WH5^Y87-lSv9X3`W?%Z&p&;Yq z5Nn?(q)26#X+exVgC*!i0m3{pZggtvZx<-{WIr zUlS7gu6Ds4P@;vRprI8g74s1>k=|d1On$}`QVRb;@?L5`|HqQLcx`14YhoGo*Z1!I zreBkP|Ncl69h$@c94I9!3KK++VE=LugAy@t`)fstEKy?a- zmF`?HseNR~t}ISy7q&wJ{%{p3aaBtk z@tUAG3e|8Dnxy^E2N^XrHCfrAa2z_Gd5ka(YI=G;|EGtC+w*yD`EA7Mf%x~sLqqVW z_#a$${&G8S@p)bJPfVO2w?8Mv##YzX!a%^%Nf+?-0(XHJWPos(uBa%irI;;<_fP*H z*^1*X4o+{DjFrH)^B-hLE9?gxXISEMsc}!GGw$&`D$K4p12#pf&&uS1IGI%*TIMj4? z2_Y-s`Rwo0In84b=qctW=SVXdSW?+wJ8(yC_Tt>$ZuU~5ZszAH+}6#Jp$RR=Z=&Rx z{@EPJud&EB!HJ=lh-tDilVXow%9z=94$1S7jVoNs7qy4$}3 z8yYFPwWS3a1%-~5HX$+5k&9qjk^<9yp;5SS0t$kRjBLbS!;P3pl3zlTUou;$NnROM zp{1muq@kg&(@#f}^yAF{XOi0L=CX{%hILc6)gRMHVueE8C!yy&n)znI@~p3M*{R`Q z=@kpQhrS5;VKwy#5C;lQ5_8!JoK8{d81GpFiojB}&2{%1ElQh>qQ>>^)Z0 zYJqzNji!C+8!n|pfbYbjE!y$OEp6?ARrJ{^8XH_ zZ==tBu>B}LopIeM%UDY+KB_VBH)VP>AV)%IRUHyohjTMS#)`ogXkt=sXu`1Kb7Z;Btz+)ocQ$qS%ZXL=A({&9>a1U7~7b&kiV87 zQlmL~U>!0x{n%%kHv|1z9h=nc{`pzBGbdYK;(=yGEsK~)yz2rg>g?aYesq%=_z^g9 zQ?b_#E;hsZnU!w&xxu{BhKvsC`k2O{1u;a-c`8XhhD z!g??7;_Te)bL-^hb_MDw&mbnXhLKp2PPk8-l3tRU*1=H1<^}Rp;V(C|Q@0xIKMtr! z_NeiiR0WU#D!YmKKL;~7spta1hklM4;?e^u#qBKRS;cX$RpZ8Ilc8_BLQvi)DJwJS zHc=dXm6nmAq-+jvcGA<+v$bV_2uVmtd@NvlP88|f!}z2kGAIeY`%U^->8Typ{bMN} zd^=TAxp$akI*YSPb7a-;fSS{s_BAye2@>#r#Z-OUB%wX>Pc-kvSzX$<@Q5t~b2)EX z&Fb39yn9f%%PT9hl9H}gyF(4TzwGYszo(?6q@vQ&)C5%nalhNFR#ZepT3VW(fgwck z+Wlz3(b?JiVnB|RCMM!V zZhyDq1`1)keVdw?D4sjU$H%v^xw-4&rJ|x@Xqes9)YQk|2MRHBa`rLO($VeTo^Q*^$xTmBi|F@P zRWUm|JJXB*DN@Y+z|L-LY&?d{>>;22h{P# zAW__`RtC@1nV3PFpqgVx*8#P(qLwL_clG#K&}y>i{`z=(YfDR8`*`7(^Yzi9si|pZ zW+s?K2n07`b7jS0Ybf>X?CdMndw6(wqxbTutfm(`V@&#OvQknkY;2_^B~}&|HxIYY zuCAA*>Xoc)Y*3H}`(-y**WtlI@5{YJ0E2{ug(GHv=~Of~-%poogS%Q=Uk4ut1Nqlr zukYgGqO3fTFOy6_NLXD|1I9C&!eadNaEpV3116~J_1@kdSOfM;Et*P7Ynz*{pk?4H zgK`au5Dh9SD)=(c>we?=`ud)qpVQINy_Zk-SnEM5bL~Wt+v;1Zu~oV%j9!r`QOkn$#5`C)0_Sh5@uI$ z`dbdoyRNdSXVRmM^|m5vu9zf6rl%{AEnY@B_k4XcU!VQQ4%&p1LQG(7?KxPYP>|c(TNV};jay9UK!^3;sOIL|gUL)V5D^g( z4A@jY&$Es|nC+1anS3QE$ockY7PBEVgaiPu4;2UsmNfitF#9E?rSb9c-w{f|gHiVQws`i6hz}0}z$Z2iPRQcq!{ha`MyPRL zl%<6QEF2sZ1neg8Dp1IhCyNax2On=TmLm=ZWG_(pm(zQ)NYVKBH_b)|(krK-6c(^B z;2|bvW@hH*`kHYpDu@jfrmgq&B4ON^t#;p>2+xZo#7mlN}XR_CakkISK#@=K|{F(c14Sbii$Rpsi~_27|-ukR#~}Rs?G#~#7fpRG=L}m zm}yX4RAiuu6egCFmZoKCSzJ>Cs(Ypb>dwl_$<67yedhh3Q0E@&L#d7vxv%D|CFZW9 zbB@}G`krp7d3%FTnA?|)F29|cC$4bVn2y)BDh!8S+EcY}M^WN4aeIj6xXsJ0R?93o zi;}sCdLKP1D&zz`^Ril~eXT=i19FN^@PsjJa6b2~e@}K=b^6uy@V*NMTfDtOCD>G2 zv(_w<*)qAV?n4mC-Q8W6=c$KJM2!t}paIxf0RIpQ_;5iwrlwT?LCq80&3PCilqKm# zT+kB{5rMt1y1L3~w}ANSmCyB&u8z(EU{B=a6jZV|@xr%8ZN* zu;<7hCU6j80nGsQ3=Itp5|ERV)6>&KoOZp^2_q{(#P3DH#I#H{4|;WcP*7K%RQ%CW z>)@u-`Pmu3{&|&^m1SkGAUE4t{(gHC9@j^NNvFm2^&C`GzkmNOINX+&l>9e0XQXPJ zBD)p&j?>R-uO4z-)x?O1rohVl{%y*6$>4t0lIESqlalhbb9(h>k<&?`mH$~PKWUvfZU22j5l+xmRU!R}93Cd7+cQ-n* z06iVued~ywo!xg%dp8%C$%zSG=dGcsQuRXl46oZWQ*rU$p`jr#Spd>jdp=hi_W_Fe zJtc*Un>#Tnsj0acU=kP5YZw_Bx3;!uXlOvUF9+9JY5`Suc5ZHOe_ujOES<--q?Zu4 zUNtW_H$E}Z)X2!l!h)KPt{HGBQ0<+6Eo3c{4S;lknH#Aogjt!IuqH~U(-N#-S*Y+G%d%Ae zH<9y8)8&bicYAsxOQBuh_NisJGh(DrpQ(u9K*QKSt-rsYRP;-j@rL{H(%7!cLF11f zKZ=SFc`M4w0Iklg($><_l9UWP(G3(j*_)!{_#R2Py}1eQ+Nh5ncmMJJIyE)*pA7ux z*Kgi{$N%~Ic-epS3;21t-TxV^b3P&k5{wg%mXwqf(E0!h@q4(b2AF~bBOxIHU>qa` zFk3Wqblp*&d`qAvgAx(KFl=dMxTxnLp__>8o1Ez+xvSIZ4W-8 zB(h&v<$I7DRM(3l-m)I%^tN6>Un<)H^=c&%=pl!#NqhS{pAv?!{z$ZEY)a7txA7|C zg3WOmtzVyQ?NAgl`EJig1;LKzosU^g1LFht#X0(PQArSE;^D=g%Do$t5~(A#eP;xqP!c zp)POy)bUn0m+8rHQ+?Eo{?#*#^HvV$lv~Sb3}aMD32jpHL$ac>GC)90&fD_7zK@{m z`ZumTYJ!5&081(L=T$0+?eXFX*#f!Kx&wQL$MAr1)*f|yr%rRHd$hMUedKFU$ce;Kwg6vWM{^%S9O@Hz}x19wn zIgJzDc@|F2(a}-+v7L>LgxJ{AR^>*WBha4xdfVz1m6XO8+?NLj_d`cN^rE(O?_H62 z1SOOg<|Xp7Z5FHEv}wEzaW7a`%hRy)8_>7>CO;g}+zxm$hLEwO|4nKEcD8WgACp3z@tTie^tT3cX^tMJYcorABGrPaBO|KjC9dO1yi_ zJZ_{{G;19ByUHUfI4-z2==b6FiElv5VtTk4e6UBQBqf>I*xvEDP(w&)X^|lg_gA)%fbZYG zx3u`^>MnvxF!<-s!NEZ=qsrywWrN*fNpbO8K|dcevjWiRa&nB6l;VPdL;ziLaykow zL(!n029Px(JREdDfM$WJsHms+v^D&_q@-lN!5%0HXm44?0Di^8Oo)%i!o~eU5>{CF ze(McTLN0)GBqulg;|B?H2mt64e14KD9h8RA(D4k3u7#$$ef#Y#ZcN^{IDgJdi9Ng; zUbx10kMfMt+`rf8bnOv(5}-dii>jCTH$B>J-bbmPW$5VW z2nhgP4S*+_T_>izpdx$hi=f_^$A7_yms+!Pc5Lv5e) z^Yf2p@Jdi%0+Ev!;69sR&+{!$504hVNBnc|v(r;UW;{s1OIMhk1^VEMLn8pHP!K>K zua1{TH=Q);rvW#sBI!J?_PISMSP;6`pTS$ssE{W4XB>BYE79 z2(m`at2HE|2s{DbNJy9j6$$kATOWb=fxL-bpb!Ff3iu!-;BtSa-S-|Ba&j=sW^KI> z5F%LmGBPqY87b-Me4t4JdWixX1h%DoCZFcyG6m!Wt<*e-go4~UE?aeCUX0@5m@a*v za?7{@oBQb!UU!x_!#p6a{rM+d9FAG4Z00}v$=vk*0CVpIisJQGA7808?Dn5UL=GsV zq%zI1K?ily$!Y~+w3#@#y1Os`w0;8#7&b$IgEQ3ITJ9bmwguA+2}n;*f5W6_kcrDh za0EIcSYlxL88PDlbcz^6uUXU6(E%CUgg{cboeQ$E^v| zlig7WukfnClW|=C{c|&??&r^HE?ZUADIh5ir?Q!YX{1xH&}nvKm^K6j0MvyZqj!DZ z-m*EYbbt#jwRkFPY2iEb5VKn)L`6k82tDT)pY7vgVj?0mNKs|wGeAap`S~vlDR&I8dNIFvIJna7pn;@V zb~;Wh)K;gGJFoJ~#tvl8au*XnK^JlV^o46IRvFCx7wEHt?ys2cvpHNoVVT^wgL*wp zIL9KD{?+$Yu-3UX2_^;Xb5?D1NC=$UkBW+?i`|Lo(cVb}$lqpndqeEZv^0D-BcLGw z#MAQ?tNQB8moH^y^alx*xf;yd3WhtZ@S)wBfY54a;KfRAoO;Mo)keOdY(WOp18}rT zUW)hcWu&Ezb}1o{SMa%cc@%7H$`%^SZN921D#4H6-#u&$#1DVxoY-}tp(%cNcnAuD zHng~yF)=YA;Bi6Z1*=xTVFlu5V_*IL_QS^*YL-XZ{0l@LLwymavxQ70#1vfM9faPKO2xJXz8r^@w zd#{k`^Zev1EghbhX9CdsWU;dUDs9{Du^b5oFcZF@zWdS&`~%P?Ks5zciU!J7Hsb-) zhm`JBuso>f>0iY&YirL>i&{}oP*nYufj}a9F{38UN)0tn~Gg7mT=jY zix|$0q|DMJ_&3VYKhCF{BDKU3!%e-Lc=ZgB{3L%>=5We`!P3Y!0`eX*Z;zE@Af-y@ z(;gsJ?S6dM{D7E&E&|F*-G49O#%OusOJ;wPm6NM}q&l%D3kI7Ph!6CSr{EdPJL9mj zvLXf{KE4O4d{=ijkT6*0i&SU>g&bU6d+S%LciVPt6xG%HOjrTfhk}5T1^R$Sb9%bY zxk`}8ORY~@ZtKBy_&jUAe*OCEm#37J6i~;UoN55O1^Qy-w{PIW zwiXse+1aO;y-L&xvE&jj+FDAAiMcsSysf9F=iyT8JiyysU2FZZjB-s&HB$CElq%-8=b-vDQ8B44>PU zCwah7;0A-^vJ-k)g09Oc^W^66r-V=EsJOT|K&}5AZEee$*OI$eTY#9r!tz??Yu~05 z9JH{bBg^N`4j!QWsi~?O_JN)YM9D4|;G-EE8D$FkGZJt!F-7+pA@Oc(Z6UXlQ^lvI zr$1cmDjU9X^Ypxpf3E|pI@5;EL!1aq5h{&0jQvOGV)~%m>6;2G6S8(*~uxm zc6(z(?BmCR{CxK3SmxBB@bGW|DKWnWhp;3kCQeOGBJ!TxJ^&7RplKx~^_4j;b#N0G z7uN`C6l?}rDlBVzD=X07q~+ytM+Oa;(?X;C?=FgOH#=6gmp|OReNkA1sN%ILsvk@n zal|j9?U~uIJfpB4k(OffceVm|bbPUGMdA6|gidm}ZnexlvDccFwOu<9UhB*IAb-y4 z?W?z)2h~~;o8RxDF)*+b6l7L+nBF-I=wAmluaS?EZ8|VQ|bd3A)24J%9L4ESri^g$ZB7%5PeW+l(?RQ&saW zO(k`Ilu(vS$>lE+XEbm&+s_8Jw#B84IDX0H<>et}O)NLL%Ist^3IKuk7N9k(KBjPt?Xj~G87zr+eXF5=hGAR8Uhi>s?5s^ z_T&bO+c@e;mzY{3<{5@kW-wZ_h|8Cd?cppvuUzwp9RdNxNYwsjp}zOuADRAJyr@=| zwF|%m-YYW=`_}RCafZLQV5gTH($LgYRaZ9#)DQ(F?KK@b5ue6StHhH{paEzNU_HxO zSzEU=u3=ze#${wY0;~m^D?l#M`0RW8SU_ymI?XYv)=1$XMXYgL|2QOU1MX(g^csT^ ziF8>=q_q7kZT(x}^CYom&4leZH2Se%-0l0Cn|X%6r@r$S3Q+;%ucn+F_xB%ixVU$8 zPi|*sB=$VNOir z!GrxLQ)VeZ)SaE3h(X^7xG`lUlC!e1GBYc3bFr|oKX{xN0Q3hqZfI~YaL7T0z{kV; zyL$mRZc9rG09ybb85?$p6iU;xh8Mg2e-@y95lCS`SLur;W}>7_%*-Uh!BGXs84%^u z^*;Y&$3H+Z0D}T72Y{->!^2aLmcMXi_4QYO66suc7W?`NZ)pH71H%I47=Rp&)1EDZ zpSE%F@c>1Wk(2Mb+yy|v0R8de@-pP>HBei@wm^OLYVMjc?3G6-+${$p8WWpI(0I+` z(9oOk?c_zaKfG2E^GRLqzCCi4cKh+BHtJ86r1K5`Dn0ct9nT&;lUeOPIEw4&)Vcx* z<3AU`g@@1?0~n?YRfaH;AAx|C#_dcseF%_ct@$W$OpA+)tLy83D=CXhNPs!VquVsF zwvK?4FU-#eN*H(z^!Gzvl;D@Xsk-ds?(Tne@HV>7+Qw$Xf$Qb!lYAa0B)FvM@$sOL z#U&*b6c?warG=})0B4puYjo!CUoP*fgEGz9Hg{5xG|0{l95w@Z0yhGOkT6m;(5!a$ z>go#Abx?9@&4vM+=(_9c>1k|i1ZIh}lvE@EH!D4T*{y|-!1ByYLviuf0TZWR3&2<5 z$rv#KrT`lNkF*av1+uARQF$osuUfi;>fkdAK5E~XuMfUYEacRFg_pj@zAr-Ls6jm= zhq`Ed2_rGFg{O2`&r`-=mXIme&O_<3N)S)eUH9cje?!ALvEQrcKhKwSYY+%oRMnIr zA?_^~*p;AAF>2O`Nl7972d^2Ey|+gQ$zsFcV9sxS;6GObuW4vVmi{VobemYf$K(0w zQM=v-2@$bFBrQKL4;mUelg|?x(oxfB(Cd0}aWRs?yS=ss4f&*~hyi zVZp`Kwfai6RE?0!!2n2G5C{;w33=TRbKTO|zzzWP3oQI7VnHH^fr&|8dATmI#dF6h zbbc9`o6qEwXl#Q}${9I^*!Z|=S4YcL=WK~DN@27%pJ9uq0Z*krhI#etrfJw2>7!8G zBI(kb`Nz|?z4gb68KTcDxX9DfA1TCdU6 z)Rd{Yd96m-Y^}L5AyDSoKYV~dfEd-)h=qWFP$t3;=>V`_xfGaZK;g4lZhIURe1?Sp z6J1?Rt=>q(%4%<^)m!@mk)R(z1SMFxs}5YCV5I0Cy8u#b5@s~S2Gr5UGkU_}mBR;v}ro3JOODzZGpZ}3nQ zoRwp`y=(ZzOaIrh$@Y}(F-SmsJ7&H@J%1PYv`Ajp2C=g?gduIiywoyd3Hfc`&1xU( z!a{HnVJR@Bl$6c@*6b2UL10al#ht!1W=+)qxgw!P&)_MZrCRf7=Kg|2r2z;8Kxy3F zw@Yol&A#_nK%)l67>EZkFf6uQIcbx`4|L5{Kp;Z~*y3Ns)&M^Wxt6fK{-D?BunJ5I zTU%RDuA0*XVHdiC9$+z>3jSEjW2k%>bYuF_g4TKvwuB$Vo7&sZH0rB%=S+&VMO$Fb zkIrp9R~d!&@Ag56uvAZ$0|o-5ybuMX_OAd9wX{3{A+Ne>waxe5C=OK16@b``e8I{f z8~Q{WEJg;*5kM=f9#QJk?w7tD9UXx*&})HSLFzw)GZEb7- z0|gSxf;9Z~@hU*nyStV~Myy6XVf@}#h{~S``Mo~S&n|PlYuTnOso$$QA;b3V^Nq3CUyau?z#eP@aD$D zkjAo%%lYF+_^i>&;^I!Ku7@fu9dTd^Fmp6lS66#pILXV)J2|l?50)1fzy3&ge*2)1 zCD7Z~hw?f9ds-TRpVjN*<4UVP!9VQnS*?CUKqEN0x_as~zheudy?KkQcAZ;_m?T7k zf!N5Nn7HkpaFDW?kU@ok7><8)zSbzMK8wKSk)Bnp=rfB?A}9_}FMu3CjLPEors^^j zr_^weLbQCt#f@w_{qa_@ZOZ?+vaMcS{FDD5qBlI){!Cf+59rs1!VN=kXS~Y#`dRVu z&Q}MTKxQ7P-+TwMbYYOV*q>)<@zk14Hltm?Hze6;F%LrojXrB94>_Z3w%OI@(&SBy z0Z0ASE1T-eX;+ZRdyPshrQ|mx0$9?=gU3G;-+xV+R#-{ySv=gd6?s#;-zn>9YDX{# zku1X@JvXd{&PL=pJp&gHFlu9@)`7IpXPWiSUzDLdYjUvGhWMGpRK`HEt>sBdG+d569 zLhj1Dp>dTJ&^%Y7v`xS9@NHyz&eW^0WAs16h&q!P1@mvuAipO1oETUM;kGOGUtRy4hu*hPTGh^J#cy&pFK>g&Q+Qz=t-5@qaF$e@~Zs6@R zjDSso{pI*)K(BG=_D6{~J=_=dSAX3o*PA~fC+Po3q&sA(t1Wm7UzRrUh+#XGjwae+ zI3w|}7=xV+2z7825@h8J;i*JysQ{|c^s(KS0WhU4J}13fGjy1D-Dv z_O4mhR}A>ifjFMWvsD^pZfe8375dD1^MM#)HF|8=h2^DfHO-6%v>1p5`#Ga6r`xr+ z>scpO{xnfBR8o)Kf9xA-{dHC&Qb}ci-54=vB{>6)jGHGK`IU*J^VX;zus3 zk0wR%RXq{*^O5sJiQlQ2zP{pi%>IZe?m45`g2t`MxrA4Rh3yv?LO?lq@Y4yzEeD8# z!p&ye#1De5DG0bsKF^y5qM@u4d7bYopbx9cj(QQed}??ir#0{&KZPGr*tQIG@qKwS z4EFnNer>;tt`at=ZzFDyh+77+2mGckTQH|~qd)KyY> zh*PEM<Q@AClcl#W8 zYalPtWHZYO1M_)d!QO`PMXDPgcLKT?Xxd#}L;!;WP7k~a{T7evs3>{55K_`bd_Flj z%S}VrP$t2;IBiY(K$uV)2ZXY*vFI0-cm5nQpE7-1woR&PnW#xB)8(8iX2z1|{93B! z`|;_{7&*W#b8}jOLRBUMaX?wE-@(AgFDWbgW{Qo1g0eqRh@RAH1GGNSP&Bdt`iaD6 zFX{ctOsr7PZWauveF)Vbkeo3HmmqoamzENBxzxIDdmZYTv>pbSvjLEaG((A@0p=vdACl_62O z7tP{B+7EA@2;sS?ynH&>r>nV@bN1(SjFZN+amfXTmQ4EP1IBAl&E_^=Lb&qy?yGln zZ%0sdi5^&VAS?F8@5uEOu7RQkAl*UZOZAqM13^?;+GkaVoG}aVY5yE8Eor_F|JJjc ztuX~@nZYE6?o|gQBqX3UTC+HTJdn6}NL3q0cutp-Y>XXpc0<%SvBd7owK|b@A(t_~ zUR2Fg@cyK>kqd>pq2Q!H$>}_MNX{{-;liuuqLBWD$4?S9m2xJ01bQ-c$ z${H9*5fmVKW3gli97QN7D3DnMG>G5#uCAqpPYo8B{zWAv;gS?-TA(vg&P1=;D$adg5-LaUqh(H37CR*g-1<~82o>T>DE+-coS^GjJ1tlaQ8JWSG8>X=rzkrRE^^T`YSzO!{^jZ+@NsfsDz4-UqC1A$E*>Tv& z`5}EUd41&B#yNFo%-hCZo3rH)AK{7Mpv)Sb(Hf=3S^}{3nu&;^e^5>A|Ck@i)*6u^ zo7VOB4rL1XW!Hq-7!UwCjjSjwZFW1U=~ue6)PrkaEhYeJlhS8X){(q3+-T>otf3X!^+#7?lQ6GBVy?UIKsvl9Ilt9Bk$z zJ#(d{rC_1}x}8CLd7|8iJuppYpoNiOBaA5wm`@mzbe_uLH+0SYF5mzYwBW$N-&5`bTvFtpV9+ z7zlN&a0-H-h5h6K*L&M+Zvx(pqCbnBX2Vlg)C{u_wVZ|j?1nc>h zFdmp}w#2^MjaWIcmQ$poG}%*L`X|sYKu3`6`y7NLd8SP`Ex;r__xd^cXUVUu+w`MM zpI@O+$5kYYuSDz#^AkfGae5Q`X3qCYpdj+CF_U9A2U>yqExqm$g2Rx3x^#p9-0X~oND~Xqw$bRs40yiu5$Hj_4 z1xeB#S$YF&~UfQTB8)78Rj}1cU`eL^Y;16%E7a zN#|R9kb#y+7j*bt-?=~1-q!ZCdnIl5Q}Rp512#Zq{IE6E47LM$9}j%#!9ZW`FMyqg z{@XhplA04U@{&%?Xc*s3_H6+zqSQC!C5EK7FCV!bpuSgL_W)TrrH!iCVzYvRtl#By z9^Jr=-ZW2$7ifwRI4%Y3dZ#z&cskQDuYSTxrYwmi6?}QCM4)60PGKea_qk+ZjgfMC zc!&lK=~r?`PN&=RtW&NW8H(J(BHy|701Dx8P8V@Ixwih_AumB_BNNavfK=5x*2>XX#6!3x3AvQLm zslSTv_LhDp6ha^b&}j23NqE^D=#8V&;>DWyHLVJ0^Y29yu26MTmE+%c#gZ}rxV<|d zxq~->dtBk`P$6)KHqAN;Y267Uz7O{1mq z#S7pGambYsE7N>tL^oYI`=kPB!-E+~*g;kCxSQT4g1W&hAp}x;`S%B@M;7(b0R8JX zsv(1-V=qz;FB!i_-a|@AVj)M3Q#pqvzh{%BrJL!=5`r$SU!^V?q z#$~SL{gZ>|jrGgdPGd~DS`I#yM7Jrw^Ldg%L1yG$Guz9lNB1aJ&8dtyDl+u$IUN1F zAo<)Vq%pr{9vB37yV$9$Vj*+Q-l~xO>|kkmSyPpA$8B7InU4NNATOvePph(H)=6MU5|g-% zd^WsbvtmFXLKOYj%jJlHpFFL|+l!&B>(oRtJ=f(;Pqziv9mOlPo*B^HDT#vzyX#-h zTHU^w+49;FlYc9@>cZ8Pzp`H?3hRCAdBcv_B9```fLPC$7lXasy=nOF9wC`Gu zChYeU+GTK|0p#_KzwU~3nVE`3F4Ewhw21@)70|U`&VrqxQZVS^6jW63w72fs)>Juy zQd0&q!Yk4<@Elw;SdxL)6e?XXgsy!Vn@JZ%)AVw~G4`V4oVhpff9)TpPRMlrag}kN z-7IhTuccsndvV!!TnOesubiD&1Li@p98#36a5=C0#<867prmR{?6s!=^cx+L6s$ok z$KT`rZ_DgHGoW+!Jq;LehM7hkA%Qb%dF+UP@TCXQS=lmC*yv#(DwMH+tY8)xC}Br4 zBu%X*Z)WpJ*$Lke-~P<%WtHFr1q}8K#-S#>uGgDO{MVHFwZYpYD2S2=s3+(Uyvvk9 zY}dRG@UJCw*e%9BjZ=a17_7ehdQ#Jh6=S( z#V{fa=ZLbY7S6y)EEnh1!(Y|(we4zmdfD?`8+)f0P<%mzj`t5HVxBh%5bnR-Z-nPs zbA=FB0{j&-8{>Wy6UbU=i>RhVP(_?#JZ^-!e>=PGWLmLD;ZX3+cl_ltw2 z|8%>Z1;VInoJi2Sdy<;J2BURVJgwO7o$IR{%V~H--?$>J`r^8736fF=jSzDg*Q}=G zN6!={nTA;=RSYT6BZ0#Rz}&l$1*hR&n_%(6WvQwbHdG1vV{lZkizi=7AP8>wmSH8E zK|#LkH`)`o*m*y+dCx}oaLUZU9vKRgv?CU0g*@C)s+!k3NS_5qqTFY{lU3{sZQ*ob zu}WeYJ}J6C@~;s^Yv=2J6V``SCi8Np@z%AUs{YLQjaueHKMey0#~uEU`|`ip*)1Pf zOf+zFXv>|9c9wGP;(jVLJPLG<-%Kw9<!O z?DZFvC2>`qv=C8;LHaE?|5y5WZ9E-*k;QSbDO7u8byOPl_4NqoL^h8;+S=N&b1;D} zGI?bSfdTrCHmJ6zdM!Du(_n*PZiHu@1bYv>7OZP%{D1;kDQ8m=m$sByI0}I=N1Jkb zy)2d?ZyjU`@)IVL-Ei-qU*jlEcf*zy{K0x>5plh;L9MZFfU~js=^+g9NKtzsMZjo* zwf=zpq{w@vsc?Q-eVd*mD}u~uvQHk?u1S=HOmxl>gv-=>$JkqVnySWv!h>J+-uIp} z`ePigw=A)i^etX*Om4vO%`f)gX#eQi&de~Bb zIGd-uX zgRmT?UW-WWmvO0t&Y%vS!*Pkg!%eUq9N4sC7yVwtByTnTdy@&`yZb*RikXe;nh%G? zM(&pc<}0mpf9w}0Y-2rH)-p@UK`RIhgX2D)rDzQVz72|otYPk#515<;{`Xk^x;6i5 zdo+sEp4ma5&jM`xj>Qy+wu0{ z`RROedKzS!y#a-O2_#-0Wi%cJ3yDRgUlZh{x%Cdn4YSLnzA|b#8S4$|9qmnc{Ulho z!LH!_a9MD+IBo8&$;Xlqc<2E5F}O6yyZ}ob9IXMFqE{#=)??XqbuKF%frvrCFV5t0 z>{?zXwgwcZx~2x$SB`FO1!ETAgn`jdEa=xLh-88z(LZt&J?*{+oM-7#O9nWg`&cC! zbg~iY59Lt4C%sP|)YD&NHEHDo85p#T>kCg__!Se-?kWZ7#pgq6-^4<4e?&VNB8a9(D4KToo0kX z82iLSHQ8M~{a?p+)pt%dF6bEvf-Gvr(DxmFS8Jb{R`ds7;a`uIU8@rJ@+j_uykGU~ z0rB;{z*=iJ2zxs?Idw&4q^E;Whn5R96O+2Ka(5Jn6##KauiXz7687yI0l(J;NJsz1pDR!=5j7V z3)^qjGGx;qNRp~D5|WaXid7C$KIXEd%JO7XfC%A%ft_8OqeN-&4}>F`=a1Rgyc{M3f*iqi;=yKe(sbYSyGS`5eXW8s&N6KJqQHz>tste zaMXc`3UWsv>VpBRsHg~hZ72u^A}lx;280P1NRi4OID)5zh#UEFf4Tj+(NSN(g#fJ( z#JVqU_%Vs1CCFm^5IWZ!xcGNbD2JDHd<`Ym4@JW#rUU)k1u!8EUjkCc9;zL0?kkFh z8lIEHzmN{Bq}*r9IX0KfUf?95AQh#spFIm>9nJw0rur%@ zEbQUI2Lg`IgM-IF;sDAMI3Mon;IIQUK!^d5?11>Ju?@5-aP%h(oF>mR7CM3bQzqB4 z$9fm}`N!-PoSK(82CAbc#qV)mkua_`vu4#TS{U4r%8`9F)BM-fc6VT%*O$F<2)uAya1_Rhf5Rf6v zC@rNUiZ1*_ZE0x<5(k`YY}Nv2;HVDB9sEUwe((;RT0^Ul82>Lq<2_m4=MM6Nef5Z5 zc+hUjo5jPh>B260ilyDennH>m8oQYkDy#LZO3h}WCVsYs zPT}@6Q&7Z8g21|zyu6)-MOAh6>k=t3G1&S`Mj9FrA4*S6rD0;SH#8(&Fa$mdaEBzN zrPrK(0cWbGO$IT~JOD!qAYvd2%q;Lz3VQjSH-LP+3BfXD@by2eLPkorVW%;dPf0nq zeehdvJ?`LJG zs&SOomRkSXrP;JHFk;AzROlO9{WDtW_VW&DOa+p3QUT+?*@Vj)#2H$wc7<<2fYZv^Y1o~m7m{iQ+@YeY?2E*Nqi@m77^p-i30?FmB z^uM?M>FRRtVZ0!|)F0z#D(#D9lZBUx!VF7kd67i*!l5Bf`U&pAX0a!j&FLI}VmH(z zY#w+=c!d?abqCSZ0v)#D%k%tj_uD_$bhEVAa5Qn=+8-GsAT+lRuen^)( zMaaB9`qD&qeLis`+n(BebotbMo;0-`T11Ava9V5pi_@7mr;805L$G!hu)HL(GCH#a zw~lXuFa0o2CP4X94&#FAn^X7o!6ICH)Ev^-%o}#82jyd)BHGsI)qh_fKQef>2ozFP z8k^aWzrgqk=*$x7bCfbvfti%Ex?b%sXOiFE*!h2y$fNDW%isQGu~Me0z|mg3BeKhu z4Q@vlkxwqN$-w+@QF&+I_49&p5 z%z-Y7NTt9e9Iy3Jro|td)JEx$SD9#ig(}EqJtvKXH9Q~GSso`Wl~*D}6>oLP7zkf( zW>EOkM9#>(aJysCuB}r@?mU8S{WSrlVVQp?>{bg{ zWU%8nq!-f#L#LH(ZCRJ3t-acDF9d)|9743Zz#}h#Hz6Va80NjwdLRK>ylaooGGO1Z z4_*D$-vG5mS(1451)5vUULh`c<*mp!5S|C_gn z;Hc`H(UKB2?Z!j~XALq4LDlcb;|Mgreg*@V;2`Rl4fJ0}50az%Ju#v$ge>`si3JU> zjn-yrz*fj{NEy5M<15Aa6*T1dBr560&6>zwW83h>SsK#N|1KYuQ2u__;?vY@IrI@{ zeHre%oWP`_{aiF;Zi*E|pvBR~Y$QwA6&d@AiuM8t1~Km)$u>ViVUY@u!)~UZ##9Ia zh4TZ@(@{~v0yq4Jby+gr#H(Zhz2(mTeZmnWh+0-}zK{83Wo!k&_cHaUhj*& z%;^!rfI|!rHxkXEZz8|hyM0Q401fY#YfyD8pIYdG<=am#=xDQ!0*`a$?$^DXKxabC z(?09v-d||QsC~f{7<|~#5`pGi4Q-+#XzZJl(M<>-t4Y=NzX$9tKH`lgpZT%9{`mWd zCJ6T++cmD#nZ}YvHzpEdIu+69y8vloBX08ZfDsaR@XH-&V~BN{Q${PE^KS#tf(BK1 z;5gEX(HFXq^94L;zBM9^PwI{5xB6u0|NDAl41!EnYpiouzIDwL_j^c{H-=pkz<$WFp$t8(oW2K}~|edy9%HSO_is}Qp`*Qjk?4ykO8q;!4mpuJ98 zAC;1RQfZG2BsD&L6(@){&}@W=n!4mfKz*^$2XRWuLRu$c{KC2d!a|`=c(pOb< zz`w`soQS*H`;Z3GwS~XM$PFDL8~5m9WQ3|W9rr5Ts%zPI5mWPSc5A(tGR&6sml9%As>p_`5HF+KD^G1+oxhpR1}@na#T zjdfTEFVZeoN0|^UrclL+6+3!UvRqc~Ng`?s^oc_sQtwh)eB^*eHkXILPDh6xrV-P7 z!=F29{jq*;usq6wa{EO(bG-|$UCSTw_Q+6@E<#0{o2%Y1s2xn2 zjTz3`D~(M`{Wv{MT>c6?jDod3vu&uzzk*P)TgJ|opdcgJX$?BIrL*tVEk-!@Vqd@v zM~JHpa<8S$7Q~KFY_~I&keH0e3%SU#RyH50x3K1 z1`bQ(hJ&uO=6b~ttF{*86-8k9;u)n_9d?akF&z8 zfFH?#j0fzbhlG?SB#@TInm)jJhq~dB&|GUl9xjo=r;vqhiyUj8@l|y= zYn%fIX)S$y{fVh5E-tRz++4sOH|Ll^=L$MYI7}&@I^|d!aeXBe&a)t-cXD<gbn)qkiON1xK+LXvSRw<%}H6T$8Re__fbJf>(W`m#<-{~7B7@Z2q zW<+RT*J(T}I`-c2r1q@g`=VV?A6u;w+*5 zEn^WUpmm`ZkH^FA7?PHpaVg2kE0bx?`Ige-c_3zHW@3VXhZ(O3!bOq9E;QCc$nU<@ zxX+!l>PzYKCxdr-1_pQp1<&&XbPHnApIptrkXrv;{bmYysp0EiU_x~u22$;C#079W z43`wh5A2QgHe+n&R}2**Il(vJj?MCB1DVV|>TP%KV|H*by8)bB9=Kh&aG@0WkgO8- zK~%kFd;J@zU4aT>qNj)CK1*8r(V~*!8;b%7Y= z1vS_G`#WizW^2cwGy46?;^NvFv6Cx*+%yd`m0bU9Za(hQ@+lqc>`Yu+ln@dU0$>me z+I&~;GLQCDAc#O|s-&psP!%ktskuE^>I>>K>6F>cwK;$wUhe|=?G8$4fcyt8DLgYY zGz1((QuByWy(b=YpyGOsx~$FZ?QsjNKchcfO4=<|KhPOcXPx@fu(o*FB;p26Ys{%_ zhE@ZAcY7!uqmwAV!b=2orfy?6&)3=+E?L{pUsr?w_<`t6U%z?fCQ#wP_Q^TbNU|OP zs2Sk-hr(@qsG2}fxV<(?1VmVxJ$9;K4eS=)R>(E})VOjbd#pB_IwIzNid$z-mP!gF zu7RBD1!D$~yhB4nCYW^|&POAWe~%@QeM+m|y<3JZ26RwFM1;b)dENms2V9({j7aV6 z-#~9Y&B@fX8m@V2YRvXMvXJPWsm^fY^JT^0doB_%LW1utUAiDRoHpNnLi)zBReAEp z@uExWb0vXBC%KtGg4)}xH*eH9F);x$H-Mw?1AvoXTTT^%L7~O>+%yNmDxtx_2hfcW zDZDcHe&!h)DDlBX=A8E1mg z;_s`6Xm!iH#q7yF2$g8p)HM^FoOKrz2#RNXY>|b%Uwmv|=aB>z=FG-&Mv#*ABY)`A z92)a&0=Y~D@QIU|f7=6gltBP%665^iY3cOmwL0Xq8@4CwQ3(P;b_uM|vitk~*6yu# zx(#RayTejW=NRZAWu--WlGlhM z;^;ech%h=?=01rz*Fi5Tq($T8*c-2Rb!Joe#D z5LSy=Cp*oQ{`z;D;!h##5e@0UbZijN)HP|0%YC$@g+W@deS3us+;`>u zb4oQ_#$RxbUJ!b37A0Ym?39*RN5pzin$G$?#9e-Nc8v9HTxBHi6Ar)o2>{{T& z2>(@&Tqr*J>xXPv4Ykvuwn6-Xy5e`UXklmyKb?|m)po@_4ytaiI(tuOn7j8B=rmFz z`VnzWYeWz+6u)zbDnA5|5OpP>CWisPz7xFp(p=t?YS@HwUnR1dZ|!K|F0Y1b9s+gH z7-B&7L5V3#Q*!1v{k`e7?wqC+_VoAhvS$a?I#Rlogm@_RX|d%m>TDrYp$hJzc8AK2 z=u!Sku+=wzF*TxRWT&r;E{sx?NyVq9N5WPdzj2OI`t8lQV0s(YtDRZ{>VrVnbng{9 zdE8Hdckcsenz`dFBSoc=ij@EM39QN4^p`T*o>j|V-2!M$;3p|z(dTTKDF*P6c_nEQ%D(>6x}eAB)R1}q zZBy~46Q4eOKriULL1#{jJu_oB03fW8r{Q`LwB>mgN(JSgZYG%(_sOJu1Og@LSYctI z!ZW|%IJ|{)3OLYFM@KE_NO~ze0}1MN4Z)#kvx3S>yKIYdcgy`lLzQDQ-7>&i=5jzl z$t0Vt3SGM2Km%x<~Syx!ur z<<<8qYiQuiBwM!aL#TnsbC9!Ir0XMstd!C2sE`lq!1Vk=8 za*#v$`8^@JqX=uMb#>ZbrX!MlPjDAG5eEDF+dT2v7BgMhXNkG~0C)nHESOjh%mRc| zL2ah^0^Ba}c+fN{9693lollk&fbM}vHW2J)k{}@h#L~z2O38cz`3=mWQc~VqT+0R0 z>eqg5OutUl7d40c4G4~B4hA)Lb$0{f1OY;Tors1^O#pm@AZlxS8$4Au2E4Fr^AL)H zC!>yjNg!;m^*Q(X6mAaJY4SX5sk}NBxzw;9pcy~_#J*vlIWRTdUk9o>C|k2R0gsDeaRQf~eIAP-lKlNJ1Lqb|iH*#D>Q#0e~!ut9( zV8c^P(Bv>18!l|8jdV(CYv^4W900y~mL~jXv%BC5+4PnGc_^?6gy8nh6t98*8N?^Su7I2AF4*BrOt@#Vki|FA zUa|4=l>lpsNZ817p1d3ca($xP>o3MyIPe0{422HUVV!_(1#_t8L*VK&3OfJI&G2@z zx?bbs_DOehVPRpN+-}wG@|W;Z7>Jth85j{^Si;Wf(JgHc+S7o&xGV77lppM$#*
NX#7A^ODb;5EPfk*f4S?hadUtlOUvfTpbpQI1OAG<8Qy1abG;Cd&vmc9 z`g6R*t(AG%RyqYx)TJgC?cf+<&t*kO0g>9;z_ihb^-lZbyMdomQ~uOmvrIO&0P~oB zA`l4F(e+JDw}&4|9^#_7q|nfwh7pcJ9UqWOCPI+^X5aukjh?w>ZubNIOxW;(kgiO9 zKa87(qFIrek zC$y!)Sk>is;lB1~c@muBXir5v5-oHa>HvGI-|k`CL8L zEKp*qsZlI0!NSHOZxz6h12NnWTL66-q!$mz5$EUTj);RXE55Msn3U85lB6@*sn{9s zdCu0>7X0X@PceX45`xzUVn$WY?<9#!NR+>SzqLV5JP;GHyZxK@p}3r!KZfh&s9`h# z2ncZX0#k~b6_R>nWF*A3bJa!dW8>q&s11@+7j12y+RBRi=Z&g3IXLihb5}tAiQ>Ri z`=Lb^_H#+bVjvM_q^E}qcym+|(mQ^D@6}sdOh(2rb_Q@cxlaOy3(OXXV82@CHvpuD zk)h$`hub~w9vDq|LBqd*-TCxGr zDUjljm=ZhmKwne@H30iMIH?uem!FqTf%%@*z&vvs?8aD_nSl=}Fcwn-577#;Ix!7< za9>l3j7_kbqlHD+^f7*n4jdi8!8JA$_M{-YU#=Tlo){(1+φ(T,L,B))/2, t)|>collect -# if visualize -# fig=plot( size=figsize, -# xlabel="t", ylabel="F", label="monte-carlo sampling", -# ribbon=sqrt.(f_mc_err/M) -# ) -# plot!(t_ni, f_ni, label="numerical integration") -# plot!(t,f_th, label="theoretical fidelity") -# savefig(fig, "1fnoise.png") -# end + f1=averagefidelity(model) + f2, f2_err=sampling(model, fidelity) + f3=1/2*(1+φ(T,L,B)) + @test isapprox(f1, f3,rtol=1e-2) + @test isapprox(f2, f3, rtol=1e-2) + println("NI:", f1) + println("MC:", f2) + println("TH:", f3) -# @test all([abs(f_mc[i]-f_th[i]) < sqrt(f_mc_err[i]) for i in 10:N ]) -# end + t=range(1e-4*T,T,10) + f_mc, f_mc_err=sampling(model, fidelity, vector=true) + t_ni, chi_ni=characteristicfunction(model.R) + f_ni= @. (1+real(chi_ni))/2 + f_th=map(T->(1+φ(T,L,B))/2, t)|>collect + if visualize + fig=plot( size=figsize, + xlabel="t", ylabel="F", label="monte-carlo sampling", + ribbon=sqrt.(f_mc_err/M) + ) + plot!(t_ni, f_ni, label="numerical integration") + plot!(t,f_th, label="theoretical fidelity") + savefig(fig, "1fnoise.png") + end + @test all([abs(f_mc[i]-f_th[i]) < sqrt(f_mc_err[i]) for i in 10:N ]) +end \ No newline at end of file