From d4192142b85c0b8543b6b435d3c690ab840f4d56 Mon Sep 17 00:00:00 2001
From: "Thomas A. Christensen II" <25492070+MillironX@users.noreply.github.com>
Date: Thu, 29 Aug 2019 18:59:13 -0700
Subject: [PATCH 1/7] Added all results files to gitignore

---
 .gitignore | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/.gitignore b/.gitignore
index b7f3332..6d5306c 100644
--- a/.gitignore
+++ b/.gitignore
@@ -198,3 +198,7 @@ venv.bak/
 .mypy_cache/
 .dmypy.json
 dmypy.json
+
+# Ignore results files
+Results.txt
+results.txt
\ No newline at end of file

From e101ff3585a14fc29f090f1cd26566c3867292b3 Mon Sep 17 00:00:00 2001
From: "Thomas A. Christensen II" <25492070+MillironX@users.noreply.github.com>
Date: Thu, 29 Aug 2019 19:56:14 -0700
Subject: [PATCH 2/7] Added Julius van der Werf's example to test agains

---
 Excel/van der Werf.xlsx | Bin 0 -> 12285 bytes
 1 file changed, 0 insertions(+), 0 deletions(-)
 create mode 100644 Excel/van der Werf.xlsx

diff --git a/Excel/van der Werf.xlsx b/Excel/van der Werf.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..2cc896006a3e737aa02ca23d007d73c7b2c7a796
GIT binary patch
literal 12285
zcmeHtWm_Bzwr=AdAh^2*cWWTHy99T4cXxM(1Pc}*xJz&g4nZ4tm*95Do|&^}XZC&W
z4>(gFs#@xOpX$}i*HVy%fJ6sC17HCF013dR?gq^R3;+m+0szneu;5z4cDBwYw$6Gg
z9`+_qIt=bM*2MXc;MBPQ@YnnQd;A}kKzZ`GTn`gs|EWZqaLPN&S_eVciY0?@sFciv
z#N;QM^s@O}4>#9M8k{65YI#>vbumvqolW&rwLTR|^zP;5LknDZ)sWG3Wh`H^eHiqn
z#Z&cB#Ncd@wvGUUrqfGVLl<iRv|*0)S{BHgUhd(ymH`<q`sT$|9hd=xPm#PNbB0il
zzHFXrQ<(tD!n{fS=^<8h!vf)pn&LU}x5&5}ECSFRy4+}E`%;mDh2B;(jA?Z)gVe13
z`WoDAWQhH|1erOVNQ-T5v@u`N$r09jvP$5BPT5yAI-hG?b96BGbL&JEB^A_cEz1N$
z$AT0-fCp=uIgDB=W}w|by<kW4F>D6B&%6{}CyW(b#dB}djYBc<L}bQiTbnS6%l8@k
ze3r`}SUiS}zgFyu9Cl$iM|^g@bKlY^1nf{U=VsiW!*9AZ14%#6Fe2CyOl_T=`T0+q
zo7X6d3#|!5?XFME8k>n{@~D5p&TK9cbaB7r5-(nTRkoKG2!O(0EVf>iiR|a=H!`mV
ziSTN%dX6U6PK*q{&i}R4|HJ0_Z$G^}QC9xdQbSKA9>WH2W<l{NLeg&UC0a>TeElVt
zk?Ug$$nij(RQM<=xIqx&e(k=G!%HB(n1ex*>kYQ@sJH06WDRZ=p=tLHE-<u|j;Ue}
z<?FquuCrIO*BRnc9@MVwakO8WO7f(KR>{Ps&qb<H#u?Oc;m|(ghhhRV0yO(&HC7C7
zE5T-kl#eSys~S0T50b{S{N_@N_K^g_xn+)~GB5_649ylQy$7vHuI_MEmCd;=sti9m
z@{oGy8C!q56v^sDd-A52%N|rBWqW(eGA=$qk$dH*S;uxhlH<|K3D@&w@N6hBil}7e
zHF*7}k?^HTHRwVB0EaLD0P?G3+^rei>>RBO?d+_6g|u>28@ntX^k@G4XP90$_#O7#
zDH_NLS{vnA7f_l+l95a`e7aQ@rE2|)Zybr_l1s}$i)aMF!`i?&_bzwX^i69e`$<WB
z^i64m9R74Ja2?~&PU-Sc)q-9|n08scX&Wdz>P0(xWA@U<3vzU5U4Q%m@YgLeEtrX)
zZy-lHuHsvVIPOw%nlyltb}$-wub@2?sx({QaLAPhGWe0jDnzOq)ODDNn(4&iTBpt$
zo=dMLsK;R)4^tQ023^^+;fhb!-3ms`QEqbJ3(g)BOTNnq#B!ewT;fH?<^@|t@f{54
z0xm-qVX3=B!!meO$N?5VN_L(RadQZb3srgH1kI7$HRQU<?4Auu1=|^85LS1KGi^OG
z<APhrYDyBj+le5=E^A*N)tS<_fBbB<W^fzLn^XFOO_NSlky-T1Otk(`<Ijaq>@aI%
z`VTtjyn}V)v5p-mVxm@j4F~YS-u_|3us!ogKa+@~O}>=uRp-B3c{j1X@>O7VIk?@D
zi+Z`K@FiA8ftqdt>lU7j5^iM>p&Zu;-Nv#q@k$-;<SU#7jjl95ywF_W8|_NY<R4kC
z_&@q$CtROIzTBD^JcepAPm)xTK(mhDy^Nm&T1dI6FNEMAo<Hppk-XD8BiMpXRydZv
z#ZUBaI|ZV`-9MwylQxPt*wRPwH9Y2jFaM0!s)jD6q3f=@1~rLnFBS~9?yfS*|AfoK
z?sJM&B=7E4MNEaGaV4?Ywctkwd||`47O3lKnbcuzNr@Ahls3AFxx!d2D%uwQ)aKW)
zC|6`xu7tGrU6?R=y6q;QCJXwOlyhbv7oY9xj5oC~I{o8(ZJ<dgo!`EP+D`CC3{h!;
z3zj)#7QYQ60TMQL6=IH2F#(J&ERZp(vn^1bgaavIkGI!@eILY+Pu4_1><N@b?jLk8
zI%Wlx#wp3N?!|xsw|GYeBYZ|^+5mi^y>#B=^=iH+nD14#L^1b{_OOI5I2Sj_AGXel
zCh73#W>UUwVXnxh)s0m_a*w_NTP-%O)lc9iRQM^aqRBtgoMMt^2ju}rDs)!{J1kza
z_;ep5YhfT?Oru7IUS2tCk+WC5d`if9=wbA!*eoSYdxz)bVfy7p-vg++&Uji~%(dq*
z{Lv?TWx?;&(EigD6Je5%AYYAC;ME*)0kB}NrudJL_3sAxukZ!-nrOb>`+xf=PaLuA
zWkM9a33(2h>~LW5Q3!P~om5J<$@2BLPv(*_gwgWZI<->b8TuA)Tf;2TeM1=0y5#3p
zaZ65wYldMg2_IPQYGPysb=o&FN*4ajT-rK^0ty0Zes@m2RoyaSvsgcqKfJ^+2)t+6
z?la<|LKukyaDl5d`!wEY-)f<N-%!>usR(&$oSj%;8a+%YKp?r1I5ID8?8U0J%#sbr
z`!t!bjeyB=Hqj3y63{`ZcBVy@fVyS!Fq(OdKQ5H?@L~9TG!?ss9j7B(09v>gZJ-om
zi%TQ8Sr1Lj=*uTW4%a~L`cwwTMoVf$PukV_f$!tB%yutbgSGR!0)MMRIH~mn&#x5<
zJQe_e{d&hg-MN#wiHWll;~y91Uw%D1-O4UY0zI_zYuyW2Oui~9suO!h`SCEPVewFr
z)MDkmg{VG$f&M1A!`q7M&Ql!h{vCJxbu}Bb{;#V}JKtGh(pQQbRST^x5<|~h<fB97
z&|9uIiM0i;h8n+~m}F5~>{LK7M$XmAuUJU&Bl*5m+~87;?<Y>lfQ7j%QjO3pIkVGW
zr8#eZj&T>BI~yXEQ`3b`8%j#GeR}|Pm!_6rdeBl((Pl3TA7-mGK&DG4^hwU}#YA*Y
z!ys{DaD@298OL^N%4U16TFYbTa%5*JP93kx9kGDPYg{&>=A!ko*d^ON4)y!PXcbV;
zLwfBi+)Fu&48Jdm1;)K(c7Iek*qCO*<@!%vP4P?bb!V;@O^(nC3TM0;sAY|e&eX(7
zblsh?@Asmt;_20n6rePwy(GWloa>2<2CcN?JYaKOSf%d*$LQp?J)cWhPc?9VzTjmr
z7c01`+geRZ2+N0st)<Z>AIXzjU{&}&^W>@0cx4kJrsR`ixB*JY`W^PbxC~bq2?74J
z_BVW>079d)OhJh}B|Y(@86kvk2p}`O>ctCZbTr^Oks3_Lol%h^wxF|zSW8HyebTa?
zZpo5e@^0Vfk*8pNnk6XiI26f5AKxYD22R!ofc@6$16>J`I11Rax*CbK2+Zx}%EhOz
zA-vY#O|T)8kCpOo87(ls4uA9b5hDEYp?9S1J!6C(hCckve!mg#%pPtP3LO@Vuj$By
zA7LLItB0m6b;dY;==XeiVlR|prWS%dpbf(8d9$Gkwf>5bfhxwPaEOMg%$ovi)xJ6f
zLGibWJgOmD)c!hoOh;u~W%1wg5h9dTzr#b(b&Wxqun4fdH>gAi0UK#>dawGRZO-GB
zXv&_FB^(kk?5ZU(AD9|tbY%Gz3^EXsS~>)hCTcsAl{pvegjX*a6mI26N>>XKV#h{6
zhUOx(=obHu@4>;cO{QirC4$MIH%u2Wj^#ls&6o#9lNo5`ssmFZ-QyiL4%&Zu0zKUc
z*i9OO(h>dmvg2f~9!5c+T|2$8fg(>1KGs&VMzf;QFGUA$><EjnuxLrJfHq(sYDe`*
zeJB#spgLGu)lKw+5PfCBZcrkZ5ToT_=&#B(#|(Lbn!S9zZsz829u9~Uug#<>kh36d
z+eIiUu|!RDgAD72`?Y6S7)KJ|M@LWEWM;_e@~vtVLMxXP0zviL&MH=?xVRaXr<D~&
z1z$$n4`UZpkbF4xvj~LQquarX_e-FWC#jey@^ODay*Pb1v%+yOxAILj*EAm)w2t3V
zE3!cNG+Y!Aqv|)!#Js{*A{x`QvsZ*lXgGbA&@S5}9w}?Yq3Gu5!OA6)KToBvEvgeN
zQo`EY)o`E|>g-Y@ew(EM%7Cay)_C8-yHTM>tc$(@drl{z`+4|hrOB7Tu0(TJ!L9Z~
zG1H)~M2lE|sK)7geP(9xk=`5qN}6vqG6(^jn1tYT{7TI_`Gf-7xGSGz*$mPI&}YPC
z$<IkFVfQp#?l1K3jAqnUXuqJP`00Vi79K!Di|2O{NWO9=ymI!OxIG{{6Rq@TU<U=2
z@or=mtid9Ky&<|o!LOw0dfz&<Jsw%#g({g1Dd%!u9Hg9MZevfbDggM!A^tCpMSC0j
z^jAC~P+Z#y9X}#0pk*r9ieDtcic5;81H^s1VDFZN_u-ADxjJcSl3<zKwsU+h8--%j
z7^~oHd#Wn38&L$b8&L*Qs6-JrII5RcB@!x{QTRmEdbm}r5cxTmf4J54$M{=VKOxIJ
zug^<7a5xFwl7lK}0t7hBrT2q_mLUlV=Ju-z3;M>!nPh`iWW9@nOw3}??dY*hp=6tT
zx9!UA;dq!t3lPC^XCDm>1aC=%kC&2Q??_0^qsGgUfCQ+_71M;_6zu#>KuNZy+Cgh3
zN*1L4LN5DvZz?!e(AuvOi|uo0QFXRpjQus0vh%v&@Uq2w@m#I@<rWAx5o5cCzffc|
zD03<6i(3wp#@1edGr_(+#n3HJDI9r2lH?KJIcL6s?sfL}F{29$r9<vD{qHCK)!F(Z
zW;mOh*qAW>asDG*9BGV3;c=pOV%-ZMI=MZvZ^luqtc_bIfyj+>lJV;sjuh2cSrXgu
z(7-9Uuau|@Kc@-V0h6Z$A&}aC#=()-?GsH>jn0!*Y?EhMOKG|yMTtMR7Z=}r^!2=S
zKA&n&_Y_Y@=uS-4>rl7=W=GGblTWAmS#eRbdCN;AN28j++1?@0JxVWq@^wVD8G;>{
zhqx25&Pyed<X)o1?1K$WX!>T#6|hPfJs{Cd_6c(UpCD)@9rN91r3!4a#9Q4!Lj^A%
z780@}q!cMXuhKU8b>;*iq-g3@^|ES;yz0%K8lmuzJ*5)@lEn*ykIcM;8(vd4)NkMV
zn%_}8?P4`DVos`2iflt`2;1hqur$l)dw-a8y4dfsoFaw|()9YQ2T`IuaCP{k+1W5v
z3}((2m{Qa*&vTe!8*hN5a6g;Zs(z|Eg1IuGr&HTZP!<<5m#gAe6-tEIvL=e~^*s&0
zc28K2X|dViTrRqvR?nRDv93bnPb&oxNyny8hLCzox8_9*<L<p!v=zjl_w9)^<BIV}
zUyhWNWj>H{Zfi!Lu<!%?!Ij7*H-`EU<tJcknxd$=*Ls&IreWXOg}$ZKZ&IN_l2X0>
z${+AtL>Ae=#8Q}|qaX1BPwCp-uy;tN;qhJWEGIGJN1)g9+5Tl%C(!Tb-d&GgO>;f1
zv@HyUvhvAnoZidR6V{4e=kwXYIp_Wi)p>Vs=i}W-TIb6-&sEjNG6rLZ&&9=l64t8E
z)vkOT!Qv7=LHZd<MD<zay3r|RltI9ohjTE2oL(lO-b6vvJ=c<xkwI`>E#iZ2m#X9L
zCHxQlP;GpjRtqge=gs7?8DGV!p_YxJU9sJG*2*$R&d%&V%7BtCn<H8yCK9Yp7^KPp
zH@mF4>KjT6F;ylzNRe+NinzL_lsl`1qY1vW3Wsmc0a-R=sv>@#pnMGqpQaFY)TdBl
zs1r5}AL%wB4qm(%U9tJvQ(`p>^fo3lczq5g-|*7y{t<Hq#Umwd5oUIRjws_T7DvA>
zamq?YauLr+n#+45x?t%Yp0F$2_nq*Fn{OTnb$U4xr5*RCnxn8;6M(MsF$#Q6i!zCW
zy}p!hNJJCdf+VrL;}&>DfWr6+HO$M6a1Gs0rld_!Tb7y0RN?m~V+7T0(HZSTSk8e}
zQ6`Ig4fFM*5v;sCm!&C-t=P)W(f!z7jr>j@CJb0dcAF>FN%d0*j-<-p4UVZ^W$S~P
zB$iqIpm9uteP?oWo@`TPX6F#-j!7MXx>|ULD*}?z;o8nt@yv;RBfG{9zD=|hHBB=I
zwa1j$+9`P$1FK<{=81XYtz-Yr*RFOzUvC9?+yRN>ZHYsC-FQ8UqG%P!XAlYtkB9{d
z499rZ7RmXdX?>EMHJ`<S41q4D){d6671K8!QA73UHcaRS<^&BDr*;02%#%fJvMUl$
z3wJ|AXXR=Barv8lV!6`~4^=1K+1r<gGSg{Y{;wuuRvDL?{a4pl?=3UhyaR?QPcQa9
z%C}W-vyvpAj6n6Z3c2|}woOLEFqJGRwoo-&f{l^hX0oT#Ehm%FZ#qy4P_Zn=O(1Ft
zn}EJ$N6K8OxU*8`8Wo1Qyu<i3QSl&m=^hs|EoatT-jZ300F`*>PNN+ooF-MYrwKe=
zshBzR-BWSXNQ0(cn8r4cF~w7K!B*6JBZRJN7OcCF@nO>W!=%DNTArqu!!|sV4F-n$
z_Vt=*`*zs-#rL76_{qCQ7{zrNZPAtRV;N$eNvx$AU<!teti*B7$_p(H9+oZ^W$mJc
z^kmV_wJ}UnOo3SIbiAXf4;uX~8DTRl^38_Lc#f<GG03P^qNG$!0hE-=6TQ5)clDOA
zjkIKJR(EK7*gbRZ(9L8sT3GqYO|CeSA8l{l3{r8hYQjX4d9;({hM}1g(!G*<>GkOS
zQjpjt@$+rosr#ewYi?!=z9Y(VKyl_BbZ6-zvj_!YD}IzC3RqXJw0|z^TO}wbU%>Bf
zCyI${hv3=eL@oUAf$<$oY=nPQkmJE9b3j%VM#A=;+lKzI`CMD%6nAng8BdggYks=%
z$8PqB%I}l%?})JF7c#r@4jK5;%OI}X<w2Z8v%^l%7CUA|eY6KN&Pn@Wyjw_LE`(=l
zXg}LN$<aGjS}`=gc{}SN)wVw469v)5ME~hKi>yYQTu)?et>RK0D;+GN{0zsfjf)r;
z@Vt^rU@+&qri&omp;XT8{qS<ztSO9J0H3G5A6K!`o`R5F8X5Kd6y9Jmhkn#@l12>d
zY8y8d{fjS+H4I2<o<14Zpc6$UNGT7s;Y)}sw%W*3mG*%PE%KUvf}-}D>BD<4ZG#C>
zSBtL$a||b&+w`L$s8LZedf&0?*Y7tvv=_V`PEm<-9_oj4HEdihLrm*|)R66R1rtX<
zctUh5v`&#cK~<CZd}r-sNN*NrZ50kV;p-h6wO}nNOneJ8N?=l?cuN}bo^qdv-x=v|
z)A$&70B$T>8(#jl11E5l`Ed1G@@KvU0FeF>ew~~>tWBJL)oC*ymSUGUP<$9C{7`Ny
zXZsS3@Dc;BJfR>Zg1}p$?xm!YWrMMT@;MzNpS{A8jfW&T9hObkucG=$Ja+6&*;oxP
z8rHvikd_rG6;d~yCkTs>NYirA+;1|wnGKe*V^WwZ#iJ8Q#H^ThJ}1uHe;e#cML7lp
zmnhfGdN!W5wamYbJ;ZA!12?5D<@&*-g7<ZWCU9iYdgD8SkKciV43s4*W*Zkwu+f!6
zrxU`E+lOg(6BQylOh@->_eNBub#{m{S(No|i@;L&9I-2@@}KX3*_cv4YNj~-4~NR~
zq3$wImb1K?u;T%DEQQ3^k7rv)F)b)G7w8`O{XW#vv^37g=#GrekwV0lGF7+KTT}^I
zG-vjWw9*ITV`!$taPesrA;wDFxY{tOwQmEsE!39{#O5M1O)z$n2{awdA%84^c;vdS
z#rALr>9goGvQ*m{T6fL%YG|z{Byo5(^)25Z1kwd4Oo~L{L(APfy3cnhLW#IYFH}B^
z9mI5p5l4@&DoBe-_ZT{vA(e9+qHi~!GCs|7G9B6mMkcZq%1pA=?2&c%f833K2(5e!
zao#?=)RUOQ^V?W`u41r~3Qww6p)hWL)?rSVimVAvkP0$G2U(AHrt)r5m9iU$($UbY
zNCafSUHT5-@RCFbb;zJnQGlE#esn<UwX~}`b)*ui*^$KX!&u5Ydd9zGjA^h%>+sg<
z!X2-Sq)l;?;p+~wj4SL5p@Vl2M+4KbA8UEU+G(I=x5I(s&Au>0FlC*RpeHipoRjL9
zGgp05zvv%^aRvoZ_Blm7v|I)lBD)hjwr*#(18<0~y=7HwTjMJ8fo?N^wblZ;cvgIj
zE#PO&=^2Hf*%<wCCe6iQmS{U15UkuHk_I?M9RZJ4>7E&YrM}VM1H590D#R_WuzIb(
zP4=qINX&Xuy$|y~p{t>Ydw<y{rvYeyVBi3<mfyH`scyfQ@IPYhJ)T{<CZ}p+Xi!Io
zn|!27&UO+cnU`Z7O>-$9O%<*WaCFe5>y*bA+6r@-k3?qoosr0{LmdeH#=IIsbFUB2
zh)1J&N;9qUD2)J?Zo<k&(eQb02W_R{ojCnOR&rstR}oQM5)JN0zuFx3UWLX3w|nyX
z5Q)jN#}((Rt{HEm_w7D<!s75RW4e>*PcEI_ta_l`Jl@MRR?f`$`J91~6g+Xaky+YL
z4|#q!H+!V(0@`HrUdZ-3DD;-lX-H+f*5c0w@t@7v%=ecVMRI{V{K&=6Hu$%;{O1sy
znHw3rS#pM)aY6R%F{w8WgmFy+Sh}3Ix+ry6SfSo+6Uk^|qNiV0gpFvp;6hiqg6F)t
z@fcES{=O}Ub92B-_3Ak~UdLGI|MU;$2974i%Fd1!wq}303Fj$Ud2nVF=(W~}7lWxJ
z!d!1y+}c)1ZCFZA@E}5}qpBri-o6Hu+ndUDeU_*J{!?cjn@kzGhETFQef;9G)ngjs
zm6mS%jIy|^{L4aOku7kmf#fn)DE1*+yF?UHHZShNvPg%Qa+i5EV`gX2$0kYFF|Fv)
zsWhG=paCjvzPS}}@2Y<>-vCO>Z7scQQlTqg^d>-;feE{m4vtdVki*2NI&_Y2l3DsE
zqi$fzxpRn9ow~EBcRBHZM9pX&xqT~TO$Edr`F_-dwAVn(tn*}>h!&%;^<c7UvQFw&
z)HMJbUpmJlUv1wTc=uvxpSH;@2$Nbw9e!$JqSuXz1AkFP(wVeS9NhQQhCPido+<f^
z*nr%CZv7nD`SPWP=v8e0*>6o7&o`ua?Uh-h{BN(!*wMi4mrNP|`Qne_=lcf@J3LMl
zU&e_SCQurOAbG``L!m;!n!<1XR(Sj<TM3KXQX)-TsawNJfsZ%7TT(Gj12LP|AsdO=
z4_WI9;~(vuB=Ky^CoIv~4;@6bDR?j}*&}T$vFd$AelVx7lwuHx^4b>cZtD4YG{}yg
z{TM?G2&fVEj?J%k>Z_Wp8p_=L{==mU!<E&>ap%lJEoP=lt7^=o0Am0YN{?t(s^fHG
zIxAPVJMawEm<@*!6J$DZl`hS#IHZ|nNyEzX)ife?({xv(Cu?mt06giE4sYOv?%`{d
z(a+;14f=<DgV9sYrFK&{`M`F;xAPrX{f)7S$J7VCLAIT<QPc}D`L16;0+E{~iVn3<
z+jPvH=%QPK#37zB9`&Uh>V17Hp<C8Y1`>q@B^U{Y+~AHem&)y*DPwC%_~j~=#!J`?
z(c5+CG^Y#h@|Csp)ed{7<!pFP5iJp@neDDH2P#}|a*AzT!{9d70U=30QNm1bk|egk
zTQqZ;J@YR^$5;~1gpE<>Op8%TQGa?b#kHRQED<Pmfk`s0fXJxdubV3$g%&@)tDE~c
zx)bEKf*O<iw9ffKQFL&rEpAfd462!)4q6i;BW^>7Xkm?>W7!n?+*Ei^oQrR-|H&8O
zEc5~9SoaHF8<DXuJHNu=w^F7b?Q1+#0yv_fR191>Hn}z?UF+dk<v)ca_UGCoQ#1{B
zToWH4lyoi-bpu@;^o`MKH3_AK2Pw=W5qq0{8ubBD+=mgPPbFKa;3fuHN_^!#eTm^O
z>OJK@A$%Uy^&{uMfwoS^7q!(Ag}qnXMXoHtU-MZf%Q)d1q3{i`NhWoRB}uw>MW;|g
zM<oea^=@;u`oJ7{q({NZN4mLH9V*fk6_OlwmYXn@G_@~t)T0as*U@+~`~_^F#$X*r
z9}z+&{`(w%r6s3&iI1Ist!YY`^<KXp%VAUA{{Co)IGqxbk?6gPrzXot!o0)+@M1Y#
z*Yc7_ML_dp&$B|2-9o8&vmO<tfKJ`n@ge|U0%6fNZuk@0zI2^vj4oUwLmEkdlfNW(
zMMDh5?t)TkUyxMWkJrHnIkMo}$mPr=Bk7@iOr8n3w9p^Rg`HYb8Ad8)$#@|wb#KH^
z4^8wi`I(l(`YdUbKE<pXeKH0saNvuHdvoWqn;uuj$jjf$$gq1^OB8pld_UBX+rda%
zDD=QMe_`qC>G^y?68&EOtQP5@@7Vn!W=(rpTQ_@*h!|^bPh<kd2j@ZiGuHR25ysT;
zxXWYT=KB|@llY?Vk->7NC6zQ7d>HR~xD%T#k0%%(1^b#BLG8D<c<F*Jz5o%#eGIc5
z8XX7IdiWxFwKms}v!-;#km|VO2sY^plf!#8f_h$a6utJBx#t*gNsi*!NlK~T&y&5~
zup_f>*77JrX@xO4Ssyiu7r1kB^ps=O_?GT^Qj^@=XfmFT&zGOl*-a#NrAJ1g2^~Y>
zZrs)QP4jT3$Yn#2S&~)*NE=Pv1(`nA6;M!VWopkugtFNil<4<sOlkNhcM8AM;$B-s
zImIkJ*OnSuP;l(o35gTasoKP^XNW*6DB<>`6bPJU?+N75;G*#iArx?$GgWu2mCUYx
z7W$%8m{(oHFF;HeCg9D$5b_?`VWiLUSjbATSP~2~1a(jdKiQ<RO<!kt<{|Yn4L$;4
z=q`#gmAv0Y-&4uWaz}UkbnhhiU8D#`tkBfLz1aksE>oDU&r9}#z)+zepD6&eo_9Gp
z&6|DghaUX`{dcYv<7*jT`>WG^_Zqjz{t>f&4YB??hWN)+>(61tA2Y4sxKa3SCY0dQ
zz*CfK+bK_!gd$y$GukV#RO5#ro3}OR#ddZKfq`o_W@%5RMT-K(ipw+?YdN7xn8brP
zb}(}mO>6_Rn2=P?nBu1rJu-+7p3Nm=Y=*<~%g6zuMO;57iqT!As`o8~Qu*UY(tZrU
zDW_H|gN<VFDjnv}m)s8d-$&8NsR{EfycJ;$J9_+dyyVK{e73`pAT2<b9CJH3%JSAo
zTrw#6l%O+AX=gC^K*~q8T<>oK%LZ>!)6wfUA73N(Yi#*7gfTL(HWD^Bu(17=)N4<T
zO7${fhpe}+0fX&M3uAdA*>d^>qrKYw<4ub~9+Sxj-EYJ)RDfARp9oh{8=jtz2p+0v
z#+hMTK+e6^Ob*#Sk`oVEZ62hlQEv)dnLoV$#4+vL)Nv5wy2G^3U8M=DF`+Vi@!f+X
z>59j0!%(>0F?#4wNq#zFmqG)$Gqkutu<8nmeZbAl<D`R2)G~DPfi^sO7-W2Gacm!A
z>}0d$4Iw^dXs@5~e_~vhBN_CfWA}f1f*$DoO$mbf*9QctxPowQ@bCEzOTm)ui&2c9
zG<-YO?)&jxKRrM`yBf#s_*=;Y_fX`kzv_hVt7K4LN1VoXMhcF0_D+mOc8(^$MDx0s
z^WVwltLVHF6{UNa&_h<CUWA9dvg(%VQX7P+f&w7AA~&=oPJS3BC6l)tuVyu9RmgeE
zj5X_x0m*YtMUQ;~IZ_gF;2{F};$_aIY*~M5(>-B(^eYWnRFE6*0iZ`ON1~G1;goR+
z>n#PB;G=gXrpU{O6fA6#^sPU@1!jKZE=1#$La+oW+Bpa>@Pfc&VVk?b>9%K-mcgG)
ziSx-4Yrz)0_gWzAL-U})rf~70NuIAYC%PZGCuqouOOEm!cnLOR7qYK5)>g7CnND0=
z^*Q$qn;uDgbp-fazTBr+w;K9iB-aL8<f`qrQKx+AYDH=jEq!P79yfl)GG!y7l)kR0
zzf$KW3;q_g&UWLqx4cJq=z6FJdD2`b*|{X?kwJG7MKMhFwx7TGe5jmQk%eOD=)Pj~
z+^9o#RWIBq^CjX9*2837I6ku+rLhPf*aiBhr}RgJzw=TCl5i*yuiA+ET79GZyEYow
z+y7S?U$yX`M|R?SyCo*XkW+|zB18}#H%I|7pxjVG{b!W{SYm^jY>`o34S9HZ!qHq@
z!{p3@=RwStm;1elBnY2J?`QmEpF+fTJ~|X6o)y-0iu*iB#}@-C5w8B6S&$;ue?v^q
zTPX`4&DhZ&N*uCWL^@ebDc>fKwMn15YN~e`XM^L^#2VLZontZW4fxQPnzmGXVp#T>
zw0(5yqFz;bBBY7}g4;6(HyOn8*~a9y+R)q7`{G(|L+g@T->S8g?5N2727SGnf&zw=
z>_$tgm@syj!+v*2idE<X!tYgeF%CHO2NttIpVJYX{oHNT^KWfXptGfZQu3mOS3@t&
zMAyqpXbAK@M$-fG60Orijks^dtZqaZ#UZeE`$_b{4c74;b_XGzEjcKsd3#5=&j*>#
z<+zpt?0B5UICV4R@2R{9T8+CCLX@w1Ld?7+THUQ+b>}?zRhE7jm)t?g9r7d10PR>F
zNQ^<5RbltDmz_*}&znDct}Z+keHrd_%vdQ_pKG4~SAD_2>0gV$zu!sz*WdN8>)&iB
zD@gy9;IAu^{{sHH*1X!pzb#Du4*Y#>>`&<2YeMkb64~$Izb?A`2?YQQ;r{^tKh|D;
z=lOlv;ZLR_#Q%MW|5$zao#prO*PkqluN}G9pYr=K?016S8+d;bXyE=q@Mk;kcj)i6
z{hv?|!hb-2uk?Rs_`QhzlVR$$>V19szb$EhXZbx}`;&#9;t!VpU-tGp;qNKXpM>uy
z{~-KJLi9WOuW7=cSO7qj3IO<zq~Uk?U;X;O!>4Kf4gOCruOJQes<OX&uSfvZSFKN^
I`}OXB0Cu&I-v9sr

literal 0
HcmV?d00001


From ccc980855e690acce7052459b6933c12eb8520b2 Mon Sep 17 00:00:00 2001
From: "Thomas A. Christensen II" <25492070+MillironX@users.noreply.github.com>
Date: Fri, 8 Nov 2019 18:33:01 -0800
Subject: [PATCH 3/7] Add Julia stuff to gitignore

---
 .gitignore                   |  30 ++++++-
 .idea/dictionaries/cclea.xml |   9 ++
 .idea/workspace.xml          | 167 +++++++++++++++++++++++++++++++++++
 Julia/install.jl             |   3 +
 4 files changed, 208 insertions(+), 1 deletion(-)
 create mode 100644 .idea/dictionaries/cclea.xml
 create mode 100644 .idea/workspace.xml
 create mode 100644 Julia/install.jl

diff --git a/.gitignore b/.gitignore
index 6d5306c..2ed722b 100644
--- a/.gitignore
+++ b/.gitignore
@@ -201,4 +201,32 @@ dmypy.json
 
 # Ignore results files
 Results.txt
-results.txt
\ No newline at end of file
+results.txt
+
+### Julia.gitignore
+# (https://github.com/github/gitignore/blob/master/Julia.gitignore)
+
+# Files generated by invoking Julia with --code-coverage
+*.jl.cov
+*.jl.*.cov
+
+# Files generated by invoking Julia with --track-allocation
+*.jl.mem
+
+# System-specific files and directories generated by the BinaryProvider and BinDeps packages
+# They contain absolute paths specific to the host computer, and so should not be committed
+deps/deps.jl
+deps/build.log
+deps/downloads/
+deps/usr/
+deps/src/
+
+# Build artifacts for creating documentation generated by the Documenter package
+docs/build/
+docs/site/
+
+# File generated by Pkg, the package manager, based on a corresponding Project.toml
+# It records a fixed state of all packages used by the project. As such, it should not be
+# committed for packages, but should be committed for applications that require a static
+# environment.
+Manifest.toml
\ No newline at end of file
diff --git a/.idea/dictionaries/cclea.xml b/.idea/dictionaries/cclea.xml
new file mode 100644
index 0000000..8c49cbc
--- /dev/null
+++ b/.idea/dictionaries/cclea.xml
@@ -0,0 +1,9 @@
+<component name="ProjectDictionaryState">
+  <dictionary name="cclea">
+    <words>
+      <w>beefblup</w>
+      <w>blup</w>
+      <w>matlab</w>
+    </words>
+  </dictionary>
+</component>
\ No newline at end of file
diff --git a/.idea/workspace.xml b/.idea/workspace.xml
new file mode 100644
index 0000000..bbe788f
--- /dev/null
+++ b/.idea/workspace.xml
@@ -0,0 +1,167 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<project version="4">
+  <component name="ChangeListManager">
+    <list default="true" id="816df031-a834-466c-aafb-fce0c72415a4" name="Default Changelist" comment="">
+      <change beforePath="$PROJECT_DIR$/.gitignore" beforeDir="false" afterPath="$PROJECT_DIR$/.gitignore" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/.idea/.gitignore" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/.gitignore" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/.idea/beefblup.iml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/beefblup.iml" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/.idea/inspectionProfiles/profiles_settings.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/inspectionProfiles/profiles_settings.xml" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/.idea/misc.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/misc.xml" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/.idea/modules.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/modules.xml" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/.idea/vcs.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/vcs.xml" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/.vscode/launch.json" beforeDir="false" afterPath="$PROJECT_DIR$/.vscode/launch.json" afterDir="false" />
+      <change beforePath="$PROJECT_DIR$/Python/beefblup.py" beforeDir="false" afterPath="$PROJECT_DIR$/Python/beefblup.py" afterDir="false" />
+    </list>
+    <option name="EXCLUDED_CONVERTED_TO_IGNORED" value="true" />
+    <option name="SHOW_DIALOG" value="false" />
+    <option name="HIGHLIGHT_CONFLICTS" value="true" />
+    <option name="HIGHLIGHT_NON_ACTIVE_CHANGELIST" value="false" />
+    <option name="LAST_RESOLUTION" value="IGNORE" />
+  </component>
+  <component name="Git.Settings">
+    <option name="RECENT_BRANCH_BY_REPOSITORY">
+      <map>
+        <entry key="$PROJECT_DIR$" value="develop" />
+      </map>
+    </option>
+    <option name="RECENT_GIT_ROOT_PATH" value="$PROJECT_DIR$" />
+  </component>
+  <component name="IgnoredFileRootStore">
+    <option name="generatedRoots">
+      <set>
+        <option value="C:\Users\cclea\source\repos\beefblup\.idea" />
+      </set>
+    </option>
+  </component>
+  <component name="ProjectId" id="1Q4ogIZMNQ0ZN4Ces5FEsodDyuh" />
+  <component name="PropertiesComponent">
+    <property name="last_opened_file_path" value="$PROJECT_DIR$/Python/beefblup.py" />
+    <property name="settings.editor.selected.configurable" value="configurable.group.language" />
+  </component>
+  <component name="RunDashboard">
+    <option name="ruleStates">
+      <list>
+        <RuleState>
+          <option name="name" value="ConfigurationTypeDashboardGroupingRule" />
+        </RuleState>
+        <RuleState>
+          <option name="name" value="StatusDashboardGroupingRule" />
+        </RuleState>
+      </list>
+    </option>
+  </component>
+  <component name="RunManager">
+    <configuration name="beefblup" type="PythonConfigurationType" factoryName="Python" nameIsGenerated="true">
+      <module name="beefblup" />
+      <option name="INTERPRETER_OPTIONS" value="" />
+      <option name="PARENT_ENVS" value="true" />
+      <envs>
+        <env name="PYTHONUNBUFFERED" value="1" />
+      </envs>
+      <option name="SDK_HOME" value="C:\tools\Anaconda3\python.exe" />
+      <option name="WORKING_DIRECTORY" value="$PROJECT_DIR$/Python" />
+      <option name="IS_MODULE_SDK" value="false" />
+      <option name="ADD_CONTENT_ROOTS" value="true" />
+      <option name="ADD_SOURCE_ROOTS" value="true" />
+      <option name="SCRIPT_NAME" value="$PROJECT_DIR$/Python/beefblup.py" />
+      <option name="PARAMETERS" value="" />
+      <option name="SHOW_COMMAND_LINE" value="false" />
+      <option name="EMULATE_TERMINAL" value="false" />
+      <option name="MODULE_MODE" value="false" />
+      <option name="REDIRECT_INPUT" value="false" />
+      <option name="INPUT_FILE" value="" />
+      <method v="2" />
+    </configuration>
+  </component>
+  <component name="SvnConfiguration">
+    <configuration />
+  </component>
+  <component name="TaskManager">
+    <task active="true" id="Default" summary="Default task">
+      <changelist id="816df031-a834-466c-aafb-fce0c72415a4" name="Default Changelist" comment="" />
+      <created>1567039339175</created>
+      <option name="number" value="Default" />
+      <option name="presentableId" value="Default" />
+      <updated>1567039339175</updated>
+    </task>
+    <task id="LOCAL-00001" summary="Added all results files to gitignore">
+      <created>1567130353481</created>
+      <option name="number" value="00001" />
+      <option name="presentableId" value="LOCAL-00001" />
+      <option name="project" value="LOCAL" />
+      <updated>1567130353481</updated>
+    </task>
+    <task id="LOCAL-00002" summary="Add new gitignore to avoid uploading junk from PyCharm">
+      <created>1567131158833</created>
+      <option name="number" value="00002" />
+      <option name="presentableId" value="LOCAL-00002" />
+      <option name="project" value="LOCAL" />
+      <updated>1567131158833</updated>
+    </task>
+    <task id="LOCAL-00003" summary="Added PyCharm's config files">
+      <created>1567131230009</created>
+      <option name="number" value="00003" />
+      <option name="presentableId" value="LOCAL-00003" />
+      <option name="project" value="LOCAL" />
+      <updated>1567131230009</updated>
+    </task>
+    <task id="LOCAL-00004" summary="Create initial python script">
+      <created>1567131261939</created>
+      <option name="number" value="00004" />
+      <option name="presentableId" value="LOCAL-00004" />
+      <option name="project" value="LOCAL" />
+      <updated>1567131261939</updated>
+    </task>
+    <task id="LOCAL-00005" summary="Added Julius van der Werf's example to test agains">
+      <created>1567133774405</created>
+      <option name="number" value="00005" />
+      <option name="presentableId" value="LOCAL-00005" />
+      <option name="project" value="LOCAL" />
+      <updated>1567133774405</updated>
+    </task>
+    <option name="localTasksCounter" value="6" />
+    <servers />
+  </component>
+  <component name="Vcs.Log.Tabs.Properties">
+    <option name="TAB_STATES">
+      <map>
+        <entry key="MAIN">
+          <value>
+            <State>
+              <option name="COLUMN_ORDER" />
+            </State>
+          </value>
+        </entry>
+      </map>
+    </option>
+  </component>
+  <component name="VcsManagerConfiguration">
+    <MESSAGE value="Added all results files to gitignore" />
+    <MESSAGE value="Add new gitignore to avoid uploading junk from PyCharm" />
+    <MESSAGE value="Added PyCharm's config files" />
+    <MESSAGE value="Create initial python script" />
+    <MESSAGE value="Added Julius van der Werf's example to test agains" />
+    <option name="LAST_COMMIT_MESSAGE" value="Added Julius van der Werf's example to test agains" />
+  </component>
+  <component name="XDebuggerManager">
+    <breakpoint-manager>
+      <breakpoints>
+        <line-breakpoint enabled="true" suspend="THREAD" type="python-line">
+          <url>file://$PROJECT_DIR$/Python/beefblup.py</url>
+          <line>44</line>
+          <option name="timeStamp" value="4" />
+        </line-breakpoint>
+        <line-breakpoint enabled="true" suspend="THREAD" type="python-line">
+          <url>file://$PROJECT_DIR$/Python/beefblup.py</url>
+          <line>59</line>
+          <option name="timeStamp" value="5" />
+        </line-breakpoint>
+        <line-breakpoint enabled="true" suspend="THREAD" type="python-line">
+          <url>file://$PROJECT_DIR$/Python/beefblup.py</url>
+          <line>71</line>
+          <option name="timeStamp" value="6" />
+        </line-breakpoint>
+      </breakpoints>
+    </breakpoint-manager>
+  </component>
+</project>
\ No newline at end of file
diff --git a/Julia/install.jl b/Julia/install.jl
new file mode 100644
index 0000000..396cdc4
--- /dev/null
+++ b/Julia/install.jl
@@ -0,0 +1,3 @@
+
+using Pkg
+Pkg.add("XLSXReader")
\ No newline at end of file

From 22b934dc49869790c280ba25b9907daecdc679d2 Mon Sep 17 00:00:00 2001
From: "Thomas A. Christensen II" <25492070+MillironX@users.noreply.github.com>
Date: Fri, 8 Nov 2019 18:46:11 -0800
Subject: [PATCH 4/7] Add excel reader dependency to install script

---
 Julia/install.jl | 11 ++++++++++-
 1 file changed, 10 insertions(+), 1 deletion(-)

diff --git a/Julia/install.jl b/Julia/install.jl
index 396cdc4..3276384 100644
--- a/Julia/install.jl
+++ b/Julia/install.jl
@@ -1,3 +1,12 @@
+# beefblup install
+# Prepares the Julia environment for using beefblup by installing the requisite
+# packages
+# Usage: julia install.jl
+# (C) 2019 Thomas A. Christensen II
+# Licensed under BSD-3-Clause License
 
+# Import the package manager
 using Pkg
-Pkg.add("XLSXReader")
\ No newline at end of file
+
+# Install requisite packages
+Pkg.add("XLSX")
\ No newline at end of file

From 6ac985fee45f3c24bc8f4ce2994e50b59906e7ee Mon Sep 17 00:00:00 2001
From: "Thomas A. Christensen II" <25492070+MillironX@users.noreply.github.com>
Date: Fri, 8 Nov 2019 19:01:46 -0800
Subject: [PATCH 5/7] Complete import section port

---
 Julia/beefblup.jl | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 Julia/beefblup.jl

diff --git a/Julia/beefblup.jl b/Julia/beefblup.jl
new file mode 100644
index 0000000..f519cf0
--- /dev/null
+++ b/Julia/beefblup.jl
@@ -0,0 +1,28 @@
+# beefblup
+# Main script for performing single-variate BLUP to find beef cattle
+# breeding values
+# Usage: julia beefblup.jl
+# (C) 2019 Thomas A. Christensen II
+# Licensed under BSD-3-Clause License
+
+# Import the required packages
+using XLSX
+
+# Display stuff
+print("beefblup v 0.0.0.1\n")
+print("(C) 2019 Thomas A. Christensen II\n")
+print("https://github.com/millironx/beefblup\n")
+print("\n")
+
+### Prompt User
+# Ask for an input spreadsheet
+print("Enter the full filename of a beefblup worksheet> ")
+path = readline(stdin)
+
+# Ask for an output text filename
+print("Enter the full filename of your desired results> ")
+savepath = readline(stdin)
+
+# Ask for heritability
+print("What is the heritability for this trait?> ")
+h2 = parse(Float64, readline(stdin))

From 44ec5ae1da4b28895a23a0eb651dbcfa70f3867a Mon Sep 17 00:00:00 2001
From: "Thomas A. Christensen II" <25492070+MillironX@users.noreply.github.com>
Date: Sat, 14 Dec 2019 09:57:48 -0800
Subject: [PATCH 6/7] Complete through "Processing and formatting"

---
 Julia/beefblup.jl | 110 ++++++++++++++++++++++++++++++++++++++++++----
 1 file changed, 101 insertions(+), 9 deletions(-)

diff --git a/Julia/beefblup.jl b/Julia/beefblup.jl
index f519cf0..547cfb4 100644
--- a/Julia/beefblup.jl
+++ b/Julia/beefblup.jl
@@ -9,20 +9,112 @@
 using XLSX
 
 # Display stuff
-print("beefblup v 0.0.0.1\n")
-print("(C) 2019 Thomas A. Christensen II\n")
-print("https://github.com/millironx/beefblup\n")
+println("beefblup v 0.0.0.1")
+println("(C) 2019 Thomas A. Christensen II")
+println("https://github.com/millironx/beefblup")
 print("\n")
 
 ### Prompt User
 # Ask for an input spreadsheet
-print("Enter the full filename of a beefblup worksheet> ")
-path = readline(stdin)
+# print("Enter the full filename of a beefblup worksheet> ")
+# path = readline(stdin)
+path = "C:\\Users\\cclea\\source\\repos\\beefblup\\Excel\\Master BLUP Worksheet.xlsx"
 
 # Ask for an output text filename
-print("Enter the full filename of your desired results> ")
-savepath = readline(stdin)
+# print("Enter the full filename of your desired results> ")
+# savepath = readline(stdin)
+savepath = "C:\\Users\\cclea\\source\\repos\\beefblup\\results.txt"
 
 # Ask for heritability
-print("What is the heritability for this trait?> ")
-h2 = parse(Float64, readline(stdin))
+# print("What is the heritability for this trait?> ")
+# h2 = parse(Float64, readline(stdin))
+h2 = 0.4
+
+### Import input filename
+print("[🐮]: Importing Excel file...")
+
+# Import data from a suitable spreadsheet
+data = XLSX.readxlsx(path)[1][:]
+
+print("Done!\n")
+
+### Process input file
+print("[🐮]: Processing and formatting data ...")
+
+# Extract the headers into a separate array
+headers = data[1,:]
+data = data[2:end,:]
+
+# Sort the array by date
+data = sortslices(data, dims=1, lt=(x,y)->isless(x[2],y[2]))
+
+# Define fields to hold id values for animals and their parents
+ids = data[:,1]
+damids = data[:,3]
+sireids = data[:,4]
+numanimals = length(ids)
+
+# Find the index values for animals and their parents
+dam = indexin(damids, ids)
+sire = indexin(sireids, ids)
+
+# Store column numbers that need to be deleted
+# Column 6 contains an intermediate Excel calculation and always need to
+# be deleted
+colstokeep = [1, 2, 3, 4, 5]
+
+# Find any columns that need to be deleted
+for i in 7:length(headers)
+    if length(unique(data[:,i])) <= 1
+        colname = headers[i]
+        print("Column '")
+        print(colname)
+        print("' does not have any unique animals and will be removed from this analysis\n")
+    else
+        push!(colstokeep, i)
+    end
+end
+
+# Delete the appropriate columns from the datasheet and the headers
+data = data[:, colstokeep]
+headers = headers[colstokeep]
+
+# Determine how many contemporary groups there are
+numgroups = ones(1, length(headers)-5)
+for i in 6:length(headers)
+    numgroups[i-5] = length(unique(data[:,i]))
+end
+
+# If there are more groups than animals, then the analysis cannot continue
+if sum(numgroups) >= numanimals
+    println("There are more contemporary groups than animals. The analysis will
+    now abort.")
+    exit()
+end
+
+# Define a "normal" animal as one of the last in the groups, provided that
+# all traits do not have null values
+normal = Array{String}(undef,1,length(headers)-5)
+for i in 6:length(headers)
+    for j in numanimals:-1:1
+        if !ismissing(data[j,i])
+            normal[i-5] = data[j,i]
+            break
+        end
+    end
+end
+
+# Print the results of the "normal" definition
+println(" ")
+println("For the purposes of this analysis, a 'normal' animal will be defined by")
+println("the following traits:")
+for i in 6:length(headers)
+    print(headers[i])
+    print(": ")
+    print(normal[i-5])
+    print("\n")
+end
+println("If no animal matching this description exists, the results may appear")
+println("outlandish, but are still as correct as the accuracy suggests")
+
+print("Done!\n")

From fc771d338a94dde9db2d8c5ec46372cf342b5a5f Mon Sep 17 00:00:00 2001
From: "Thomas A. Christensen II" <25492070+MillironX@users.noreply.github.com>
Date: Mon, 3 Feb 2020 18:11:55 -0800
Subject: [PATCH 7/7] Complete port of beefblup to Julia

---
 .idea/dictionaries/cclea.xml |   9 --
 .idea/workspace.xml          | 167 ---------------------------
 Excel/van der Werf2.xlsx     | Bin 0 -> 17061 bytes
 Julia/beefblup.jl            | 217 +++++++++++++++++++++++++++++++----
 Julia/install.jl             |   3 +-
 5 files changed, 194 insertions(+), 202 deletions(-)
 delete mode 100644 .idea/dictionaries/cclea.xml
 delete mode 100644 .idea/workspace.xml
 create mode 100644 Excel/van der Werf2.xlsx

diff --git a/.idea/dictionaries/cclea.xml b/.idea/dictionaries/cclea.xml
deleted file mode 100644
index 8c49cbc..0000000
--- a/.idea/dictionaries/cclea.xml
+++ /dev/null
@@ -1,9 +0,0 @@
-<component name="ProjectDictionaryState">
-  <dictionary name="cclea">
-    <words>
-      <w>beefblup</w>
-      <w>blup</w>
-      <w>matlab</w>
-    </words>
-  </dictionary>
-</component>
\ No newline at end of file
diff --git a/.idea/workspace.xml b/.idea/workspace.xml
deleted file mode 100644
index bbe788f..0000000
--- a/.idea/workspace.xml
+++ /dev/null
@@ -1,167 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<project version="4">
-  <component name="ChangeListManager">
-    <list default="true" id="816df031-a834-466c-aafb-fce0c72415a4" name="Default Changelist" comment="">
-      <change beforePath="$PROJECT_DIR$/.gitignore" beforeDir="false" afterPath="$PROJECT_DIR$/.gitignore" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/.idea/.gitignore" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/.gitignore" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/.idea/beefblup.iml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/beefblup.iml" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/.idea/inspectionProfiles/profiles_settings.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/inspectionProfiles/profiles_settings.xml" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/.idea/misc.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/misc.xml" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/.idea/modules.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/modules.xml" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/.idea/vcs.xml" beforeDir="false" afterPath="$PROJECT_DIR$/.idea/vcs.xml" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/.vscode/launch.json" beforeDir="false" afterPath="$PROJECT_DIR$/.vscode/launch.json" afterDir="false" />
-      <change beforePath="$PROJECT_DIR$/Python/beefblup.py" beforeDir="false" afterPath="$PROJECT_DIR$/Python/beefblup.py" afterDir="false" />
-    </list>
-    <option name="EXCLUDED_CONVERTED_TO_IGNORED" value="true" />
-    <option name="SHOW_DIALOG" value="false" />
-    <option name="HIGHLIGHT_CONFLICTS" value="true" />
-    <option name="HIGHLIGHT_NON_ACTIVE_CHANGELIST" value="false" />
-    <option name="LAST_RESOLUTION" value="IGNORE" />
-  </component>
-  <component name="Git.Settings">
-    <option name="RECENT_BRANCH_BY_REPOSITORY">
-      <map>
-        <entry key="$PROJECT_DIR$" value="develop" />
-      </map>
-    </option>
-    <option name="RECENT_GIT_ROOT_PATH" value="$PROJECT_DIR$" />
-  </component>
-  <component name="IgnoredFileRootStore">
-    <option name="generatedRoots">
-      <set>
-        <option value="C:\Users\cclea\source\repos\beefblup\.idea" />
-      </set>
-    </option>
-  </component>
-  <component name="ProjectId" id="1Q4ogIZMNQ0ZN4Ces5FEsodDyuh" />
-  <component name="PropertiesComponent">
-    <property name="last_opened_file_path" value="$PROJECT_DIR$/Python/beefblup.py" />
-    <property name="settings.editor.selected.configurable" value="configurable.group.language" />
-  </component>
-  <component name="RunDashboard">
-    <option name="ruleStates">
-      <list>
-        <RuleState>
-          <option name="name" value="ConfigurationTypeDashboardGroupingRule" />
-        </RuleState>
-        <RuleState>
-          <option name="name" value="StatusDashboardGroupingRule" />
-        </RuleState>
-      </list>
-    </option>
-  </component>
-  <component name="RunManager">
-    <configuration name="beefblup" type="PythonConfigurationType" factoryName="Python" nameIsGenerated="true">
-      <module name="beefblup" />
-      <option name="INTERPRETER_OPTIONS" value="" />
-      <option name="PARENT_ENVS" value="true" />
-      <envs>
-        <env name="PYTHONUNBUFFERED" value="1" />
-      </envs>
-      <option name="SDK_HOME" value="C:\tools\Anaconda3\python.exe" />
-      <option name="WORKING_DIRECTORY" value="$PROJECT_DIR$/Python" />
-      <option name="IS_MODULE_SDK" value="false" />
-      <option name="ADD_CONTENT_ROOTS" value="true" />
-      <option name="ADD_SOURCE_ROOTS" value="true" />
-      <option name="SCRIPT_NAME" value="$PROJECT_DIR$/Python/beefblup.py" />
-      <option name="PARAMETERS" value="" />
-      <option name="SHOW_COMMAND_LINE" value="false" />
-      <option name="EMULATE_TERMINAL" value="false" />
-      <option name="MODULE_MODE" value="false" />
-      <option name="REDIRECT_INPUT" value="false" />
-      <option name="INPUT_FILE" value="" />
-      <method v="2" />
-    </configuration>
-  </component>
-  <component name="SvnConfiguration">
-    <configuration />
-  </component>
-  <component name="TaskManager">
-    <task active="true" id="Default" summary="Default task">
-      <changelist id="816df031-a834-466c-aafb-fce0c72415a4" name="Default Changelist" comment="" />
-      <created>1567039339175</created>
-      <option name="number" value="Default" />
-      <option name="presentableId" value="Default" />
-      <updated>1567039339175</updated>
-    </task>
-    <task id="LOCAL-00001" summary="Added all results files to gitignore">
-      <created>1567130353481</created>
-      <option name="number" value="00001" />
-      <option name="presentableId" value="LOCAL-00001" />
-      <option name="project" value="LOCAL" />
-      <updated>1567130353481</updated>
-    </task>
-    <task id="LOCAL-00002" summary="Add new gitignore to avoid uploading junk from PyCharm">
-      <created>1567131158833</created>
-      <option name="number" value="00002" />
-      <option name="presentableId" value="LOCAL-00002" />
-      <option name="project" value="LOCAL" />
-      <updated>1567131158833</updated>
-    </task>
-    <task id="LOCAL-00003" summary="Added PyCharm's config files">
-      <created>1567131230009</created>
-      <option name="number" value="00003" />
-      <option name="presentableId" value="LOCAL-00003" />
-      <option name="project" value="LOCAL" />
-      <updated>1567131230009</updated>
-    </task>
-    <task id="LOCAL-00004" summary="Create initial python script">
-      <created>1567131261939</created>
-      <option name="number" value="00004" />
-      <option name="presentableId" value="LOCAL-00004" />
-      <option name="project" value="LOCAL" />
-      <updated>1567131261939</updated>
-    </task>
-    <task id="LOCAL-00005" summary="Added Julius van der Werf's example to test agains">
-      <created>1567133774405</created>
-      <option name="number" value="00005" />
-      <option name="presentableId" value="LOCAL-00005" />
-      <option name="project" value="LOCAL" />
-      <updated>1567133774405</updated>
-    </task>
-    <option name="localTasksCounter" value="6" />
-    <servers />
-  </component>
-  <component name="Vcs.Log.Tabs.Properties">
-    <option name="TAB_STATES">
-      <map>
-        <entry key="MAIN">
-          <value>
-            <State>
-              <option name="COLUMN_ORDER" />
-            </State>
-          </value>
-        </entry>
-      </map>
-    </option>
-  </component>
-  <component name="VcsManagerConfiguration">
-    <MESSAGE value="Added all results files to gitignore" />
-    <MESSAGE value="Add new gitignore to avoid uploading junk from PyCharm" />
-    <MESSAGE value="Added PyCharm's config files" />
-    <MESSAGE value="Create initial python script" />
-    <MESSAGE value="Added Julius van der Werf's example to test agains" />
-    <option name="LAST_COMMIT_MESSAGE" value="Added Julius van der Werf's example to test agains" />
-  </component>
-  <component name="XDebuggerManager">
-    <breakpoint-manager>
-      <breakpoints>
-        <line-breakpoint enabled="true" suspend="THREAD" type="python-line">
-          <url>file://$PROJECT_DIR$/Python/beefblup.py</url>
-          <line>44</line>
-          <option name="timeStamp" value="4" />
-        </line-breakpoint>
-        <line-breakpoint enabled="true" suspend="THREAD" type="python-line">
-          <url>file://$PROJECT_DIR$/Python/beefblup.py</url>
-          <line>59</line>
-          <option name="timeStamp" value="5" />
-        </line-breakpoint>
-        <line-breakpoint enabled="true" suspend="THREAD" type="python-line">
-          <url>file://$PROJECT_DIR$/Python/beefblup.py</url>
-          <line>71</line>
-          <option name="timeStamp" value="6" />
-        </line-breakpoint>
-      </breakpoints>
-    </breakpoint-manager>
-  </component>
-</project>
\ No newline at end of file
diff --git a/Excel/van der Werf2.xlsx b/Excel/van der Werf2.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..e7884259d4fbad99a212e383d2c813d83960595d
GIT binary patch
literal 17061
zcmeIZbyQqSvp0-GNN^7l+})kv!7aE9?t=$+g1b8e2*KTg2Mb9E?hNh;?hfxja&FGU
zIp6oL`{%pv9oAxcYFGVsbyaos-rd!TGSD!X5O5F(5D*X~5OXT^-vS{aAX;D{ATS^h
zptMD8ZJbPOoB*%g>`WYW8C`+a#JMm~)bAjm!2JLF`5$b7(u6+wE@o8mvsc%ms|?b!
z)h`h^c5*+V(<%vmnu`Bi45*%)Twrs4AwtbAvO_7qoa(U~v5hBJuZ3l`FRY~!Z@@te
zi)$0e9-iho0O*rce)}at|J(DzugFlT5w^AxN0|YN`t?sZER38K65!_jO>?S|pY5_w
zSYaQGG8Gzme#L&d2Oj|ctt?dE5HB7MFFGs7fDS7>nc;Sx@8BI2&4;T`-NjEF<)sEl
zdGoEvyp_H|=X8E2z|i{9NjO7X8(%rx&*VGyb*WyKZ;nT5$5sV?g0Fs52C^sO;`}qa
zR)MxS*9230sjk|FI4^CD92pNnpZfD=HP{Wm&><+Xvf>uCG%3f@lV?Kk@M+0QiIMLw
ztL>trL(E{zz4dHJ`i>z)XK4^7Z9K4d$Qr%}oQDwEj*cT;kL@0i@Po*o9`8H2Eue@9
z=nP<H*SoS;%wEUu@(9i+r#3Jy)*kHpK=<0mRm%j6Eatr(efguwn2J=sIGiqs{V8k+
zGG~8{%QX@N#N8b<gyMgU-qi~V6*%zdodplhr{K{Ga4@lUWMaJk|G%U5fACQKTh~iR
z<}AC2G2!>GI_8B3yo_cdStR1;BZT)pe<50Yt$;3d95Tb<b@f&Mg<tuz37M1gvx43e
z)9O6eJXEA+#3b=2{uzcq9rX^YquP(~Xc)YaWPHVouh2M8PR1_-C1l-Hwj-lS3qEQH
zXUF-(uoh?0fA$w>Y#DK=OgUtB8p(TJN!Cc)$A{-_`cj||Q(ISf>L)&G|L{iqT|5H=
z#aN0WEpqKXNuLPo1!Hc<kW%r>8ux+GgS^}?mN5mwC4K(N?G1dvru@PeT6D<uH~lr&
z&B34vYIozhhSiZ!qAzYeTo%vQ{JJWQhz45wcPJRxDIC&{FPwAEvPpTa;!pU%miXte
zaZ0%xvw((x5CG4CC}0j(YbF<42P;EcTdVsCuT&Lio5ziM>pgM<-|d39PB>-^OTI}9
zR96M<c!xSUe8Md8<pHr572oXf1;x~Z#%U&7Y3r(b`MpioQ<tts+kmx_ov<czB9dlF
z<+0`Z6g#^(my@7bgg0`pW)nXe5Nd8Jh;3#k`P|yK8C@InO$ER~WNK{IC3sn<<rYeX
z_uy?w!6$_sV9RJ6dZ-*{EGmz^{Zh{Y>^tno9+l5VHS|V=E1}hj!J-RFa&mxLdf$P*
z8J>1(C}9<@vZ(d_MbJF81_-)4iSxB!<h!N`?1LV!Uha9$#Zw};nXInCU?1ZS{Xs6j
zA9Jb<4BQx{ab^%e`S%MK(exNOTqINN(H$1A%D;XkQ{~~WUcQCbypW{7gMFsi9vSpf
zD^JZ3VIF?Cg^WuZkfmNr;u7z~F?$BED`6eq)CvgHNq!yYkON<7%y(!Sfa1kF1Xo(W
zGnHk-812xutR+Ru&whygmcL8X<lLJ*6UKl%lDwrW&bT!GT~`30y`*uCN8gMCI?!rz
zHQUbtDiHP(3SiJJ9?`Ij*r$c>S~FDRdZg-mm7Z&kFpt2w?{CfhBX!F9Nl|_OknO~k
zKYO~mR?ptk6`?{b99oU_*p1_BM}A+%Xfts=sq&7ao}X;SYj8p$Q7@mY<-BDgezLSM
zsEG|5YFi+Y@*#HVL#QpeUnU;)p+C8ufLZ3oLdVnHUf8aHQWQ7!l&5ni-UdzGSJVEb
zuJ6|RZU!fo1{IGKx~#9t`Be#bqvA4Ml|r;EuWW1SrjZ$6wv|{Z9I+mDG|ZV_Nyumh
z*g^5}@7{=SdLg|fGZqE}@+!7qUpnsQ5j&Cf?#E%`*(`HQSYAY-!Gz*-7wxzg@g1@N
z(4O<rpyHdA;g9>|HPQJem6>8`yRKI?1ppGvzt3H3`K3dRm47%OwyP$4<G;OTM$(J-
zR>F_Sp^Utebh2_H%Iu?oF@Db0l*D@Zb{$Ts`qr4rv=YGo{BAQiZ+_(hBTr)@rB6F5
z9R8eHH=ec!%88u^#uE@P*UDT+@NRRRno9>$QzGrMe5uT#_vsL=wntt9W%ak4p8<2!
z?7a<7lAvD~(~j5^maqtI6LFPTti%*>QQ!b7UeBV~-hn#s_uELQhz9oL`BS_aTUrU~
z4rJiJvwX57Ffkag(7gTEQ5grHunhyA%!I%r^Em_pBzRQ*n$P||HvgL2Ai+<hVBY`P
zpVGKK%Wmdp;%9-k0mE(ftX_&i_NK$iDZq4ZU%LcuSwnbjuhl&(W!|r!Vr;5eUUi-k
zg*JnHT*}VLUpzO%GL}O0FLgFCvVz^~>FXy8`D89*ok;-;4Lh|lsotz^8M{)TpC%Ym
zXcz$1HEWyyY(_Dd#9m;UyC`EX#^{IDbe^E0oJ0Hvl#xLWVxdvYU}Zm{ga+cUoali&
ztL9=$b^*Sd;n=mO*sS|QpJB!P+9=ibwW(s!R!y$@(@qHnMdGhCz8&-@;Z}3twPgsw
ziFRZ37GbS&Yrbs+U`QC1WIy9{_UEZfVsvP*q<-d3yD-)JeXxeb_O9bg&D4g_zs68o
zI8;JC3<LxmDg*>BnB#Y5cQiLKadKpO_+q(F@fpchw#(d@LG2Z_cO=!9If;lQIOs1w
zQR_gB4;1G)o2eaHBvK-a$r(;a<&0b2T@m%#Lm^ix7U3{6?eTaIR1u=Iw@_78_ppb*
z;uX<q4#3Q~Bo7%GNASAq0-nUgt?4~?Wsy2S@LGK4g?D$^XD7L7RsWe={S!&Qgn8yh
zx?i#B?Q<%<2}=BoPAZa|(Bv-dCgmEW9z8`*no!l&(W<fTaOG}7xMOq_^q|;ivu+<5
zyEI3PFe9em>zX&UU*Bpk0}PX0k57Q3hf>v$;mOodhtZvwG@(MS8u@SCnOMGG>9Pp`
z+zR&dx++LIF;^mANv*2H#WLFl)O>Tak#Da!RUNFzNIPo(BEhB_)__K*56t&rq$tyN
z`%#yAAjMX~YV>h)kpWR+Guy3wL*poG?kCTBsR^%5!T9KjRA1_VZf|N@O&iYzPGqn!
z43=q)GBvW~mtnX*BIVw$Tpw>S*W{g6HHW=5yxyrnmhPdL&Y7dY7<Ecd-yCFa(?DMy
zcW)x#0cP(3f`BjYRdWa=6N#^u_>>PrZlk!UPoA`|=pMsB0r@Q(;vt9FPAw3`DciM6
zdhi2dn!&zSYh%A^niCRwD^GZzr7IN$=~-wU680Dsi<ob_G<t_f?Yp~{v*o*}=B`d1
zNnzh<TJf;Eq&0FhzwpvV(b%i+eXR#EahXs8Js%nlrVz*BDpNFRvB$^&lDAuLE4!o2
z#$gQw+3Zrf5UTbf6yfqf_GX1PkWm0yN*&1263ovzhBL1>gTpl0!ki0J9AV!T<sKtr
zF{;1ZQVYAyG3>xnd?!6D@S47{b6zDSQ3Un_mIgu%%Bi^ttjGsKl23rZapNGJTr5wd
zR<lW^a9Y_wSSxI`B)k`!ICEWa>-o@DCp%M{(r`AR6b1Wk+}?H@3JTE*4L+^C1~XYa
zLfN)dmp9QB(b4T*cg}p{d$t+)UUoGt!zTl+-g~@&><t%}GlV&$aa5?vK+>P>sEvO1
zEM0Yy_P*zxx|;#8d}6tVXtx=B3aEwbcvd%U-3HTQCWIR(NNy}X9}>@H<~Ww(My~_K
zd$Lis510tKeOl>Y?fPJyg4&gixXPEA9OBO%*JJiiaW^S3uDH!LCUF8xwVu@04pHeS
z`*S*pYUQdJDEX~J->`+7>A5GBRZL>gB!y}3%9;(ObGRST9+QG<HuKdqIB7CB?WYPE
zB?fEJ-7Rpk?7mtYcjk%=!Vly!3BbC$$uKs|ux&PTT<)rFir1^>*edf5pV$rPL}Xb|
z9FOesfwkO>)_A4Mwh~v<L!Fn<v-w3k=8eNv!%tF;JN36!22#Q6BE~+dG-)b=Rfvns
zZ67gZ=@LbxojWNCWbC^|7-5sQcW&~{8qYu@IqNf29Uu<ZD~v_)a&<pj`Tg$a`WvV7
zH7nDS1RoMVO|H?`FFi>{*jJPE?nc8NKkb4)Wllym@xf<ZKZM-+ZaP|Ge^x$YEoNgT
zEF(#%qhU<x5PZ}J=cB*iN_X8Xsjiep!-5>%2bmCdcm+MnxCldz?`)#a6@e4ih4r#U
zyQ}|smnd<wJT(bPqboE$^v|6|!HOZsX5ZmpUw~sV{`L94=eiCY)HZo|2#C+b_e(bq
zxz5Sl1Zcwa@c$vrZEN<2<8xuQ<6H_ob9A}pSc#^XUmUcKpCdQUOdzPM-&Rs%V~uOU
z$AF^XIZ>v5pPwvZE08eq5*n%XIQj{B?T;73RQ*$AWozW=*3w$eNa2#Vtpx>VS>EnP
zP6s2cDejUfPdnq10Bwqg0vQogDdeL`K33e+?4Am*5+cw{p4ePGrN5Q|WqUiI0ly;j
zPC;LYS?46Zkm3Q+V)q~f#eV!`%I&v68PWTylPnv%mw+%}J_VaDU%3pIEbd&--%!!R
zi<N|I8!1uR$D^o4VTmPH1Sx`gLA|(&BByGlt6C%^a8r4gkYwi2Ad7{MXxU@rj2a!?
z+x&v!W&@{z342(LQfv)QQ`F|&9c!bkzNf~p<Kd4E%MoIP04<Mv0Cb^F@5$DUR(t(O
z0i-#*f8vMwDc-F_n-~Km#mn)WX7xSQKJ57+fNo7AVR3Ze<h$2Cl|e6_H7$xitq`UW
z)aeS&G%YaOntX=|(C(U)+0j#MIJQz0lXCbN#28p->C!lZW!$+Lg)#puK)5xIW>6^x
zsbpJOMOK59YfUR+msJqL7fP92YWZu=Gle09>W|^nJd52Ricth~+aPpG{g1C{kfc@5
zD+K*+KahphGqb)=)YT8wz*jzXHSGQ>TYpXWZk&slDNDfPcK^pwaJzud@#aMrpt`Y+
zR>lUNLPceFJQ{F!bAvMvXusY6dBF8!jOw7XyZ!p2FS-5hfcK<wc@~SQ&Fk>+M?B7g
z*U5%LG~o=0fG}mBB(!S3e935!GTguq`RV{tD6^Ycq&w~<+NN{iZr>LuJ#FHjoo_04
zIza>)pJ7}0+pVUXUK})%N2OLsRKd;~ML6TS@Gcgo_3iK5Wy#LPA2o(HhYrPB?J`Q2
zLY!@|y;EOSo{p?ESw{*(5B<R1F{095B^p6k(kvRXHYvcmEL$0Ryo*{95Hd<3>Yz`d
z%vdXG7}D2iLi~2-uzwy{(N$<QF5qcQW&qZM**m0Y?Q@?wqtdoAj~ELFVOzNIDyxG}
zha_b=6S<gs7|qe85&c`4b>85U=fdrX&sLDHh;+L-<76B*M;gO%*<uBpry>>k9cN_Y
zzI1z2B9n;6x&%n!ct%h2i3x}jC|0w~Haw~CyfG#H2)k;TmOvG9X)-`q)e@1~`U1zv
zzcSophQEHQu0NEGkN2o3ai$qp#VO)5u1ABQqsEW{Ti-_GusW%JBH^}lDczR=)sqZ;
zNRzl?s|gy1WCS{svx5X+rJ1e0zbiI%DB8mNi|1l<(z@Jh8L!<lqmbnmIiS{FtcH)$
zOu}w5$2GT0ZABtznkBnq?|SOm(Rtg}^y&lV1$NqC@H{Q?h))|%`%#sw{P_)n5D=f?
z%=w34-Rg*CmS|b;CZtcLbD}_FN~pDB#IHv742D)yUAufEa)ES&gN@ccxJuwnr#9IT
z^Q(Ds_JZEZ-8^gdlU-b?<Aj^4qu%)W-Bq#asGeYj$$(Yrk=Eyv(-UFK)D};_Z<Kq7
zn^_7iRcmY{3A=r;J<TF6UN9}g5%A1~Af+a%MoY*6((^Qq6#Cf&GKLj<N+BxNndqTs
zTB0U%pEAN^PhPvSQNA;JAN+<6D|@KyO8&?-I&xIrtg*Byt>`Ii+=VNRPNZnE^oLC?
zfxY>%u`Tb-*DjjLaMaVIxCXMO_)2eZl>o>hn0jV!^`=ubhMhEq6@Mn@Xi3<wAu<E8
zuw2(p7sXrGf`w;<gG>n$HjJ<eYExSx$`J=rCEVlLic%pJ4Vl=8qn%Wyo9x{z-&hp4
zioa(di*TxmWFBGm$62E1>rcAU{M?ipJjSZfXxNDF!1gl|1<gvFl<K1&C8f$xH=oT#
zoh3p8Eg8Gj1;!?B*Q6_4BiWcXPOi#FXFRDan{yX~Bs`qzV6g;Vop||ga4fMY9tqtH
z00y5#B=%u~Tp*pgFRGx{O4>`h7wPt>PJCZnSv$xqg643QvgBX*Evb~--4^#O5SEfp
z6Lhw|h>UK9=H1{zd#|CvL<b)g>Kh*5@Ux%AFTD~gcJ0DtS^t~)WJ}lxPeK$KZ@8j!
zZi;bMCr4=c_hALP7q|-3X&pIRjQlCZ(5I~mb6hXRzd6ELtebu4q5V1L6#pZbZxyNJ
zP;{)CHs9t(p24Bqim?$HecVmDWogVS9J+&<A^SV4oMwxBS6EGr5~!Aq9^sk780R_g
zjRd#AK{=Jsm(1^4Z(h=GNoSs4ew%F>H-&fc<9E07;Vw|#R1~pICZoO_!T*xLsUJQY
zuNjH3(85E-aOX{94L>J6#gOpapdIyffN~C6eMz7*u3FzsrOwYcv?z=Eu}V6~qg$7d
zItD}H&K4EDlZ?A7YYhE?XyM_qfbTeUOP9-SI@6x^duYU&S9RasX#$-s15N7$s9{>=
z^M<x3cmwsywD*wQ=PHK@`1f1MkdSA_Z4|e-5bGQov=J;ROuX|n3*i%``3f8GZ{FPy
z(;4Zn(RdlRL7Z7O*WZC(H2=AX<T=x#vjToqjRC)EfEQ=)Uo{+^+^kI;@82)RG(a{W
zZd5OT{+;k*b8a}rBquqRJpm+aQxDWvD6f2qZ<OD{g@=n~du}~~6O6wCN1Nt>i~Zp}
zByQ{91KEt@CXvhrFqAkFNyR%NqX|e6HI=KRh5^xye4_GYRn#o;{_i(mS|0OW-88(P
zfR`Yue`^j;ES8yDHqjioG2WA8r?c>d5RF+ys;4R_XU6qJ9!+P$6vIq!^GYS!h*mMY
za<PYpn4kpZDpUJaTUS~+P2Qz`RZLfSPkLkwVg;Te&>sOu9z#}jlf78oFdy(^saCN%
zDX|SwM?npX;mf4c#{6tok1x!ZWD?k_jbN&JuAd_Bzo@K+rtRm@7rh^<Q~A<5M1$y8
zXy36M5c5LZ<^yr(Trqu-5RXENANfnYB+NLeHFE?)*8ypm_4m5=!FalKmPx89D$eFT
zIcQVnE?eBN3-~!^uC><H#z3``rtK_1PG51UFqO&nre6^YHjE}zQEn(Y&$Eo(tAk2t
z6<FgpGSlT(9X}&GgNY*s7ZhbAWV#F;%?L|5w=mZl_n2;`xR|$W{KMkd-^&iOS8tMa
ze$L*Axr!>k4s=@E-)N61SHJCUzp>DrNQ|MfDv~aAztLn#9SN&`8!Ihpj5%lBUy;Oz
zK~==@I*6WzX8x65s??EpFCJe`h)Am}8WqKy<In_=KU3FyF%LinrKCYoA{)9cDdYL!
zi&38NpaZLo3A&s)$qQv(v-e#xZtHlox?v)nIB%=q;{*0ej3&}ehh9fb9c65BDSE!U
zN*|f?5gGp>2zq7<PkGd9&Txu*%BoXP6on|S0Xxi78wDb40abN13ji<@K_JP=;Vd=|
zi1UH5Ck}8b!t*MYEOICu8gmGz&+_sasJqUSh9ei>&&5y-iGXry62rWqh(|r1rcl{&
zTomr_#<FA2erzbpLirxT;B@X+W^tAYeFtydF+ACSN{&%J-(V*aB>T~Q;j{WT&a^v!
z5l3X|w|mFzo|>&Q7y@Is)aV5ihrQh$KZMG<Cfo|Q^G@wY_@jWISSkWBTiMo<A~DO(
z5-e!jhPfeXtb^2ynH*;LXRYITE&ZZi$06>edh^SFFh9TTnlD)NbS}Alub`|Qf{C<X
z+po)*l9tZ#db)vQR&;T4nl3H3eJ{CUz2u7D$?Vo)T>=BW`N*}=VP(^%QESan=y>lI
zPvFH}&&qI{Gl`RMZ6O{%i>jHm<gGzWzB#-350DGw5Rqecu*w}LvL{#exgWmV!S_r#
zs;cx};g4w&eOKJklCvvMU!=C-!0b4V$3)|#Idm3o8jVYu(+3iF?DI*udXcv*0iH~o
zfyeqkv%yw%lE(eA7Aklr2J?3|FgI{8F;;PMu&^<END9`&t@j&f5Q~JNF5GVZId9P4
z34gIy4WY$Nfgq8fYottCcB~{EYZH*bi^+&Oo69c40OjE=!e>q<(I!=}!M%zEa=gL1
za$Ts{#XEmjs2V^~I`u>V^sVpdl9%t4Oq-z>!V2_8+mB)D_Ifhq?*yf6wB18*W}}f+
zw>-M_;q;8BgYff1g^j2XvFN1XFn;Vk$u|8WLm4=!4AV#k9D`u0TNFmT4R&0X<KyqD
zV>LTgY910TgJS%=$db?Wnz2rugI??Gpj{?Y-79zZGDd#-4Z<nvg27_38Wjqi>jW^|
zDN}8o29sv?y<K{zE@C?7H=o6W+(VpyvzOpB@h5Wdp0m}Tg&JcA1DAVKG5u~lyup0e
z(6)WfgXs;hxa*pBn=X5;@B20#5eue})I>BPyjU;|4@Ki-s-`~pWy$T?&9)!V-L~Gs
z84%HX+UT=!o?<rG(D)vSOU;}kBJG5mv>ZJXvF^D9*W{5F6X{lpbm~|bN*-au@X>HZ
zJ3r`yshw3p(5IKvkb=R(GmVBr=2OP~j;32pi-^ssoIH-_>?IPGQw$m#)x%*c@?F?L
zHDl#GE~A!IGe<^u{@Md?!z6^+Hr+e(2}+F1<mwXeGSzFu`jqZ8yR;Vz{2=1RlOmUW
zl5W(eST3E%N|SDmE7#A`cWi2Es5<nBpI9&SEo)1Db;c#^xsk6JkEEVP%5yH6<BM7;
zRI;yuU884l#}!|FNgU`N=~h?7sov8w|8~{W(csnlu0pI>L!M>($Rmu_e9EXANWp<p
z&|o3EA#Q8tpmzMq^mQR;Z>ck;LX>6vv{E}~FnW!ei}#47%e>o9Q?nRi0_0cRawrDo
z&Vo<G7~hqQr<PFh+HTTgJQ0GyHUJWFh!Fc{*?QPhWf7cxu9NqfEhQ8-1+xC&v>Ls&
z_u%|dUofh9m+2x#X6|#2r1l!^kaN$-iO5sf=4m+O#LM+O*GWlE*KnlzRN;@)5SM7k
zXMFY*lTb}dx8BX5pD=sysenk7F<v{`h&N0h2fiu<#V&<Dl{dpT*VeS-UWTO;z)`+@
zQ(>5|meQEHdh7c-nLi-;oh||Hba!bI)d;HclQdFD7RDP5dlZ-~R_GOvkzFY$_UdQF
z3-MVmLi|0;=0`U~ENUiSxDB>K;!_64g4EFb+r}>f=xE_x0q3ti@C`bI0`&0{zc4Tx
zD=#4DiG`7p8t<n^Pud{S1;yT!7zarjXYdclAh=Y;9H0`hE4%CPW`=sX(vozXP3XFD
zENXC4UqeI~_#(z#jP$DRtRv>{A=yooDkmb%ZYr1(a<TI?H|~ds7aNqmYL{1#G37zW
zr)ieV!XK6a>IE-L_7Un<aD-7tZaGzRLK!Bkg~=-j#XwsKgZ7jz97(+b>#2w7mL}F%
z$?cGtpX80#SFQ4Ef2(#?44q7q?t@>2d?gaA?dV`l!!suWKl5W`C@uZiMTW*PLORhQ
zPcdj>Y16B9yRQ>I=&DkKLHuV8i|2Mcy*_&s)xiNt`!2ze(uqlL@OjUgwe#)Du(!0(
zH1^n;Ti|Q+J6LVK{(%8(rWc*KwJ<`UcQ&(8y$c8fJ^gDv)7y{^_K^}Wz@g--K~?vT
zXU*4!q37HSyU8VYfJJ2U&eB(#PGLXOlsi<aSTo2{O9qf{I#FInc!1_2725osyxxF+
z+3l-mZ<;lNWu{i5f9cv|6<*3|!b9<K#8uC{>wA=%L^ABmOfw@%P&&L<#jBJAzRr|{
zYU8Ma)#fh6*%o*d4wlz=)gM!r26^RhoZ3EqAjovBKMH1`+S^rTy#lQ?Ua?qx^Tn3w
z>xUz92zqknswQZfgEz!27l^_dzu-sOVCwpkIlnfKf=WA0X9_xq-Oiv;|Fh<Zrf)*K
z=v~e8Q;TrNNYHIfk)Z_z=ccWQBr&}zFlH%L3{Fw`c~@ed(0<0IK+cor7`**Y^SI2J
zs{}U+$CvU&N|fK{R8<QK5z_|?c``Bv3ZvNf^;qtRSSc0QLt+P_eGwr@Fez`**Znqj
zm6T6I@RTTM1J#L2!RN5&rtoaGtutn{dl>2>ObjbZWMul%YzRY-Iatr@E@N8g>-(4d
zA0XzIa*l@k_%crKFd}Z@{w=f0^xwtQgV%#b93UXD!84Jut&yUGt(_xyS?CSW(8j>R
zn$ZPlI-3cy#W!>RY(^c%tC3`(nS@oQRzlFg@UwpUl#)xlMy5axB^;u;)mDqVnURyB
zvR|%3IyQa#H239MRD?ph5eAI#E54{q;#bweF)>Ly-0-U}h?&C2)LYGQ?6Md5GDv8Y
zjy_-b6BkDGdWw1N*|w#UTnLhzec?_E&n;Rx6C3W>6e9Hw8+px-FSq3Y@Uyae8o0Bm
zutL*_gnzJe{y9>rz09in4m&C5=LS+q^wU6)_x|AB@uq0Bn2#sNM4b2NRWtq_y7>4H
zz%66@-DPt24U&ui?O^qKb^hl-d@h@*Lqz)yMjr(uT<wT6$(p9?7x-3BYqp|4+q35Z
zeDb!6<Ro4YzD6<YcR;_i3mq%NWD-9MMcdQpWAgFj?K7=O?RT)Oyw%}qa=y;8&dH6`
zU@VnrNacI|kvutS{Bu>d@wQX#yM<fh-HFZ1#mgEkCfZSx=qN_HSsON`>g4#l-Uh|$
z*Y&d1W%2m~rWwN=>G<D$IJL$$`92i#l9tUx7#>Mhvpli*?A=Z@Xwr>AC2tI}<U&K-
zQfFTBS}sy4a0r1#dncE=FL>}hm5ryH2JIb|oNM;^YP`ZaYPf{Pcp0SdQK*8z$!laZ
z#M-zgdl=bn>!}2WQZVLhlr0LLyNw$TA#9{C(>cI!M}tDyq8V6_GZ|n$fa_8q!Ey$|
z(E2^zayZjE>V~*){ovcPP^Q*nWn}wogduEM@TJ2(rVUgaG2i+H@TDDRzTl4l81ry2
zIKRuT4zj%xj5$1*&Gq*Bf?ydE?Oh@m&UL9XkibwAeB~Z{#;Rfs6(tM;_!SIL(SGZW
zY!6m1_#XQ+zU?Y1N+<+y8hpvN?ZOn<9xN4nFLesC-9QD$9<T>|X|Lr-9@$<T#$57V
zYPWfB9W^=xg5e_hI_Y$ywQ0*^uQgdz*0kaI9O*hkK&pH9*@BNyFaqpM+LOJCHw~bE
zlJZc99el9Tv-9%<B;KMh=3C$&Ee>na$o9H0<_qj#^Ng3MJBBP^Le2z%%~PS<oFqDn
z3OUnof5&0T7i*~P6i~Ht81q!#hhM39y5-Q`%~^Kt$;louC-cMOw(>tuhHGSF+M2hn
zFSC@=mVSo-KpLo5&BGF-4}lwRK69I{ryWT`r|!HGV*;r~s$o`9KYwi-rB@$VnGX}T
zFF4!>W=<5N4^2{Z+vQ20Lrgi^tU(5|iq9xRiM#rMIm?~{)w(r@-?@YCnPck{UZTxY
zxlh+y*1Dl>hzI=ae@uo&!x8%zvN#&fzj0TGn1HbzkE?;I?-?i*!DfHBrBpY?zmYY&
z|CL+38|x4E8?(PaRJ&FGkpE!Ib#we1S-1PmzmY||k^YUW(k=IgY@+>kAw7lOVkYT)
z+!z#esHE})a|wTU#=Qc>8KH(+MGQ9j71$Rf%J&z%a+87h5_Dl2m`KH3l1v~Gb6P*r
z76q^gy38riWO6CSDy>1NMMYo--4e7~gP3PZLuzt`PwMl)a-<xIza3!kHwPfarh$FG
zhw;GeAMzB;y}oGz3VJ2?yW2oqpi;N^Z;eIYfA}2|VC}H0=zoO-Fl6jC*yq{za>Get
zI_ST&;}o;}qh22=V>sM+`eQ|~(A0ms{67>a4k`IJMgG|T?~2qjE;pGgw(@{O>+fa)
zq1ubG7tY_?;g?D7V-_h781|9okJ{zp=dGatuX-TOU_ID8ZENLa$v|xhx)2RiBun#p
zVCb0I@17+Mhmm572Ru&$lGOfXNR!&Q5Q<@ROO3x6@53I5_bcp!OhHMnfPYyL52)Ae
zHvH}}QMeoSpTq}yKjJK-H`HV=%$csIP}cq(3EC3c==M_oBx2kD3n!Yr_&+(Z>{b6s
z#ITq9ClSY9_aDSx27aUx2-9BWFV05nr87$Jlj-_Aevm4Z_<2-wy<YteZW*zm@B<5(
zvik3?@cjyEu*Cz1z?=`E6g4bOnj|JuXfO?aExc#Ex%aS9ttOh;gB$or-~8nVqhA3H
z_RW6~AM(j#;zKrhOnk^CV2ek+(*5nQ9bDi;CV5PJ$Rm%54_V|f@gaviCO$;)W8y>p
zcuaiA9*>C+x#JgcO$(ekz$0#zzT<BBnE&S5BSS|bUC4D)^6q?S9k3vD?lmdinW-YO
zOfNjJCo|CJalRA`+(G1ZjAlb=XY<j~eY!n&DbV<V^W!?-G>5{o@{6wRw3+~`+gFq1
zePM9Pd1mWKgRsW3?xxGbu4j#ICA`VZbJ%=?J>hwLh66n^x3}U^!h&vx?~}5}2F*OU
z1lo7iZg$jw0U+6Pr|oV1?Fe1H=#9knI5tX!Y{5*K2yvQ_ko~Wkam~@0gy@pdekts|
z{?sbDB%TgQW<&=+3bF?>R}!ejLFI}$3T<;J8|T|KgX`$wxw4>(q-};xBKzZzKq)H$
z25Bj=w+d17U$wV}SDsRW)86-~)J(-34@ZLjqWBQD-0w4u`58=zb7b{eWwSfnbl);U
z5VGfy&)Ni{h>(~qq8xL-hX?I<If}$L?zRO+fA}WV?WS}Z3MR35pfTySq%==liD{xD
zEuA%j?2lWoq#SinP1?;7+85C!*%rx{y_K-VMbB>-56T+mzv*s?D<Ovj?T0hkA7`FQ
zIr0-rYqGx;0<&U($PG0U0>J&7rifZEk?Ib}r3c>xUUxZPCfT83C5}2+O!S;0+4bUZ
zO=p=zDG3ZJ%4pAyB&W4yqJ#`qv5)bf3q|OA=8A!m8kn_<@n9udXn>+T#6gVxEx<lB
z+xyaj@gRiy*!Jzh)^?s#xL2yA`MS$(+PbcaZF%~BPUZ_Oxz<LT&R}P6bB>esrYA3E
zE^@|e)wd`;%{PM8{{NW9(*GMwpKDTnWB<3r{<{k3w%c_2c=7ruWVRK|&vGq^@(I;u
zxjQbrMcKI?a{t-yg&|mbB)V0^btZQ{G=Ng}{~Xg6Sn&#ljf&0wC}{%8UfG9UpyP;%
zJv4F$WB5%Ue2LZC<?#KX*R%(<Aqe+gRXxUI`%7RDJkIWYFP`oD7u)6PXQN+EANElF
zi*FAD?ZJl9EyT}%;d)K{V0kcNNq0CY=|hl&qJX3O9*u$#9Fh+WLF!-ZiLO%YzqKG4
zM!mOvNp~MO-tS{86y+~&Jv5|Wcku*A?Y*aH-5YLG?n80Z2=^G~20Uz0B^Mg<@OH=9
zLRZa$xZS7AFUqsZ-cKfs=j3gAOxHwGdjPUm$buvh$t4x+=61>OoQd*3@!@+uKh&En
z*b9sDmS&zzWKUjgit}zRaZ1iwqq@l9-eic-U;X~B%i$}P;zILY+onPOW`OncX8ND0
z-L|P4E?zb@<I~wxx4rfxlNEpkkCy`g8~kf)?Bkqqu9cqarg7wkmeb&!z>nX1O`?f1
z9i5DABX*0VvQaaSSq0rZ5}uV63ibkjl*eT$P6OAXigKfP(vpz3G(iSOGMu@JdPXm3
zVe^H4ED2p_4M`+3%Vy5PQQ424q_(sL#!i1IC7-jCAt+lqD&3<c9C{;keG83}u>B=d
zs$s-u^dr-r+U4C3?Brbm#tZtir833|I-xx6WRI)&gN7N`A5*rM*w-i<;#!6JGOzCT
z<G<O4+6)rG=S@HdfTonz3OlIUx0+V%V?YZ7y^ZbBL^XM})mNe%t%*eAH~KG!-az{k
zIJTD8;a@pjcARt>bz7!kaRejQZMsx@tnM0hyWZ`bb}Wpp^!?NX3S91FXw93=YF^cp
zo|?JK%I)zqC-V`quwHmH&GuH$^yb=bjA4^WL`@22jSCXZC93dYu{gAxZfNcd=J7dN
zK1CzzwepD_za?!VE8vX6O&z0ra@BOYx7sS$_Li!ef`4gS3M`!U#ZrOWBQAEOj;dxR
zejv+u)pe(+^faM{3Xez6gSSLr=NH>o%fVl4?cXeY*1$CHSOGUq?M_WI1y%=7=Z)$a
z$6vy9nmcsAS{J3KEe|kPmTFhUzW_-}S{EMx^}bs&@q+IQuvVBiyI59SvOZSlkz5=Q
zGpiwU8kKBxY@&2C7R0$Rkq6Qkt_}#KPPC!lh|BLd9e*pIIhuADl=R6`(T>1vbX54z
z=f%+^RNzeO{ziQ#(7F*PSE#RuK`dgnk!o~JZ>*t<+y``1M#k^+F=c#yHg79`yWZhz
z<Ad5D&7FG3?E5~im7Z4rY9*GQhW%=dm!1~?YQbynzz=xcTWPX-#_GNWSQY$w_LKaq
zMCSe%-<h#i^Uf?)gUh(yK$6GwJ^DtLY_dmBqk`VFst1SmuIh&1*DskKth&2kq{9GR
zFw&Zy<8J)Dn(=IGw~ltE5{EZ(VMXZ+9{W_icc%rtPuekj(B+dtI7~+xnDZ;1?#k|4
zyx~nZNFyTOhP2K9n9mcx>C)zaw#7ra;dT8g8oh<cVNZy_XSE~SdLE#ZyU}Q&SrIP_
zYCj$yICk6jZe>b%eOvNbY}h6G?s)ixts0T%G09KU<LJJdOOWj8z<y;tf8mX-2c@Rl
zJ{m(@)tW$3V(C&pF!-21dDWxpbM`d()mDwJ*@$P`O=3$iFTd4?^4^T&(p1rU)052v
zf%wL&vy;6nk6Yh8L+NfejvtL2RB&%HwByGz&vo8xfSZ2cCVoY;$Ht-H;HeCNPCjdC
z=Y|0=8#Ai?W7ABHqT(4~FpK1{Vf8e~2#zm%A%F7dnlW{0m+{M%TGx?Ri~Y&=0Z(=>
z3sVNo?tyCsHBT>9<5krBD5clJOu1PlExu30GD6|%K?`leB)yNdm4FUD?fEHdRO!;R
zttw-?^T!|yhP{4@x)#yinZ5a$uMVT_!8^mojU0p<ne;B&HWh306(=UE=k>)VAWJ?o
zZHhX+22}=-8VFp0R5aIEX6fF;n|U-3{?VGM{bgYCDdKRZJ%_fHj};04-ojM7E)`U^
z8-d5urV7#of{IL5U4dLBJG^xfvX*D=H#8SpezCnF@2w6BE3e)Z9@yfp@}U&HxD)Kj
zb8ojg4`B=Cyqr2Qnaf)uZ%Cq+fr}rx_B!?7IX@+sKpr%k*7f+>By;&H&RSvO+Bp%N
zSr0bZJSrfu@nU53&U(qU>Qqg>#(TTAlZ4%ZPbvI!&k3udn7<Xg!bb<*pC<dY=X75x
z{kz!pSMBtllGlgI>9^7Sh@H%+Z}<H7P)}_}+)-mc=!xyqo<JrUUj+css}BlnZ5jRj
z7lCHUH>MwEgbI{qX$}`NgOss}zu?)zPg;Cr?;Xd6p>o2O+<Vm}`|QfSv2cLh@SDOc
zil6ug?unrSOlRq;9~L4>f-!x`6TMGVlB!f7`?2_xw{oWn&%gR!hSSKaiSkdQi?Ic7
zUuW-toSB{W*EwTlgy<6@&%gAuq8mv{1tjbdwg)S(e|h&)+Do+*0PfpAd*7%1tg#(@
zB?#Pm@Fvpz{=1QZwUMZ~frZWe7Q9Xj2;9|YL3^+ZIjeengQQz6s6U@>LMJfV1Q1DR
zi!IF4f?{T24%f)YecZe?yRvW%H!^+4t098qubLbqeXrygqdplbm#cl|!`%Ojuwf}k
zSI2reQan)A(jg8Bl)ec)D%7r}W%P<>R3H-6NOZt|PITp49K~Q*9aIr+I0)2(P%Y9#
zC7%4$<O$7uROb0<5-s|r4V*12b&R&9BC*TZyO+B~!hZaU0Z^%3S6=M8m+<Ea`X~Qd
zY!KO_fvE|;;tk$jL<85x|JS0fk*$NteQ&|Xz5cf(9vlUpaY`~>%$R`-uy>+gJ<@C6
zI0n^=QU&-ycZ4l#zuKKJOiv(h+F3}i*DjOylpSaU3<!{??rnX$f&ofXs9?guxDhv^
z=g)rF<N;i=t>X2)*H4U#$%V#cA7o@gpI0nr=eHd9t*5{pjD4TeBBf*i6lLLT@l8wp
z#Pc45OM3EqQ{Mb?(1A-mECCv84~GV2ipfXV%l?;QlJ6Ixt1pIB;ZX@hVL4t|+2OuD
z2U!VEt`GXj<MC{b5It`gW!X((aAiW?ZKZ)<v1xefY>(5mhH|iVxm2V*&SV0@f={6E
zS1I??%&%!tEQj0U+{X@NjD)p#>0YAhFZrQ<C=9RJdCpcE)g^dzb}z><&!ZGn=7aWj
zV`Rj;A^s4ZY4?cY8!7sLr^8CTlsu`s_~z<mS^t4io7@6mte5=`ooqRND7uS3j3(5A
zpAb#Km%hVHi1G>?asM2A22xK@L&5Qg1}=a>{d+tb*xCJWID+Hw_a`Gx*cQb6EN~C{
z^2M_`e4aVQXMUxIuhfq#4Ity{&E!59<y4c0l*VpPM%NFIO}qb$T=j6h6qA}G-~}AV
z4EHF8uH|CF!r)usoF=+X&FOk$!NxtQJ79T95#_rq0q|5#C%`av@P(CxDHW4RP*ctY
z=CFO#=cyd&p2b_{+%vI$ZnnfZlZ*`YwWPWw$&MI-V=VaxCiRDIRfQKwWqE|>`}Quf
z&q~J|5?X44&XX<+YIqE-3(vhP7ZWqW!+zB31FG}#7!xxZEUhAgIpDWGyTVW`!2Klp
zURfJyk5~88Vmu%}<*Ad8D^NZ69Eb{+A$?5AhY?Z*2O5j0Q+TB*)N>ud0Fe`Cof2fk
zb2eagCe9=YjkEEYL?6mviNJp23-qleC*>$#cOTEe7v=+bZjhfXuj2rhUaG<+l?P$7
zac69x%1KwCndhr!S1Sa)NjJgQpb6u`3t0IrL8LJOTh=QQ<GHlT;LGu&c4q$DmE*3H
zLw6-_#tU6DHj0JY>f5J(g)bx&1GqBd-_O1O<*vWJ|K$|CqRf8+{O8fzzXb2U+rVSu
zZwGN73qC%4_(${vyr1{`z{+1o5g!Zx=NYL#q7V?_h!4X5A19_B<2*hm^9O10*?%eV
zk5e;`Q65*I|3PU+dq8<ykNz0oaSh!c0A9QYfIk&=k3}CB-24%JMfgkfarw<-gvTW)
ze-K=W9uWSuAmuU2<L%ZzC?nu!vj3W=|6}v@F~;Kuxjz^!To0Q3-yZ8813rE>`UBX+
z{Q&qMPfCxa|MOwvk1Pa40ngvM^w*=uWAXo-8~!fN#s4?)-!p}x3@kY8As~>!f2rV+
Kb1HcM>;C{s!H#VJ

literal 0
HcmV?d00001

diff --git a/Julia/beefblup.jl b/Julia/beefblup.jl
index 547cfb4..16c0c6e 100644
--- a/Julia/beefblup.jl
+++ b/Julia/beefblup.jl
@@ -7,6 +7,9 @@
 
 # Import the required packages
 using XLSX
+using LinearAlgebra
+using Dates
+using Gtk
 
 # Display stuff
 println("beefblup v 0.0.0.1")
@@ -16,19 +19,23 @@ print("\n")
 
 ### Prompt User
 # Ask for an input spreadsheet
-# print("Enter the full filename of a beefblup worksheet> ")
-# path = readline(stdin)
-path = "C:\\Users\\cclea\\source\\repos\\beefblup\\Excel\\Master BLUP Worksheet.xlsx"
+path = open_dialog_native(
+    "Select a beefblup worksheet",
+    GtkNullContainer(),
+    ("*.xlsx", GtkFileFilter("*.xlsx", name="beefblup worksheet"))
+)
 
 # Ask for an output text filename
-# print("Enter the full filename of your desired results> ")
-# savepath = readline(stdin)
-savepath = "C:\\Users\\cclea\\source\\repos\\beefblup\\results.txt"
+savepath = save_dialog_native(
+    "Save your beefblup results",
+    GtkNullContainer(),
+    (GtkFileFilter("*.txt", name="Results file"),
+    "*.txt")
+)
 
 # Ask for heritability
-# print("What is the heritability for this trait?> ")
-# h2 = parse(Float64, readline(stdin))
-h2 = 0.4
+print("What is the heritability for this trait?> ")
+h2 = parse(Float64, readline(stdin))
 
 ### Import input filename
 print("[🐮]: Importing Excel file...")
@@ -39,7 +46,7 @@ data = XLSX.readxlsx(path)[1][:]
 print("Done!\n")
 
 ### Process input file
-print("[🐮]: Processing and formatting data ...")
+print("[🐮]: Processing and formatting data...")
 
 # Extract the headers into a separate array
 headers = data[1,:]
@@ -49,9 +56,9 @@ data = data[2:end,:]
 data = sortslices(data, dims=1, lt=(x,y)->isless(x[2],y[2]))
 
 # Define fields to hold id values for animals and their parents
-ids = data[:,1]
-damids = data[:,3]
-sireids = data[:,4]
+ids = string.(data[:,1])
+damids = string.(data[:,3])
+sireids = string.(data[:,4])
 numanimals = length(ids)
 
 # Find the index values for animals and their parents
@@ -98,23 +105,183 @@ normal = Array{String}(undef,1,length(headers)-5)
 for i in 6:length(headers)
     for j in numanimals:-1:1
         if !ismissing(data[j,i])
-            normal[i-5] = data[j,i]
+            normal[i-5] = string(data[j,i])
             break
         end
     end
 end
 
-# Print the results of the "normal" definition
-println(" ")
-println("For the purposes of this analysis, a 'normal' animal will be defined by")
-println("the following traits:")
-for i in 6:length(headers)
-    print(headers[i])
-    print(": ")
-    print(normal[i-5])
-    print("\n")
+print("Done!\n")
+
+### Create the fixed-effect matrix
+print("[🐮]: Creating the fixed-effect matrix...")
+
+# Form the fixed-effect matrix
+X = zeros(Int8, numanimals, floor(Int,sum(numgroups))-length(numgroups)+1)
+X[:,1] = ones(Int8, 1, numanimals)
+
+# Create an external counter that will increment through both loops
+counter = 2
+
+# Store the traits in a string array
+adjustedtraits =
+Array{String}(undef,floor(Int,sum(numgroups))-length(numgroups))
+# Iterate through each group
+for i in 1:length(normal)
+    # Find the traits that are present in this trait
+    localdata = string.(data[:,i+5])
+    traits = unique(localdata)
+    # Remove the normal version from the analysis
+    effecttraits = traits[findall(x -> x != normal[i], traits)]
+    # Iterate inside of the group
+    for j in 1:length(effecttraits)
+            matchedindex = findall(x -> x != effecttraits[j], localdata)
+            X[matchedindex, counter] .= 1
+            # Add this trait to the string
+            adjustedtraits[counter - 1] = traits[j]
+        # Increment the big counter
+        global counter = counter + 1
+    end
 end
-println("If no animal matching this description exists, the results may appear")
-println("outlandish, but are still as correct as the accuracy suggests")
+
+print("Done!\n")
+
+### Additive relationship matrix
+print("[🐮]: Creating additive relationship matrix...")
+
+# Create an empty matrix for the additive relationship matrix
+A = zeros(numanimals, numanimals)
+
+# Create the additive relationship matrix by the FORTRAN method presented by
+# Henderson
+for i in 1:numanimals
+    if !isnothing(dam[i]) && !isnothing(sire[i])
+        for j in 1:(i-1)
+            A[j,i] = 0.5*(A[j,sire[i]] + A[j,dam[i]])
+            A[i,j] = A[j,i]
+        end
+        A[i,i] = 1 + 0.5*A[sire[i], dam[i]]
+    elseif !isnothing(dam[i]) && isnothing(sire[i])
+        for j in 1:(i-1)
+            A[j,i] = 0.5*A[j,dam[i]]
+            A[i,j] = A[j,i]
+        end
+        A[i,i] = 1
+    elseif isnothing(dam[i]) && !isnothing(sire[i])
+        for j in 1:(i-1)
+            A[j,i] = 0.5*A[j,sire[i]]
+            A[i,j] = A[j,i]
+        end
+        A[i,i] = 1
+    else
+        for j in 1:(i-1)
+            A[j,i] = 0
+            A[i,j] = 0
+        end
+        A[i,i] = 1
+    end
+end
+
+print("Done!\n")
+
+### Perform BLUP
+print("[🐮]: Solving the mixed-model equations...")
+
+# Extract the observed data
+Y = convert(Array{Float64}, data[:,5])
+
+# The random effects matrix
+Z = Matrix{Int}(I, numanimals, numanimals)
+
+# Remove items where there is no data
+nullobs = findall(isnothing, Y)
+Z[nullobs, nullobs] .= 0
+
+# Calculate heritability
+λ = (1-h2)/h2
+
+# Use the mixed-model equations
+MME = [X'*X X'*Z; Z'*X (Z'*Z)+(inv(A).*λ)]
+MMY = [X'*Y; Z'*Y]
+solutions = MME\MMY
+
+# Find the accuracies
+diaginv = diag(inv(MME))
+reliability = ones(Float64, length(diaginv)) - diaginv.*λ
+
+print("Done!\n")
+
+### Output the results
+print("[🐮]: Saving results...")
+
+# Find how many traits we found BLUE for
+numgroups = numgroups .- 1
+
+# Start printing results to output
+fileID = open(savepath, "w")
+write(fileID, "beefblup Results Report\n")
+write(fileID, "Produced using beefblup for Julia (")
+write(fileID, "https://github.com/millironx/beefblup")
+write(fileID, ")\n\n")
+write(fileID, "Input:\t")
+write(fileID, path)
+write(fileID, "\nAnalysis performed:\t")
+write(fileID, string(Dates.today()))
+write(fileID, "\nTrait examined:\t")
+write(fileID, headers[5])
+write(fileID, "\n\n")
+
+# Print base population stats
+write(fileID, "Base Population:\n")
+for i in 1:length(numgroups)
+    write(fileID, "\t")
+    write(fileID, headers[i+5])
+    write(fileID, ":\t")
+    write(fileID, normal[i])
+    write(fileID, "\n")
+end
+write(fileID, "\tMean ")
+write(fileID, headers[5])
+write(fileID, ":\t")
+write(fileID, string(solutions[1]))
+write(fileID, "\n\n")
+
+# Contemporary group adjustments
+counter = 2
+write(fileID, "Contemporary Group Effects:\n")
+for i in 1:length(numgroups)
+    write(fileID, "\t")
+    write(fileID, headers[i+5])
+    write(fileID, "\tEffect\tReliability\n")
+    for j in 1:numgroups[i]
+        write(fileID, "\t")
+        write(fileID, adjustedtraits[counter - 1])
+        write(fileID, "\t")
+        write(fileID, string(solutions[counter]))
+        write(fileID, "\t")
+        write(fileID, string(reliability[counter]))
+        write(fileID, "\n")
+
+        global counter = counter + 1
+    end
+    write(fileID, "\n")
+end
+write(fileID, "\n")
+
+# Expected breeding values
+write(fileID, "Expected Breeding Values:\n")
+write(fileID, "\tID\tEBV\tReliability\n")
+for i in 1:numanimals
+    write(fileID, "\t")
+    write(fileID, ids[i])
+    write(fileID, "\t")
+    write(fileID, string(solutions[i+counter-1]))
+    write(fileID, "\t")
+    write(fileID, string(reliability[i+counter-1]))
+    write(fileID, "\n")
+end
+
+write(fileID, "\n - END REPORT -")
+close(fileID)
 
 print("Done!\n")
diff --git a/Julia/install.jl b/Julia/install.jl
index 3276384..1a3698d 100644
--- a/Julia/install.jl
+++ b/Julia/install.jl
@@ -9,4 +9,5 @@
 using Pkg
 
 # Install requisite packages
-Pkg.add("XLSX")
\ No newline at end of file
+Pkg.add("XLSX")
+Pkg.add("Gtk")
\ No newline at end of file