From 2ff8fb487f7dbb6af3e2285f72f323fed70c6bcf Mon Sep 17 00:00:00 2001
From: "Andrew E. Torda" <torda@zbh.uni-hamburg.de>
Date: Wed, 4 May 2022 20:44:06 +0200
Subject: [PATCH] Behaviour changed. If final_temp is not specified, it is set
 to the same as ini_temp. This is consistent with the Uebung description.

---
 .gitignore                            |   1 -
 Makefile                              |   2 --
 README_version => docs/README_version |   0
 docs/build_notes                      |   6 ++++++
 docs/doc.go                           |   7 ++++---
 docs/uebung_ackley.docx               | Bin 92316 -> 92325 bytes
 examples/example1                     |   2 +-
 examples/example2                     |   2 +-
 examples/nofinal_temp                 |   8 ++++++++
 mc_work/dorun.go                      |   2 +-
 mc_work/plot.go                       |   4 ++--
 mc_work/rdprm.go                      |  11 +++++++++--
 12 files changed, 32 insertions(+), 13 deletions(-)
 rename README_version => docs/README_version (100%)
 create mode 100644 docs/build_notes
 create mode 100644 examples/nofinal_temp

diff --git a/.gitignore b/.gitignore
index ab0ce1b..1c4034c 100644
--- a/.gitignore
+++ b/.gitignore
@@ -1,4 +1,3 @@
 test_tmp*
 ackley_mc
 ackley_mc.exe
-
diff --git a/Makefile b/Makefile
index eb745ac..e3b5e3f 100644
--- a/Makefile
+++ b/Makefile
@@ -1,8 +1,6 @@
 # Use the normal go build system for everything, but we want to automate some simple
 # commands.
 
-# The annoying syntax for debugging is,
-# dlv debug --build-flags -tags=no_gfx -- examples/example1
 .POSIX:
 
 GO=go
diff --git a/README_version b/docs/README_version
similarity index 100%
rename from README_version
rename to docs/README_version
diff --git a/docs/build_notes b/docs/build_notes
new file mode 100644
index 0000000..4bd4222
--- /dev/null
+++ b/docs/build_notes
@@ -0,0 +1,6 @@
+We only have one build tag.
+  -tags no_gfx
+turns off graphics. You need this if you want to cross compile for macs or windows.
+
+The annoying syntax for debugging this version is
+    dlv debug --build-flags -tags=no_gfx -- examples/example1
diff --git a/docs/doc.go b/docs/doc.go
index dc60785..9ca73c1 100644
--- a/docs/doc.go
+++ b/docs/doc.go
@@ -13,8 +13,8 @@ If compiled with the graphical interface, it will run the program, load some def
 Input parameters
 
 These are defined in rdprm.go. The defaults look like
-	ini_temp  20
-	final_temp  1
+	ini_temp  1
+	final_temp ""
 	n_step  1000
 	n_equil  0
 	n_run  1
@@ -39,10 +39,11 @@ and they have the meaning
 	foutname  ""
 	fpltname  ""
 	xpltname  ""
-	adaptstep   if empty, do not adjust the step size. If anything is set here, use adaptive step sizes.
 	verbose  "" if empty, do not be too verbose
 
+If a final_temp is not specified, set the final temp to be the same as the initial temperature and do a run at constant temperature.
 If final_temp is lower than initial temperature, we will do simulated annealing.
 
+Adaptive step sizes have been removed. They might come back later, but they made the main loop hard to read and the aim is just to provide a model for teaching.
 */
 package ackley_doc
diff --git a/docs/uebung_ackley.docx b/docs/uebung_ackley.docx
index f90b40cf1d35f149a1637193abcd42112f50055d..b782f3021c4297a6f75fc4d4bb9cbeed630f07b5 100644
GIT binary patch
delta 12789
zcmbPpk#*@s)(xN7_`deNZ<Uo`V7NPxfkBXA`ok1P;rcf*XZs#|6z!A#V69T}QD@t#
zBlTW~LUg|b&&ph!zU<;%Yi9Z26t6|yDMBitT4xU(_>pj+!M|;N_07|ENAG=O{KNNG
zrSF8wQTaeAFV)cNRn-@Qm5e7(6pTL>uQvJT*Z1?7=G$zt_*hYKJ1~C=U*f{UvHCH0
z%If|XR^5F)L1pQ}dhc4fO}Edrsm{HB_wQr-n!R#6?nJWc_;60(n0(NxD9`$|NiW;o
zyhl2b4wGjeC{CL9p|a}nhXtpf7oShOYubNz-Le&3JGzqIZBS|acg&<_iKn{qbDOP;
zpB*dV<yUrG{8}e-&Zgx}ZB-_Mn{U7WR>qmBF2KZTS~`77or&h^?wGgr0!DsHEHQ5b
zjQkg0QJ7OIrJ6f^QmvF~?ew7cCSKOBMBZpw^ha<^EtX?9X}w)>kvHSdhOTvu1unt)
z&wj@3Vw}}Ab@IZv9|cVAAFfPuTCm$ac*AbS@|o**TZ-1*jy(EMq55CdZ<|ea6V9jG
z-&UV-e_iy8b@x;KSI%GZb915en<w>;C!g8Kab=g`XGV89HE;b}RqYQ#vI2IV>9mn!
zdSLkBbJ=T;+{~=sF}Gd#oE<mt{=RvGb>EXe2Cku*hO<19e!EDEt*U)Chh@P**W>Dr
z!A#+~FKddJZt*dmPSkC7^;<ms;tQ+RQf1loOCz6_JXxV2yPhE|_a*b|GgCNJ@?Pon
z*DqT1bVu38$XS<EeYbV>iQFvonHgO)dqGR!mWp|bHNu}mLjpd}bARjmQ0V%MIWmsg
z#w#{WR6n2m?C>PF2b#?3AJPOr%{W){fa6SaLz}Q6(+_bG*7_r=)k04#uCdhmJm=fA
zO0wMihqU8d!;TD}`wc--zdSDcR$mfaDlj=n$F#lv+l`Rw37gllPdK>qeZF;9Qs0z^
zNA8LF-k7uQPs#RYu`4}Fe&~tW%U|l{U&AM-Y4N*1?$u_$@9$o}+t&DF@;%<U?DKO^
zl!#hAn7h$VbZyFPBj;}{7U8b)FJmL}HM<IV*gRci8+5uN4?S|fwp01~llAp;)*fP=
z`TwiF;m-}fE^MoxZTrRU?<e~diKFxU4<u~iGMthxURG@+7wKFZ(R%mXKWn~k@AN8T
z-@cr7>w4DWERAPB?j5?+^k+x9@7kmV-cL5qx^$r?#>i=JLj|X&rRQVbhpVU596vl|
zqG4~XpOAjU$rWA-Th-r8Y^{_mUZ|0Bg?*}8wtlRT<t5IUrx)_&N7XZ=f6WQLdd&1>
zAImze;@EiZ#N>eR^WRVQOe~Q!c5CQ1$vn7p;a5@FW8a%%%%ZMF>wDg3<nNy0=Dkwl
zq1FUl+y58O-`2dn%RDb`S?aYJwi}YC|KH;k^oVuh`6G9D*xygM`eNF}*Z0$poizI&
z!n64TPeIE$=9?cvV=jnPU1qyGuYT3HW6MHa!(V>&34T9Y_u|<f)w(bKiBu=JOli2y
z{)zF{^wZkQSPv{ssE-!ezbNJG4Q9v7UG<kwJ-;J;>w-&z*M^tQr53W^mWkd_dvNqi
zpm5GrDaPoqFF7}^2q*lmZCZBn)CS#%6G;bJRBRGhYD@~|F+?f0w|$&A<7M#SH@<IX
z)vt;a4&FbXOQl@o-`CZ*X0PQnuY9(i#q0?4e%{_m(=2A5$}oAg+O&PS)~k#qy{q{0
z?frN&T~2nV+^LUo`cfGA=t!K(8-~39qOrn0`b!!s4DT#lo+s=dl6>p^*PbYwz%O<F
zySM*ZZXD+&CB^<D+O1G?&4r(3$*ce0y3v+-=I+$`4+6SN*cbk~wYS6LsZFBITmR!(
z>y>KN0^+{t<!R_R2CZG`JgLw+$sprRfQV#oNc7CE^lwk&Y*+#=R;xX*V}58L`_RC6
zet^3E?tMb%Syxzi|I(P-aX7HI=>Mw^7HNL3^1B$0-jP=N!aw(f=%1;x*!!xUe_ZwM
z>XKPW(=04CuGTL+-dioPO=pXl{Xg5Dl!Gb9;-lxhFH^8yS8T|#*1CA!^~lO4%myoF
zH``8c++S^c{15ZR#tRNL?-QlnCR`9r3iD#n+`=J$L7#2T?t~KpSs5NTgk9gRP&G7v
z;@j27$$41#qFU^y4_#*zqQAMYE<VYodE^1xtsR=N8{2Lbq+Cg^-@YzN`LE>VEz|P6
zTBax;`nG=O%2zV$dzW8YTbOUNC7aPL@K;xAZ{fFH#wXq~NI3O0Z3z8vZB3fJ;U=%J
zEA06fqBh6s7lwJ3?kd<J-n8;`67#>)<(o1?lpaMLsAW^tS^07T`$VZN8D}4}TP|Hx
zu=}Rf!skC199|Y8DEc5W*h8Z}<DuQew_h(6PkYfQ=d{nxWsM2f?*)<E*)H~M*ElP?
zXS6=@N!q()rN(#52Wr^|BP{C$%xc5j@2r`9?7{~3@4_aBKjruv|8~Bux9J7j-6?@z
z8V)(VvV1x9tnoDaejUH|ecL&TyS>AVX9}fGQ@?ff?Yc;p%$}#5vzYsG+Y&=oDb_z}
ze{e+0fra^>tnc)Sz%=c7O`qF8SS;v#T)BO&A#d77cg}6DKc^iwSSddJ2QOQaEPMS2
zeXk8~HRSh2xhJr^H2A1-<{I+^F1Pu|zeh&S{PK^n%(me6k7-56gTA{gO5ea>_Ele?
zNZ|gav#YL0t2s0)a-3C|CSahjB52C8Q)Y&T>rHsJn$Jv-3J*K|@$UtrH?Qub?q{2(
zqaw^G-nebnqphsEd!@G|M7~{l`Q#0k*_!^_-GcAF>b}(UEhXoEw`#*CRY~FGMLn~h
zNHli({pOuk^2Okp+tf_8bC>io`9F8rvM`7LnYh(Y>i5}0GS}r<zHDk>efVyh_Owp1
zCvzCBe)yEt?|+)0cH-<YPv4&9yCRJ}|1x#l-KTb9694MRf7Nk(vwrLSI_dsq&IPkx
zr7xPUHpTH-d-X5+sy$-;wR_jjojX@$o|}HrI$C0V-@R1X1Ezn!#QzVJ_Tq{B`Tq3I
zrbXL-emA*MxX5#tVx!LkB}Lb}bLTPgeiONRJn3}zkIyrePg$nc@11wm&30MP!HjDw
zPmA5J(mU-aUVPnKx^!pp#kmt|^#08e|7vTcVcod%2dlvSkC!6N^|<2?t?tMRa8p$5
zHUF8<_QvR3`>)IO|F+BI+~>Ymq^=*=`i{SU<C@4LHGdb)3N)yaRDW~Vj6LvAXIxUn
zro3%4FR)k6kDV9(lYi$NulD--*9)&!Uo^Zbb$rVKlj&x`&o8b0u=n=vPo;9*=g+(<
zZM~|m_`kI6kACC(mV1Y)^5dEAEAEI3eXx6ubk$bA6Aw?mE_HZ5ZL88t)$2TUx`%hp
z>6vf8_oT}i>*I>LvEIVB?_Cp`czU^9^?lh;lgDQq{(jr^B$N5?^m9pee}B){(|%GP
z{wn{~t7*DVo`)CJ?!WyfsA`4Zmd?QH*oCXVAO3e`U7p;mZ#ktwdk#)muYIhI@Bg;e
z-RomFNqt;qwOC=t3WwS|yAKH850R;{E=ro9bMK9*!qp?UKgWN2c<!mH@`TlABYb}B
z`SQ(wkB^(siAPacx5ZW&TYME3+%w^_{o#svx*Bup8>esHr(?lvFJ3Y6gKNFs3Gc&N
z$Dd!F@4x@gPthITvvPMfZrGH>`Hm^^N~D<I=54MoJ|7jk#kWT1K^)7P)z_O`zsjCS
zKlms?r#8`I)#uIUS=@VWH#^(BSMbf#I$)AM?N~{MPnGi-#nZQCk3YNernf6rzy153
zkm7TPJ)c}W^z3K-5`jmNGw#+tNfF(pBX->J#CgkU^Ue#)>6(W#nV&z}+qdIm$8pA8
zig)(x{pzIqA*8L%ZfpAv#l4$0t}vJRy5*jt&1T)I@29S4-J0rfZq|Iki{iH;wyapH
z@xOmsN?g(Ai_`bAsha%T_f2Zgf?2hdTfc5!kQg8KtNi%mJAdCA?z+WQzhU?OS5Er#
z>vwOC`*qYUesS--$ZK<3moxt=KXPi*+Nj41^OBB5o(unavLa>mR)4QqVbii(*Q&Gj
zzL1%GI$i#%s?OS!r#I(^$={68d%AS{tRVUI-1Tvvs=`?Frp}%lx^BPt!@1FK>+35E
zA6^s?)6rc!FHoq8PuY-H{7u`Qqxo0rKlNX4F!7W;llR@RtgYaak4n7Rg*k~g^rbue
zQa%VZzT%0PD{w+D)$-@{%Tvx&^rQzp@lOBr`g((%!A~iENul`*buH@Hej9gN)~Q#0
zh(4;uI=5Djb*<WFeY<Ddci%NX^te{9=A6vh+4laz_fxNJ7F|%*mv8_5Tg9*I$Ch2%
z?pE)-^qcA)Gb!KD)F$?>w+8dXL&J65k9|8MKaq9X5}QKBVq1=%yJvi^yt&^tf98(P
zCshpxKiM8uO@6j-F}GQlWUav3&5Ly(8cAKa5cz4wj<`v6ER5EVvbOB4s=gf<>Kwl5
z_V=T?N_@o?$AvG3$qD~Du`l&r66b@=Wmz+8&l*{!)E}D`_rh;S^KG?{zdy6LWuGnZ
zRQ$F)|B&J5qtm>cYUi)7uSnVGuY8fI|LJLqDcTo*U76$eK|*Rx!a8-~x?~fn*TEnE
z=9=92<T{U~m-nCR#>I|Dg!yeBZjUZoZpfc!ZTRVrY_E8nMyA7I`Deu*=Pz#4`u=n6
zlHi5El6$pQEwY(e|3YGxeb=!C>npwp%saGZWv#d40oe&I%r3mD)7#HJKjZm6HM&rY
z&v9?-MS=SC9jA|(Z9bNr@Kl-0>xk$V2D|R%r}^id7CNaEw&(G`G{%EW2d6$(3G~gJ
zqqJyS<k2hNoxV*Jk9ckH<gBUfWdZfLkX0dz?@!;lZCY{Em1lFm{;sX8Z;B6Iy6VQQ
zy!orHZcg2r^^*7Rst-LYuK762-3eTIEz~zmde>I%Rm&tNWnVSmt=!_SoSwEn|Lyv%
zIl0pfd99{g-In<7+nR50dh%~yDzCn0v`n)$b^4BVk8dhOtd+X3Aucz<_Orv*ROv}u
zpJiPQ@Hsh~O`LaurkGsq)^HQwi(l^5Gln=WnYSkYim30^wLUxcU9t|?d%M_m`9%ru
zwg3^1$AbR5WQwXDJzV}nRw(Gp$@ezuu}6EJ<XGRguWR1+-tw}9)B^Wq%nx}tP1^o|
zee#<R!BvMUk9-n&Frk0Cg|C#FU`10-yp>VSivQ2QMJbss)Kp!4!@4m<MUXkGz_epq
z&cXVtn^(5>8^~!G-(wJAch}Bu_E)H~W=Y#8?rVGBedotnN0`OVT(EQsVN!qGQ0(K{
z_NaBG1N)rUw^h!rQ?aYRQ046UGJEmsjVpWB^B&!|U+x?4?^|K-a-9vUUPZ|Mm$u`b
zm^*brT|m81<<Eqzwl8*;JlQcx)3eYoD%@boe=f1z^+9S+Og)pQS$Xf7)4b!o<KidA
zx49a98QJftDitnlUS28lKz~CZ_a31Z=Y?ycZRC!6{^EVic#xYxqC%kKlKX=&v$L~P
z-gvMbQT?MIHD$`?n*KJ<2Wt1uzp?%N+S4K6Nc|b*xnI=Jx9UW(ye@Jqn;*PVXzg<b
zYnRlzkY%%`hMcVbr?EC-nVr(g+s|@(Pk!H!&%0T<;GEsXKZ-}3^PGf!#;iM(qdq0h
z<C3ZVdp_TfDJMnx|2&#-%(f~=ciBq4_MIy%9QiU-b}Tj#U>C^c&stP>Xdz>szNo*G
zn#Y{bKOPr->Yr}ev~pd)-WFBW_~oJ}{ZC~0tn2-A!r1U@+>aW!+WJRE?2}%%2u@&-
z=dh5QcKo_Y*<`ab&s}eB`tnlf?fl&fkHtzHKK%KxkL&71cFy^;)ostuFuJ82J^RIX
zi~HKXl4nn*tm4UDe*EgPo|_x3em|c$b(Zbt6EB$DPtCQPbAHamX<HxtZp(d=+vdFe
z?Z$-)IZvlsU#q<@mwC0h;=6LZ)r|U=SKmCzx)xAjVEt<MhxHY|rbY;^Kk!v+zwZtX
zHPu$j7g6Vm-`tRsQ~UJfbo08kB{q+Yc~Z+hOt{|7_3}djpXYM#G-<&{v+e#`yt#Xx
z>7d-ZHpPXS?IqJyW~-z&`iR#QOjLTjAca#i<e(<Y6PEMRM@?q9uN1H7w8{D3k^Jz@
z!ull1XDh^bLKd%&zQ)2Yq04&y?{2yIdn0}X-D5w0+k1LNi8hbr+QYi{PQDcW;=WJ2
z*zBf+u!9)8(qiW)Ol;dZm(&>L3rXJB{uq1E`}ewCTV*o3j~7(;S?ma@v25u4V40Bp
zIqzv#&bkxl)f+z7GcYjx|NkcPY~CXefx4+;)6Ztq3nw)nY*6#ywBK#nW9H=+b#~)M
z={sF-x9z$u@XX}>hW7^iX8xN(bDeY8B_>&1TxTNFbNAi+RqOwq?z()QvHOOr@}rFt
zgjao8`e;j`-JgnmQYWM>|Ll;;_$1KEJT-v-nd&4>wdrdNRXC@bB+c@?suJ(wm2|A8
z>2=Arsvm|@9sczp>=nsnr_84W3ue29F|V>Yr!S=}byj6L%X71Y8BbH*sPNunpT2RG
zo_2r?Bjei_k1hW_4xN?l*SYYurYd{C4Hxfxcb*eM8<nC2=4<+Yl2p!~=(XsY*X9+G
zhuGF}zMDQ_>RTRx)7k$QJuEx6xGwvN*_jn-{8H|3=Uv_$CCff@dVSaL?c%(ri|)SL
za{VS>@69>O7S{cGw#*=f@5~R6cVDbh4AMOhcORSI^UEgg#);Qi0?Px<I@?bd7$3d0
zZ}&Iz<hL@j-!@Hu)RpnXJ6?0%!*s`Ek2qGmEzD1{4qdVK+tl0T%fq)i%nfci_VD5i
zl^~r1DIEQryvhqr7Ji)7ViLEcq&_jEwdK#exk9TvG#;`3i=RDJ=!o-<JTZCAGZ$UX
zY?oe<EV3rKeD~vvH?)f1`bGT>Oi#JKGQnLtaP=A&(a;;o`SsrGBohq$ll_~oP87Y)
zXmk9*B%ZERf0pUI^5&bgmj2l?$+6;4dg$3j_i9gEQ<q<5Ew8ZI(~dpHHfZHb+u-VY
zPv=8=xu>_?Zrfm6-StDWShCScXs;0K@jttoPjsFa5?h?C;8A^G+pcZG%`?w1R;zTL
z>vCQ%e(<pI26g*=f0@PHR!k|W-Ip>??32xtQxd+J^PT1dmA}mSJfBtA|8&9j;@iUe
zzO|?uv@AMYedouX3aRS64K_*=b|J|NEnQN3uGQB}i?uy1GkNn5w&!!Fo%^l2|N6YE
zzf|}C-}>j)d)@u=_xFYO>ny+Zd(Y>49G~{6Zf&3OW2Yk9W2Gm0#%hKeHV6cTxQK4r
zQ~U7y+(+`3)@GM)Jh0o-XSn`O$`ZAHskyVxTnQ;Zu_|d^hvMHU-6E2$b4rfSx?gXv
z%I&D2=*c(H;8lI(yDNq&FCIp*KjQpR&2jJWr<DsYET1K`H1*ZX5Thdb?x{+R5feE+
zzjbHK6mn?uZIn#8Ag<EZV)a;fvPNk@8*9_4B;A0YMRRW-|F+?$Pv^U1FS)~(-d=KL
z@xtkCA<k)C9~h@<I*0Bv`S~ut^-C!q)6$fK({=n;^DPaMtj&M!S}%Y1j2zFZbo=F9
zI*)7iZJzdO+nkqM?jCho5whm{v1><GZH-vJdfKgNz6#mCVUyH2Cuf+fH@05&T|{M}
z0Q=khOK$nEn#!y$w2`HJcmC`0PwS5b9I#r_abjk(-(ko8g)CDR+%CJ&`70sDN69j4
zp~u(76aQX0Sdx4C?Ji5Zc~OfP-Rc|6W-PQkdZA>QMAj0v{Zsc}V~n}5JaKcvJZqE7
zowF{h{#nSUvhaeG_U=8}g&jSsH`?ar-YJ=UXJZfh(nu$f|6Jv1)dkDF_q5+KJ}}o;
zgF$riqekXPofzTW$rD`^gpNzhxcaXDR#|?C*@kaZOK+CWsXKX~Xi1qw#Lu(Q(l>Iu
z>MvfHA+BZ5CGq;}TB`|@jhEjnu`@ZzyDi`1TSLvHCKDfT^@nO{vdn(zMrH1P=}~OQ
z&({BUEhw9Pv87orG^f)1c-i|^*}KcupDEmKz3F!N+oJcY4i!EVn8)gT@C)C;4>t^J
zRe3^>Mlxzd`d^3$U4A3s-Vz;et?S2xl;8bM3_M-$WyiVPT_ktoquhH25zdz?cIZ^N
z7To)8%f$O?{cUBDmDcb4qb8TRPP}t(%X_O_lNI8#&fR`_hViRGznV<a`}GS{cP`iL
zn7`$iq}S?(YR14u4`+Gv9p`IIS6!N;Z;-jn_{Jv*^AzjJh4*)@`H`{x_RL+E1HV6a
zn|XX%6yv(H=Dqby&+g_Q&QwgD#&gFc;m=gh#!VIWhj(qbcYB@7QOoPMb!#6U-)H{n
z*_973H*U-G=iPGUR`1%zdELKu*Y5xFoHujjR_U1be}8^_t6k50M|S?4d!g)?o@A}r
zzvq?zUC*D!`KlsQSDiSUplCOB`RpGd8dG>wmMu*C)m)%gH~sSR6js)H4SwfIGrTtZ
zC=|$Q^=D92R@QRfsCj0Vb@PU#dHW`wU-FMn?BR;5nf=F;ou9>UoPDA@+u2I6Izsw&
zDbuBo#&*&!d=?K>*Z+-6Gdua{o%We1$0yr=PviW+c&EHtZLU;8VC}^h`OoSWc`O#=
z^vh4b^3ipHVSd2}<u3_rsVh3<=DRtY>F?o{_$pcOC~DcA(<@qUsh{Q+TWQu`x#zn}
z$E)e?2G>t%erx^S{$vJ6i%pt&OtnM7?rgES`4ddG8$VosmRHB{?AeF=r(dvT6y01c
zU&mB0AoBb7ekKNnv)l{}LJSNH<@rS^`f2(3C3*QJsl|E~xj9#+-OihBAkpw%{jKbT
z_rc7TVx}vvESp&szV<frj$<sVe<<8;dGa>B_>5PWVAa+bmLG5bf39!a{ag9;9_GMV
zT>ID#_VU{t+Ss9afBnJt@w-oQRoqgLTkKP?iA{ac!TP_?esABc<keyjBeBT#^aCp`
zmAK|B5|WN93cVyQnfD~AZt*z7E5_5MvT()$rxg-ybq}VO2XXpr*~V_CROcXd#iB;W
zty@%U%_in^3Vsufp4>2OTCG@I!5v;3V8N7P%hPO^Rb`sRI5*VxLXOCZsH@6nk5AOI
z*e>DD^?S)vi%a!e{G)o6zn<#+${piz=xvgci1RVKRkkS}d!H929TbxM<z}4aH|?dq
z6wjhRwJok;t5r8-)|e)8&g}3s>y_?`XpZ@~;@bQzzs`MGacRB&TDFsIvf2H+Tc2bY
z<x0hxs+3Lr74~3`VBV<}j;==1?gvY~QmtigM9hlykv{QoGlNV-{k1pRDFz#^{xiL`
zEmX?tfxze2cQ&rs^rtxK0so0*&Re9MeSZ8tb4(~rdxv#JNXEG)X&xCnUWGS|Vl#MZ
zB29Wi&o=CMop*d8^X*j+Lm~@hZn2t4STNg56|-H~*m-%iMxt=|i;{_snwN4HS}W~5
zx1zG!dgEq^lYg&EJ!}`+`6i?OOZJxMc`9P>UO(XqyBqM_cisN9Gud5B=1$>pKVPjU
zQgwQb>o!}P=!bDL8NcUqAB{Z~U8M6=@A|yih6%;rBA<#*wDj!0;G)#p6L6p_rqQpd
zXz`X)n@_B@jeeN(Ys2kx^70OvTlIAxFFF0|dB|JyxIe|bHTBzPyqCLk%xjA(=czd3
z%kh<(uBR2g%wNSEKV^4D=Z4bjc}|^*Wov%>DK6Gye7L6e$4u-0fq5(qv#oEe3hz0#
zJj*jw`qckp|GPIHXp%nl`BUwOUxs@BS~s*fR{Z|4<@P`M#GNnhr7K!hUzk1pqZgxC
z{T9QIA)qvVoEtQe22JCsc_~Oq{AAv31A(^ptiSl0-tnnc_6k^M%-V56*W(s@g2pW+
zB|QaA{rvk2!_<qIBQ3A8drKAB@ztN}DJiagFmwCjTLzlb>MBHk$cehB`Q%LhH~;?Y
zi-*4V9923vQ>5gjeErIg`Qi8X-+hy;C8VfpDyPuW`Mpc$@0BP~<rEJuOSAAfg1WtT
zZfGl6PMXcDEb?}#Uw}x9W8v~GDpw>T-m8mQ=ltPv=lh${d91pi>;R{d=Cx^ZQEnB6
zDP8L>F5>QDYe;Ahyy4ob?{;=V_bTRDJWP%uvX@FEUY>cDcGx3ux=Fpfn^d`Hcj6L1
ztB-k2-ma`Xc}t$%IJEelQs39;b#4q{H#UV&GpuGwGgxgnqjNi>U}L!RgU3&j4#^bG
zypS{Jipb6{#k!A9)nECs`@-W(`*)ehnY`Yvqf_4d;$YgbYx{Qfmo5FhSH*JqA+c48
zzQz{<xh5N3En&Bt!q%M0d$U&UKz(2O*FAf=s;AxBA9cI7l1F!g|BtITKd$NdQ=BC6
z#lp8e!{%4Zhq|*?LVS{U^lUt42^SuocXm0Oi|m1(3-k25XKh>|(Xr*~TZ=X6kyqz@
z3Xy5wdZ^8|`oZx8KL)jpo>hXXt4dwZiAT>#x~#HG`O@j8S!O$w?Js)d?hoVGzvWD5
zy?1i_V%u%Hh41bB^b$f@OQW7V%RabgNrdC_Cn|M!xK@75W)pO*teEgHu8`@SywK6u
zuIMU_H}BTiUQhlL`aAUL+vMsH2ZNBr-nCckX2jI>9ZH(Vbl>yMF}Z^h|7M&!R?eTi
zf5Rc`(^jwd{IU%Hamah#^g4kx>6gFWjZNOQTzXD?`ZnuNAv2yOir7ryw#zIx<PyH2
z>vKrjd*N#KbBf7+|1&Px*gVhsZL$55`)cVcJEr$k`9GbYr!MmU>bLWnj!QQ9JUM#$
z!PCWB2PFg5y1$kP*{%7{A9ybBy>R2L>dn8nr(d*X6ychDf7k1FW(Ec;p~=V8HP~g3
z3ucKo@=yM%u56#RPqL`4@5j5p(^AjP(wfU<SMay^e5r=IlGo3o2X!@b=6#KImp-FB
zhtrS!U0>GszaP%Zzc)+FnZkT8rd`^TZF0DVTD?=RL45MtcZDgpd&D&NFs5_AQJR#H
z`r_Mx5B_a23$(ASZ*CG(n34G4df?p2fwn1|S2gZh>o6g-x8!2qmorbQ-g-FBExFlu
zW#!~wZ#_&cH|{#nx<rO8d7*^ahj3lN!>=+`&ok(3^c2sWHcQrEZ(Qd=&HAi~f)$}G
z*E;>wVwBYCm1dPLnxw?`&g)~r`}(k-x33&NRA0Z^CgYlpPRnkml5C?~ulTDeZ@Ip2
zRM2+Yc`}?=bJjsCu62_`!tOS1U3xGf$mv>qbc5;XQ{T-l<?px}mzDc?5f6L9<S*a1
zCPz;EdA)Pr^#!xU0xL~CAD*6(s>UNy!*M4ovrp`L{*4+hhk$zKFC4i)Hfl)zwPl`o
z_}#WSTOMxHvUCmAQMc7ubAJQ-Kl_==4ShaS=D2X4Ds{akE<Gn9mh-OA%S_MMkR7tt
zQOY}O?z-&=Oj70F@}+dkZD;R0H;PYkuKC+$b??cq(vvnS5l+ifRQALcE&QO=+>=mk
zv7kVGeZ$_SfH%7JOK#uEtIvLWL+&x(LkHD^*LYRKeuy1;;p=|+-HW6*&%%C)&AfAK
z_pUj8wHC{Aoh-Jveoy$cY<k@TU#-hNw#hP=CDt1S%g>AqHh)p}d#+OV%v1UNY119=
zJ(O~OKA$UbMbSiAPs=mPXEJy;Z1&skp&G09ZT<W5uIYA@FPX|~^{jGuU4PoJet$aq
zg*ukdgI6~PXErnhXNvVcxcxT7zu%hkYQzQAl68M)ujiRqlzT;4<j5=aD?-8?3Zhp-
z)~>WJS+>32&o1V^uS~k$!|&I2&c1WsXIcH-`?2;`(}EfId0Oe8i3nb;;yi!umfnZG
zRpAF;Sg?s3tGdj&{_1H(n#FwSGhzx~b~ZBBPJW{$>t4EViu?T)`}D7VeZQyTeyNas
z(6Y_1rL`|^KWAF9X3C!<Ube6N)U~euXH99Z)4f^yqauCk_Ycuk_WQs5`<cX7F!iUY
zS=JAEhGu{5?f%+~(}cL5GyBy_u`w`M2v2`+%P3hNnSa|tpe{b)fAhlUrE8k4a;|9V
zoV>Lx<Cc=&#`9ZTZo8-iD!jCKm-PQPpXSN%8A_Rf3LYzdK6`Ygz~apA@BX2d6)G_S
zVg9<$Pi<pcaxT*G_17!E%I5PbXDg_z<Br)h%V({v;I-ec?f+M7pUt79x}lSC#fqa#
zPA$LYb}G$Cs<XbvSxaYKgmV@{cu#o5A@%rvtuyl5&KjhLHmk;7nsn;#iLGxByt=Ys
z<<_TPmjoNG$v8JNL3^c(TF%x*mM{N(5&C0V7Cmvw)7xEFf}&*#gPQo)h)Lh-h>QNZ
z`4+>SFJ5x>TRM-gdh@fa?{C-JcM|%1LhB3tIZfs~b>Z*`2>r%)ss2}ngl%xihFaEb
zyLGcd72e&bDt)!{zZvU0V^&AMdH0qnzbff*_!PpgdvBS##fyc<4_F;MZ&|#^L1fd7
zJ)YgSkC!e<*?2@+hc9%A$)xR;Pvx(?U7g(H@^jXz%!U3XCk^w@r-U1Y+U<<I-t#z>
z^PR&>9);h`Z|DBvbquNhGRY>hepgp>DX;A#>*T8m=9NbtG^+R;Exy$BzU_kyAD_vb
z!i-)K<!M6Ae;>^Y;hp1U^gwuOXkD;{Zu5lSlWLBdp4@U;Qdjli-YoT!iI=)h7cKm8
z<jSV2<@-b@g{G@4-80ShWT490s5@8ghioyu+?H^a@%)Zgxm#8Kw5+mR%5_R|>4G`+
zPd0rG`g&ad>)mbdrrNOYPwKlaV?KL&q}}&v@lSSy$H|z*3Z3cCP*;h3Tx)TE+OxVP
zO>udQt8ZKr(D41+vs&z8rN#-ad~Z%M@dcaw9bS|b>}&k{L}6b>0pB*EFz<7x=B-bd
zKhtNy&*B%pYmY2@(9N6YaZ5|=`uiI$8cyk+JaYBw0hYntBC9G~oV&T49~b1t-<Fvf
z@<vVjSjt(IpQrryAK(uSJG^d%-if{YiX(%>&)eyJnG;-Xn%?xbce(P4B!l{_2SKqp
zFTCyva)vud+n(K$p&u4nEbwb(X<6jGO-(CgaxH6XvoCaBV?X--Qi{ou1c3)z6n4!i
zd1vf#zo|~We)Fy$9O<1+ZndfZn%M1&8m^rBB>uwp`|q|C7Uz_N2S03&mB?+HWVU@u
z@v;0@Oj3vQjorev4;^G^O)lkj^|HME(qgqU=jZGBn|6ynkxyt_eQ<Ys^G^X|PT7x#
z_uX_@Q?EO@cBZ_x^s*C9@^)KyG&`3hY`v^9=eOhQl}`N~Z(Hg$Ugb4%ANhXC#8k<6
z%GupP`y3`Eyirp-!rOG>W6y~*y^}X`eq<8qeeR~>nkO?`;7Lo_rKeAr3I)4*<Q6FS
zFeR4+ym4JA=^5CwNqoW;vttZrbI<RQH*x=W$5hdd#ni~{c>kfw8(;r&r-&U(POkfN
zW8TtZIlnU{mQ3X2l(FXNsW(cLNSS!iZ$?#J#u+CQ=Y9ESoTmJmnLcHv<+kK?b8lMA
zn7OwvamVN0$c{jX^kW^z5@(#9AlSKZP0!7~VD;i&BhI4d8x3!BetK}c!=-G-%&o4D
zViKF@s+7Hu+C1Uby4RBaCr|%Xm|NYMKV|0c#GanAXNpdZhQ*&}&OYj9E$MT=e%;31
zDUI`Io-Y%#T9@0CA=c`bHEm|IQ*vv{%;!g~xEG5x-HY_SVja}W^VTc&K=Z1p9tB&U
zi16oK^)6mwG3}4yveJHi&YsUtL?Uyh+giwMYgm!saqfZE3&y#ZxZ1of9{985`cvI!
zuNU9zYnEPVu(a0ZMUb=Y!}{WFxltwTlXL5<H-6LMbDFEN_NjW**OHB~dQ~f><P_8|
zY}4CkQzW}QxgudvYquA_#{3QYKA7mQY_BiRJmI}BV7lbaqrR8#b)Md8ywz*hy%S|>
zZR?pnv;W-~vs?CA)|p(xzm3IPjUT)dbDwqlvfqzA^Pb6Ou2{bHy~y4#)wAlCu-fWg
zcrP%csiNNVHT#Zvzn3mr9}&4pR?|-AWBL*^j;cxjZeHB=SnuP!+^3T34~Y~PuXEg+
z;25*~l~2U<AFcYYbSDOD96!1*N^U+!_}a&+q6TM7KFpfX^ZiD=RI$LB3qN`1s6=E|
zSzOay`*vI3`gp~l^VSo-TW{*xsNmZEby1!6ul45MNiT%&8OU)>x+&n8EY&A$E_Ewe
zMDcq6@s~FL)>|+6yMNER|LowF)^-^eMgvB!C_U4MGg%oJ*i@&7xiV@n-kx6X$|zRf
zP;h(tOXV*ce{oGs;gJx|k<2I(PL8y)vifEB(&Oi!r&EvYRb6E1H+f;&J&7oO)kLX>
zPbydZ{b=&@N%r132UVI1XDwK`&28tmr{9jt+pT9i+2kYmRJ2+~`6T;l9ht1+W1AY+
zJ%1=8IOD#6%q0{5^fzigo%7D6S-AYtS~`)B=W_jo{YJCeG|pR|WA>S5aNVrnVYc;_
zjzgmPIt;nCbDU3Ic+Gb>gst$4!}r$?3z!Au6%AeoOf#CYnmf!?VYjeAdZ=^JE=5_%
zIk$N^uT0#N=Mg4*@_kzfw|Wr&2NP$`{}$ewjXnD>rO$QvwdsUi;R7Bco1=zl*?pNi
zRC$h^^k#Rf-*Z*0F=C=o)AGqCP8&t1YP<}HzL2sZ^MB^;xBz*Px90u1+#=gwPb(LH
z@nY&xOYPg2%x;OEb<f<Ppph74I(N~G*E7_9%B0`E+F0WHSpK={iBD4)yoIKJ4V!t>
z=<2<wTW&{9-Z0$0aQ}>KC-<}Eg1(E7=5Tf7cAJT;s-5P;(f^=6_W0Z-H-)d~?EEo>
zVFBBNwiQp)0+W77^Vr@=-G1Pp?&VcSG}BxkJea#7?X>BK31{RS<ThrkYhm(=Exjma
zWYPZ9>;kJ&kn+=<hu2d47r*(r$K^-ovhJMrw)4)_cXT_oANq-GFyq`h|H;qFlzrzD
zkA-^)@yX4fJmpa48Oe`dUVN<Qcz;q-puTcJNBf`7xYV69UeAk2-V&X+`hiKhNZXs+
zO<%O~`)$m2+CI);?%JWgc+rM}gmUAy7q+e+TtfskJoSD_>L|u<IC!@riT$0&<)h-8
zIwb#fP1C50N?)-?**iHf>PnJI_v!q$e3rIEzb^fcK{W^d1b$Q9Ic4hEMPlad^}?HU
zdjdK*?>Ky(=3kSovH$DY*||UU6C<S-OKNv*zQ3a4flH6O2*2rRQ7fKj29?w7Lic>$
z^=sKCy^VZ3XH9>2i~Ur|UWSktx5P4*aqGEl&(}0IW89!UMfPysvYgmSMfbzf-!kq$
zbZhJ1v}>-{V*5((2deQkg!T#u8ulKUds?{u{q9$HY)`K(;QvvrT>f3_&hL4@f|%aC
zic0!)b=}PU<+}fB<7+wQ#rlTL+}rE?VZB95-2xk~eHwB;ExGb@IreVS<eD+rdfCm)
z<8vZ}tKal+Z#!Oo=4X{x;g{E)wZ$nv_u0=p5wb2z-~H648S~ZM%ag?!wnwzCI^{23
z*kjfu%TfRTef<9wTyJDwRBy;!$5!%r#uCoctrIFccG~~mr*_#y%vbQt&Dg$&!shOG
zZV8-;J$Jr5#=KLz?#*?+;%(vY<=kCO=NJeE8=q*s{YBrpF8lI5uFxF<i{GjhJM9Sx
zwUX?L-BRr#E~0$vcUxGsl<btY%Dszbc@=$LvEs=Ip}H4O>(!sDX64uUwcWYDY8T_@
zErHd++wCiJ3t4BrUm2zwEPP@WX9&x+R@R_%o)3OMS4-7>FS2xE^l`_w*chp&J&&zw
zUi4^xdi^*xCRe<6eb3~1ADGi*>`&||KUOJk^jK#7<EX#(sdH>@)*k42oPOr}DGOJ1
zotUOg!rx?)lX_}bOVrngtJ|2}{rETI%LdVM>$TN>Hw<bM*v~ni&fN21?xan&cGpkN
zay<T{_p<1sc)`b)%W6KnG7s@(+rIRz$@}w?m)>hOp8O~%yZ7LtAC>DK*q$|e{%cp@
zoD!K!CqIPjznxXorF-`7*CxNe*H)i=ytO1QoO{>e^2FuM>*_i$s=as>_fo6gae1BP
zwUlX5zh=K?o%<^_a`M^COuKbrCGPs?{Mcj9HtrEDDcwE6WA=nMu~Rzh{dpI){&_xC
zqeS*T%Q8Whq6$&xz)9Nqt0pVYxKt#QG>ebnK#|P$9V~WST3f4`L)Kh04LHwn=fI~q
zkAnVaq-k~?(V5rEyz4JNXrhjpA;6n?`Yl^V?&(qKi~`%my&1!VStT?wRNA)}#WAWg
z*5@BG;Awk5`J4QzJKmQz`L57%{VlLuQL#7ka$aPT!eZaWM;3p0KYiO4p0={hpZC>0
zy5Xt5TZw;L!d$JdGhCS6m@X`u;<f(RlTt%oca?$@Tc1>}O!WEo;DpNVeO|4eE}4mW
z-ppd6I}S|V;@u)S*>sVoMrjf+YdQ1#3vn^iwmo=JAG~`}TQ_e!W3S}y#jFa`Wh~b=
zE_~Vfq+!Y;_gVR~%^5T`o&+?wxi5Bf>5pTH-qFq4`{{I~#vSP$8E=o=+grTKJS^07
zMO*LH@A4PF-~8<T`AWdcsxT&xX<ic-9Fo~IMb24%=lk>e0ypHoclNik_ANd8ZtIo+
z*6vGQ_hwivl>1rF_HC;89DVzfdwiRO!z90S=R8zrzc;tI-FW}arA$9J*zbP6K5lx%
zja9$WQw+0pYVUn~7M98K@~Z{k&D;}h-5E3P^a&latIA0C;y;}k&tmk!+4q)c$*J2P
zK4u!Kn(A}i)JQ2j+r4;Nd00J*So^mB>fdDVsl0hN<p$%eo%8D6ykUgYpRD(}%C0}#
zZkWiZ%4l{l@3sNYUj7^Q61&WcCN0%way>dB=|!S5^O4+sX6?(J2GK$)rrzJaE`D<?
zc-!1(mme$t{o)cdy(fuLvi|M{B^LwJE!Q4&%{AB%{qx`48x^lFxc4qp;1mwakzX8L
zBzHy5nnT!2`jpAFaskaR)_JqtwD)X&wbfbp=Ng^$mU7Coq+V;Ab8T`gtnB!&9Th$G
zPwVnmPdxm)xU8=`xp<m?{PB2p)4QYDGp?A)vC3~=v`O%mSGk7wqT>!Z3GI&c=UC($
zDz@$8D>771Ye?G5&aV}u{&dfr^Q`ANKAzkerrmKR<!7jNesgZp#|zf_lNugWb}qP`
z6S%iZ*8il3T$lbW6>aW=>@N3g-q`G@oghE?prKGStM0m$f}0x67ChL`S=c|JZ4cvI
z9=S)K8nmL~S^iEwW|MrBo5kmC^hV}@)%89p&Nl_V>?l>NihOxv(l?3k2KuKr{!>_9
z%h-7OklhPgfqf;{=L^V5?TY@>els&~+joPAvu~c9JD2`v;jV3SnVor`DKBEH=kH&)
zsXWv2UR8x>T5<ig!u_QSGSogFeWZ53K}U7lktAPb*>)qLNAuP_xLv~VH@hrie)d1V
zC3CB!1n;e#6X);s>zArk@fOa<&6D;|H%Vcf$c5OpF#UWAqd%j;boo?9eQETG1_n4-
zH$5|zQB@j!lol?}Fn9XOR7MqP#0VDyg8)L7fq|iW`ioRXTP6qX$%WeD)9un2wHPN&
z&rD;~kVf?HAw$0iZHK0>Nn^BQ<evUIjZt43(VBxygCNvdPd80xv}a=Ro?e#DXw8^0
zeN#H4BjfbxU(y*>rNN~OBLjmfGrI9v(@ioMRiwdbosogzIWxMOHKrG5Fe<aTq%blt
GFaQ7{3Q(E=

delta 12864
zcmZ2_k#){R)(xN7xZXPpMny?5Fg%$s`45|9{pq~d$0YBGzp&=nlG$`aOmmWus*Lpl
zmz9MRm#}~3Y|?Xoaw5gIs6#O|a(7AK0*7$d?RQHqFMpwzUo>yYf1dZ^F9S9`dAx;>
z&!e#IkN*Aoxr%ek<L!USRjKF{Wft6>ynczCQlo3G`HubX{(eX-T=_jf<ne)Ha=a%>
z^Su|{yKY{;|KHW!XOAqqc1zu{r@=+Q*KA?w-Qq_UbH4nZ_*6jNJ$zrbq146Dn$*Dk
z{^8o!mHzUkmenhqob*G~C~KkyQ=PGl&7&35CvBd$D{!;qy#t5mEV#8@=1kJ3kP8#|
z)fRm{Uc1*bqku*6LV<$0rzZPksc$<&c%5Dawm7U}Q+;dUrCWcBw|}uoMv#=OyJu$5
zrT)cN3e@+%d|;x&=(A4ab0N?4H~Zcvf7`&P-&oL*T>i2uZx8FNO9sl#cON|D=zOqp
zPSg4u?7<RWYqNy?HyUew*|k0Gpi1pOqxyN>cAoS6{cj4d>Hc?8Ic=6q%d=WHn*&X3
z6Cbx#w){RG)a<FQ{Aq`OeW&~V)9-G#6df<PId$F{mQzg$;va+m)ap;Hv3HGpl_2P+
zRKKX3BXr(O6}yVWgFiGfydIpA74~pTxS(LL#p8I({EHWYrD|7q1o$0)D<LI0_lVY<
zzGw#R@-xY$iC!#9*4ySbui!C?`TZwSCYjxR=i=h7Ck6|*%(-(-KT+hKy?(tw-RT!r
zIxBXbUtAWwWGipJ0pEhDY$28>7C+}Ke*QpU)7*UDy!XtfK3rn%RIm}Kun-b^;9ek7
zCuDzevc&UkZZ^li<Sj6c3aUPRUh&Mz$pXeT+<q&yE;lafJu^?Db(hfRz_Su(O(NU(
zWg9o@o%wgT)PK^Jl_@h8*PKw%-t+Czv-&f?FHK1JF|YG-?4GX<S57K_TURZA=Z5m$
ze`j~s)iqljF>m|&V4`l+8fNZI#`F7+#-6lIZFE!f(Bqc-=o_#<n5S6UXOTcOL*xk$
zmmA72?=JBE^4Wg+>j!dc>;CaOm%rY<E~o6m?;n@#CkPr{cNQ>c<}zFtFZ^wy<!nCQ
zqbm;C)qDN7=J{4Pd5_ZTQvq`PSz)27xn}a8E;jw~k=JW$UFiMf^2|w@b`~s&*WL44
zlDd+%xt7Q-exD&AY;yR_VvVA|r9q*BTb_5=Oj3|mHsfJsseW)N<MA;*FCE#1p-HbE
ztYR*3u3d7nInsACbK~k)G2gE{MY6qnRV5cCBx2noz)~M{B#>S6>cNPh1N+jaY%*Q7
zuXOPb+Y<%n3tp)739kA3#P53bmrcKy)V#ZP?igF9$ufSEN1w`HYq+hGuX0(le#gge
zuDwPfUq9OL`fd8_o0#@castO%_AN`Z!!C%_U1p7qo%VI6*_uV#FLy1|OkS(JYh}%P
z*1PgExgA>sm((h}|6KoLQ>N+`h8+@*^xZ#l_n5WVD9A;<IqheezK)NlmpP=#*DxhH
z@kKXRj#`Rg>MgeDQnv%s_-v)kesbKX>!`?>WR&w=MS;h_QSd8glML@9_l7Wshi+vG
zX{Vg7mK28+Z(X=j^K+}iriUNa-@lsoSt9TBnRt$z3EY*^d0keImUz}*ep#Y(U1*-p
z$;&SjqI@fM7iMKDdqt{D{nyNV^|;!^qjv(1Fvk4rT&@-or^5KiXvd?OyS!~DUOgId
z{c4o$#3>u?X7C?B&K;K}6~+69_h60n>JLBjwuRcI-d0Gv64%Xh(Xl!p%VtJ*SI9fz
zwC;~dYOh7#2Dffn_<c?41{TlKW%YcXT5qQ`Ts*N)%V_Q?<xItSH?LM-WN>(T)41u&
z%kGxU&X$?aTC9Gyd~K~b_q^}%gX1cl5mR;tPK(z)C9yFnKUzd1DPHkmy-jD#0mF|r
z9#?vguih7aIqTqVm*VwdzvaT8Cv4lYWt;s!+lWm9MpM^opEKt@!I!M(b7AJpuQq9^
z_1?b)0>Y~;r#J538T;{H^J3<ujFs=4<J>G3a`#29WIR&f{^7vqMgA*!8+n$pv3J-c
zzLVXm_KfT9F&o3eLmymr^vX)bEKn^uy>5nr#jPSPo=XnXO_q8%Cnxw^t$G{w<Y7qV
zF{8D~vxHQiHf;L4(j+wYap}pBQw7V;_XMsz;Bl!wLW(cF=IbqkB69|w6@dZUxNhFM
z`^LTL+C{1A3A|R-+upievs$Gm{;k+?PwCG4PX9iwDUlTI-1RKD{{Tm4mfaJlrEL1Q
z0=`YqS!ocH_gUr39+@rnQ?ynw*+}#7vg%z@7q5Q2W7=Wk3(gMazdNR`*}y%kx5I7T
zaRsY|`GTth_SH|iVPPTmSdR71q)3i1FL@WA`6hzfe=o~gvH0bkCy5bpZ+W+DG?WsJ
zE$Plz3oYqC?DRVG)zc3tCtpt7aO1<B-wNG3GWPmJEi(CZr1tBM>~KTQW1pDLHXYga
zB-NzQ`31jXXK2F*^)J>k)_>e1b+$D;;eT4#q<p$`eYNo${h8f@_4}v4elFn{BCmJy
zwf_-2KKqZ#*A>#I%v^trqe;O`N&IluUG}vHeoZ-hd85(GyKLKHJ#PQ;O;!#n4_TDn
z!C?5+UBgo5e$kmz({~CyI2h0}Ex=oVMZm{1=yuBVSqoRLS~FWM|M96j%h}S-r6IAq
z-+Vg|)UrhDgYSYn!FSdzIBGt*o=@rYu000nQvEY0onQ`*_t<2%x+CYDTk%HLfIgM0
zA-*ejtY30m#HqwEuX9QB*%_>%uV&UvEjw2c5&G!h{2PxZoQX30d~2#@kG#m2O--y1
z>&|s8S@m;fN&{cL+TL@j5nM&5BSV58uavu{vapWr$hvspl#6$(X3q83cC>!q<eyf*
zG1%#Ln(Kp^D=L!WisnvFn7pcjtxh~VJY3xLeC6hQCr?#et$OauF#X5l`~QNBRrMop
z{(Ktbcrf>??Kf#pldkP2{H1mTWCZ#zS9fT(X)TR@vZ>>6`Z|r}SH2cSZniqy+5B+I
zO5b<;3v(wm>Rns+obQ!U@j<_geM0S*7wtRDE?lp+LCl`{QLoP2O&dKw2t~SXcxc7N
zYj`{DxkyBN>au@7uJ2br_ATarbXMx&*9ZRnRrXGMCwX_zigb~tdFq_IufKJjGXKzY
z1OHQM?_1tYdi~{j=G^-DQ!5SrES(<q-8I)LV{2w(=Ca&Z=aza;tp2?E)B81I^CPaD
z?pl$r`oDDApY8RB%2~EAIQMq#hqukILq+#}Ub--Rc7u-8-g_)1S+Ux)<`tiR5G;G{
zXUU&GZ+%xK9h)<I%NJ|Tjej;>>zK3u>y3|DjP;?Pd!KwLD&*)raNPd?YRw~>l{Jd>
zk=rA7iO>7%we7_3x9>}NZ_3Qb$~_)G@%ga|!&P<ln=iO;irnMpUG}FweHGLHPY0jP
z?k;iL^Tt<`LsKh2{+rd?6KDVJ)tBg7V|?`VMzJNgQmX&`Hv0E*y|cg;Ywtq>u~Xv~
zbbjABhef<_az!-XRnAY3CUVTv%%68)p8Dz4>Jy$W`TbUX!gHSLkNN+K-6EsRxa`*d
zule)x{(6PjM_*QpYG`K*edO2V*fytLDf8OCO+Rb4vX(Pz>+VQreEn!^LSSw0hsX^#
zqP^oK`d$e=H$Hgqhxj(e|C|3?>{`We`CjBIZfDcd=I;}p>a(Taa=N_yXuI-nd+C+m
z=KP+0qIv4sU!Dt(Y&O{a)8V4k22SDs;Vn<}m1M3j7UZ4#w)Mc-tv|On{kbiBs%Bo>
zkMHN}E5%+lEKgZ+Dk<ed{;t#YT!vSBSK3|p{%~g9Q?uz}hac?|I`Q(Yy33c=eP?&B
z->UK7dAF4Q&%Nhvzy2j8^*?_7>-Y<2@7?43D`fK~ck`ai{@~BI1xvR2hE#vQGV|X3
zx~n(zYr8LAmsqTRZR^I5TaN$dyAdgO^6<=s8G9b<d|%yVSwFk)$*uU*EpK(Mv;F60
zVEF(4P4wTq+ZICgqKE#+h|fOqYg3B!Y=(b#crWBNZ@F>wuc!AGT^AwNu5FLy_+D7w
zv3~RVj`f?#cX%{!PWp3kN!ca6JnkzKUApE-o;#l3c4dO^OQlKIc@Il$w=KSrKOr<c
zC9xp9s^qh<O|X6Wuk$n1>y4WGdeSycZ~gLk&e1o!1lu{q7ypS|T=_Y;_;c5uopUFd
z#!u++tkEfboR;D#_jtwZrB~VugAN+0O<PmwtX2D1uipKsjQ+Pt&c^mU>!q2c{9Qaa
zzqyt@+pDtLUQbZ`Fw2))Ms-<_ME`MpoqfL4>h<r)-4@*W?k6jx-cQbceCk(yed@ja
z+y4Fjd3C0{^9)1JX(mw{dRjOas&Ll2vG}IF>Du%3Ph)`RlQ|PP+SDd>{#bBgS@#kB
ze_8u#G#AZz6+SVg*gAEy&+#weQk!c+A3MF688dkY_nKy(D+-mTv-Pa2eO`C(TNwJ6
z<<;}jr8{4>v8T?O9K`Y2G&S{_>QNKdrP=js8};)o!=HprGGNZ0Ss<up{OGgFWRr>F
zDVrv!EaOk`%`owpz*FgzHpj!y&hX&LtaFo+^U~)x`R$w9WEa-*?^Sl)UhBVcA6xEw
z*`@bgTm8=EKi?*)oSmx{i%Rf*ylDGKIX&gZ+`z|cjO|!C4(4~w4c=rlH_wZ8FYh{O
z-Nlc;)mOLO4e(d3p7P_$tMA8`8NLxcu+FFAW5y2G#&Z{5T0iFE(kYzzF8;Z`@$3)Q
z2Y+sPbAQ|7x!ET|#hahndNZj-6u76a(U+(V*vsq~{=o2o(c4Q$<`+9Q8BCUMJlMAD
zRsOnrappdkH}?MdR#!c{wNPrK^}Ln)np{KLyuWdCrxy#<&%0`LH|%)PT(+4TDrc}b
zzdO8}SE%rRBkQZV9ZqkiZ$H`KIx%Z;;F0k6(|ue_LT$q|b>=M-n1254?Z?+DGM@74
za9mxr%WuE^hU3mLQ@8A1*ln}Q@Y;{=oaJ|?q$^e(iZ^@G$J)PMnklo+{iu%57l{iW
z|CVrlkjY(L6o2Pcq1NL1e__(1p@nK0R`u6r>|c3q!FrvEN7%)5Hb0iqD3G_<=&-0S
z@WbI{!Ou=ld?&4$7nq-Ddgc9-U5T^fwBPSx3bIpL>=JvB=fUP3bLu&mr)3IleAwsy
zT9SoZWL=!PZm4IgOJ~*ERX&Ri66Jq}T{F}9`EOd;>5b*He3oCnxOQG>z^?jSUA5OG
zR`(9NUH-7<uBGqxw=dH!Z~cBI`;hgjn^~6E@;2v2?Yg<~ZfuzD>5Cz={Z1S52VZ+t
za_8{cb?I+~(hj!lx>+2x+jw@a{_@+!Z|}}DZ;5m}ena-oHrr_P*?c8y%&qfwmoB;a
zRmS($zqHGSQM1gr^}cbOE%3GdU3X3A@TU6P?n}6IwwC^mi{0pR`|Qlu5^t5y?lulD
zk?P%;!f~7ZRK=CqjO;%Cha4=|wk-1B$6dQqphqO-_kO*4x7+-GH#YFSn3x$}+5G2Z
z^oM52W7}0;UwVBk=G2zkIcat~DrUugoT7Tc?5WZ+p}eDOW?H|s$ck!vXqL=%u1lrH
z;JC@k>Gew=a88&vr#boae5>rt6^c`*=(zk4bzT~kWFQq&dYLz}p!8G5Nug`ZmB*w7
zj$XW~_3MGO#@@M%YJoHQ4r@5=`#7sZuUq|vjG|WK<c|j}y;=hfxoEv?D6XA-@{iZ)
z=kZ~?(<NkA*IQ;tHhUlc;5L14-{aV-V<(Dvmgn-FIk^90eQrliH_z8*x5FknVy8>W
zmI$@)dmuUgl=l4NzIJn+uEm^SU6LfVbCO2Is?un|kUya+W~Tlg*=nql|H<etKk&{r
zV@h56<pn95McMDAa_zYA?s3SGAL1ugu*T`OBrjaEv({F0=AW4#XFOo#kjMxSTI#(Z
ztXy~3lP4KWS5*Im+j`ZXirnXZnB`&GL-n%ri>x?0wtti_T%upk&+KqfD`L9*lkZ10
zt~$P!NJtOp?LW2nng6Q8|CX(daF%C%Wi}&M|Hr%1JNGuMO?Y`z`0RI1WjkLFslwRS
zn<`Thot|_t8}vx3Xilh0s_Wb;^WrSyF<s|vL77<x6vgLWR%QOd?($KSvEJy=j8`G1
zE#ixMoZHt_2v1P>D4yhb(P#athW*EWRpi{R{hb=2|7r`@AC(@iqL(oXcn*}WxKb_l
zk+aw1Kn=rnN#z5K`vjQz6}m-59lwb@>UJ;Zxa0TmMZvo^sqcr3Z05;7(~3&{%Ch(F
z{4aAR7**Atm%62V=DU|>a<+wxw5@dg<cZCii=RF9kc-T;oZ!Ap{)tu61^1n=<0W_x
z+%YwJtJWP+&b2r-q~u}%^X<ImM>Csh&R*h<-#te?ZAt3dAEk^3gBD!OmpvLh?eMof
zjmZmcc%6Q@YjNN8?Pi`o#f^R#Hh6ifw#Cn0H$PtCx2Z^+itT^Jq}OZfV+%hn(0U|V
zSX;lFx9Q1Lb9OJT;wxL{1XTSp;&tZJh?e{M$o8_q*~@pIx6~S|)z0%_EMZI86sBl*
zGkNpzy}u(}3L19BIjwlYQ0}?yoJzN8OveFNts{;nWfZPV3*tB-*?XBI(ZD}dqFpld
z+Y~Ol&WM=QCk@|UO0K*7>UvY-%aWYyvR0z+jpZcklgp)Vl>Pk|f8h5YVTbNH?I)S_
zEZ25Nzdd==`HOqt0>jzwdKe~jMkrcMi8=GYjQPa!!pTlG_dDNr#osQ!y7s9@yue(Z
zza6RKErIpx8a`Ma`8H+mRPIgt6z0c8TcoBs-u0IFG4H)g@SE`PsaZ^-8_tM7Y+1<J
zd52~GlTA^F9go*bRhUk{T6U|gwRPv!YpIp<td9l7|6B27(<)uARnPjp;&kq^hi)$j
zDc5oijV=29Z{aI<zK8!f7n!xco%r&KYwnHO^^beakGJ^=mHG)vwR`Pjo!@v;BV6XF
z&0(2pU!_@R#7dThw7NWG-hF__<k3l<|MymzD;n}D9NnjNKj_Tii}jh}%vYLz*R6|I
zJF@SfK=4AL>-SE);)!DZJcsLBifYi$i(#Evs|pU>s}o<?60)~tv$N&#VD$%+rik!2
z+0K12HC%K4+9Uxx{)pJkE1JIL*HyTM>`;<_p15sk&??!#&P%R3+-Nkc+;n=5+63;<
z(>JZ2T~<@dXPUN&MV8s|*42ly^=q@MQw;=n9@wUG(evlUm);+K@Z6kiZ~x!@p@QOz
zT{~0uJy|5#I@4&@DT}2RQ)?r9y7pVmu=%~^O!%C(o9Dj0VKzK(SG|yb=Ycu1E*hWT
zTD-%}ck<>ROrK+~@Rj}9Tetk|*5ZOs@4x=`UOWHqt=He5KevAs9JcUVagnyD<Gh)D
zbL*{-sGR3Dyc)&b&9G=u$ZOU8n+3kF`QffNzsrPc=lwfLQxlp!MEj&(Jl?Cd+d0hm
z*=Lpq%V+3L2t2T)`rglv$GL4*aWtsdx*b}OTaoNOQ6c|Y(Yf-L`i2M1>Juj(Vldxg
ztSMI2{6q2K>Yg({yN+~h-_EV_`R*zQLx~qs3!e!))`ygp@Ora^PA~Nref!g@yyLBY
z&iOa%UK|XXc%vz3zS}>gMMn17tN02RYrPKoX7BZ(l#AQzq)yBxnK0SrpjrE-e_j-A
zUhKg3b?5xau3=B_+?nmYbDFW-Z9O%osVRGHP27$}PX1c_G-&n~&TCQg6}@LVEVexR
z=0#p+tiO0+-Iag!+N<_XUMn**FxQ|y`~R%^w+5m%YSa0wc-d{;jw!O<ozi?S#a=^v
zG2gV5mesqm^HSnV4Ax~<tebxN`1xrW4Tc`lX9CVH%Gh;AB5MiD{;B4tE&nX=Kb&;X
zq&#Dlkm*wIAI(+`>WnIBWtE|iBqaYXwA{vJUh8USaMoEZ((A(p=6Zq6Yekh4BO>R$
z3}Dbw`;p9%vN6N=%rwWKzDX5bjgn>Q3A5iiUOn(;-^JpOA8yHczBSyCb=q6{#$Tqz
z3yt(c>S7YkTsME^a4*IxJp3Gg_!GX<<`p?iQ<521FZ8ndDY#-e$8uAtD(2;(x9@)V
zRDXYQ!ME8LHCSR--TYM7ced<)y{WD6_DrqaCz&t5yM6UmXm9dcfqANJlU~XSAKu6N
zXT6YfY+=}|hEvB@MLlIIW{jK`@?z~i3!bJGYf?q$x7__0ETq#Ga^%mZXrbemn}nM_
zFZ|14vf}%9{{WTp8<LUAHY=hG=YHGq*s|<v#G&#<#_m^~@2)C{us%Mude)^+!r3e8
zmu{KW{%=l7@xBfwnFo%)vqF|;S~y?uQ8#!Jd}qSTo;KOV;h!(gyDe4yxbb)D;xoyo
zBUx+A?wiao{+}r9=RG0Ppy0B?WA)mO(y+t3%;IjZt2k>pe|yxgsk3YKL(W<rzV`U$
zD%<o0PuE<%VSFL@KhOL8zd!eVeZJ=EgW?zG=g<FNZ}EcR`*HQ<ck>y2)*b&9x8L}i
z+tu8AQ5vgKCLTQ;RI{oq{Nu_g8k!SVG|gJY|L9|-jZ%?g*o20*K$oCersW~AjWG{p
zJ0?tI6cyAL+Pr$Nqeb$*J;LXg{A>QT?}$OqG+UWh!c&9KsS1D3;VMx2XlOp|wMdxy
zg&MPt7mvclADPu(nE$xs^38Vt3)!vn<4bu&A3T5Mt;@Q6*Nl*O^DphL7rx87N~!v9
z+`V??WqI*geT&aM_EX?J8vN~O?U7|muURf;WcIUhURpgX>O`%H-VEI(yUsq|dtNzY
zzfDu+Pj;6dOV6KIa#3iA=f2x|Y5rsNyr1n(J5{beZDG$YWn91F(KDvA|EDKLGKy}t
zldof{|JmpLXEPH6!$fWd1|bFphVuNP6#caP{F1!<lGI|oirkzl(@y5yHV|-l&-#m>
z={=unWv_s>!=oK1bTe-8Cup2fQqoh<JpJ$Yq^P?mcvGH+$Zp}W>^D9<Q}F$}ex1*J
zE}2?Cq*%=DK5!gYT2`~4ZGZh+HSGtP%@sXH2Tw}PsqbL>KKuRscU>BShYAuqzUef-
zUNND9Hz+xUW#KW+M34IhT~jY88Jb5JXiREJV&YzyB=M^;exIkR!Nps8KRSLeW(GX|
zP$;1rwqnC6K8t1x*HaUck8`YRJM&OC=BwjFjtO54cz&&Td36PAN!VA1&4CKpq5UTE
z?z<k`Ns`fgGkIG5V~@*rnWlYLyw$Jy7btRWZ|)4>k^H^#+XTh0^G<iNdZ)aQKD)v)
za@qb2gN{$%1;rxQPCKyd<JA_`6eYErCV3i3yhTSBN7Y|?Wqo<E*Zw*asjQCwQ%{9G
zm?N0yb$+4J%-+@qR!cXnEijvw#(tTt=*O8(=F_t(TY1=~o!)=z^0s=eRCbx9I@$VZ
zq0{R9c>XY2KN47bOhEWw+~$ofH+9~K$%)R|{D6&F+Jl*+x`9=VS-S3U#8sY--CU(I
zRvS#D{I+WMIqyEe#{Bz+d-6F3E3?8XQPrGBTyyn|G+cN2mbj~2mo&{UR4)@#-uXXD
z@lSx#blKjz-fwrdZap;DU|mA}Ro2qDC%>|f?OEn=EAM#0mp#v&7H&Tim2&#<X32L`
ztslHI3%s+qM)i2myKlFOZP{|p^|{3f&YZl3!&fEf1PAMy6RIy%G}OIT+jQT#v03!o
zGXC9f|9oVe67>7Aq2I1G`pdJ<XMQ|y`r+^0q<z%|!rB+7sc4+{seiImM7Q}u{VK8G
zn!tApRS#W@wdPc9yRq$^b!*QqUXJv)$IsvS*<Ym0aOcj!u)RqpcA?IzW(C}n{;z+G
zi91i@{OSJ(UmxG`L;R2s+rzhyFX#QPZ+W_?dVcGRH_rM0z^S|5;2v|*MkWS^Zf*t!
zL0H00%}YT_*SCvi8wj+$x4vASP?FF6VWZC*pSwrq-btPGjZs8V+V?`DO5&UC>yKRI
zUB^Co^{+z>Oa5;8on8L$-M<q-?;WNrYg+GM@U79^Ld@{#Z_yurv#q0;{H2d1PvP<l
zd*D)UvLn9!Z_#x3mg6xJi)_0eSe;UhYdj($>A0fMN#atn+eV?NH)kZc*^cJ98#e9?
zFy#4rK>ytmp&1#u-+woP()Gg+8|8JwG`1h(vuL$&JT)QtG)G?BnFr=O-ZULzRFF=T
zC<&ZAJFxi0o>vXanjNmK=F2qeDN)?<(;(=Es%~5TM77vUMG`9~`>mP%z)2`iz1+z#
z<;~(-&TUt!)!fB&CN9i3S^6wx=D$716i(e|uSi)d#Weefxv<!T&KXN4?^xC$dP1aM
z>vv%C`<Q;O->-tASJ-`=uEx)l-cu0l!YAsvR;=lD%c*6DZU)Sn>yvfjVUz&-v}>=6
zB){lPuiv`7J*fVa*p1_Vz3Qh;n-{OnR@X0F*l}!8KX?86$m&BmV%wzoCq47Ivml*$
zC1WDH!SN;Ohu2G<?qV#AUfAbcJ<aW`rvK5x4DFl?oBL}^p7AV57FIV8JQNlFB4?r_
z=Oy1o)=E3iu8eH9-e}zW<X>^?hu}`1Z6?1gZw24JR`k8Gka@LM{ff^@*XI{MGj+UF
zX5F@Uxve<Y-lf-+qvhtEEm&{C_&s0vXlz$>k<PKZk=5&tH<bQPee!gorKk4=7p2yo
zfCGFx|H!Z&waUz0K5=&Srsjj+Z`N$8oYRnWWnX!ZbZm%yde;6+kEV-$`1`hI+ufX5
z3)we0M`rCZK7aIr635@;MTgS^p6W#@<t*OXec_zI?W?uv6M~l4Tja{c8HwAS>R~9T
zWtaXuDe!oR-^wD#kMXVbo}4;67TRdXv8|U6{ixqz6q}cR`9b0D`i>&QJ=#oLzgPX8
zu9(Ou%Eg|s>un`71H(VT$&1xB*!c@vZn-|VG5M0ZvVD==xq`aBFYo?NOFcJBYc7{v
z$=~Ahr5fr=UO$Uo)YZ(H#~b4=eMWf>r(gTKzO3(mU!0YHZ<d%dh524gyR>K9WD^aw
z`e{<h`;OhZ`}o8*N!_3d=5t~<TvZO9y70E~;e4K#3%yNNG&hMUoG|$Cdd0buE9^D}
zUUj$|rqHp<?bTLg+sQHeWi?k!*1I*)s`SzV*_CSQVd6U`OgWsOBXXGS-zwFl2-~$!
zBpX!CUMigp+RVJ{y1t>sQ-570ruAJ1tS+8x;XdM0@3Pr1B<RS6SWUTif6r^xOIJnw
zc%Q$tJ!?|^!y?@iU4eO@&P-e+`KI-Om1k<{p~z?{9yQgJzgtzV>1aRCU2^}!jy+fA
zt?OQE!!WyL|NZ)Dq7&=qzNqHaN)nvvK3QO2%;uh!MlOEt?YCwvKKoVLY=V%$6^2?T
zX|q3*3~v|<Sk9>rIPEDd*E=aV<Y>{Q!^eK{{=V_Q*^=qS373<eQx@C_PpRE_kn`Hq
z3%V2MiG{2`^!CO@xsy*{^Bw9H(ym;vbp7JHa<^_BtBam<;FkP~$Ik2G#m_i5sg_u~
z9lEV5&ZTeE!O_>w%vJu0>xMSR(R-iVcb$E^?ONYk!*?8ZT`D(rEj&`sb-zhs@5_o?
zvftO;y~FihIC-t@`TcX)rYCP``|wco?qQews)y|NmKdGAVtg!H>Tzh=D`n&CS8R9R
zEeLio4F37{;2Wb3>zKzmFZ@~!LUx|$TO#v`!&r*>&6>+!ivo7}yota6USzV}<V&XV
zT0N^AUe}*C+@H^Wp^hcAdv$tnW<!00{%IYl!@sYuoF`W<v})r5?-wz@XReP?dU|Vx
zyGF-V{}rlg8yZ4Fg10VvcVgzff0jQBO3e%B>}<Io{XDOz-gxG}()#kRFCzR{f2cj&
zGhyShFwZu-*(GtG{~g_N-#)Ind!~@mx$Li>s?tp6OP>)___DLHxZ(Nfguau(ZsqeO
zW`?^rz4+Os7gqCg<%W69nc0l@P2*23&JVdi{nfAXxIOP@9qFI)^3AJbUSHN0v&w3#
z+HVfxzx&c<>973_Z!WA?-OVrme$6MjI`8+-r{A-Wf5H$xGd_8-+keJpEA8!8+KkhL
z>eKi$esQodFuW2%suhY;OG+~H(!s5rEfKl-w=H=7)gAcTeqg?S@BOzd-t9N~@@22(
zW?0=xKN=>qa9fIp*X1__|Nm;-JQsa+;a!6T8L4(poM%nZ*!BLN{ne`v1dAOd=0=^H
zn#*gT8{zq)Dyw$)`Sy*b0zwz-g+hBCt=xM2P*wfk=ko8%gd3N#m>M#0X{nVOZZ!)s
zn|#qjL^?*f@TimzV~k;(GG~8{?JB|g8)k=_uJTSw@zhfKu6%u)^XkQk!P%!@E?F#b
zZOH?l4O2tZgttXaWGb((;{GSKf34@SuH4?0SLEDf!UdQYCpk-9+<evkG~0t^>h+a>
zE_&Ns3thjx_)YkpcYPl}C48LqVT$A7$i~(NrmLay3Jcwj|1)0H{9gOR-Mx#JF%=)b
zw)*mk=Iw$z&T|}MvftaQ{LaEHq4JdTsqG(I*yW0R)(HQLF;Hh&H09|tU!A=9*C%dr
z70R7p@kZ%r$9=Jof(5s)AD8KTWx9J=tKCFZ_xe4ylkZMj*?sctETg_@qAkqP=?qup
zU%a@I@3H9D;^{v;Plw8_$$lYtUvf43?vu(*94d8_EGLTDn9mR|C|K~YM>%_m$iyie
zf2M!mk@q0v5##(7t3J)XlquSBC;W%@B8%&a&zDS^@%@V5lV4rpp;w(=cCJia{X9uq
zW%YTNt~saTp44khTqb7zW%j9UVYd>NSUs32cK4e18UH0w!K;?c@$@=WsCZi=wEmCN
z`)ym)EBWeNW!73|ndz^td1t$SMdA8I7JEus1??BjotV=5{o$XjTfQj^KKUd%;goKx
zOXkb;knJ7ooIAraqCG$CU^umpVdCwBN7?QzXZY;e%=KZak(TA^s`|Z?YSWBbpPgT5
zzEy1YBz?0L2P3_8*4^LO8PIfYvV@Ia!gC+pm{n($rGl<;NVR9P-+F$+W#jBM5|gKR
z{+T-8&QT^Sns>Kvsp7kz=W-2o=iA17nB!eN%lPoOW0xgEEN9r87Od1WOHqBlp+$Sa
z*?D16Tf|r2DxUCbWocRDy-h-^>!ornD{Heav|eLZeQ_$K<Vb?RgC$3*tY6-G9<Yz~
zXJ7i3KOE_u42$;})hjm1vo@^!^oO<T^51<2E!`$aN?6qO8^^ZqQ97&paLz~RD5vBf
zx#vz~wOU#z6s*ZhoGf&%%<BBBfJ1-2&i<3f`j36e!>kYQ9wykS_&B-$*^qBG;dOld
z>XiLw`>myyod}Y*+q$FKxg=ogW|LX$iP9^b`a9mb)ZMwuB<}GbH8RcVjzN#?`HnMp
zRAgrh=zn?q@lHp|zS$~jN6KsM&3=76tijv+{;_)8oS75P`2YFCs<Js)*=W0M@s)pf
z_@Axg<Wwl#Xvpp;b9V8%;CO}28j|H_A2NBqiL6gs<3GVQQEc{v3mX?-$Pji-l#~AJ
zGE>~;Zk(O_qnYw7$Fz<_PN;U9Ju}(*nXHj?yK(Uq&YtAPja;4bsyhuAdx|_YzTctZ
zmm<A-Zi?SKsobFZbJdS!E#Jp(@20YC@treb?Ufm8%4Xe~n7zg+ZTGr#8*csi4|3iu
zZM>K!=hB&9|HZh!#@-|>&#w6Mkr(Ril{ZblHQVqMedau;e0*!K<k`n3PM$F<Rxo)t
z!RGrOm+-_ZP8Z(3>iqR7czVt5pzvo4@^phj^er8mqrc`_1W9@Sn;>xOu{RI1g`KO|
z?Nolgc3p)miR|Kxd$tR_qYg`GNX@ADvTK)XP3@DnCOZPTitDc)XPz!~*zo`C&AFEs
zS<i62y?R4<Va@){Dj{`BrTXHxpDx<I=)=hZrrgY*Pba>%k-f}ptP!JIwd6<37ty7q
zhO+<ePF=X-hf}Qm4>42U>TccZ=dP?eQKfuaQ*55s5&18Ni{6&aSZ%od=oe<KE3+HR
zwp*OK6IZtA`B$U&8lS5pf7hojobIU4-M4AO^wwXe8J~&83I1Xb&)OgINp!W=w!=@`
zmH+u%QFA=M&)%A+{Fdmg&1-Eght3GfJ12R?Tz$stpt+K-b?XGXXZ_=5x?Q+&{q*At
z|E=QPqr`NqpY2)w$%_8JOAU{=uiU<)aWadxZvU}eOSj&=#`&eZ`HC>tmA`!Jh4fVl
zZYJ){yfEwE>$8_wezlnQ94;~Yve1zC_)2O1w;Sbu=S}F@T=^-!+~aqBQ0#yH=@X+F
zMYc=2Fd8s&nd;B-=w@YLIHfW@*p*SEKKhgW<`i|#3(Gd&eOYq(<<<?Z7mjSQ-PBlc
zd-_Y|FB^YxO-<pE5YCazC=yPNw6eQ@@w?~oQ_uC+o4nFBJnbRGZ(5L?VJpM#>GA0K
z!f(efKbgGz>n53_id~Z=__vy#+oHeS{@;g6(+=^HF3I;d9!nIOMJzoSF#CC>)rGG*
zGm^XNH%WUjO-nhcv!7FFmHK8$&JFv#B33vvOMD7HHHjfz@_41gN1@bPY<d;B?P5k-
z9yPNrSjt{}r0CDC!wjbsw!M2DC(BXB@JP@ppxtQdLYdeNs~98aHdw5dnsk@<rAL9E
zxk`k4(Qd_=r7G*CgLLO+*dClLa{SLjokWqRf2P^K7OCGj-St-|^KuRKV4;aMXAOCm
zw>YZQy*$v-e@iLAIK-`>bXJRn(P1G~$t-R$ezV8+xyQL*H7EUf7SxxJ`$)g7XHtS=
z;hfVaO52pzxeBEhPs-DtQ|9`Z|GDakPay){Qd9qiX=X&~<?H${b~+`SApJ}FZ)S*r
zeJ@9kn6Z+=Hs$(T0W0=xPI0Ofdlz@dbCzwP<7w6<EDfd%>IZbi{|Ftvo$%t|Iww}X
z%t$F!KSOP%_BSo8y4N>48&xwz&S)=I6^VGYs(Yh`d-0T9;aG(v)1Yr+dA{FGyyl<R
zZkZpk$)h~uiOrVdcQ<Vl;&W!q4i5P3I&ZG~v)?jHYql;v_~4^uy`RJ^w>g6P{c`6&
z%$vjbak_ffjuZ2Lgl_BJ`L5H{ZTHh_9KM_9EL8ZvZNc7Cca_ik=54gN)sXmZ^^6Ri
zm@9X6j=!0g_(f?K#}o_U_mVn+CRZNZtw{QL&g1e?;awd?pG5Twzija<yy?EpZPC^j
zQ!krn&6k|_P?~Lqi2tIc6%Id^*4Ipb5>&?Ks-O2kt>m`EqlSPjjn9<le|Y^ZQvOcP
z?Rw5r*Er69T4edRx0Acm)n-xUikW>E5)IrxZvHas(Y)_h&F7YOmYxueY0H+^x?EZt
zu`qk)hM9VMB<}88c}Bt^BUIyd+f|>NdsLq8TYBzRYh3Hr<jS*ax>xO%dA)C~kBrmi
z`dNxf$x@xSr>oy9zp``pvfzjBKL2)mcP@0t_t>vXS)#6PW%+P*-OT;fy8mkZZCS3b
z5?j6I-Z6na+HDuUPhI?=`_Vy{h-()+FKEwR)L^+YxyjG^*^xOB$<;9%Z1*LfyB_Rz
zDBl0gY~{<J_ndT@Ruoz-nK{qx&+LK{^D>2!uFw^$>i?7ox>hW{(fHy2zpwqP4!mpJ
z$3FX6)&aBWo<h%O3K(&7d{4iBy=0b4l7?Gk-E_e(N6rdvH=9y0{qtwp>%Nb=<bUsN
z*17+>-2c!b>CWPb!Nx~gZ|{lU_hYSZwRXrur=D%I&v1Se3VoTPRd(^6Vx&P+<a_D3
zx0xjuq*A{2WoezPKf8F*B;DvA?)v|W;vRi__B`d(-?gixS(f+Sn!W4q$?ERpBRT$C
zvu4k7Nb_Bwv{6rILXof2dF$;O*`>=hdhb~AX}x)LNh4h5`a|dNn80}X^^dZ4zo-tK
z^QA?8!=oSVXR5>J{8rRHS7yEL>i^)QC-rypM-={%{=8Sq?;<Dn^aXF~TVij|kVrnS
z^XUDnmg5H`@Atd$M;C6tyz6-j?*qrv+vTJD9=D#KGUe^Zoql3JWIvmGh0UlvYB_(~
z<AbZ`2bfCUnYK;%&V1oZ?>SvRe(dO*mtgqGCj7&^1!v1^?`rs2wUrpzUOm5U)|@2o
z;_|&qg6+R%eX7u%d;RL9cP}Duyj-%nUOs7y=dNAhySy%~l<za0xpdmUXH}ld>mshD
z6i3bKU#)ugV)N$9p5es``6{t@^^ROpJ5tS_`sjXebAjUjN>wS@+spzI`}7a5mgt%k
zQyC;He{jmQ#{$AP1e`8!oBK%MhsKI4?|1_?x?OjOJzQRMd26@x$Hg<2Y6xxKA;$OO
zn?1wyq)0}Q?Lyv+X~L{Izw$SAY_E-DRA;RBKWxC$_TJ-9yViF-Pxg#eyZmZ8ggrcX
z8h5lFxy_U-rF!qwKikaRWjQ`~-t1XF-=b{O&GJYW>4LcMG_M~W6FGX;OffsHSvm9m
zabclK5x34oMh2?UA5x~xpPzbZ*`}foPqeuCOn+!h=<#1VcarELU(K&YtV=&Eeqi7j
z6S7W3qCRN*lCTr{8(C%<l(~94EM(!j#gJ?*aFTE3lZGwX;aeC~HAMoN9Q}WGT;3~w
zaO%Gl$JKv|rb%rV^NXr3EU2%xm3X~N?&{-fQoo*mdHeC<&f2UszQ)>|D(9#81ST_9
zx_wXV`19+ruu9o+XQ8B*4hc^`WiOVwa&b?Q*Op~*8E*gTS#Mi$*iDg7Ik@bA_j27f
z*Rl`I%@6cijvwD}vz4uMv+b_u>+emmxbfmy`jeVVQEhvk+gj!dgv{HZ(s|Q$j!N;l
z%Y~DSzS%ic2XBk^`>qraD3Nt*?Sk;Q<J-$WS6}*2aYoIMdr#r44}Ik!``M@M-u7R;
zOg7H!$iBnL4(8Vn_s?a73@)=)N?hOee7kibqbg(lY0&ub9`P6U9D9yWd26M~Gt2Se
zf)fi}c_KXfj<!YJbWxq8AbI`W^>WRlT6e#N%uSP*;!B_Y=-Z90M;U7`Tw8OaBxn;;
z@20J0yG5(l*u~eTxOOCDWvo3YYMq#{{nKCTjUTTq?3Z$O6j57q!_IB{6Pp#bWddp;
zW+9nt-YW#Xs4u@`)*Dulel@#Y?NelAT(OOtiRrcQJe3sY$DaiLPv5dl>m|4U)suqr
zMMTTiwRTL(d-U<RIqTi->jo=xZMf{hTvL^|OnnzHZHgT8#zt8-OW_)ZM>l_2oj6*f
z^@HJP6~A1Fiu%(%bI!A#7kV_gcx{-#$&{a=+WF17Ngpq0>kHL4KB(+oa62b(Z<VaS
zr-xjZ{w)<v?t|<uw`|_n?5Le!?)gBXb8VB@^^}&J1-u$R-nZ;<_i(<~;K#&YSm_|T
zR=??=XY#y_k64-1euuqbNSWH1lOf$Cw_VRgUt2c%k=$dsXStTw<ezNGS17o7gMBaC
zq3N%F2_JHHF7N&?V14#z;x{R_l{4qW-2B;NA8~zKal(lNJ0+tZ)fFXE_j;fH`0Qb-
z(d^%w`+l!wb~T=<FJ}IeEk$+XQk`oimkKqKo$g!~t(ncYp#Ci1>*DkGTek%JAKpFn
z@@r4e`28-=WpWSwSnx!Dx@`*MM5JvF({H3O`m^N9GcZi{lb4z9p2{d7jXu%900*z9
z7o{?)s-X|q!sQvxsX*r?Qu9)#*QGK_N+X8B7#IW)3K<v})=mGI%4o|NI$b}F(TRDv
zCFAt<mW*Q4yV4j1q|rP72#wdLZ%t#gW9D;VoUY-*C<}2hqQM858bK(xn{JoRXwUT1
zYkGY;qcvmf^xf%<j*LCi|E4plN`ng{Mg|5|W(Ee(z%xQy{B+w4Mipsr)?j2{c+QOO
Sc)95{8H~zo)+vk(3=9A^&yuJB

diff --git a/examples/example1 b/examples/example1
index 8ebf089..752d27d 100644
--- a/examples/example1
+++ b/examples/example1
@@ -4,5 +4,5 @@ final_temp 0.1
 n_step 100000
 x_ini -1.5,0.5,1
 x_delta 0.5
-seed 1637
+seed 1699
 foutname ex1.csv # This will store all the information from a run.
\ No newline at end of file
diff --git a/examples/example2 b/examples/example2
index f24c2db..c444546 100644
--- a/examples/example2
+++ b/examples/example2
@@ -4,5 +4,5 @@ final_temp 1.1
 n_step 5000000
 x_ini 4,4
 x_delta 0.5
-seed 1637
+seed 1699
 foutname ex2.csv # This will store all the information from a run.
\ No newline at end of file
diff --git a/examples/nofinal_temp b/examples/nofinal_temp
new file mode 100644
index 0000000..63414be
--- /dev/null
+++ b/examples/nofinal_temp
@@ -0,0 +1,8 @@
+# This should do reasonable sampling of the minimum near zero
+ini_temp 0.1
+n_step 100000
+x_ini -1.5,0.5,1
+x_delta 0.5
+seed 1699
+foutname ex1.csv # This will store all the information from a run.
+verbose woooo
\ No newline at end of file
diff --git a/mc_work/dorun.go b/mc_work/dorun.go
index 044b4d1..3ae3121 100644
--- a/mc_work/dorun.go
+++ b/mc_work/dorun.go
@@ -185,7 +185,7 @@ func saveStep(cprm *cprm, n int, tmprtr float64, x []float32, fTrial float64) {
 		cprm.plotf = append(cprm.plotf, fprev)            // copy the last coordinates
 		cprm.plotf = append(cprm.plotf, fTrial)           // and append them
 		cprm.plotTmprtr = append(cprm.plotTmprtr, tmprtr) // then the new ones
-		cprm.plotTmprtr = append(cprm.plotTmprtr, tmprtr) // then the new ones
+		cprm.plotTmprtr = append(cprm.plotTmprtr, tmprtr) // and again
 	}
 	if cprm.xplotWrt != nil { // The trajectory
 		cprm.plotXtrj = append(cprm.plotXtrj, xprev...)
diff --git a/mc_work/plot.go b/mc_work/plot.go
index d9cf7de..c775d44 100644
--- a/mc_work/plot.go
+++ b/mc_work/plot.go
@@ -5,7 +5,7 @@
 // What I have learnt about the go-chart package...
 //  - The default tick positions are terrible.
 //  - If I give it a slice of ticks, it works fine for the main axis, but
-//    something terrible happens to the scaling of the axes. I could get
+//    something terrible happens to the scaling of the axes. I did get
 //    it to work by: Make a child type from continuousRange. For this type,
 //    define the GetTicks function. go-chart then calls this function, at
 //    apparently the right place, and the tick marks come out nicely.
@@ -61,7 +61,7 @@ type range2 struct {
 
 // GetTicks is called by go-chart. It has to be a method that acts on a range.
 // This means we define a range2 type and define this method on range2. We
-// do not actually use the arguments that go-charts offers (style, formatter, ..)
+// do not use the arguments that go-charts offers (style, formatter, ..)
 func (rng *range2) GetTicks(r chart.Renderer, cs chart.Style, vf chart.ValueFormatter) []chart.Tick {
 	tmp := []float64{rng.Min, rng.Max}
 	if a, err := axticks.Tickpos(tmp); err != nil {
diff --git a/mc_work/rdprm.go b/mc_work/rdprm.go
index 1be0d54..a74e1a5 100644
--- a/mc_work/rdprm.go
+++ b/mc_work/rdprm.go
@@ -20,7 +20,7 @@ var cmdDflt = []struct {
 	v   string
 }{
 	{"ini_temp", "20"},
-	{"final_temp", "1"},
+	{"final_temp", ""}, // no default. Use ini_temp if not set
 	{"n_step", "1000"},
 	{"n_equil", "0"},
 	{"n_run", "1"},
@@ -84,7 +84,8 @@ func digest(prmMap map[string]string, mcPrm *McPrm) error {
 	mcPrm.xPltName = prmMap["xpltname"]
 	Seed = geti(prmMap["seed"])
 	if err != nil { // Only returns the first error encountered
-		return err
+		e := errors.New("Parsing input parameters gave: " + err.Error())
+		return e
 	}
 	if prmMap["verbose"] != "" {
 		mcPrm.verbose = true
@@ -138,6 +139,12 @@ func RdPrm(fp io.Reader, mcPrm *McPrm) error {
 			prmMap[key] = ss[1]
 		}
 	}
+	if prmMap["final_temp"] == "" {
+		prmMap["final_temp"] = prmMap["ini_temp"]
+		if prmMap["verbose"] != "" {
+			fmt.Println ("setting final temp to ini_temp")
+		}
+	}
 	if err := scn.Err(); err != nil {
 		return errors.New(spew())
 	}
-- 
GitLab