From 11e965be4fd5a8e45605fbdbc4213d58d8da936c Mon Sep 17 00:00:00 2001
From: NunoSempere <nuno.sempere@protonmail.com>
Date: Sun, 16 Jul 2023 22:32:03 +0200
Subject: [PATCH] rename sampler functions, elaborate on README, etc.

---
 README.md                                    |  25 +++++++++---
 examples/04_sample_from_cdf_simple/example   | Bin 22192 -> 22144 bytes
 examples/04_sample_from_cdf_simple/example.c |   2 +-
 examples/04_sample_from_cdf_simple/makefile  |   2 +-
 examples/05_sample_from_cdf_beta/example     | Bin 14204 -> 22176 bytes
 examples/05_sample_from_cdf_beta/example.c   |   2 +-
 examples/05_sample_from_cdf_beta/makefile    |   2 +-
 squiggle.c                                   |   2 +-
 squiggle.h                                   |  40 ++++++++++---------
 9 files changed, 46 insertions(+), 29 deletions(-)

diff --git a/README.md b/README.md
index a13492e..9924e15 100644
--- a/README.md
+++ b/README.md
@@ -1,15 +1,30 @@
 # Squiggle.c
 
-
+A self-contained C99 library that provides a subset of [Squiggle](https://www.squiggle-language.com/)'s functionality in C.
 
 ## Why C?
 
 - Because it is fast
-- Because it can be made faster if need be, e.g., with a multi-threading library like OpenMP
-- Because if you can implement something in C, you can implement it anywhere else
-- Because it will last long
-- Because the language is honest
 - Because I enjoy it
+- Because C is honest
+- Because it will last long
+- Because it can fit in my head
+- Because if you can implement something in C, you can implement it anywhere else
+- Because it can be made faster if need be, e.g., with a multi-threading library like OpenMP, or by adding more algorithmic complexity
+
+## The core scheme
+
+Have some basic building blocks, like , and return samplers. Use previous samplers to . Then use the final sampler to produce an array of samples.
+
+## Getting started
+
+You can follow some example usage in the examples/ folder
+
+1. In the first example, we define a small model, and draw one sample from it
+2. In the second example, we define a small model, and return many samples
+3. In the third example, we use a gcc extension—nested functions—to rewrite the code from point 2. in a more linear way.
+4. In the fourth example, we define some simple cdfs, and we draw samples from those cdfs. We see that this approach is slower than using the built-in samplers, e.g., the normal sampler.
+5. In the fifth example, we define the cdf for the beta distribution, and we draw samples from it. 
 
 ## To do list
 
diff --git a/examples/04_sample_from_cdf_simple/example b/examples/04_sample_from_cdf_simple/example
index d280840081b9440e145852f4e89daf1d8b6c54c3..ab05dacc039e0c9e66970ef8d345233a110a471d 100755
GIT binary patch
delta 3838
zcmZ`+4Nz3q6@GVrV4niJD=d#iU;#0LHtve9XcgF1l(!KhAn4i}aVEx$&BP%MmL#<S
z$C$dKOd)5+Hfb7hl1bu>wp#2&TVqyZQ|({Yaca#^4KYo!{zUnSO8n{Tci+3u5@&j6
z-aYqx=ey_Ld+xnw@9x{nI`^{f9QVLDNiQKoNb#PsH)K^dG_~ggzZ}dolpwj;k|MQI
zm==Fqp0F5tf;KV59TlvoU`6|OvugH#J0BFC50Y$c$twTn#_pFIAAjkwkB^)`w5#FV
z4<eT~(?)BlhfkH^OVNkgq}Su<kaZG!oj$cH><q8TWeZT$mcw4uwU8}`dTjxAoKD$l
z*-jk`*jKaPQICBm1yi~ll0~vgHc6ttq~urYJ~4<@AIJqYgfG^!17oI2AEC{$pX$0@
z`a<X%a@QeZHffNGQhg0~zJ*$)L7}&KIm;rY@-?8SA~+61l@cP*_4nZ2)RQ{jJ4Ud$
zNG>rH9*D8@RjM-0Z2)>bGRi7;SI1igo(hf8VuwO?j<od7QK*%E0(2wL?ont8?FG6A
zXzwW0LH`2!G0=UZ&@@Vy6-txS(&MAh45|b=S58Yyj6!KsIJ0l{g#9TmxQZR$Md*jL
ziF0VO>%$pnd{re$PG5~A4Lb{NLT|oDUwn_gfO6b>I>JqAS8y2dRl6z^1(MXf^bf8{
zh*v8P7%blymiJ&uK4d`Ix~P1YLl2<+4o(92nx=uj+u%PB<DaCY6+U8rRAa}IJ-k!w
zA{cKNj5gl=$CeHgK5D@0IsAw*(3N~R2P{pCCyU`e`SyP_EfPDvcB<3&oW4J5>}ql=
z2HAXxZ5U3z0`pqu)wtdVBMqT{U0MII>YN-u$o&jESJ@|GfMMtN;>iTqYK(@)mDip?
zQ}Y~~D@k8>BMcYNP`m8sxtKlD=25rVU*YX7@itlOcltW{*zXDn{vn+o{4G6}na4KJ
z*35#k4S=K#Q<~Mym5scq#y)RLQ@2$vgOZHg!CKZ}$oK4v>HW-L!R<Ugr-zQhr@nbu
z#7U?BM0@(M)0O@S6=zM&J`76o%iEeZ<c!+if->Y>VSkI(W|gpKXj|5Os~!Wk*10}`
zbtlJ5Qh24FA+XXvUl>3)JtlRr^Dxhb84>-^M3dak&rEpvw0L+~T~}~9chDS(-3d2E
zVqc@nedP4yrVx0^spR(*6gklx9qNbCqpoa34>gbSFKNAJeD-v(u%v16y4c0K*mZpl
z$pYHx`96D)Zg~PMpT3hlnGMnM>^k{_n|yIA?1yMi_N1{-=Zb9_xJl=-3s^7R$*#7f
z<x-VW$SzWg(lX}AH%5n7ytFsln{8{aM7V4v!|>7K6}`{U4P{FHGhkjiY#tiERCV*^
zoNcM#d&cO#&W#wpj^uD(img3-*H=f5@O6EJ&gPWOTm))z2dYa4oxYI1UrW{4EqHVd
zJKc39w&Xe((H%_ywi0LQ(d06!@D}9XQV{tBz?u>pj~*ZY6JVHfI?)fk4~{u}!<fg*
ziaw9D_IcpclsNdJR_aWg+|zW^J2|Bk)Z~+WlXL&hEC;U9VBRb$&#Pi*>505&+1oUb
zw~F1QCHX;lBa$S%q@SR!{0=rk=>;*io_=2N@Yq$zbw6^$7L1bJH`OtSoLJXozeRI=
zk0br>_Klb40^-x$N`LeDGgs<b^dsl<3vr)##|Oz-xSg$|mkVEIUJ4cMbyPTYCt<o)
zG?&ex^5U7SpVk#uV+(f|7f#$Z8E0!*BzA@$uxsI_aO`%o+I9a%z)S1O?jVegq_6Kv
z@wjw~TU;*$cCazDBhW(E1J8E!exhk147?9A3fbAOX*tLh-Jfb2|MlrSuW7B2y^zmA
zCLn(S>H195_HjGpUm-*P#??o^)sQQBe@LG15|Gb9_I|Euy^#Ho2O#5LXj;M|Js|ad
zscDxvkkm96@=ogoO)G)i2RRE;x(o-9p(~np95OKgJJO?T5O&CV$a=^)RhCa-TWMYS
zcm!ul`GghGbMOcUX30XM4PQO*gb6}&&USnbEE(xdFC6Q|mj!H!j-U;B0AGOEbcIF9
zJAFLcL&52bI#$D@6G4y`X1KOnmSp5SZw+S@Znjy5Gjir<xaMVuA6M>m4R{0aikztz
zjy(nJHZcDvY&)<V1be2yaAl-kU^T#waZJr9wETKZhGW54(x%Iv5#V!1aTDN{Q2vbZ
zEJ?Fb@1f=y6UyF1_;(8yA321`*I+AV2uy`+0@g~qW-MZD^wo@N)b|HlJg;FNrVAF2
zL=)c&+7(nRBbVv-!G)}_kIk%3#R+rB<{P4$B9+zd(g{X4BoMu*^kSut)zhI$KkKFo
zXlIgK<#$KLuQG<$%Pfji`P?oMJN{narJkOyTEY@^vdYI^qYG7!Vka7^J1G_(OPz8F
zIYLQz9SkY37KJ@l8M7=3E8cq3f#F84t<qS@MGfIvCeeOyozV7&6`UWgxusJ2G~B{e
z3eQt8-|#$-Rgt3f1hgo{=KaiDDBen@OTI@IoTqx4HebOkBJ=%T^BrY!+65;;J>aOc
zdwx577u6_U^O<Tg|DiWLl1dFV|BH@kJxM(^K2}KkYueqd#Rh8b&kCBoz>jC<rUi=k
zb)187nA3pA;F#xWKOL=6<m19B(WwQ>-JLSWIL$fJt_n^d+rW<r^DTkLg=h2hrc!vJ
z!W6wO<+ZmKTgvC4Um^51Z`5+ywNSCIffehpWU?pe6xf}D-7lhR-pOuawTgGAuqN<p
z#*2A!_vj<NBdq(14eKL@klyE3`7_fqZK|CV;cqWJGU9n^&R-L_Bo3?(>*PhoiRb)Y
zv0n2AXVGA-5-AkcR(?59S`8sh6!>bfjJfd6ONGGu#iLozAHe5QgQ|F!32O*HP`sQq
zgtUe>!TPAM7V-m!^06VLr}VxtVRhjF!AnGp*rQuth1G5GO0IgtYSzCYcSH%EG(C&x
znz`b%z<YfLe~a+>8SrAhcSV|D1d0uYc@~|>%{*#>zMjahSYK$-MQdMR$!f+<P*1(z
z-EBAMrp%;MKwsnkDNsy#AGsU+azda9iZoQn=7p}OSVP6#G-0NTS#+wQ0~fiW(JxCP
z#1h3CE2y*4h4gW%(eE}xZ^{MaUh0?4bX7``r8lUvDL}1FF4?@s7_F&7Hm_=wPBvBC
kO<Ym@d*J0wayM(Z+U{i)^45rf6;K1T*TtWzT&&XfUx5VA2mk;8

delta 4001
zcmZ`+3vg7`89rwp$zI52cS&~dZeE+jXb99~69|cxB)ds~O$lTJ2_?lsY{iCZYdQu3
zQ*FWk+7-~?Kch0OouF-%sn}rb47FMprLA@>#!j`U#co7&A)t|nB1&(+bMLt{rPF`r
zp6{OTf1Gpv^SXEN03SKPhx{Iv#dHZ9z$W*UJ*I3{jN0ZsZ&oOaq{Slg@V26GouLXZ
z`)?ex81-Ycm8+hZw2Dcq*yw)V!2hpvLOLf%v98PB^gy`u)=kUftG$lQV?#cFVf(Fa
zy2Co#D>}mQr0Gp<)ypjUt+kZ@ipH%fC05$#J!>|(Z6WHl`6+Dk^PlOOaJ!BAY|Z>-
zy*6eG@t+g7_uHRAjAA-y$)-K_yKTK2`)(&S?Ow&gtjq?Ll@6rk(^#6nD2bXI&)s_-
zK<jokA?q7+HXv6v_9-dp!FCIJC#($C%087=eF0%*VHu+L5HooWqso#p#Zj3qlOLun
z=}mklJ(sQqGNi45$u^@R5;2R;rmIdFK&N}i<4~Qw(k308h4LLLwK_7KgVR_$EdzT8
z*r91`I^7BO39!S{ScMKa)W8g*S-l*YrghSB_+%&<&crm<MNUPfg-S-o*ff?7MzTg@
z%GR{|-NlZo2FA{5SI?t4?zgggoCtjGlt<IFp($6vMJSF_t>=jWyqv{?as`_ilZum9
z)Q?|Le|>sA(h=QVeFC|U?ymcF0b`vjf9Eb`EY*!=2hAx%^9yKF%M2<_hq^}ub^QP;
z&x_@d>P6S!u-9-{5AB~=%8DTIKSbj%q_zsL_-R-?XjrTi{<mCa4;t(;!FC!g9L>l3
z{(os&>IjOK_>rM#{G-&55p&xF-#C?e671ZYT^}W$hp*`DI2>Bbl2cLFoWvo)O}W<E
zy@)vF`bHv^q~o({mK=kV=6z}*W1kFR$Zz5$y4(I&0oPsV^G0{t?-lPm5`Bu6bOi?y
zg}o>d_WO1F>PKizRvuqZo3jea)`Maj=XOT7)isLOX#8YfMs#Q0VpLKsm(cy2(BvQO
zb@Y5zb-|fD(VcG|MofLa(EpcR6)*KWr(AC5t28US-1h`bQXhP&Y3E$4?T?{y&edgq
zj8<ir@bA(Ov#(jZ2HefA?~h@~f6rwsa>F4cObqhRrx@DP)6p7#7wR&oD-)^H2&S_s
zp%=>}^7YC@WL4|#>LnsTXG{E2q@yMN30!=nW#^<l3XfDh1-%7DLA*}S@(OsgmR-=Z
zEE2qv)_QY&O|Zg{Dif{o)2;DOjn>NPKJOY{PZQn{cT%2j4u6-H_*#`y7scSZ?EC4b
zzS5Z+a^wuLFX*VRfRE5wUxVe#0F|qS{4Huz*UdQkx#4hcfFAY*e762Nq{~(^g&1uC
z>HR)^pw2D02bL4hx`Nx05uqr!7~OWpzcU^7-iv;UscUC^c)<wvq+bM6ZT*qU!EW=*
zU}>MzOa8J2YharCC5nkdu3%UXw=x?40wLW~E>CNTEwu$oxZ@Rst;99^aO!3%4HV>0
zs7U=7$mSB8NT2AytKiUQ=QjFQVEK$!J~R5UOV#_4(ccK0<`Rb(RjY1^ReOR?2j-;J
z!ZfvQG$-dMw;VZ7<9W4IoLA3Zq8sym#Q#FcytVu^HRM+-x8krv6z3tjGk<`m=u&>1
zucf;SzCCmKEXFo@psSA3`GRuCxChg?+de^agWuzSqrJgg<tv~>H#gIB!HTR_y|(&o
z*IkRTmxN<IX~A8*i*^-0$`uL~9dMMo^dOaVvS=Y+O2x$s_!xB+H}E0aTU>bc?Q^iW
zR<*?65G(9_q$3jlvNO8-+M7VHY%RNlG&)kQ(RYexIqwfWKo5r2(fQE8z{o!}EsQ1;
zz}3Lvq^9|?!iUE-P26!q?`m2va0Iv=m;^oqEIgrUBSH^+1sHn|CldU6fj0<$U@tHU
z+zuQ$scFN&N#N7KvG+9%SM}G}$OoEsMi41YbK|%fJf&$Rzyz=s=str0z}SbHb_AH5
zgdT@r;m6PeyMZy_7}ZtH<qy$@s$4!vJF2eg>Ujqd^_6K9*?>Mg-Qbf3rI*5H7al)`
zj6+Qq!42an0XJ7?7!&?a;|UX=r}7>O%**A^QuVy01A~aD39H4K?p>C)O#fZhNM_+S
zo8?5NzbVt*m?`g35%ne{iy?_PDRmLtW^ftcDyDI}z|97?KyuhDY#3YzxEBN$%`CJ$
zIV01tXeMd%6t4pl@zrUzW3Zi1`SWx61-cI9v(!2Ns<J~!|9)vDS`H)g9hi#wl8XX6
zzzx#A`AhkBIy=7s<+au8yia2u=1D7&i4HLrw69vVIL^{D)r)y(lrLyVSIkd{VQQ~y
z@YG&qxZvE-B~1J3f;>ux>nix;^fBK3<f*Ul#N@?{rge!?dwtL|TjouCItb~Zt@Ul#
zXs^`=d4fK!{|=8)q%n&U78~_8Y!Ej(#gr0qgi|7VScO$l#A{Vp3aAEp%s@ssy|l8K
zY&Nw=me5GVXLYecL{S3_8~94GG7qbz@yI$JMyo1?8+}%l#pnl(D*6#`ye|-%{`CWQ
zF=rNhh;>tClgfvwrKuudep;E)`=w2i_P{1g`<wcapQ2_pV1BcjmapqCUaY_N=Ko~O
z)V@u7nu9z<2b=pngT)4G&ett;-J%M#f4Et#6nDRy1(-W#dKxJ1@LA%c(Pz}nb@!s~
z^@h~Caa0O1S6Z&%Ey#<UXG6=_N!^QMIz`N|jQN68J0(BnH3)Ov2Pm>w<z;$#y(lZY
zq@`OdIf%IzH-HzTU{{P{6YX29-nau=k0rnyS7`5-(I(|QnY;ZFVK1gUWVC|gM3;lo
zK8Z_Hh`CFCFZvdEO==UxhIYLHY+PuS^HMuW8<vy~+$L-0UyeCq4|b8rkXvq$Ao3|B
zSMnn=L*{954R|-p`JdUVlJ<i#PrXJMR%7txtXGb8Avz+YMLLL43(=1+;Bvv^M<i=m
zkF<!%DL<y$;a&qTl2u!25N5KrNc)gj^^kuU!0s_4=Kn-m6yl>sNQ>;tqoxJ_4UI=t
z`(f0v4J}!fZ%dDmoB;FaIw|>)pb<$-xRCP{X;F(j=#swdiVfD>21O=6x~1h|?h$`k
z;qa2ymd2)aRpnLX)uY~}6IL81X-~CvbDpFdV-=pq?M4g}{j?8k!ua!JCeEhuScQ_3
zY>EQyHHvxI^iWrO&E=!aJknycuYCX)RndwH8d~9|?iC?VAwB|xn26(bMTNdDbz06;
zHI1*RP|OWiOM#W2DdvA%m;xO&<n9R3V24{V50Dt`?WnnY#K!QeBV>feJ2c)+<DE6i
VaErl(sAyG<GA952axm85zX6$0(!>A&

diff --git a/examples/04_sample_from_cdf_simple/example.c b/examples/04_sample_from_cdf_simple/example.c
index 77eee02..76bd474 100644
--- a/examples/04_sample_from_cdf_simple/example.c
+++ b/examples/04_sample_from_cdf_simple/example.c
@@ -54,7 +54,7 @@ void test_and_time_sampler_float(char* cdf_name, float cdf_float(float), uint32_
     printf("\nGetting some samples from %s:\n", cdf_name);
     clock_t begin = clock();
     for (int i = 0; i < NUM_SAMPLES; i++) {
-        struct box sample = sampler_float_cdf(cdf_float, seed);
+        struct box sample = sampler_cdf_float(cdf_float, seed);
         if (sample.empty) {
             printf("Error in sampler function for %s", cdf_name);
         } else {
diff --git a/examples/04_sample_from_cdf_simple/makefile b/examples/04_sample_from_cdf_simple/makefile
index 69adb6d..ac5766b 100644
--- a/examples/04_sample_from_cdf_simple/makefile
+++ b/examples/04_sample_from_cdf_simple/makefile
@@ -19,7 +19,7 @@ DEPENDENCIES=$(MATH)
 
 ## Flags
 DEBUG= #'-g'
-STANDARD=-std=gnu99 ## allows for nested functions.
+STANDARD=-std=c99 ## gnu99 allows for nested functions.
 EXTENSIONS= #-fnested-functions
 WARNINGS=-Wall
 OPTIMIZED=-O3#-Ofast
diff --git a/examples/05_sample_from_cdf_beta/example b/examples/05_sample_from_cdf_beta/example
index 481e2d085132be87ea531397543dec0e4a616ece..071bd691fe55410009320000e66bcf4a724cc15d 100755
GIT binary patch
literal 22176
zcmeHPdw5&LmA|&+#7S&jxx~dhBdF39l*m>dGzp}!jO|=g1q^Y(2_}kS%dw~}8(T_D
z7t&D22_g^D-7e6=eoMFg+Uzdvre9y73q`g=a6$`Cc+@34OcIKO6L6IPCXn3yoq5Pt
zibeZ%yZgs}p3l~q`OP_J&YU@SX6{_+ext%unVprT$YfH!qEM;&a~vrZ6+<_O3`nUm
zU&+Vk4CPW~63S+V@%mCuP|MR`b}G#kxCxZxiYPM<y_my_q|!p7B$p{YJcZXts#tlM
z<Z>u0%IUQGe9kAStWUOA!t0eZT1yT<lGG?iRi%qqwNkn?l@6r&MamIOa@9hvTF6PN
z2sufmKglQbKVSGWo!U@G$;LhvC0(vo$feVA!7nMbL223Bp@zSvy=o!X<l=hL=`ju~
zlFIgOf}C)*;P`AJhxumFUZ!z~3V%vU<4vi!Ik4vHxy8*5vzi0JaObSf`B%@ndakph
z)j6BhlU>vX#nh4;meX)V2@`H)V=-Nnss5TdOTJTk`JLW5Q~v((lcVqd=>7}T%cKr!
zlXxhRyc}XGQ~o^s6A$nIy~uHl=t2MUEkrY=AEQEC$|)?bXTYZcx216M!W|&v>3J*z
z{y+x&whXvG1O8kF{O^F<@SmP)0E}m+Hv|5|4EQ}6@YW3Y{TcA9f!pw(o~i+iXXlf^
zFHdKtGUBm?!tZLuTUE2v+u(2auMKpB{OvVMuWN1%`fGe^n*EB`ySAk@=<Nvk+CyG1
zXC0T(-`S=#ul2RG_>|_>wMxgG?V(1>H7clSR5}7dYO!9l=xu8c1VfG9`lhd;J`h4X
zXrlgW1b7>Lfo7$?xwZalN{g=<lxl+p39q-l)8}mr1bxkcuak|WS=iX(LrabAe!sHB
zQ&o1IceZo3bKY2f&RBNASoX`w?A&B-j&m-#DVtqEW_;%0dLgPPlQG+wP(z<|713uF
z>asxT8m6%bUl_p3k*gdLSl726#N{?wc~6uhtUO8iSd_Q0vRN^c5~M6>?Eocxl%@v)
z`7}GIs3(l;HzoClE}-_Zl?#NtTp#4TTtXK;k|ySViI0e&$t9zaBreN?H*%cj3hJj!
zLk3(dk*q3az^6)P&=CV(V8BNWxG<5`HVJ$(N}@{)Z!_S=d9Blc%Q>C|HX3lYUNbF<
za-RX`^B}L-X~3y%nI1Oa(nnPHm;tA`PNrQ3JWqn4yA3!7G@14oa7z-!B4EH(1Kw}I
zWsFeufB_$GJ}{`5$WC}*!UGc?nDD@a2mWRrIH_KGTI)KVt8LBs!(wbVx}%|NeL(Bl
zo7>0Q(C5ttoX}^^0a0ft#NR|UiTGQ(uJ2@JTJRI`gUK>2^oe*+vP>J1MEs}8GA;Ot
z`1g`!TIdt;ZzapLz$fDOCd;(2C*pS|%e0^;;=W{=HZ+O&P02ECWD;?AvP=toBJN6-
zX@O6~FH4qbVNb*hWEs7E{VUYla)tF5<wa@bd1>WYY2}$|<xA4ag=yvdv~o^b`3#N&
zjP`txRz8tdemkxFdRqBqV|lr@eO?GY*S0U3cd4Spb|L>XjZtg|T!PA5d1VnQpH*TQ
zO)c__7U_+>Pk7}0*gI!+-4(sPPsE4oL{;SADlKwE>pC`4UE|rFGpUfY=RbfuVl;9e
zwCaUxV8W@LC`ORcw&#44$b!*OAwuyhlt(B|(TCI-8)1n8-Hd{N1>>K0Gm$IKYLOAG
z=e@;R&!|bu+NT{n8=4LUg+f8DKGewkCgmw!H!h-TCH$r3TGyh)wa99bw?lc_)<qwp
z5c~LyuE!c+;=Y^{C}**L8Rg^aQ7P<MraDF&E7bRom8*`}NwxN4&_A6lS68Cl_xCrt
zRmWar4r})wuf|o8P^~lXxP+o=<*SKu)q2TtUU;EF+iITaV3usMYv6Q6PpSoX+vg|T
zYHeRat$iz@u7alGL3MKkY$qXm=3+&GZe4XecT%l8pca3uIu5E^uEhAL#m}i*E<*<G
z@5BG^Vl^<Cg!&<*MfxK9w4PyGceGIr_Nm=;SOgV=u0xILhC!`o$fj+5l9|z<&AFD@
zrFMUx#{J&C)Q#Kv;M7ep0|rnJ(cgk*yM`VWG}udFc#YkYPz&npHhg;QMN}WMTVpfe
zS>&zA!|-FMEVt*LbgL`>*eJAw?AHQqu$RP6eyZzm7o-^*1ROoFH#o=@m)nZ~!Z~Qj
z5&J%>Fz{yd4`M)=n5|J$g98b*?nSWn#x{^V41kqPH|#^c1fMnbT)2A0Z{g}1yOs9@
zwJY`l=qFj|K!XW&B?e|?6eU7?Vvi>^+Aw5Lozr*`{Iqr6--Do&hp4&y4ZLsyO|04r
zkB#Dg1l$`IJ&GDMpWgsV;&3C8!VpP#)ZiObVqRG}z#N;4QG|P;AgazN#jveF*sa<J
zFV@CX$Lnsj=2aMYxUrmK0F3}2VL?#SONOgk{vD$R=^^6ihv?g?<M=w#I&b%dh#bd0
z8q%ZC`wad^MQ0Js2xKUK7BSebZkdTGzbl%xXg`#xo0gyy>7_WQVLf;{p<W3it~>(Y
z9i{+&ZC6(K7HHjvAlnPA2xP}mRMZ`XhhGE0OhbfYFjq$74t9Bo_#!LSEx$lJFmYA^
zOpMnW>Y){O)<SD-gMGe1cWqL)dTY+-A;x<mns^2?7MmBe?REBQS^$4di<m2V$E@F+
zLQadk`2<nSU9ln}`M(E=z4x)Md(p?JyH<<5`Wq6#6%*AWJv0lwWz(kg$9oVFjZS%e
zpjDlE`Uy$l{}2EKz3UKqK8oJ{9`XZp#l<5ehxt(0!<vUZ0cMyq3WnUk?8&)_>{$hS
zR_q}&V9u&|B~--j2WD*l($mY(F3mu6n3yurj6%{CRf}I;U!)cn{umn9tFbFT(e>dn
zX}=aZO5Tyy%u2T=igR-~sShMJFlw|$+q%t84PXk%V^c_Y3Z@UT;kLfm+hAq>r}<cQ
z43KxolHT>_W0-R;hbh_~bPH+YVzqc+-DKAxOu-cI)_f~!<6mGwfbylM-8Ce07A90h
zKGh;0#m>iq-yIFl(RyMwXx{IM9C@0YgbOXl1Yyj216bBBhJ?Fr4=keKI1t|rqZl9f
zQM~}wu4wGL5{*$kYwuQ9)P0Kw>mn&wUKKgxuBnQgg2)#EEY~70Ymw7T(X#)t;-Hco
z7+VJ6*c%_|dR63Etdm;gJx}DF*cog^U};s?a8~$u>}gb!CzveMJ@k1v_lnVzntEjv
zS8C1-sHd$N8}Zw2r-*jfxo>x`a(m+^Aq%~(s3-ESC-Nav!P+C$-x~hDSg*+bt&=0Y
zJxB6gy^(>QSOFtLIj(`o<i19>pW8Z_jNm*2Anz8S5)$BHrbH&=I?I^~WBLTY>cDLH
zSt2Cx2lR%C_VlI3C;U^y+T2e{FqVZsAy#H%TKkb+W&yYUAk7>-@dDR@$ibesa$E<f
z?mz?)M3eN21DG#8_Pki+Lt05z<lG6RXkb6hF{)!wbz=6r=`fh;?78b_quFwM9@#k$
ze54|_8-kdrSl0mHvW1!+3;<gO`eLh5Km2UtSZBokhxG*)dR_mxeIP3_m#--7$F7Ju
z1RbgEMR<$-Cie&V;fLUWfoPe#m-mP|K_2)Y`X<jF3x2>|a|_zZ;2QJ#MX<zdC(m5O
z%z<Yfgj|Ye%;XuqIIduBtcA_wslM1k@D4vI=M##*vHtDDhVXy0e>+a`{;fkh8T;o>
z_U{^I4*K^?$o*ygi+^$w*Yz#nAKt^~C(O56<P-GK%0|l^y8xAF!eciBFxzR@FbU9B
z6Ts>ANU!QRfYC$<ARg<m60XYm6;4)ht<k`nY2d+F_y{nJLH_xmY1gD$kY}gf{UsO2
zc_|l1Zmb#R)5TE(wHZ8=x4sMpU_2=VikKO2(LX{k#YK53A%MQtn(c@H>uvGD1Y9%+
z{KLOu^F{ZE>W*!gbhlUPU8C82LmNqGkx}Yv&%4=L&k?hhb(&odsuf`rg)3`&Xm`8s
z;BvZFC|iiqemJ&v03Ib!TQCqFh>Ua%eOmhAK2vC(Ru31}njhL;oweYV4VP>6hq|KK
z3|X*m!+F}nl-(lt_0%--N}nrA+2LNpe2Ly;#w(cAbvkRq<z4S)^98r7H=AAWP(tIn
z8D?f_3wpyR`V2Dy>#zG(_bpYC=V_(*z+K~xytG_INKyA8b;tc!nYNd4PfWpD1P8!^
z!|cB1h*@>KsAVM#Q;ZcWK-7v^3FaTb%Bk)fAeXRgAxcql36($&l1qM99vOm9+#l{W
zg<RwlY{+W<Nr6|~3tm!pd=AYE*YFDWg1sB2z&ljxR+iR(ir&@P^RC~1?Nb;tx?HhW
zj)ne&pmwX)r2uI2j5|VPz=!yyi8t51brPTb(pDMw_8R+HzQDnQJ4r|RQ&6%Bkxy4Z
zf#)go@L9O>_u(V%$Rrw#NR-X`6o22w;W-%lTpteS)Agox?L$9KZJ=#~=NGu~>0fu;
z6ZuG=esVaE>X9$gYd#spZ7Oslh+_*!*ySkZr`v9qD=43T$@TDgrf%Ou!hcn_yFSfQ
zAJDPwM4!}dKXxzFGdJ@@vUk74BLnJP{1RS&7<(eei}6afTiwXrf9Yu+KKCAXAX?Dv
zOHSbd1Ze&yCXG(!rR?&Gzy>`tm^`~y*gYAIs#Sjxg3UB)7&}~L5T33U<3^&dRDI|f
z{YABy&E~uk<2)UVxFu6JABQoen@*^kjv}Kv-U42DH|+Q|H?o#)AYM|PFRG73)r#It
zQMH?%kgCoD>La~sIe_pzt`lgl4(&NG^av_yzQKMmKG6ahh8wOtIEf}tRtqa&+^D+c
zO;T0-th8m6o1i*y(@}ebEJ@%*hZ}(QV8AnIFJE=ODjN@VW82;DI>c{vx7o=>$$hHY
z{X7(6uS)*Cjqd*0x{SN9JheEn?jAI9SS`-K5WQVjfY#<sM@}t1sc!ifj1z6(ZJ7P#
zmbvL?P&T$}#s>Zg7>6zMPHg>fn-Tl&3FEFg-op2_s^bXO1vGUu{v{R*;(i;*{^5HK
z^MRAx0Jo|4n@F<z48;lwzslv0#D{^zLTE*<PuQWJS%&)o89zUVCFHeo<~0g|apGsV
zgwMytadq=z%o(_0C}R9@#^y=@<cX<3$P=6AF@PdWj&7d*GB|2vsCDC}2`3NToDai`
zkYnM%`Pf84qu!0k*##!r%0s@{KAjFOuz|-k9ehqLen#CKpnj8CyI33`04VybIVW%z
zAi{CQUe!4;7K%S2M(3c62DSS^0Lkd9!-0j2M?@gAn_mUY9;jv<%<0`13)?78VG%9i
z9*WYB*N4^O7u3ynGn)o6B5W+_SV_cF?IF1EW!R{0=^4Aeli}0!2_%Qt(d}WzIQ>ts
zq{S(2?D^b(G(HNgJ5}ck>izGL{_YrkQZ=lh)nO?9%6RR*1EuVCj|D9v)D>lSpNL!?
z2djwOS8(1MJNj-)*sfw7E5$*c2>&LUA%=fBR-Ot^4gWAU3^~IB$?iw^BtPDMKRGi#
zil!rz;)`&CWRMp3D}RA}9hZOVG~|;v=VHI^zSXmR(VQ>hUJlFjE6HX0Q=E&=vIAb`
zx#A36(Xkg|xW0SNm_jpLhsAvAi6lIc_sS#hxb^98YF&L;w_XpwOE+nTmBnB&^|ue6
z5&K};Ba4%EfD>uL0}~#Y@W6xzCOk0Vfe8=%A9x@O2lD<-UrSrF-&wDeUS$geZ8JM;
zj+qTtC{@9A{`L;Pt+kOLU&uDI(H7`f2r^Zv3AFfa9c}(#h?nad8og`$A)ho^$~LyP
zlW?##WUKcz*N2;ZA%6q*VpEs+L!m%$t*xV#6fpxjY>n-$Enr_bRjFuiha`;QRqeLM
zaIiiUXbsx9M&8cG#jDHyHUB`O>P|&@v<%o{#U84XHsD)~_|gJPrH}OvpZ%)#^92jH
zl^oe|R|OugQyQZBjg*rK+t_c&JW-S>B{8S77=Imih68KYGM9npuv%3LHN}V3{iQqa
zznSt+Lj;@aZ)`v{rD)bgrfZHNgU5Qxo;P&82)p53Z|Zs}=pN8o(1(X~y$Sc6F;Myv
z<=|Vo-UizFM_s2gr=oXs{WxeXC_P8%2fY~dK5R?x1KkDsDCi*QZqTB*uI~knfewK>
zu=PF#nmerPRy@s(fnEigi+jUr&^@3HprfFjpc)Q4_JZDr^E<N{`f!p<U;3{F?F4-o
zbSG#o&cO#kDd7Z_rL3>ttYmdsvo6Y;oVyb?6OKo-$`JZmOdREt+1DWL$j-(8!>{Z5
z2;mlMrKRY4b=vyejmqK+7hXBXemTJ;Ux)wI@EsE*5Q*>*{%e5`fyz{l%r^Y5hA$r^
z+uRmwS9XhK$ejJ`REyPZ$-RD>X@9mQX0Eb~n0=s?RJYmyQ~_W=_jn27dyeV)BcMtJ
z>G`T<r+K~QVaPrPS)x15E4f<tH1>VpT=d!gwyrNC86-S)Qs`c>8dw_vl@{x^?7OF0
za&MTXp*SCPKS13RqV9q0a!b(;Q<cTG-CSXDY|HtorDXFYk7a(><SI*P_R^`A`IVLu
zx5ZIuv6Weh$}HA0ORhGJ`FIm_IpB|_qzm?L&dy$9q3_Hfk?($h`cbrnuXC_O7{!+v
zamL05Wx{?7%m(ag#>M+E1k3}hmoUv@%l-+vP{I5BDe4r&<TVC<E8^`cV1tssAN=qU
z`Qs|ojiBx|V;ywaP}WfQbHwos(#d?j#9}i|%HmG0m`3(BgJ(Bl``?hKMA$}P9>n=v
zhCmPKw}Ev6s{$pzKadTN?J!j&eOpexb(3$GSgK9mL?iRv%(-O-=T=U;7UFar-F;lw
zuO(f~w{A;Om+5*-)NEP}!6L*Nf~#^`yQSaku?!$Q`Vk&c3PFlr8?L9j-_!N(q&&vC
z%d~{#zXSPxjP`)JY+92+RkmBeO+werV|-O(S>Au>e;n8C9JED=u8(g68#G{q)AjK<
za0TbM^xF&Wm#T66qH({J2FG&>Yx5NJIH1S)eQi=XktRGa;eiPcOn6|z1OHz=Ab+PN
zf0squ6-w!}gTo3QweUon6(l?`VW|_}6QXnhC&}Ln)BOS^v3*nUMK+e?@4_w?Wm*5h
zS-q8@$7pRwl56rVk-;~1S?U+{xD{iGzR*Vr&m&l({TU^E<0B;rF{S!O4%0S<Qniq$
z{RXAUoGAE`5KEdcP~!KAjI=u{`0!|fCE32*qe?z}IgRxX%`xb#D36Mc>7p#{!nTL;
zKQHR><v3QB{SeLl%}DZh+G!d~#CVqqS}o{mLE8k~DCkZ>9}{%9piw~w1RWA|L{P=}
zW4$yRKC38yKe*(&>lWG^%h!a1p|EX%bFQ;wmMhFkt_`zGoF#Let}8hHS=xA=30FbU
zIF6>@V$xTR)00UtDYfI`X65d2@tkCw8+DoR8&4VEsrY0?UI(dot|H?j6`!KWxJt#R
zDl&dk@o9>T^HdyHy&NAyg$Y-`ZCu=<SjQPhlcFlNad9Xa=f{I41cRJc49&vJU@kEt
zJa57;&E@a$jkPA+MO4o=B9gx_IZqjDP0Ht!zgtbk&sWB~K1{e9l=Dd{KW4*fBf|40
z%!1Z&aXDL#$<VZ9#!b!`!gDmwWn)lE@qH#r?2G`E=GXD~X96ewjcI<SX`ArfY4}2c
z%XI-omdNZZgZ*OKRJ;s7*KeDFlYY6*FeUi-2gZM{|9=X6Jo_{C>q(|(x>6%=@==y)
zJUs_kW#KvFdl+~YX8ZK<Jtp|&@0+FWk67m1c%!Gpm^IFchu1NoXVdCwl^Qs;f&X&=
z+h`z=<a$@aaE=*2JsI#-z->lGk^eg5w<^&uazNf!-vN7Uv|6Xd;Woy9uA+YydP>9s
z|Ab)q6>ux`q{lOTx19P_Di(Y>Jo_`~e@*aLi%Z}dq5n$4v(KFehBN4yjDAtO8%4Xl
zLeF^uw+g&i;4>I53I<Xv^@xRBo`cQLz)xp0X_pYwM(3gApL9Qmf!hr2$?T>KdL9B^
zkbQw7*Ez}m1H!505{}S04W*wkJyyjwgO}xbQZ$497X`l};?U2<l%pB=UqQVU^FUEr
z`~-!bYT*p&=Z_?Q30KU{#gNUyx=7<SBH~K!m+}N275;A%_Rx=lsNM8<n<e;%#D?t=
z!9P#nyD#GaJNH6%iNNKcv9m_t9)=5H18EfeJ4JjBin<VR@_)ID1JeKN1%I2^;JzmK
z@5`Y7JA!|;7%z6#2bu3>;QtwLiidPRJPq6?+A@+q2>ors!Sc+pK;R=^;N0w7580!@
zFE=U>d0p^tOq)+sbbP${H$0pP1ncPui?a>y_HdRc_3fdKP#9k=saM8c2I391c=0NZ
zpuYpJ*=T6>u5E5z<7@Udgj(A>yuNU!Qr}9CdhiU!Ilts;*Z55Is*Zry*WT{C%Zmp*
z?RVh?9ljR7w;|lpau-AlB`-KbDXfM^BVKRi&F-ZY-ijN_>0KnLybX%Ce8mm!rB&CZ
z*05KX0KyA#ycL?zpq1aOc$au?EOUFjH&#|ItElnTxXU~hq?Nr)q`o7}tVY@R>XFh?
z`Bo8nJqW#KWDKGAjC|IsRK$ZPZ)0<-FXYAZE%s&<gZ!73kiP?+4>ov1cpB!F0&7}3
z6|cX+hlgS9RLK<Z7o-@Z(K^+(ws^zA0G?^Jkm)J!R2gHQtAHou4T8(p%$1gSL#<vH
zy|HDCZT12s`52B~<zlG$vUC)GGYq|T#ZXJ0@RGs&r7S3<yhWy&y(h-NKZn>mTD?u^
z0KIgk>P84Q1cKghhrdDLfm0uD=?J&*7t0t#SyRa&qv#PYjY;JZt_ig&orv6~Kx1gm
zZ0_y=e<(=<;$@!puIcCyc8e=SG4gQlv)9ac4X?2oYmOpAG$$NfT_pi{QDAr>&X@>=
zPD`LO6mIvk*Wskrc)fTfj|_TxSI@XF;Ne~(7p1-)Csmgrrep{j7#1jsv*WIokZ%oW
zsGZX$nZp}<{OxUu6VG1#&b7g?v#q_g&EFon%K)tj2k?rkK!bqXWmU66zO@QtZ}N3C
zDb9wwf@p%%P&==o$HI6Nn^N$iuHE14BLR_XYYr(+HpEV3ooicBq)dmuUU7!}ohZ^w
z1MXIKwK)AvVs2?_fHqm;I`|~RC1f5=`C0-{#>Ei`iW3uD3nr~`rry5_qU8!(CG0(U
zlJ~Qc%JUad&3I*=C@l+g)FI{NeXpc8)ESdJ=a3lf`;79nSU)I9YUjn|r(8cJKo7%>
z^71}fQk!HHG+q8y;CQ^-D(#o&N|MU<#qmn=oJy+oA;a&rI3n-MB^@H?B9V=hmHn6J
zKLKFWW!b*Gf0tCATY-fop-vW+JCUO^At^8K^Cfi%PHDg7lXL^>>5NI1<#~XlJB7UT
zv9w>n%64SPr&3;?8%P=zjIw<xFZ=(HkS`N_@_a#3dG1E?B$F=xJz(S$sUMp^o&*)!
zF#&R0lLcfR#lKPBCgda?F)|~cre66`n!G&MNYXzezdKD{o_|O>^k?MhKFer-a+4)C
zTj8_=PHDeX`*fPTJWr8Sp5IFT^y}vtAuq>Y-e*YKM8+9Xdi%da9kEDxc|IfQF@r2+
z)7w7`%qTC<aU{)66Oi#E)xM1ke%`>YLBiyHynLTxI=?JS`X|V!26=gIt%(~#$w#^=
z$+DzpfKlv<X)k%cAm5KE9Yfb2B`GKU_#eQi|4I29i3_)ri#etEGA1c6=_Cjn<&*cJ
zrCdt3BXvkQS<ZtnVX}RBzhN!q9O?YjmN6BG_DhT)@>2dwX}FX(QW-l&*3wd5S}hLT
u>9&&EkoL&=hpMIRT!Mc-MhgqljI{WX?ebc>w;NynX>r4NVVXc%QTcCMduX=+

literal 14204
zcmeHOeRP!7nSTQe1ZpzKmjOQ*t*n3&f=gX=Ynf0+hm|=?TX6$LlMp5lnoMG5a1q%E
zi$fX4DD9d`x2>XEdvaWka;)bx?zS~RS|Ufgsd7BAQq3Z(??}3v0&Zjd$o_uMz4x6v
zlVovE|J{3-x%c`0xX=A~-wgMMR^BqEpujP0jPoUjQ0JtepwzeHY`<PWja%aQoS-w!
znd*!O#)nT(DD2h9^d(xg>BnilvFJIj(DaG|C1?7}G@t2>hwN!D+k$$i`i|=qXnMy9
z2!e9-)!;j>s}j8F16rx+jocV5cYdL;JLeZFIny6fXvg#>zJ>L5>u*@RuztgW`nuSb
ztqZm;xnaQ#i%a8;rB))ftm)|E6a%1thX}N(_}KU#3)XR%8yo*UjsJWe{wErLZ65wT
zji*Ppo*Om(hCKZJ8h=wBzDVPj=HXw~_*?Su0gb;c58tBkU&+J!SS)B2_}G3<)cLtH
z55I`aXf^oQ{1L%A&c-}^sm9ljg1=kv<DJfVLE$v@w^QR+X?#?uj`J;ze`rQfF!yQv
zcQpP3y|$8titw1mmtGkZ#H|w){uzzmGFRa@E41@-jo&s^;pgc*9MkxJnx^o-0L*q0
zpDC#RP1Az}oj=q3|BW&66E6u8?W~e-$KiHB{2Gln`6<%)XEeT3^Uu`y_L87rO#NJ=
z@&B&zit0E^G`?r1!Z5#V%QQZw`B`UdU(xuk%Yp)L){CqL_<&&YzA7uR`f4iUiR$J=
zWkYpc%!$XEn(JbTTBjk}5RWF<H8kE6b?P@(H#Af`^^F^y=+-8uwt9Ws(VWq(bqS}Y
zzOm*TPD6D)dZ)HI8g=4#H$zPvIw4R?$y!I-tBG%^tgVaH*)Tv@qYHJhn)T5{HTrv^
z&GBeu&4${_^^IGdcy&WleYBZl;vB(~F%gX?D)YGoSsmL@nW$@sR$8*o*2d=ern=h1
zpI`4ZlcB6K(O9|2*#ZYDV~x#-LGX<Yc@QGd<m%?;>U%5WTN(@_WNRUr<8p0%V|4-%
zZft0b=?qsoNOuGGwi`IAQ6n`sSL>-U`#8a?G;Xr)XQH>)>1@DP;)Cgx#^q3g&!zZW
zBpAoxo^uI4+|y@i9J|lqbA^Vdq32%sSsM3wz0EjWmMpP}TCiC6>?Yyk$H(je7h^mf
zpBeb%x63U=0AudsY@KCQ($jZ4?rsZjc*)pf;Te-@y%yf&inKnBXWkf>X|9Dgv641m
z;Tel*gBISLYeHL8B^~vfJP=W6;TfB0lPx@DOe?YQdMgxCz{1;koonIE+)$v*!gGF1
zTVml^7p4U*JZs3bu*TDO)|+XoEPmFiX%!aU)Fo+E7T&JSO%~pqndEP>@T_suwpw_5
z{p_&tc5Uvl@T^hO_E~tlCih!-vt}rBz`~m~N_@M8XWg52$ij2&n$}_A0|tq{)53F|
zo7QdN%^pEokH&Le^4Ie@{}q2+^t%<QU*8o@oeJ+dJ+SKbm1*C*a{^8{UG$!GH@oPD
zroWHx^ZUqubagoOAK{(-EZ$oG>!Cq>MF#M-+O74+Lau*Wr+;ndIe$yhV^;>iv9NpF
zc>g!AhV?x1Hx>CX$y^ATr}!nK?_t!j0ADLoo#9qr{K^1EM=|=kjE;mm&g&%`PW>|6
znmqr>&`>yaY-Q@zm8my6Am_L%kjoA-rILN5c8nt$ZAI#ca45;so}?1H3SyqZq)TdA
zDM+^w%Iw4%WFUjtINu5coU{Il!$3y{&-p7N1O90*`PaUL985p1Xvfd`?-DxFLlV>C
z-#ZTe`FDQ`r}9$seUVz%T?_77VQ*=D8P>tpruzu;9t#b)J5hV#RI)?jXzN7R!Ky=)
z2t#3+p&p?#Ln_X-5K9Yo*P!EnpadcaRJhV@y<pDlN5N#)11|!pWIK>B=_`b9=ltuE
zCj@Ru9;D_DvjgR3_Yi3^FUo2c77)x(RXE*7YGx74FcO_l)KPLj@W`F??qP_4ZgNWc
zKG4y(?$Lbg-S2`ALaiiXGLSmx#AJjp87xQ&>`;Qqty+-1Tcrg_JnPTR`42_FD{30w
zq-EK=i?wrbp0j$^zs8<q&i|<^SR3npP0=dR_Ig`|R>9u=J@zImKB{0C<kev$>4geH
z`S*J3Dy@ON`=Zw1Il)SshJYD#OELi-81PlW6T&Mse7lCj8V+lCIbc2IK@8JLe(xOY
zMkV+kehVjGhP#(S@9Hc_3H?m3LcG8e^SgieFqQnGC(8m$@?@kJ?_RFFSORipGwMHA
zZr2~d`?%<*f^+`)hyjz2g6mF;PNL39lGcP=g*xk>-_2Y<1l`CVa*cIpW~AacN%sk9
z-EmK#2iVISl)!Hwuo?od`q#dSoWBMg$4$O-b17awKg>PWnN@I}qcBo&1kS>{H4*NP
zzOk@~SPpxBh~pnzHlpJmlQ`R}nDU)v9KR6Kg&c?dBm;mmNlvl~5%NW?N<L@1AqOS0
zoGjrok#CXSy#vQH#qWT)+R%94@bBjS!(eX_e2b{aRYj>XF?!dCMwOd)C0UzUH^-0Q
zrmkL#lWxw{ZnAfu#jS$||LxL|-F#8-Euw;MQl-fe7tGawiRZ#0Ucuf%5;*JRz0hWs
zX^BX*B-y)jwHsBQ8xKmhcH>*Gw?fNOR^wkNjoevUxtL95ugOW>MWBqf4NoS`+J<i~
zk!_gfunm)%c^Y11uK5_F>;YT{oH~Gy&oT$_b21us08jhKSnv+uZd?xa0ql@E$^ksi
zlT=5<6A;6Br4%L$Nlhn3LAs4l22MFInM+6z04AK5keII&zWEp8B-wih(Ei=4@U1c%
zlEpU#bWCWdf0R6h2J9r2(2;b=4X0NRWVq%t&x)~_?^iRahkFbU@nN=l>B}?1(-Jy4
znyo(D)<G;B*y@F?zD$_5zN~F+)_P7*X|)D=0IX)n5spG9-0hT4g*wBj)!oq32|eAz
z6BSPV-C#Ji1$VWWZx&R6lM%a*Dd`?a_hv5-A;7pfgmN7HJXq7AKDjgSeg+j#2{zv=
zB~X|;)&ZpBt^=FQJ=IqN$ShIugifov=A3csi||XydYYcpns}7BQy}ckb`>xt2sn$W
zG1mXVw*D3l8}~zeAtim{P%3mX9Xh#laPl>%i`M(0Cc~>gaN*Qp9%7-BXWoJ~m$B%q
zbSniNFw|*f1u(PHZBJhjSks|y$x47edVXDM`*jQ0WLD1XcGoKHhro(Yq%S<L`&?L^
zc14pQbWZMzEun*);Dgzu9YDT92axZ~B>~Adj!?(#$AhhxC2qgrw67i2uBPT!iq#!Y
zo8B_y)QnSJ4AS8jT%7HxLJoO(O#3!^3cjNh;PU96f+ZE+G51}JCB3<(S!Ga?E=)zv
zr$grt_Z6ih=NBE;kbAKW_4y1WB9ZfOAk{shKKlCrEP843!z-pKhg&uS)W`tzNpV+G
z#J?sq=&_YCKiWwuT?!6v3+GK-r5eG75;~vFeIIk*=PsQ}jV(ow&X7~v-i0G4DRR>N
z7%UhSZ3hdC7geb~dx>ocVV8J4!h=23<Nj8$M|wxx+9Ul`+VA-vfd6{W!-E$h;W|L&
z9&@*7Jp;I6@JWPM$qFm-(Ivl;mE^IauFl!ajw%RNKKp^+T_P%4LZ=Cyc5ejexy+R3
zx(p`Y`so2lyhseXmuLljR50KcYEV*&)K`Q!Dp~{~vpBfZhE5{b5M?6RTg+qERB-Vs
zs3s2&Uq`Gfap;6j7=^g8hI-wnr&#VXT(i%B?0aRJa<pG>CNA=wjttz`?QdNSR=t;C
z3e;ZWo_U)SQqvhe+&2b`doWzU!#mWEAPVFFI>S4vHKkrjVt3XX8wyw+@hG6r(a`BM
z{zY8ct<2gzZS0<=-P2jS1F#FkF1NtLtp@Ih#gvNR%rOp*p?^pX(%Ryj`H?q!O;$R@
zyHunvcdy)zM*;?hUDk^`g*_6RBJ1)8lbKH(C~1DAU<xpYoKP>eu5_qN$F?aOTbGHg
zi?MZiaYT9~wyvDmfJ1B^Q!3IyV=DQ^UThMAq+=bhU3lgPdPWZc(+5{EEk*B4;cUlz
zZ=-jggX@g!pk{1OhuXE9_hj8{H*U7m&GuY3In~_8C(3o<OOW#KUV?9rajc$#a4K>@
z`>`^8=YdH2((<9#X#9C~j|v^i@Ju}OFwP8lsN$a%kZ~^wmZJSxMf)v9z~Wdx=3a|8
zT=eUCZ<4F<euOp8Whys>(-$)hxGgakZ)(8RX<#MFaVNLM|AHUndO^5ik$#Ev4PgDd
zU(>WFK*QzC_3u3__<vRSP(MN)lsNf}zpz^jY#Zx;v=fInU34`Aga>RQYE9ikW<nEA
za&d<*26eJdXTBx`#Uy@Ss%-jZ8b7f`Vz+CF^k;A}@RSJe49)xd$xK5m;u>Bz;|$@k
z*p6Q=&)IwPk0Gk=&A)`OdU5rPjE3EtUGHY@-JAE}WvIS4_mb)z*txG(PI{7(6t@my
z;a2s=O26_-Cq+TJt&cbS?Wmp1-{MP6p}fD*!yeGng2>M4&GKwc%FaQ=yVT!$Eif`6
z+>2QiMIZXTy*T+YIDwPiUmEJiTfGN#z75VLj<2xA>=~bqvRkb4TdYqe18+4c@yOh|
zMTss2?)a$DmSmSQ>Px^8|3U%%drKhe-;Mp8-fF*4#+h1$@!bYyP+BQ0g*vi;oB4tk
zxPsbCwe|pz{^Tr0%4)YwxtV7_%}MAwaO+ebQXVsZEy?ErjPiV-U)^P2*ro7p9#2A)
z$1__8I`z!<6m(n9Y&`%~?QJVIwYOeI!)ouR-jkPCd(&V<=f0iPQEKm?Cn-_6zndtv
zrxYd+k(y450=(n{l=~>s$4Km*=2Q_UMDH|jXkg@><_**@QtkVRaG9OvD6szBKhw0w
zK~p=;KMVdTh4*$EoJ5^C#=5ObJz!qyxixxkNS7ce?3PlKc}afHOL7Rd)B;~1qhU$r
zvQS#RBtMIHr8>!+#?g}eQ%_PQc@M;}z?H&eC#h*wC24u>n>mFY2)Xe#OuQA|FFM)|
zpvy+9u?Sebh_zx~4jY&`u#XgWzZZ7HMJDXkG8z{454{TU!tTJz*J1A?b(GqA#FJEE
ze+y#pR4GjECpDcE1?e_InKL-J7&b3*`6o#7L2!`B6Rq_i;>tZ?l6yc;3zr0Z86r#z
zQox+kWR1|MP-wE<n**ikbF5&Lqd6q{GxKv*3A;H{LYmz=Q)a-o|Lsiq$jkjK2&<AC
zBcoxtf8;~tNG8?G{TK0)UFUuWsiWln0Z&rpemcaEd!;bBhtzaZ6r|e-Wg7W9cr(6)
zKb71Lq#ka?FfE8|y2`K;(@9dynR10*o5+MdQzDbHXG#K`NNbXP{+Xh2qns&CVx8Y&
z-8xg~V&+M8l!;<}reL+K0x6tot3YS17KP<%sS<j2hnRs$O$%0xJcTK;cH&8}p?2a2
zCz=^c$Y|J%t^32U8QY7h)-x6+b(9%f=}D>?tKluKHK`OPSCN`dih^_-q0G;8<Qou+
zcZ@Cv1$p*;P{^G#NihEyRa0Rfqc4M*t4g@}BCG`QXU|<DkUd5pr}k2<JwRsf7&Xki
zRmi09IB+xE0<9Xdj?psZu{uT<0~qZXT|%0*!!1LxA|;f^F-q!)$LM9=<b5>3Ox`;(
z8a8=f{J>c7mfI_+YCU;#NgZYKuJa_-<oyU@m^`I0Sw?C)DGJhUgfhrk_L({ez&xae
z^~;U~qBAEcUUSN8sAKNWsh5vqsntE%!Y0f^uh&k`m&vMMM0$?~BUzd>qNPWe-Opy}
z4P-cL39+#8@5|+Xv$FN(`#>$fV@2xq^3^L+AC%u-p89xoxOKq>pue_q-Y4N*{e>=n
zAQe9R$=Jl~MQ{3N@1Xq3)Y+A(|6HDWuY73Q+u>bD3&J=4cFX(WU0sD=zq5Q@`JLr=
zl~*3Etu6J>et_&^?<o2L{Y#H&f(5+ot4_hz$px2QG@)=G%n;ABp{H+XXd!u)PvDj!
zc=q}HBpQD1{HBeje>oV)&nJs0`&)`9KRjkdamhntL&bsCaiQY55BioDm+c-OE?%;0
z!mY)@Sn-nb;<EDMxyy<J%Zf{u6;EDPTqyo*fcRYKvF$Ri%vm_?n~5q8Cp!~khI*Kn
zqml!m2}I*?;|s6)@K^siC-5!soU<5hqRI?XFf=sheZn{%oJ^0vrwH?IpXn6%H-vGU
za59GkroK%R2_F#nIN?cz_X&K2Fn?x}*&*=15Wa|Tlfcgr=6Rf{68LGt#e`Q0{3KyN
z;h@0(NSM=*DHHhHgz=Wb$po_SRKk-5-b*^CN1pD6=KMWjyd8Bi1Gvi%&Dlnn=Xa)0
z;AX=33jyL6cr#({fSFE#*AwQGY37i?YYFqIA#*_BuM+0<A+t~5zaq>M$m|e!8R5$b
zHwk<bVf@V=@e6!C;n{>&2|S;0fN)UYs|a67xJ=+#gs&nT0E~WwHn22!U7#)&m=h1o
zowH#cexfHot}Ct4yaP-1AgFr6I}8pe0~?&_9gANi3Ga}R`U>i??v0&+r3a3j`SiB%
z^EXy4dy;dn+lA>jJy>??p?gDo*MfGTxX43a&FA~efLOXh{n9Xi-|#UbvBpH8rn<gn
zOMP`Bx?z%I$c?q6Rwn{;Y6Ergo6xbUrS&bHbZaz`sEch3#2Xu;0rh*sc%Zhqu>nFi
zO>#ob&CrgN>z^V9YPZB{5_OF+ZB^7<omg=@ey$jZH$`I!Wx)HDqc?DzH~lMuh~cl5
zZTarRcW~d{S2;&Pd`t_fftjBjFEaL_q5Dx8_qTvXdN8*ZN=p~wU;OSZbsIO<M|DY%
z&ubUv&$V*_{tJ-BY!Kg{$<NQC%mi{FUK|I~n)SfI+x1|armq_xl#$toPmEJAH%+!g
zK~Q>AF1xk7`K=u{6E<$DY>5d$>A9)0g>lxhaZ_b8{0hqbHov#e#r{7uY6eOZ(X9!m
zbVD_MNnE-<9(PKc8wr%w#1o61(&(nj+UDwps8cGxP%f3<EtfVmM+sKSFM$=wjLBfe
z_$76;xe3bejm7V6P<^7=Dcw{Z-y~3VLtTwiy0H-sG)L>J)u^c+V0=?N(VVDWKdRsV
z2f{k=vgF2S8S{TfDs&yneZ<nI(H>?4MyypMZ|XUr8LM=C>drCi!qjO!WJxpfW__H{
zjAlK6MOwc71h8CNM&7Kqb$Wr7XaPev_D#PBvK$+_Szigg;LQ3WofiiqXZ)vne$09{
z|L3PuFKFLv5UnN`MhoT(O3Jfm>z)_Xy>y;E1!e6UKd7Ga#{V8I-!o51^cq3kn|;He
z$!7cCK0hd=_W6o~mlL}gy9UR5ur1%E<-4@JbzqA~zC7Qq+Wz-y`Cct=?3g~^|K)$8
z{Hi}u{@YsK%zwWBeGJZO#($FsmKj;h`Y`{O%KRUzzcd7`!^j(z+c2PQ!)Eq)^M9=l
zX+}dgc1<1b0^PRXUREU%?PW@^>;nFCe(96Z-@B+v1bY`rH|#6l|3|@LWNl_6|AdXi
JIL{sDe*lX-8+QNz

diff --git a/examples/05_sample_from_cdf_beta/example.c b/examples/05_sample_from_cdf_beta/example.c
index 03198a5..7fe0fe3 100644
--- a/examples/05_sample_from_cdf_beta/example.c
+++ b/examples/05_sample_from_cdf_beta/example.c
@@ -136,7 +136,7 @@ void test_and_time_sampler_box(char* cdf_name, struct box cdf_box(float), uint32
     printf("\nGetting some samples from %s:\n", cdf_name);
     clock_t begin = clock();
     for (int i = 0; i < NUM_SAMPLES; i++) {
-        struct box sample = sampler_box_cdf(cdf_box, seed);
+        struct box sample = sampler_cdf_box(cdf_box, seed);
         if (sample.empty) {
             printf("Error in sampler function for %s", cdf_name);
         } else {
diff --git a/examples/05_sample_from_cdf_beta/makefile b/examples/05_sample_from_cdf_beta/makefile
index 24c2c1c..30cebab 100644
--- a/examples/05_sample_from_cdf_beta/makefile
+++ b/examples/05_sample_from_cdf_beta/makefile
@@ -19,7 +19,7 @@ DEPENDENCIES=$(MATH)
 
 ## Flags
 DEBUG= #'-g'
-STANDARD=-std=gnu99 ## allows for nested functions.
+STANDARD=-std=c99 ## gnu99 allows for nested functions.
 EXTENSIONS= #-fnested-functions
 WARNINGS=-Wall
 OPTIMIZED=-O3#-Ofast
diff --git a/squiggle.c b/squiggle.c
index c68a14a..7d020e2 100644
--- a/squiggle.c
+++ b/squiggle.c
@@ -22,7 +22,7 @@
     } while (0)
 
 // PI constant
-const float PI = M_PI; // 3.14159265358979323846;
+const float PI = 3.14159265358979323846; // M_PI in gcc gnu99
 
 // Pseudo Random number generator
 
diff --git a/squiggle.h b/squiggle.h
index 682a8c8..e75b64b 100644
--- a/squiggle.h
+++ b/squiggle.h
@@ -4,29 +4,15 @@
 // uint32_t header
 #include <stdint.h>
 
-// Macros
-#define EXIT_ON_ERROR 0
-#define MAX_ERROR_LENGTH 500
-#define PROCESS_ERROR(...)                                                              \
-    do {                                                                                \
-        if (EXIT_ON_ERROR) {                                                            \
-            printf("@, in %s (%d)", __FILE__, __LINE__);                                \
-            exit(1);                                                                    \
-        } else {                                                                        \
-            char error_msg[MAX_ERROR_LENGTH];                                           \
-            snprintf(error_msg, MAX_ERROR_LENGTH, "@, in %s (%d)", __FILE__, __LINE__); \
-            struct box error = { .empty = 1, .error_msg = error_msg };                  \
-            return error;                                                               \
-        }                                                                               \
-    } while (0)
-
 // Pseudo Random number generator
 uint32_t xorshift32(uint32_t* seed);
 
-// Distribution & sampling functions
+// Basic distribution sampling functions
 float rand_0_to_1(uint32_t* seed);
 float rand_float(float max, uint32_t* seed);
 float unit_normal(uint32_t* seed);
+
+// Composite distribution sampling functions
 float random_uniform(float from, float to, uint32_t* seed);
 float random_normal(float mean, float sigma, uint32_t* seed);
 float random_lognormal(float logmean, float logsigma, uint32_t* seed);
@@ -46,12 +32,28 @@ struct box {
     char* error_msg;
 };
 
+// Macros to handle errors
+#define EXIT_ON_ERROR 0
+#define MAX_ERROR_LENGTH 500
+#define PROCESS_ERROR(...)                                                              \
+    do {                                                                                \
+        if (EXIT_ON_ERROR) {                                                            \
+            printf("@, in %s (%d)", __FILE__, __LINE__);                                \
+            exit(1);                                                                    \
+        } else {                                                                        \
+            char error_msg[MAX_ERROR_LENGTH];                                           \
+            snprintf(error_msg, MAX_ERROR_LENGTH, "@, in %s (%d)", __FILE__, __LINE__); \
+            struct box error = { .empty = 1, .error_msg = error_msg };                  \
+            return error;                                                               \
+        }                                                                               \
+    } while (0)
+
 // Inverse cdf
 struct box inverse_cdf_float(float cdf(float), float p);
 struct box inverse_cdf_box(struct box cdf_box(float), float p);
 
 // Samplers from cdf
-struct box sampler_box_cdf(struct box cdf(float), uint32_t* seed);
-struct box sampler_float_cdf(float cdf(float), uint32_t* seed);
+struct box sampler_cdf_float(float cdf(float), uint32_t* seed);
+struct box sampler_cdf_box(struct box cdf(float), uint32_t* seed);
 
 #endif