From 7a2015e3e0b21813debbd6752ce2a77bd370b4ec Mon Sep 17 00:00:00 2001
From: NunoSempere <nuno.sempere@protonmail.com>
Date: Sun, 16 Jul 2023 21:32:17 +0200
Subject: [PATCH] divide simple and complex examples into two different
 examples.

---
 examples/04_sample_from_cdf/example           | Bin 22456 -> 0 bytes
 examples/04_sample_from_cdf_simple/example    | Bin 0 -> 22192 bytes
 examples/04_sample_from_cdf_simple/example.c  | 102 ++++++++++++++++++
 .../makefile                                  |   0
 examples/05_sample_from_cdf_beta/example      | Bin 0 -> 22176 bytes
 .../example.c                                 |  91 +---------------
 examples/05_sample_from_cdf_beta/makefile     |  57 ++++++++++
 7 files changed, 163 insertions(+), 87 deletions(-)
 delete mode 100755 examples/04_sample_from_cdf/example
 create mode 100755 examples/04_sample_from_cdf_simple/example
 create mode 100644 examples/04_sample_from_cdf_simple/example.c
 rename examples/{04_sample_from_cdf => 04_sample_from_cdf_simple}/makefile (100%)
 create mode 100755 examples/05_sample_from_cdf_beta/example
 rename examples/{04_sample_from_cdf => 05_sample_from_cdf_beta}/example.c (68%)
 create mode 100644 examples/05_sample_from_cdf_beta/makefile

diff --git a/examples/04_sample_from_cdf/example b/examples/04_sample_from_cdf/example
deleted file mode 100755
index 03e23d9d29fee8c498194405321ac2c71bbcb16b..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 22456
zcmeHP3wYDzx&PX<P@v5}ouKswnS?E}k+$3{DkCWc{?#BTMT=BIX;Yfnn{5JR$043t
zFmf4AH`&|;H@C5~oufNb$1cva1uS@3L8ctgSwUw)Q7~1kDCYd$?~?u#njN3*oafo|
z9KNUN`+e{C-oE#Hzw1Bw-!;cQH$5$lsbpebVN`0mKwzA*Vqk&H0619@%fbI~Y%Cjw
zat6nVdZz%W73olVBFvI_6Cm;BQDz8unLslIm4yU}FICDKDQXmC7Eva?Ov=h~GW1>|
z=oD1-aeW1%o+ZI;!v$VJBOm#uV5X=qn4buHrC-7lDn2R_pHuQVC7*&!(MzcONj{-Z
zk@ROWw4#oZmH%ffnXg>(CBqU)uOR85wCF8=_1~j!q2%i>6mpVbHaQ$gL8b3T@R46W
zYf{~Olhl{0A9|!86_n~(ew}Z{^r`uE)syOc4T084twqx(O`qy$X>?5HCXijEgJNpl
z4U4H;LWGlU<b5$(mZ^Tj%aadp81|7bb6NVji+|F0Fz>c$CzTA+Ni>v*Uv)8+DL)1u
zq7m)C8##d^xZnF?J;7AzDO7xJh|0gFz<)IbeiZOle3H`xKnx|nAqD>DDe&7;;FqSr
zdsE;Kroev@cq=~1sSJpr?DwR=Kb8Vdy)u-XmK69cDe%*Qx8jqW76LJp{N2D`kxWcw
z#BV)?>vZNREt~JD_BMN0`da+n=Cb*->lz!pWtA)Hyv*ZSS>M>;Y4KM!`#l~(IwYaD
zwTacOtgNrEWOa=zS<4;G{u;{FFmE$HtxYwo#n(VutE5&>Q?swZU*oB&y&d&FKWhs3
zTc8lSs%|H&r>4?Z$ExZYt8QoYm34qr+hAbxc&b_}JvF|D$~xcI$XF6DjjgYQx|(LM
zm(6pR7SHxfc1)f!n4K}0y=pM)m=@1ZjpwE~rovt67<pp+W+FF?4adl)v6uy{3AOZ}
z0jtJ$8Zc?V&>TR0({}&2F+8(apTq@uc{n>Q%Nuw(lhIO3Jb_DmBWT7}^IRGvua>j`
zY_!jpLnA&ff!@UFiQj)Q@u#yZ$Ot60zNmS9Ax%C=MRNU6c#Fhmk<myBugZtyI!|*A
z`9q~%1D@srmBI#mwn74)GT`M_g4e2j1j(m)O{H>4Ps^!n60cfrz#HeqS_5A72{AMo
z@O-`JQW$GB;7PYi>kN3gl=G^44S3}vs@Q74qXXh;n*lGkAiQq70gn!iryT~oIgVxQ
zaRXj6;CCAE)c-1l4ES7y1nf58jq^s20WU+rz;!-*-UH`7aNYywJ@Eg<1O3|AGkV*J
zEPZ3<f6Zb{Zx8v?V?BD?uB<NZg4nbxfR4r{j0e!hF`}=fnrP&WSS+@cmubO{M)t?c
zv`|MQ9q}?P(9y`R;$_-cL?e&H%d}8OBlpM4v_MBAcgM@LFh?VI#LKiGM<bQ-GA+c>
z$c^zbEx^%;D_*7zOEgj#FVjXP8u@&@Obc@~k{d76f*g&c#mgEmzxSzZFTC;eSI8SB
zjQa!SoTPGQQu!>7^G5m)lFIKSmETG#zm`<~b5i+*r1I{h^3zG>CzHzmo>cxxQu+Hy
z<!>3wt}=b|v<uMb`sSJMU(VR#MZb3zYb<O*KcVJgQF9zM;j`#Upr#Ad+k`qUP(S8i
z<}pIe1LHFkB;on!1JFN(2!rs4Q28tiBb@Z$Zavr;-a;CJ$3gLgplB!d);RV>fxU~c
zE4Vc?`zgN)T^YU~|4M_qYIwis!IR;Ih?};)v_L<vuhCA1o}tKwbTgRTSN8Vn+QiUK
z!u4y*GT#ESU%NFF-U&g%ioYOoVdq%V@t@p^li}lMW3i=$q1#*)U2^;szEc|9zf2Dv
z)!W`YwXn>+Idjj2#Fw*xgPA)ZFE*z5TE^br3TZxM^v#(o2+SDpkHJtKOL+{{5wTuv
z+&b99Ae%$b&*1daCKI^wQ$2V}?|65X-Z5a()As25KlP7>fDf(_0<vPgqWy}W#`C(F
zRLufkTCBItJV04J_?ACg-#GJS6v7{U0_Ro3#66iWpqzFawadsKS&d3*56RQ)-Q-Ly
zzhiZFzh=qYb{RSDwhrQ>dQ|S!^1D~(pfcxQPzlOa!!&zGSIt)N_G?SJsSqDOG+&Y9
zqv(WckcS_FMNzFh(}O}(Tf!ktVNnou)zmn&ardAReEwRb*#~OoXeZw+(d^-Vt^6av
z_xek;Tj7V3duv>p9RYFjh<^8pg^1CpW|?paxAxY<(Gtya`PoZ(pZ}>^-<UDM&OLLF
zO$TLRs9(#y(^eF9Y2{t8<qa4IN%@Dg^+AxKWtTq!UdWDV_UHPwihWxCN1A=VwqYW|
zS<8P;+wgg0pnnfOFNDkBLSm`|lOF5}?$JB?tnHy1t)WY6r^`KvI8=C`Mq7JG@94Ga
z8~?z~sMe=k%k9$Ie@esT?p@T3rCsP`;R4=5v^TfTrlUo<)ix&#FSEI$T5g5SivMm~
z9@YD8mhfXA!w57Rc^LjCDobox{Vwg+zt+G{g6djet8E40HqgRd;O1mtuy=$T1X2ht
zvE=~?=RlD?JOxz<#9!jx6>Pw8HfxQpHS|QaikCpz88#C?41kqfHta*b0RPKuS#b4|
z=i%xyn?<w)wM%vZfoQXLLIxC41c*?nAZ$nY+PFk3x(uQ-N5TSq<FxSrXygG>m*WtH
zcc5a~F7(F$KBqvvcBWg{n==}a*x^P3V{F9X)*AjoCBF?`xwVHo7Lf{L-4GDc#yQb#
zOAvO;_P~qfVa@)UODj7J1CP{{Pz*o`@TYhXly#Ef+J=8a_<*~YD7wM=mS#V(ilj~}
zyBLvU-$Pw`403nlGaws_XhtAI__K(?Zf(N^Oj>QBw3&Np6x=fprC=w;Id$v)Gf{0K
zjF@;7zB@tz{OXRhz|D}l2SK(AQW41ZW2mS&1`odq1UC&4j?TO_6tVN+F5`==)HXaR
zVt&$-d5lHM4fT);JIf)pyxLY|kX;^^t=*Ej2W+AxQpIjeiDI7KTwz;C6Tt6i8ZHcd
z^}0X4kDMNS{kH@$mlj<_Am@t!;bjA{m<MePxytq6;Xe`!ModT#cF?T&hE*TY9qGX2
zF*rp)Ym8d@F>&E@00;zq;Q_RK0Ihum{5>?{B0=KA{4DL^>I*})>!E|&#_fV3Kft68
zvmT+jJ@|U~xesHpC43EmMav?yX)?VCW*X+pMbJw#5MBNzOf&;%Kq2PQ!>jYO{4qD8
zg{!sjL_~C7u^Qic@ECbVTJzQG22Fu-SYr{^TmTSuCYYjc++-sKm_oAo6cQML>4R)o
z+7%9fG(P`o_8#&MS<<=sB1EECzjmVKNb6>4`8}(K7aqWbP4RBYc>>BJ<9QIEJob#M
zjCg(oY-#WlJ@{dGI2PXaP+*GQ5w=3|UU%>)UwJXkGampMX1yM)!gaVXa#ievMHC$S
zB4j_O13jvLh3djkc%s7c_6_QrJgDy#YzR~78$YP;6I_rO_4S-i(B~5Rey6Zus^|J{
zDGa%8c5j~MA%ho}2G6?6N`vpiL<bOy_28fN;2CbX*eFqlg<HTJ{`3Rdl01ubQ4hZB
z4!#{ei>)avEv0SAiSX;FCeN#F$!ZwrzH)$XOK3;32;yno<KeaRHX1>$3fFC}WiC(T
z0?3El!jL<7+#Nj4MR0v2`u@Hj8{*HU?EfJ-z-WI15Gm~cHEaui<sAF5)gH8eT$262
z!Ofu2{s$mFmHmsrdyf62?_B#g4i9#A9L*{04EA({b2;eEEbIvm$9{tPXX9{kqM+#k
zIG&42aDaxpGB_ObQs#Xyrc2Vx9ZYH#VpswO&We3Ye9J-l(7$<H=U&o_*dOB{Qf_0C
z{$MA!Kv=(@X4;NOZsES*{*E^?3-?joz92>%Eow{lVSRDivcr*nT6vdbz6Wp8=%G1T
zvmeqNSS@Zm0-_3A*6PVnTVliZYf0vlpd%4s%li~gVsi!rra6e|KnL(;vnw2e;J#-Q
z_fK*EVPdT49<mK&CDy(rV{|Cu4naf8I}`SA@CW(f3)~diRfrxDEuuz{2mV3dWZS|o
zMulr$#uaD^*JP}o2}?3;<e7&-#yxWyd<mY(AkUCvLB`!!4x7nSUE%M7w(k#W{vGUZ
zYTsNg@c*uT@4hG6cN$l}Dck3Yx9|5L<L$ffGq*3Y6e_|`;a}fQF`th;vq%qqjF`6Y
z-pUBCK_yhUZ5cph*l3qk1Jp(nkfUwEPR+g#y@?P&JXT=muq^W`*nzP|+sO>-cyiUH
z+_oh-8*#ghuriKGEjQant@|v4;{s|QI!O$Utnl-vgoh47Yzhx$uP%lG=uZlPZJ^^W
z>IY*&2xKRO0NPrfVM7F1mPU#|8h!x(`hLUbi}usnmQ7fsHqVW<4Wx?$KrpHY2dJ$b
z$J6zWqZxYI8NQxs7KBj<t}O4M)4-nni)pM-b__~;;n?yXc$AR(jGjPG@KjsxC(hG*
zO#W$l6<k=J@zCalX)|74dxc(gpe>ZnF*Ek8y+Gea*$o(V>tG7CJ%l{#DhyG!uhTH!
zqxBe#xx?Ddq^-T8?RdJ_owRkPgPseaeq9eU)ASjgfp?52D&Ah#Ev}nOgD=ub@qw$%
z6+F0DMo3=!0d32D*i&pS7M>V^wHpqA1xNU~^k{}=e@Rb^4$haF1uH<v!jFCB`>@l{
z_VkcTcy<g*A#w?o0QQqhUMLCn!Y8iNyG;H<@(DJuW#3}(ifhI}ZOiANUb;q9xMu8H
zI|AOJQahVp^$A*6VavX5^R-Wtqfo-xB7E>&1hq@E%m;!F8@O;l26TvDimfcy#$jUi
zOIl@Io6BsUiUk@T+)6S^o`#TRh<q9Wx$dXY!e`;iR{}>}!C}-J!4RMIDgH*|?g@ll
zg{S*+XuRodd(e*e*U~o2{Ui>G-K$QxgCE95_xEK}J@Un|vX2LFe*xJD;_%nr;lokR
zPfKy{Aa(#^{w3GL=c$_g8ZrM}&2IZ7P5VX++jO)^YxiQeM=f)4PsCd{9orddUF336
ze*}9-`%96ze0RT2xPR;!5k7aHup?T~?8|Ghi2%&W;M{0rcEYYW4`h%d;+&5aeukx9
zwZvWmV;%Jx`VONE%%in@9K5=m+C#fzFKPLFHW!uX=g}a<5msA&0>(J+c}KhF7&4mu
z4dBO|!nW~u!pL$ua2?bfFKG{lv^kykgtYeih|sZ5d$?080V44A!grvr0{ZOedIS}z
zueM!=|Ik2&;eeVAC(-1|YheY98_+hqPNMRkRkjQW6Er&xY2`=Bk{Vbd41hiuup9bv
zG{<43+~1C4N_XJ_aW$~XMlOo)>$Ub5ArO0e^6yf*y(!`T_1m(w{OGE$L&*^>Kj*v9
zu__nZrag$9mfx>!_#yg<4m?)F5u^6k5H@(qz!75x$Z_C#9JjJ`;5mFe;lNWb&H<YJ
zDAon2Iu`j7j|F)@0nFaMyAAV!gWLeOY4@3kv;8c^3NarR{6`~wz=RJ%i&~$sLp{3)
z=O-0Eb{sFrYr8n45EvqU`U=E+oFCEF&%&I61J*W94`-~O3j}%MK48cb>!)!9MVRW{
z?AV_{QKmwz9fwTZ)U*pabT2}Vhl8Nw6AATtJ0fQXi0J48{yN)ey6C}y2h()Jb6Wmx
zZM~1$O=j)jaex4z=(A*2!oe~em+aCUJ%ge6Ga_^xQqiEbe;bH+^i|-3N5vx|klQUr
z0k;RDIR$rm`(@HLic?rbOSqe&^rO`QE&os2`a8KzhtMOuFX>uP##8wLxNtFS)HZYs
zj_-K*)P4-kzEyObP8p|<fg~wTaS<Zs{$r5=NZqPA{-oXaF3E2X(|@XlHMBbPMqVCD
z@7oZ{kH<V{5ut@4byU9#SGI`U)40VAj}Iq=?K0jlCoW=T_}9`5(f8C~c_cj5_tWt8
z;4^HQ`F>eU@<Z+S<1^!BP#qi=nTg9{1GhY{d<OnkaeWT{r_X>te)!<)H_qF)xHr$7
zJsxK{EYmN?m+4Oy>zgM{1$vSD%Cr0odh#ud?~b!3|G2^<az1qjqwe6lCBe5{vC*&V
zZCzNmUJD$jLz-b_F<4CP?Sg0IKA85%@<|eq=hJx)ocF+a51jYFKaB^{aKY<st*mdV
z^E#@S^UGFWgLOiS)jpy6N><vi%G=!HwKmoes?u+rP-FGATn#XimHFzu)|Mu3gI|=Z
zs%t!f2478MbG@g)Q^=^?az~)D*;^eiH=wezE?!*W^;gDqfRU6mH2SSom338tx=O#d
z8oRBL^Spk)uVJONrI8eJLtCsh&5iXSzj`E_)7%VB7%i%rtu=v$D!;F>!73ywNB<2u
zFx&4Dwh{HZS<8!ml(R2ddIw|QD+Xc~v-wM3xaO{!A8bF3FF{Z`C1J0G%(=a@@6+q`
zCuUr|so?0AHFNOzKc3dbbI4#DaaU3qs6m{=k&i!k&bM+U?-CFn(Mp}LEB}bL*SYn+
zn<)PbSg^UdYps+OO1so_4Sko?3Pkx|VlgN3wXerwwScXF>j?J7V)W>6#~ZPj9p!CD
zW3g?3I{|+Q7zPXh+K<IzhXA|Zj>W7vBh+F`QwF#Luof^c5{umfcrV~~KplJU6M$KL
zsE2$z;GKZwfLj5#0X_;?djj$ScK{v)ycdVfQ-sG^J39mAld)I<Al;{w0JZ{_19s!K
z>QTTf-1fc+xD7Y_CHT(BA;5CLJUl*W0=yS+Cw-R$2^Xw9rG4!tmey)XyEJ=v)>gEM
z@OYrheUVQTCBt#!!<Y@9srXo4i^VLUV`j@-bKZ5@sMT5P*sM#go;byJ1+bjI0&sxn
z2q6>6_2Y92ct-p(5zQui2H^Y44Y=*dXP_UxByq(51U^>aCKFCKTfF8hR)XyFR0RGz
zXghP6Ep6%b=H86-XGfYXE_2p(qfC3#&Ebqv^QnwVz`0bn92vJnlmfwjwlW%aZO0bx
zVTfUKNY3@<tr@G$+rYaWyacyq+$zMnM)B_{>)^@$Yb-XCc#!a5j#U6x1Y8p#=9(>=
z((fE;&bna~)wQB-(-3t!>K;U0FY10P`M!}}V$R!QDm7a-XUs9%H)URLE?7UzZ7ymX
zUTSux&mU<nnrkj_neB7U)?#yBvDs2=&eBIQlGO`YWr)vzRMwXv>&uAUBJ|rN$Q12a
zV9r}_a+$4d8K&GcbKY#OnXKkMsQ~>p=pM^40@noGqrmZ)l5{jiLcmRzbmX6Hz?}fD
zlW@A(n*JYf-5lZTCs4N?W9%9OeK&B=19wQ#cY_{%L$cmP-BygjSB-U$Wu-F4RT;+R
zIBFB`>v?9YX;_-*+&QBtE^VM$h_U%|qDLYe%?~Glo60dT5zq}>9>#1b^pStQk&fQj
zVww~0w-V|%7xmjb^Fq_VLP?Q}cU-Zd<K~XK7VO&~>oSb@Ye^RDTA%JR=e3!xGlw!v
z%fXn3_(o9A9o1~^Mg;dDEV~hwAqssOcZY}{^Fn)^A2HNsnn(QK0e?3lsVAd&RINc2
zR0wpfly%b}UFl#}v>%Gxm~W;)7bO~(Zs1xCIKr3UQw{uf;Jt+4ekmQoFS_uHvtcNw
zaGk$|9Gc_g)J4+Hr}G{-?}76kIPZafG7qTVRjJ=m(eg?u8PadLC}j&E7FXoytHzYd
zCklwKaPg#mzjdiB<2sNhe65Tp^*gT1WLedJ@M)}(P}{LCQc`QkA(_EdCQse69!Dyk
z9MWLi)$>IA2}<||NkWor%(+0IY1yT;PV&?GND0@EJmEV@Jl!h|RQPV0QFezU9qxyD
zQu@`NO3~rVO}u?j$5A-+@N$uqhqDteE4#2Y;N@Ym{v%md?U3sJ*FgMt*+~)$XfuW6
zlyISh%Oz}*aGiu(CEPCI;}V7>?2)in!c!75JT|9fk<cokz<&1s>i2*1X3xIbYG1q}
z(BKbPXE>%h3MLf>c&Tvh<N`;*R7c^J0{=hDHepoc4dHnDJtgfklamRX*770o8H~=8
z$;pIhS8;Ewqu+0;cu&L+XKGv|;<K2Fmqh#srs6CSKa#0<PQ;I5D((~Um<3dS87fSe
z4Xi`r&8%REel#(SIfuksm=oVLNlqpNhni=Mbtd{Hae)yd@+NixQ@_VI)|#*jKEw~D
zMl2mM#>9VjYOFD_3*)~#O~hZsh8iy>oMzPgl1Pu)aiI|-@+Ni(vkZw>v*#cWO;Z`{
zz1Sc|<Y?YYN2i>U6Bs{R14r}iQ1lZd{dRm~$(X(<@jXfSt0i8o6N-n-PBYjqr%~qN
z_;dY6UneCy)w-kPe3R3k>;GQ?Ka~Bc+Vux6XEb~8DkXUk4kc$FuO4$w{~qCT&h6hf
zCB0hTl=KgI=G=Irr)Zcp&WQ&r#(z4krWPeZKr84!NB)uPsapREI6je<O70Z+WhwBD
zz+2B1!!~ky3)_CBs8r|0?_(jdvUB71`<(t<L64=thopRqT=4ll2iOk-Z-M;e_<uvv
zTjeE_x^p;@LVirr7sv%y-8oDoeEPZbh8gV|Vv1oCIo`sSOTC@4iS!E}>X%c}{(Olq
zP9f(eNnap0JnFt_X$tzf6!?u&&N?}+)&1JTz?0vS=aU`4TcvGA@S7C!_W_@qelb()
zrQ&&o)58x^uXG1{o9GSmEspv;otE_K<qcCM-h?<G%Fh>ayoL43c=Jm7FQlNqTGB6+
zQ=*EWm%(Sj{F9e7e?B7ZQSbdx{dKdHQ(zSe`F#?y%O$=?_OIHv-7fLHvLBlS8(RlF
z=}nH?O_F{qx`xuj0<Z@qz9~u0cND!%5b}E@WPc&?PU$>;mk0dg9500$z&<5k8qV+6
zP*3+u)L)Ml3S!k?M<o5iaRUFUlyxG7{IinYd4(Y4cW20ELf}ySdlB#ySIPbv54=_C
zGJ=_04(7ej3#_`sYm$0<z98`Yo($P}pqHi?#jq+VXL-^**23`#q`c53=^vMIK3S0r
zK1A|0RMDdyM-$#D;wWHM&HfgD0ACTPVuLRs@%ZaKc=Jbtm);0c-RN0a*SG?YhN}II
z%`KkFKr5?iq(@SC2IMFzm|i$E5xp72=c#ONu3Y26gQ(^;cppe*z1LG6sIOlG7DLGc
zN`C^W;Sr0+GxsLf{5hUEH<ZvjP!egYnWtpQ4X*j6vlDChD^q~Ni+enCbSXhExruq^
zxfc|>+@1w<=PsI4<|%U(yXTNp{_c{hmH@XJW#gMpoKE!;6M7p7z13t8GZ}9!!8=wG
z5r6BQF7ip3r>3s4((l28IR569M1Fc<ive>L^Ltz1pN49W9}o6CO2mrBR_5_mSK>ik
zGgWet;*~81ZWu_ljrHmYAhj;xEiZ$_Q<?E`8y8yO@i*d$AzmgkNH%!}e}c##5z<>`
z3^i9NUy7I3(0gDEwXh}m(IWy34G^!50dK+!Z|e9fZ4CTVKvCV&=&416@P3@+$C{qf
z1rSv2Yw!eGywyy^UR9vJB~UNkiDM82Qw+UM0VrQ{lSnhP=Y^-GTm6l!6$7KzSL2^D
zS#*_8Je#HN_3&QutY~SGZIlD>TtXp~zd6Tacs<XcJrsRXr+gr~uv8(?4HAbJ{S2~D
zc-H$`{efmLf6Y%~jmLx64XJ@e?;#rU<wC-Z<io@_{v^sW1SOAZ#IS)t1dQPkGjp`8
zsrOf|0Q5HtSgUe)F_E{qi8=7N+3Q%@5O6d#H#T{j{c8-U6#*aK{^hHdC|7amB!A^f
z#>s0dTWXo3dQAgV2<UGXHT1+BPoEPC9@I5^>ne#s=9=pK%)vv#fvjU?BZ`!1@m4X1
z-`k2JO>UrW<fG5wt(DVRZ8fB+l8_<hHNm0sP*qv)gD}C4;9w5Ssr8sA)r|cQLE2Ve
ztB5_hNa~!gpxhS7+QIuVg&2<vUC}81BCK1K6twcv;6HUAq%gF<HS(8Beg$ocP(mzQ
z38@q%Bfq*|P|%t{O=qgabPJ03MOCA+f4kJLpwh4I-;(XG1iv_g3c2dOLP3iZOg56A
z6~DSC^8rUHm40=vp`g0o0|`$g!<g<sf$q)}zq<cWke-WDQuZr41=pgE?&4Hg-IFM|
zm2(gJSlKVpY%?<CQ^l|DQxptILZx5vtM)%6`HLl;x_41<ooq1i6HhY#Bft@fl8?=)
zND|6zvqYZT|KCH6k>4u$6y)C_h;!54S0n|0mc*~_dlXc3N`xva_;?b(x))NA8>g57
z6`z819x&Rk_RR{8QI3`}6rU3LOcKAkXHrnzb1VAf@v~d<tNxGQM->Q;ZZgi0lJ&oU
zIwDd0>fTC0`YwzyCF?(eN+Z9z-%>D?#Gv9w$$SeLdQFHN_VIf#^`6OOdR12NJ@9G<
zes#ZJBoDrd4kSFOvVv!UrPx*a)xDZ}Z>Vw%%^#E$pYr1`fTQ-u`4<osW6oIuh6ITJ
zSNsYN1GAAoey(;37NtiCP<*OP{X;mVU!8j_PC=1OPr8gLSL&DH60cJH>fEm?6u%Lw
z*ikeDola3w?~kQ(Jn2yOs9`YzIAy!w5dZb-1pH=F{3yMOHW!(p_z&M8C@x82NGh`b
E1~()T5&!@I

diff --git a/examples/04_sample_from_cdf_simple/example b/examples/04_sample_from_cdf_simple/example
new file mode 100755
index 0000000000000000000000000000000000000000..37f99c995c2b72506a88f0c7a7631c0a13973801
GIT binary patch
literal 22192
zcmeHPeRNaDl^<I+7_jx^22$q}HBq*v4YI{Y3P~I!W6X0azy;GlAqZj1SR31N^Z}+P
zC8Wlrf|$}hXPfl2r={s`n{0Pmx=nX!n>4Z^gwTALe3gdH#z1H!Oj2qRN=b<JcjrTX
zdPr$ccK7W5;W?J?+~3?gbLY;TdGqE;@16zT`oe+&rjmu-#yEGGK-{w6%wi?Mn48UE
zW%xgpUCYJ+F2ctxi_P>%p(s$aM9Nt}Nw0!3Iq1~_W{N5ciIQHv)IC9zD9UVtC%t0I
zO5RKt$o?v-$}{a%i+Y%e+Q<P&iW>E($X&yWnR`(-9W?nx=@CqNjZ&{s>M6>ko}$X1
z<P-YLk^VGOJIW~8`G3aDdMl)!nbt{pMX3!+%T7jHK1zFyQg5wG)MKV^2$(6V+WS28
zq^qUGhZ8l-OJ#fc#vvj7sVI#%rK<Mex|?QIwYN-f4|c}7r+3e}Y5GmGoRN@o1}`VO
zs11s#g^QQda3lzmZscPzS@Kj~-gJ7=y4!dC?7~02F|z#&zj^O>`-{J=>Yz4Bh7#$k
zAts*k6q+O>`hPES0wa1j@LUJceCbCh_*jm@qj_*&9()pTJ3i+0RRFo_Z_a~%ClB72
z2XD-Se=!gKbRPVZ!0q^$Qxkw(_Wwg3d_x|*Ef3z42fqop9UpW0C4lS9An}OD4hp}U
zn6IH}k*_5X4y+GGqJePJqWSHi&OnoYU3-A}eCs<xoxVuaACCHbLNupfduTn2+!c<t
zGUQrWAdF9US1XGIJE^H=*_5v<9PEs?`kLGBM0qgEx?<4?T0sNNcM|Gr^#|Kob9<=y
zPS)XX2c^<ZgPPCR-0k<Z20Q)j!7q`8RK2vY!;j`#!+`)>=xwN-@0;PAF?TfknbE9s
zb|yb7lbh+B1;-R(&WM>;jNCXj9&?grQ4#(Z^8Yle3iwCNSvN}KU;hTC(|9(7Jk8_#
zIA6?$B_F#oh$~4EI|rT;>9>Qbb!BoeSVl9lus~o!-ooWcCVB<w7qU{6Q&Q`{T6c%!
z^+uQLyTW%%yo3TDN#Tk=vP9st7ErHM8ZzKCS5-<H@KObWo-yFY_&;yJr72viP2%Ih
z$u4oc%YYl#t!@KOqoL9k1J19<vf!%*+_=u~GT_v<N)H-v<s&NFZNTMLfU7)gz$dCo
zK=&AM3}`0pHQ?3^ibdXlYX<zF0ar0X#e)XC+#o-e9rM7L2gW=w=7BK}{9k(Dtaj}M
zz2|g^zN7ee^BB{66Vbx-px(2;WFI#xJ$nY=k@O8SLA0ri<lCrZBy}vEPVeG8Z3sqE
zhci5F1V&Q*8J;!(BdMQec-pv(q`sHoX(KR_`g(?^4Zuk1-V9F{{*lyO8J;fqBPoA|
zrwjc^>hl?%F7P8MPll%r(MZab;b|i@lKMo3rwji`s$B8t?Q6GDZ|fNEFZdc0KikAl
zH}N-^_-jo36cb-&;)_lEMH~kh?RnS4zh&Z2n)ugD{L4n()1>#!K7m2h`)Yo2En~}<
z{mNadvE*8eBT7=D<R4LzTqK}zoYXu-(9Os{K_i;%f&(G{O(7pZ;ZrP$vDf2I>hXc(
zy`&UB1&IfQ#4=L9GlTwApo<7yE*g058ss1OAe~O`$G?X7{#|<f56Q2h&Yt0d*xBTd
zfwpRYNIU^mGW8hbz0-!y>e{CgdkH(Mttp;@3eRd^OeBXOP5N^{SGf{*cvkP@@#}i4
zA%1v`9zU-4ym6+n$=g?aIks?mU)guD!-G$bKt4U?UvFjXon2smhm77=`~s1==c7|F
zb#_u7Q)fbYNSnF^?XrE6LjDmhKl>pf(>~DSXY~HJ=jr|DEqcKL{qP6T$*3S8D=0}1
zweoRR`ZT|`)KD>teR{dxQ&US>J$^D;s_&?&2a{~L2>-Oe!~?}O;0x}6Zl-<GPi;h@
zw1;@D@*p`)tLooadRDU)pT3$LbVonwQ8@~q)~XI}EJI=0^C*PmhH+YD|Gw5;&^@cI
zJV;C?9%=m`2amwn{|xcuoA0O7Bif4MAA%XtR&q*9$_Rw^wYEC7sawfMYCZ8ZkrjUn
z;fDD8di+nx-(rIICSo)7{-m9Z@y3ty5gviqDu}geor4(d-(JO-XZ2o$3XSif6oVtU
z9OO|v4aKfR@)V5WXp+jQf1cZ23C|b2eX|!*1<M=a7d=f4@pqtdIe_JQ{AE3Ufvbp#
zK>bU8AJrsp!@q|3Q}B}>f7=^>Gg*VB1C};WIG#>UMX@m)m!dlFwDUY1DH3;pCqMAe
zTy+Qe%(L2ahi8q)mzob<)ay!k<EOmwbHZXLK&txd!{0ZI-;VL|f&SxVu7UVqf3lpD
zp<>ryd_2a7`nO{|86ji_K~9yU5E_u-ro_i%2#arnF(UqVjF%CSUxgIPpuklT2bm$H
z_TZli-saQPF2-`o@hc(9Z8WtXAK(@U>kp&Wf&Nsv>rnh~|FL4%Au2l*|2X*tMTgF6
z<z7c=GIfEXX=U*vsFb3CqDZSeqB#+XpFawr)sB*lGtg|E1L3@~I0!kaBFWyP;Sj+N
zpw58+T-lESsNjK^+=GgTpSm<2d0$>QtF1nmvkzqDAR3zT{uR<8=*WtJtoZ<ckRM)!
z1PnyWVnjrbs1xLY|E6zB9m(hZEL^h_?c{My(Z(8BQsf}dEaB$BGrxykmS>8{GvpPB
zaW}4j&E%<l$?cFGeoRf-(ebAKU8wlq?cdwyMgM+`cJlVmlj&cGn}hz1`>XqxT7@Q(
zKfu4?z1mbaM(f%O%k=nrh-n)it)k>66ru^QqX<Be1501B70?a~z{!sIfL3`3qlpke
zJgy#~c|QDctUL%2415s{Jh|$5U~q0(4=BRQI48C8QU~?!ubCW|WlfHf<THN~9=adZ
z=J8PJ##$JF@uU#&ax>tf6Hv?wfzqrHKwnoBIS>K1RjD};P4?j5@W1o*rT3h+^L{M4
zef8;{^M%^BIP&okJ$|10+JCB0?>}Co7hK@i2hD~sO2CyX`nAf#+JVE%X|7Oq3b?1?
z*cF5DD1rLi!PsE@OwZ8!?sEq$(b;-4T)3j>fxgCqxvy-#PH%ppCsD|exd%31raws8
z`!MUapn2+h0(rL2m7whKfMLBv?=c(8$Msw&*nC~jsY1T|>OBL6{CbB<Xk53!%mRJx
zK<q8!6^i%QbGzpY4e{sbO7X6z$rFERxr~sC-WRl;+p*8+s}-J@fNK#P01J-t>-_N|
zt@1^^U}SW?)NHr{By9XLy5f*lIj9{NB$x2)6mSW032`8Y$tBO$#fRV%&$<1UsEd5!
z(N;A5J%d*~b6?VSevEca!Zo76Gk5>y3Gfbay=+nQ`{><jN9kw#Zhf2_1+K|alHC6;
zg4(0m76GsUz#a`5$RU1d;Vtp(7$;Ug(^baP*W~zsU!BN<yQq%3$5F`|L_W=ca_{5l
z;ZtzsZ)3+j@o_X7@dRJ>DgF-Pp$G)`xy}uj(R|Z;4xk_JY{sH*^Zpz=je{Fbd*gpf
zPd+<bO6ADcrkmb7k9`~JMi3{z^bVhna(!BbU7WariuIRV51;4j_Oqn?k-FXUeu4J2
zG<JgMlhzx+j*5Eb;hxC!?i7wpsduSsMEOw{u2e6k>Y-+cM-e{vo~}f+pxf89b1hrh
zUQW@;(iVM3(G8UtEj_-~Q2{a3BjTLT6<ii8;EmM|Tlz&*&`zU<vBNBb@?@<lsa5WC
zYY#k`eo?F9tGOt|I8P=wJ1Vqor%}0k>s#8^6Ub<l$AC{s;vkrwKTsD=V6a}&oG)tM
zN@xoPwkEXR?Ih?tq<w2Zs{;`Gvg<8w!(^=zLyw@M=35+B<3Cy;!&akhIEfZdUJ5H<
z+<EQ3*Qu(ir<5({iLE@QRlbUbkCP>9V2LmQ?ZJR2(O#M6d{s3b?QPXIA9TH-_Yc|i
z9rtq=wdgZTxzDuT=TM=3?|$;{D%uV6`#4+FQL0sqZ1^%7IjU8aEk+v~%F)_v9XYM)
ztajgbF;05lYKNV6e7|@HGs~4|`xz?xCEZK>ngDNK*)g2s;vQibM329oeD#m%^h&X$
z#9v*L>fn2A^yfIP3*^Yur+F;MEiaI#hwnA42TpPW+@@`}kY?{iiWO3RRp=j24FgHO
zgjUq`2|Lt_%k<Isxe=C-*Sa~S5Xccf!_{Ivu1aa!=Ak26)f6w6hcmX-10YY_282AZ
zZ8ire!qn)NreB6clM1!ogdX3g_4W!m3@<{Chl7yg3ki*SFCyn*2+>v^`t6R%djGH;
z8+a_!oxjnlp47Glso!MQ!#oZU02F<;;+41wli|2>zvdhq4aFalp!0}|2CerS05Z|H
znokZ9kBC5Sx0nUo9#qXGxYK*DmbOux!Xmnadnro)yfLO#J+E!Mo7;2*Bf`g$ZnI=O
zt#|=0oB$iO`}#-ccP4z=-h<}w2E^%L-Z=dUL`-q2RSt=@|3vCMYTczdpVzj(P4)LC
z=|2_28oD|RrC!O^?whET$1V?AM5rsFBKMQHpGD-lK$6%0AuDXx@Q%50V=KeIjaG=^
ze;ef|!c)UPNPZ4_hAlI{-w=y@uKj*yWgNsql=!$*4Q{~<+VXzouh8Em^dHBT-7tB?
z{i^47Z(mL1Q@EGIW%`xOW%~W)df)UGz{|YTF7kUg?Jdmj{)?9ARM%0to_gaW-uT;f
z@i#r`$*=1@`*7WQEq01-(hOG?gT>U}eejIj2h$!|o|FL`OJg1w^T3z~#yl|QfiVw^
zd0@-~V;&guz?cXAAA6tx-%+I0(77QHjs)zjp|Jgih`louwKx0Qn`7<%XrKkp(<d$r
zM5DpZ_4Y`pBVdpCJG$Bf5qoPm)B*9EC$a_MFf@amvM6kCjdeChgP~4YqjJHA)D!Cr
zMtz;3Fuq2bF(=eY*7&3L8(Qtb$ju-VSyPZK=n8a31>fA#>VqyU@9<UoT#We0T`_++
z(30UrwVWkOw&2SQY#wt&8=k%8p4RW|K1V##8A*rcIh|GbLp8zm>)Qj)X6C+;dIw`F
zZ)lmu*io&)je4q%YEQd&ZC^_HC(`MAa6^92W~xc1M8Q>-TiyVJXJ;(&pXqc3?xOd+
zo=&?#6QC<VcMk!_Fzh~-PSckwcO6fsyGRanE9emDgP`^k>GU2d2Ym*#@eSP0Laz&S
zHR*$PfewLg1>JKpo!$+42J}(T!QZFT^t*#SZ>H0yiBG1}Hatt;l1ir?pb1d=-GJ>h
z^g$cnPN(;S4xK?cp77go1~mtCEodX?An1djl<2$elnTDIlofQ_3a%<0U$P4}5{@Uf
z>_POgiUjJ$<H&?D2R`)e-QBOH(`N{<+Ul(ppVcO9EZM^5U3v4TW;(7TnDkfU(+1yh
zMFNoskK(fy_z<W{Bty?8D&Wg&BuLn9d>p_&Eiux66dyOR8HDLpTfkbv@C1oHgHIdq
zF2d`rw)+e3o@gyuJV^(e#OE=T{X&*~t+38ovD4CEwf7Y*uvXq*{5fm&wsBtToSyLw
zR(Ii|iPkyw)@qNnvfgU1wN}(xZMD`CeG<D7`m^AlEwDc=>$<;iiM3*z#bdSi6cwJZ
zR?O$ZR4Be<+yZ%qSitugc%puF11kl_eG8uVA6O+Y`i?mz>eFsuO~3{S(5?2uAEQeP
zCUN=4P*#GNx`oP+NPZC5EMP|nAo+uk_e|ntr%`quW&dF;LtS?8C5WrX5Pwr4!xmU=
zJ%tOccFVYeiB_A38YKH}gv?&V<Ri#aBFqh}3Gq3LBWMz|3D_234WKj*Un_(McUl%?
zd|gMr_K>d^S{p52K_hcK+}X7TXV*`<73wq(wj$ndrMh6(wnC4!qQ~-CYof@q7K#-f
zZeabSu=QY(*E(3_2W>zp#3}xUNdHti-J8+Jxc68VlKwZLe-NWTSX4Wy%}^DZ5OABU
zYw4(5!)R9YAM@-{%)6OrixS0W6U{vXMtB`Q?Z8WLlocR=`=ud=Uv%LYcV{lA@HQVs
zJ=Lhk_&sc@ax9H`V9Wz!9vJh$m<Pr@@Zla%zvEKB+oEj`B{OXoFe{amUx1NUzd4$7
zy(qvpw0Kg#7o+<IN+p6Ye7B7!^}DgFC9lfg{U9A8Xg99yl+?BTVVP0CtvV>nacjpD
z-Iq{Oza67JAf+jSWd2Jv*0@AKv>l?<DD`RILP`Br3tu<lsZtuK@HUxIb|<9V1W8r<
zYOky0@RdH^KQzao4<w(I9m8&#mn*w)>%#ejEXNnaIIsF4oBJOlncr!fYOJPBE|Oc)
zMoHI7+9l~0Nq0%QThcv}CL|q{bV$-Ok}^EO|0qfQzHs6E`8V4um#>R;Mq~E5&RNdt
z>8=>(T$^W9JF90oUDE{qms_D<T~_39H2p4<?laBFg5OxG-=7=H==Ys>=aetb#Cf*-
zIHuw|8z0ZqJjljNn2L{Vd;(K(m5oniDt@!^Nle9gHjY`Z#>Y@#!R)u^#I4MhV;n6^
zWA>alD#3SH%*ldaQ0t1Z%tFG|MnvQ->@uc)k8dos;BG|cH_1knepNdq^Sf4KiG^Lx
z)cTZ-f1Kr-9~RtIs`VsW9;;!a5fOO{RzX`%T&<R)I<zbmu|M)#>`_##bGh)#AykOL
zs>XM*jEUAG!q=K`I!_>6t^Z2LgG_<JPPsht?iA$k&qi6#kg5J2UeBd|eh7Fj`)U2o
zh5rk$XEG~QhXUYL$|e6RUO45FaeM}N0lRb@UzYOfx~1w)^30`i^&YS1(s4hEg&V(`
zzGS|X!+{^8ak*a9tgdUH;J9jSl-A_Im*m0M0=JK<XL5IQc^lhuoq*K+IsI(K&MuA5
zom@Vnp2_UagVT>*Y^cW^&wIH%<~0Q*k{Y&sdFpvZ%G=}x_ZFeXt|Gkf(skfuo_Y$<
zFY>cnwmTr}DU<l1+z?br{0fds#z3a!!RgmRIhHE6Qr2&i8whny7X(iJF|Pv~fZGl2
zsVtq(lowvX)b&i^e<$ldBj=Yo8+?f54eKGD_fYz&lpm5CE_E*TYmN&U<Nw1_K5?Bu
z1G4Z~9{G2r{2AHbUr70<Ky6qzDopEJNXlc0=1KX#9M?q}FB@D$iJ!ZHUM2B8a=g?&
z>TiHk{pNU^A?3G7d%AdSjD1GpUD9Fv{0wCc5?6z%7TCocm&yjxCgmGteDZTely3k|
z{-5I#!pi@fq<pp92=cQ=h-}YO|F@*PU0wqCnIYuApGW@Zz$s45e&`2o7p7$WKPv05
zmhqs@E>B8)&r~7K&j*ox7q|?Cta{4PK|a4^p;}Kqp>XN{%;h%F&7Za42}E-^8i~g6
zU9@I4`t}cBw8MwDZ*&IetsE^O-}?5@I)A&bB^nAxeEwKBYYx$)8a%{t&Z)l1m0O73
z&Jpza!(sm>A0DfOH{qQe{*HjJCDzfg2`UE82gzubsNorj&sV?HvuJ^D!Qwi4`AD{G
z3-i^jT<lrYFh9G5zpDfg-hks<pvxNcx~0sw(7UA8<Ml16uV1#H$=Bqm^)8@V`TIkf
zBQb6@c;nkf+-~(E5qdcYy;)=wnSqyx;5{eVz=yp|MLgvS`&!#W{-_U+mH69Fvi0eu
zDF);-%ok|!NBu19?`#o@;^im?8JIz(p$_%biu#cCK9y17S<LuUj8|Iii-z!&3$JY%
z6~kLo=rt||0EU=fv|=cQ3FfC{h#uH2-pc~r?AOTnS|TA|8;r&)X3US<d<{#W)e`LV
z#Ug<gCZec0))9$yh&Rk+w`&+b3MKi{m~5G7Ac7c!1EYA<=gS%f<+^BybtA;vg00b+
zGhpw!NJQGsm4o7OCJmY#1^%iTpW%HrqfJt%$Ptl`9~&DK0JlmEZ^ju_p;+n&c1L63
z0Dmt|c8Sl2m-47tKyT~G`3fH4DRPt1+I69BbjSG0ob0-lHz3DDl6q)-y6KE;>WKQ+
zF=r$i7PL*}@Ft%?xQjXQa5vyw-x+gug+pC|aCDOaS{DoAomasY33+N8rbqqj85eK!
zN7|UPWm6}b5HuPVCG_kV&w8^MAIidkc0Vb|TvvOPIr*qNk#(*Qfu&3&(9E3CKsQ)g
zi69-~lgJrplWR{~3u;rGs6#A5LPO=zl)od0%7i+Cf;q7Oc3@G<T7^G4Lzfq9ov=q0
zN!{Nns?K9%F_%?&qSxa?M;}UG-3KdbN0~9Ha}R~lKG3MY0@oQziaG?F`LC|03TQya
zsITs;6}2ltNzMAV13!l|6lK3UXHrzPug;mw_WPkv_b5tV-KQ%$M9xJb8!4;w)p<}5
z7<F0ctNVFH)j1YKc#_pAmfeLMoe?R0b>FXOrIb|mD>+5!;Xa)?DPElqD7r=JD<3QS
zCCvJeA)hLJb&jBDLJF$(mA>l#15&?M%Bk}PMb$YR>64CG|9im5C#rsI21SxoZl5G5
zY|Svpd>bF5zFq1mdd4V>ys3KjLzBKb=g82%pufkYug*gh9r_FUw9hu$ulBEs&NiK|
zDEn2lPnh)8`HG_IJXgt^=g*T;UyXm}9Hv@mtR>?NNsXU^o<)X4l)gHzQS=QZD5<&q
zqb7ZIuA``$hpHNL`zKLA@57SIPUieay<gERuXsiO2wlyfug;0|8h}VjPH>~V!Y`tL
zVpp}V&Xd&pGL>U!{!vnTY8-zAjQp3;UqVv2y__c~#g{QDeMQGX*{GknFLeto)sCt|
z=_#Jx2S=D{U)_J$+(N=EPi+|!or4f2!zELs^woW{Dp2}Hs$xgU5Ouo+w^km&)9oJJ
qV=8;p`a{LacA+8u@1_e2Qne}IRJ)>-?)7r(4=xfCSDF+|Ec-8~ZoV%7

literal 0
HcmV?d00001

diff --git a/examples/04_sample_from_cdf_simple/example.c b/examples/04_sample_from_cdf_simple/example.c
new file mode 100644
index 0000000..ab6ac61
--- /dev/null
+++ b/examples/04_sample_from_cdf_simple/example.c
@@ -0,0 +1,102 @@
+#include <math.h>
+#include <stdint.h>
+#include <stdio.h>
+#include <stdlib.h>
+#include <time.h>
+#include "../../squiggle.h"
+
+#define NUM_SAMPLES 1000000
+
+// Example cdf
+float cdf_uniform_0_1(float x)
+{
+    if (x < 0) {
+        return 0;
+    } else if (x > 1) {
+        return 1;
+    } else {
+        return x;
+    }
+}
+
+float cdf_squared_0_1(float x)
+{
+    if (x < 0) {
+        return 0;
+    } else if (x > 1) {
+        return 1;
+    } else {
+        return x * x;
+    }
+}
+
+float cdf_normal_0_1(float x)
+{
+    float mean = 0;
+    float std = 1;
+    return 0.5 * (1 + erf((x - mean) / (std * sqrt(2)))); // erf from math.h
+}
+
+// Some testers
+void test_inverse_cdf_float(char* cdf_name, float cdf_float(float))
+{
+    struct box result = inverse_cdf_float(cdf_float, 0.5);
+    if (result.empty) {
+        printf("Inverse for %s not calculated\n", cdf_name);
+        exit(1);
+    } else {
+        printf("Inverse of %s at %f is: %f\n", cdf_name, 0.5, result.content);
+    }
+}
+
+void test_and_time_sampler_float(char* cdf_name, float cdf_float(float), uint32_t* seed)
+{
+    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);
+        if (sample.empty) {
+            printf("Error in sampler function for %s", cdf_name);
+        } else {
+            // printf("%f\n", sample.content);
+        }
+    }
+    clock_t end = clock();
+    float time_spent = (float)(end - begin) / CLOCKS_PER_SEC;
+    printf("Time spent: %f\n", time_spent);
+}
+
+int main()
+{
+    // Test inverse cdf float
+    test_inverse_cdf_float("cdf_uniform_0_1", cdf_uniform_0_1);
+    test_inverse_cdf_float("cdf_squared_0_1", cdf_squared_0_1);
+    test_inverse_cdf_float("cdf_normal_0_1", cdf_normal_0_1);
+
+    // Testing samplers
+    // set randomness seed
+    uint32_t* seed = malloc(sizeof(uint32_t));
+    *seed = 1000; // xorshift can't start with 0
+
+    // Test float sampler
+    test_and_time_sampler_float("cdf_uniform_0_1", cdf_uniform_0_1, seed);
+    test_and_time_sampler_float("cdf_squared_0_1", cdf_squared_0_1, seed);
+    test_and_time_sampler_float("cdf_normal_0_1", cdf_normal_0_1, seed);
+
+    // Get some normal samples using a previous approach
+    printf("\nGetting some samples from unit_normal\n");
+
+    clock_t begin_2 = clock();
+
+    for (int i = 0; i < NUM_SAMPLES; i++) {
+        float normal_sample = unit_normal(seed);
+        // printf("%f\n", normal_sample);
+    }
+
+    clock_t end_2 = clock();
+    float time_spent_2 = (float)(end_2 - begin_2) / CLOCKS_PER_SEC;
+    printf("Time spent: %f\n", time_spent_2);
+
+    free(seed);
+    return 0;
+}
diff --git a/examples/04_sample_from_cdf/makefile b/examples/04_sample_from_cdf_simple/makefile
similarity index 100%
rename from examples/04_sample_from_cdf/makefile
rename to examples/04_sample_from_cdf_simple/makefile
diff --git a/examples/05_sample_from_cdf_beta/example b/examples/05_sample_from_cdf_beta/example
new file mode 100755
index 0000000000000000000000000000000000000000..1ad5fec3996dfb9477adc1f8ea7c9f36f073abee
GIT binary patch
literal 22176
zcmeHPdw5&LmA|&+#3Z(^T;gJ$QK`}jN@Ob!ngr5V#&)i$0w!_52_}kS%dw~}8(T^o
z7Sd41NkksFyIr7#{Ytm}+U#!IO~1ZE7m93$;Di>O@Tg08m?RVlC*UdpOdz@YJM)mQ
z6r1+zcK46{JfE#I^P6+doH=vu%-p%s{d$?FJToIhk;$ZdMWIsVXE;(UDh6*78IWRS
zo|22t>B<$#B$Ukz<MqXypq8hB%v72!a1$uW6;NgjdKrfmNu`BENp8IK@CCd^QpL*4
zB$q{5QBJ3w7jZsGWqq=}B3`eg(Q0x4lB7mCsw!T{s+Hm;skA@MFH(+RlB*PQl|oKZ
zMaW4i{YgHd|9Qfn>C}cgN;dYXDCu(5LN1+_3Vunc4NA-23f28B?Nti71{c?pPLFd~
zkyN&K3*>~W1;?ioIm}B%d*h8mRQOX;8gEL4je#}S&M9oHo7or$hTCVh&%1W!wR4<p
z&CXe@p6sGFD5e(Qw48<`N|<mX8;fb8O!Xf-AN;Om!`;96*<Gua*lW8o9(jApQmKR5
zBpym6FNc`Ql)n)F#KZf4A95TcddUA=6VdU~Pf#H)<rJ3J#=)lox216M!kr*v>3M7%
z{DE=s+sDEE<KWMZgMSgY4gcw>3cy%)ddI<kI1YaAIC%3o`2FMH*8;cUKRs0f7|YHl
zfM1!;Ol8Dl6NTTkinpR_iMP(*>R%gZ3;A2CmfX<T9Q0TD)-?JRuXk-zbI{us^0kJ%
zUd}ouqrbgHX<X}TYVs+K&1;pmyIMo_l&e=zQ?IlIg4ANIXwlo!8VH8!y|oQrLwz8G
zcF;ua*9h>|`vQ$hZDVuo*OVq-BPi7d4H8~&ZM)A~9|-yy19y{+q*>V5<U>pKt$x3<
z*i%t*gLjs5mUHfCe)ec~{%H2g$?TkDZnkp{xha!fL1uhr;d&veD3dYUnNUNYbQRHO
z2I?|E=^Cc72wxJw%8{)c6<EhNAH?N0S$S8KBdk10`AC$vva(q*lM<v1Xl(-}eUzpJ
z0=YChsHi85>o+C!hc2e}GL?&kyj&mTyj(;VJ(4Eoe~Ax^pvfkqkt8n512=P=<_hYk
zOoIkoERn1#X27RNX3${+o@c;E47f0n)iwxxGD@OL3~w>u#(AyXfXg|a1U4CPwq7$W
zin7gs^Ldb0bQy4JTc(E%xbzX#J!ZgZu9IoE0nd>j=pF-(0Zpd82HcWFu?QG&)qwXI
za2X?1-EY9hnhy*rCbAPAnDD@a2PQl);eo%K2TrM1oY6W?WNX{A{<sL+jm~H&Q}5S0
z_GS07HuSl(04MYrvq9A93h_5kO(OoLuIpW_ObdP@ekfU{g+3AQPL^pSl8FB-S*8U)
z5&vGYObdM?{>@~W7WhQ`zGRse_C)-yWSJK9MBJAw(}pGyza?3wjZ7l$PL^q*PsClx
zGA;0l_~(*kTG$iuJXuC>U;7I6wp3yLMR{Rbd2U*HW?FehTKV#{a(-GlH?5qNRz8d4
z0HZzcr<G5pmETG$zm`^h$yi>l?VKBe&$XQk=U$;GvE9f&MPn4(1(%@mHeOkP%4d`q
zMpKJCtwnlb?-3q(FZT91U3W$A=oRtdI$04pv`UK{)jEz3S5|p;W=+Z`?YR%&ju?$>
zgH}C%4NN%Qg<=F5ZD-auh|C`e<s%foLV1MZ1^S>meG@EEpqo+fuVDOhmlC<^oE90@
zy5C)-b&r^|jQ!f7bD?QakS`Qu>x1>oZ&IG(b<;wsR>EIeu5~O-T#u|4c`KBoZD05S
z3bBvQ>UyjWChpHViE;+(mr*{x9+kqLWvXMOzD#}Zc&X}$ol>hm0{!!;QgtQTeeXcM
zTXpP1=7@IRiAr1*3Dr8|&dVvPR=$!bRjrpV=Y{9%wC&~@4ra+_y9Q2I^pu)+k9}Uk
ztycFY)ao}A>MCd|98kAJz;+6<XD?F}=+;%ov!~RWgKFVNs^gHl^=gcdTKKHG^>fIe
z{r&j=eXI&5lTaUov`BAczt%lu>x|Z`!CtkK4vU~-z;(D@-8i6i58AZtPcSp;wAt4)
zyVTC_)41QakGgStFPyp=X21aIA^MxwXxGrAygGX^46m|#5^7$J-G)z(y@2XNc57@p
zJd3;)c^G~im8JIVQ*L$TpX!B{ko|g~b@rmzsZVqr?t(O9gMgzu_Bsc-;!=A7KsW~t
zIbz>O6$ajn{y_`~6SLK8YOp_{*1Q1Lp4di`hXJsX>4tsC7vZzYo()&8_#IqbWw-L4
zpmxPR0DUA29cVD2uEfBsjG{znckIVWjW!G!RA<#+3O{Y1`wt-K<Y8(qcOx&HL=&s_
z!DA!%9|rfvg^!{J&F9vEk~rK*BtJwF9yRznm6%sn_A|$3V-(?DD2S@li!p2~5O%Bf
z!;951)$y8Jt$GCp9;q*-7(gSyhglF*^^oD})_=#SL3)rl`XKt2>Nv5Ew9egg2_na_
zpN8}p^gfOM5z$#hGXfdPpFs@vsat1Y%I}C~EIa^Z>gL5LMS3XCX;=@PNvK!Dh^vpn
zcSk6IU)`M%z7<;cBgpnaD+1YZ3>7uU;Ne#RFw+p>7|fN?xPx6@BEHB<b?YzD4osYx
z2NUDfhI(j)oz>7<U1y(X&|RI>t=^XPS%~qTh$f!KjK$^!ZD)<Wk`}<<&?4rF-Z}HP
zr;*blZ#+&Eb62c@NbVm%V()&W>t6IR>aNxzul$xoaK%KmNH@(wZ`!mA`r_S)h<c~I
zKG3R8J^h%Z@P8Npg5Gr)Js&}D{{Z=Zy5iyylEZu`>|xEro&YmU8U;gcV)kU+LiVhJ
zJuCK-88By6yc{ZG_X9Jwf5n+)XqRRnI!sKNXhtCEimHXLtS?Xt^M3-3>($uRAM5&1
ziL_sf93$^YYi6cf6UDhXoYV^v>mM;%qix?|rv@;E<gh6ud;z8pvf=jL*jr#_{-^m^
zb@Y>W$daD*7h#yQu7oMtUUUm-(;~I7f8Au)VNAgk@77!^YU7`0L4fiVXWUgJa}Fj{
zL_X0XAI2`ig5Mbp&(^wQHfTQJi5z{3oP-N4%LHM}di_||E`x--W-lzF;5ZoH1EUxp
z_)$F{)vjplI}(jiJ!|haSJZv02kRm!SY8o1>#nMZoQB9304&!cFKLl8OwqFcDlSnK
za$sy3gk!INsOuGxXRuCck#{|jw_|6q6@jG{9YY!66S1dIO`c$~Q1{Sh;oPf6PHF1Z
zQCz86)1jWWW^BZ7zk?#$UE{vPy~^#4pMos(x}u)QJD$h~Oa*I?RDW~m`(nK!`?pVy
z^mHH1b@fF0yJLBb3}(6dBa?gU*?w;OWHN&D^n<*Uhe}9*hnW(YjO#3G3XJI${Hg=9
z;irj^z#q{YCfePb8lUh_0c-OBDZyCg|CCsnjcM&idYA>=`a?8xbjR~t2P22N-^_9y
zq`HF<L=a8VD-L45^w@J^kq>AkS&?-Yl%jzHG{>lp0o94w@0KHAs<CIUpM_>i?Kx!U
zT=0>K*d7RCrea+KfXfzYx-kH38R(6zM*Yw;^`o5;`ybX9T<A4@WA=fp#9Y22zYn`2
z<`8tGx(DGc_M6-v<cA-E0|ugH?q1#_>I8Y<f9RVWdo1`qcg?M6XB^j<*Dr)6W;=Q2
zQf3Z3^C09>JYy!$@WpWjb7M7ZCQtRo7Jzr?2|1ro{EhZ+KQ@H_oBi8$n)h!V+8MWh
z?qvV2W9FcLzk=Lf*T4A3r*K{03jU$Je15`wt3^IWAFXV(%(076i6%UDGXS%lb`6sN
zZ8rg&W{>o!j)NFYgaG2P1}ovJtY70~71tUKyqN|bjQNiM!x-dV1e$hDs(Cqf>fK*+
zaa@>kab(9DVLn|PRZu&QhjP}JzyORVg+Ku_11|a}2&T9wCnW^X*J`sJ5n#PNK9GQm
zW`lp|*KEG%{6O8c1C#E~a=l|DlW%Av2`w^0eeHfHQ|mrz)-ukp>p`_5jG}O5bvNy9
z_a9nL*9v9xQ91y}R`<iB1Zwm9!~K!rj=@igKiF>y&DCn*!fNwFJ1aBhzr68Et@dz7
zG?O9o_iwyVdzi9Y<-VSpMqcT4MJYSfW0)_|d(3!wlRD01Y`n7LolL&qcJySj>m5pH
zT(`i?3~hc-_++nPMqvGQ-{!uxBJv!q6z{vM+>sZTiwG&`Jgn}zA1l+&67Go$uol4q
zu;2*0uQ_T~9WQ7Z3Bwd)#R?F$Vpf9r2eERh`}@fyESrx~lw3k3kVE8>-<L)P;S=`<
z`%EDh`2-uXs()7C75DrX)m@)K^TIW}!aaZA#tYybDs?JLYCl2mYV0{T?!5j<j2T_7
z*el0Ee@0NdRqGM}w0XuIAu`}Y{L;jm?cP3#&wgpEjC*I5{TyH5;K44^QTim5tU~0|
z6_Dq75<PqduKYvzs5>%=Mk5ksvp&V&w{Un4#$MM4L%DRlX&w8~kJB4z+u->nZhZRI
zo$y3H)Tf;q%AtDXOZ2LbM{t`8-3a2?f?;+!iuvjGJLC$==U;L?d_G>c?<L{CsoNc&
zWT+45*mj~%YNsE&7wVatc_P`nU*VAf^)7xnuRnr4k>iDUIoqvn;_koV3=g0CPB;)P
z==SBO@c;recQcblCv#GEc?DpDp6N`UT`TOKj7HU}zW~8T8a0d^t}+NuQwwn;(OaxO
z^tAqhTF7Q|UWsv@21eYHsasCKnBvVR)y>C{Q5|mr&%XzD{DvD@O*arPs?HbGN1|$3
z&*rGwNl!>s=Rx(69<>xe_+HmZv{!@n92j~86*XUHzYL#ffegb9R}P#+lP9Z%6)<i@
z-TDTpDtt!TGQv$z9k}VJK1!A(aH7KvKzlIYX|$KCI$x2EhdQzC?sFaHx4Jv*<f7z0
zRqcEZ3b9uu|K3h_|7=~xU09A<m{@l&8abjC=3au{uFFGfbEhGv7M@bKei!3J8+aRL
zf4OCD_&Jn~?wYZIe;meP%iM*nA8s>Z-#KaAHOHIyzE*V{#kzo|j>W&kVnN(*1355s
zpJ6_5k{jSQ^?nmccAlkJA>miJ{L%OjkXQ(<$n^<3)U(TQKOp1h7qEo9R?56aAuvY#
z3>ESDxG=76S%f(QHw*=gAI{iP4uCu{1qgX!%UlLfgvrs((O&{bl?=5`+%)0jp_B7r
zcoA|e95^4FNNChM5jnfTL|b{tH`=Gs!38$(n5Kizs)bLhTLRQ?GHW-B0|Wp?pEc_w
z?gB(OuGptK`$t3ZN5tqHkkO!aJ_sNgeKk0+knxBJWOnnbfY}4pjDtD7^D<!@#VIVJ
zCEP<%`qBEZTKK%W<sN3!07it3B^@h?c&a`O7rq1=)vev5*LN~}8a{^P&^o$395+t?
z6D(<QiW_@A_aBRoKx>!kd|ti(UDDqfqfe@aHMBYm#$O(*-M68X-R`lVMTEMd?Cukh
ztK(o5k^3^vTVu!GNeSCktYgJE$P?k;Kr_V9uSUyL;HjY>#)cqgSRmQ`2%qG~+V3Z4
z#z)a~WKw(~PLK@J;(q0?kgwtLPo9B%^5$IZ*WI^yb}pR#1>DPFnSMFBOn-uN(V2F@
z%RE<|r7Jr2d<@rj_gPbDy6cFTPd$-@C-QD-<ZZV;?G3G?7wgt*;dkgJ&9Jf<ET;bU
z!ZTtYOnYQ;(hhJUO?Y6!0}~#Y@W6xzCOk0Vf&T*!WZ*#F-|lN_Y4kg5mEvn`fuL<h
zo6Rw!?kc4sxX$0&=C?K16XXloX4KmPZ3{rAC{=+bzpbss9}MwwZC$;0jX&g*CQI4+
z=2j98HivArzQ)>cqc7yI!(MF4Vt*(U2(GoYHIpJ{K%1?;wYdrG3#KS#t*wxRF}$kP
zRv!-5h62q&8`sF&*|cbN$-m|vOjO*ZD36u^Tcp@S719QLixFR1V5#`Ao}qJJ)qXL5
z!H%M%yEc^J@j9hJs^3I8nXrxhy37+rnNkvSIt%gFc2_vCb}e%mc#fzQ#ZXgtL_JX4
zb^lVzKLrtNuD`w!)s&(cmzu6Sjtm~_DSKbn^#bgMcfX<Q#h`mZt3e+g)b$43bH+gF
zOOyj|>Us-k`=4~3&YTL~*7Xyh)u8kor4RHn&~4b3ZUfy7`Y7lC=pN95xUTO5je!n=
zI<WOV4Vpcq>sCC?je%YRnvHwIO3=Nab)X}l?VuVCJNAKY!}*=r41GAsr7!(ggSLY{
z4B7>njdSn;P)ayKWhvwCrAkJ-HRIBp$=O}7nQ%OsRR+=5LgFZ$%)SO`M|KYWAAU{O
zhY7b>%Pj>rs#DizZ&DUrvf%33_A3b{`5OGMhVPglfk=dh@Lvsl5LBj8WOm?xHGKIX
z+2*!bJ2IOrgXYX<rdX_QOZJUZO$Rb9F>{4w*z5x>r@GY!paKB<xyOqT-+Nrw9|2X$
zNY7U-UFP+chavkIWQlf}S8}!PsqFi{+32(VEnQzoGDvvpq|m)$C9oC($}QF%nfFYw
zWZyJZLvbGJet^2iMco6LrIvzSrV5K~r@74H*pc;BOVO4|9?QIr$rYC3%q3GS^U5tn
zZi}PbVk@x}lvu1KmTYY*^YLcra=;%;NEhtgl9{>KLf@G|BH#T0^&@BtU*}+nFp4iT
z;*5<A%7pz6m<`yqjEnbS5SRy84`G_cmibe3p^W$W6Vxe)$?FXKR>a#izy>6LANb)T
z^2ar(8%Eu$#yaS-p{$|q7l`BOq?7r4vBhSZl);@`HkIsa1kWDC_P-%diLgz;Jc#o-
z41pfdZvkruRsl+We;^Yc+hr<C`nHsO>n7hWwp5zFfkx)JnR81F&MlvMJ;docy61$h
zUr)N2Z{3!H4%3a6sM)j{f(3{(1XuagR!g7RW9dhD^dUT=6oM4LHe65lysPV-NqLNO
zhiNg%e;e|B80~&@$<zjes!X?l8-%W<qkI*kS>Au>|1qxH*=UOrT_4{9HekRAr|aX#
zz!jY1(r+)gUn<7%i^ly@930Cjtj*KV<A5IH_q9poM4Irxga;-(FyVm-5Bz`ifc%}7
z{9P7pS16^^HV!Lz)WQ>OR*>+(gr)W`a(b~S%in|1{Q@PieN*s7HkRb?!Y&hKS^xey
zy_ulLXl+N5Yw~W9!8dkU>J#<26=R9M&_@Z+BUqyS86|w<BP9tjrSfJD(>8`urI4rn
z1||B+HYI#Xh$T%JDDiC~Bkhg~K0I1rNwzQdsFDv~PGkK;a|}8s$|Is<x+qJ#u<c>|
z&xv|`IgXWOKSXnXH<J9FcACZ_G2X?3RtmaW&=x^A3ECy-V}kAxG%9Gnpo4-A3#u4@
zte0lPrxoSz2N&OP!vdRQ`I>Mr6t>NG&T$sabcI>TwQ*LFvuKXfbrr`yO&hN>;VLK?
z!_o9xO!~@kdNL^{rFu-<tlTpuo|TMqqb?JE<0<1i6`!og>mU`+R%Cpn;uk0~u2S(S
zij3b>e5xYjJQc@PFUQAFVZzmK8xyxE)-lG>q^OE*OdLwa`0=0#!64@qL$mNQn2U@E
z&ztZ|bNPFGW335y5tXxyh~&>tzSqE5Yf?U&{M~9QevvZP^<l!@pqx)q`7s+-8WEm1
zVHUKGiObn?REDM{Gj4K55uT%YE)#=NjPEm1VrK-PG{26;KLa@NZ%XqsP1}U;Ny8Ti
zT&@c!vP5QQ80;6*rs8Gz`F`61ob=0ehAF|vKQjLF{r@xIW7$7mzn)-vrYTkOCLd*)
z#?o_;Rpy^JzDIy(V75;m-{XQ`u4hvBhb(h`ywOu)%o^v#!)ut(Gimj-N)4Ra!2cP5
zEi@2Fa=j~JILC~io^kM1z->lGk-wYqTb1Y+I3Vw<?}R-z<^1^E!T8Tt^e;kBkyzj#
z7c9R9ZiSxoc&6``Q@@JEf-i^Xz&QF}75tUr61Yz2znbvO^XGw~ar8_^zo^|!qTL>$
z=R$#71zsrd=?oVI11Xex#6m95!RC#_PiHe}mk`rN=b_}E^m!l*+-7J`W;c(c=ON&E
znHMW^os;}OAe>q*;s~A7Q2II3V^wU^d0Cz(MaR+qg5Xy~9QwJKa%>#_mr-xUJW!A}
zj|GLEO5qIY=MN=*Iakci#gNUwx=7<SEaFP;mvRIi75;A#_Rx=lsNM8<n<@AQ#fI$>
z!9Q2vdp^ejcJ77jVu8y+V`q)PJq#DZ22wBhyF`2rh`JDP@_(s|1JeKN1%Hd!;JzyO
zw~eF!+k$_!7%z6#2bu4U!~b*O6c6cscnY{pv}GiJ6#83)gXNiHp1_Ac&$-#T9<s-P
zUujez^19&Pls2EJ==fOiZ+JKp2-eaQ7H13I?cpp^YFk5Xp)kH$Qmc%<48$91^5Rt-
zL4O-wvr*UVUEA2a#@Fbr3pKa4d41t_rM8(K_23zdb6(N4uCbZuRUHAZueH^;!HWkx
ztsC%y4qubsTNiF>+5izl$qUX<3ajDKh}T=b)V-w4TXs_^y^AE3w@&evuDHp)q~eCu
z8uscEKzJdJw@ecnw9=)Dcd_T@61T^Db9wo)vMO(tyTnsQTG`7)YTLrhYLtzy9w{!C
zZxx}}gV1Y6MiF|?$fvzZMLc-&);Bi$LS8)IVsAz<$bVT0`P<O>V4XLFr(s?xu%@|P
z@%rn0co^17l}rJDL5e{dty67tlQ$d;;F(qvnV#}al~Lw73V1@^Ah>*uTxpRv)a-T9
z8(T)%X3bZUkKyQ5E{2*fOGoiH!_Zq-47KD5FB!~V%7Q}5TVxv9dtwayvx&W}+1r2)
z&`W13ZiZl8Am|OZ`Rf!OIJM!Xwr~@Fv5Y~KHI*DPiXQROm{cC&nozUSj>v5Y)Q4uz
z;_eRchmtfPUgl}<nzlA!x41$SBM<jJd(Dj3@EV)Z<|r~mbHc%u6%v3K1%?;mjEYd`
zGzHp2;Z{F;9ZqVE*Na#3$e^cp^^Exf9_}S_QR?e)Qgs<(N`|0;VZNd`+cq?Xd}}~M
zt(-Q<9NyUDZ*5VWc=qadt__BrEv?Nh{?^b21GFX_z$>l-bpmpiRLl(d)+&s>!PnNH
zIO{e9(FCWVR$fDoh4CmhrQk(ftH04l0wUMa7*d>Uh@Hqf*EXX_nKpl|;tcuQQKXp$
z+|BH2arzs?+|p17ZL-95@JWVC$UK_zH3guIiz5&eCnmTiOj=`1y?+x#%N4dt*n9FM
z?`I{I=P#m~@ya|=S{CT2L(0qhUP*1JGbVY?Au-zb8Re_7eo&Is&Wp)UxqeE39)=s`
z<$biIHpwVxy8Lax@p!jc+Aq(QB$e%p<CWw&l~n6PhTm&(MBbN6I!Mk%A{!|y`!COb
z0>G%tvVD2~E~z}X0t-t*oh&MMAxCFIQeNKYOX?7u(tgP&=|<Gk8Ivr_^8iV^guL{z
zv|qr=PGrcZQeK`LNE#K4vVAEp`~Q%TFA;q5d_huq?nd$?lP>=~VB{02ADcg(1Qpvc
z0diZD1!Nw@zfs;M<Rl$7G9#a+Uind)ygb)P(!U_TCrw_Se@Hs`7v$+a%V>XclO;A=
z;j{xzX}?taRGPdzPmxrf-%9@U>*r}9FUMcrXGq#W#u-w2`@csWu}FD&J|pRIgDhp!
z+dl%#C@;@(B+X6}kntnczJ&~a-oUOw!sLCte4k=Ezbs4oXUM1qd3kQFi5o)6N4hD=
zvZQB$QS6FoFL}No-;XIBL)RZADJT8-AHb;pN%@<J3%8StIHmY9CMhrJBnTVjllP&;
zTuQbhbx1i`&VevtvVD2KVJ+qy>HO4|G3AN&i;N)hQvOS6xRf_i89PSS;$mJ}Ee_o2
twvyVA_Q?5%s-^8*f`2|n3k%ZpwD^(j@>;sL8(aP<al?2?nm}4n`ESK%X`TQ8

literal 0
HcmV?d00001

diff --git a/examples/04_sample_from_cdf/example.c b/examples/05_sample_from_cdf_beta/example.c
similarity index 68%
rename from examples/04_sample_from_cdf/example.c
rename to examples/05_sample_from_cdf_beta/example.c
index c99986b..87797a6 100644
--- a/examples/04_sample_from_cdf/example.c
+++ b/examples/05_sample_from_cdf_beta/example.c
@@ -1,4 +1,4 @@
-#include <math.h> // erf, sqrt
+#include <math.h>
 #include <stdint.h>
 #include <stdio.h>
 #include <stdlib.h>
@@ -9,36 +9,7 @@
 #define STOP_BETA 1.0e-8
 #define TINY_BETA 1.0e-30
 
-// Example cdf
-float cdf_uniform_0_1(float x)
-{
-    if (x < 0) {
-        return 0;
-    } else if (x > 1) {
-        return 1;
-    } else {
-        return x;
-    }
-}
-
-float cdf_squared_0_1(float x)
-{
-    if (x < 0) {
-        return 0;
-    } else if (x > 1) {
-        return 1;
-    } else {
-        return x * x;
-    }
-}
-
-float cdf_normal_0_1(float x)
-{
-    float mean = 0;
-    float std = 1;
-    return 0.5 * (1 + erf((x - mean) / (std * sqrt(2)))); // erf from math.h
-}
-
+// Incomplete beta function
 struct box incbeta(float a, float b, float x)
 {
     // Descended from <https://github.com/codeplea/incbeta/blob/master/incbeta.c>,
@@ -149,16 +120,6 @@ struct box cdf_beta(float x)
 }
 
 // Some testers
-void test_inverse_cdf_float(char* cdf_name, float cdf_float(float))
-{
-    struct box result = inverse_cdf_float(cdf_float, 0.5);
-    if (result.empty) {
-        printf("Inverse for %s not calculated\n", cdf_name);
-        exit(1);
-    } else {
-        printf("Inverse of %s at %f is: %f\n", cdf_name, 0.5, result.content);
-    }
-}
 void test_inverse_cdf_box(char* cdf_name, struct box cdf_box(float))
 {
     struct box result = inverse_cdf_box(cdf_box, 0.5);
@@ -170,23 +131,6 @@ void test_inverse_cdf_box(char* cdf_name, struct box cdf_box(float))
     }
 }
 
-void test_and_time_sampler_float(char* cdf_name, float cdf_float(float), uint32_t* seed)
-{
-    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);
-        if (sample.empty) {
-            printf("Error in sampler function for %s", cdf_name);
-        } else {
-            // printf("%f\n", sample.content);
-        }
-    }
-    clock_t end = clock();
-    float time_spent = (float)(end - begin) / CLOCKS_PER_SEC;
-    printf("Time spent: %f\n", time_spent);
-}
-
 void test_and_time_sampler_box(char* cdf_name, struct box cdf_box(float), uint32_t* seed)
 {
     printf("\nGetting some samples from %s:\n", cdf_name);
@@ -206,44 +150,17 @@ void test_and_time_sampler_box(char* cdf_name, struct box cdf_box(float), uint32
 
 int main()
 {
-    // Test inverse cdf float
-    test_inverse_cdf_float("cdf_uniform_0_1", cdf_uniform_0_1);
-    test_inverse_cdf_float("cdf_squared_0_1", cdf_squared_0_1);
-    test_inverse_cdf_float("cdf_normal_0_1", cdf_normal_0_1);
-
     // Test inverse cdf box
     test_inverse_cdf_box("cdf_beta", cdf_beta);
 
-    // Testing samplers
-    // set randomness seed
+    // Test box sampler
     uint32_t* seed = malloc(sizeof(uint32_t));
     *seed = 1000; // xorshift can't start with 0
-
-    // Test float sampler
-    test_and_time_sampler_float("cdf_uniform_0_1", cdf_uniform_0_1, seed);
-    test_and_time_sampler_float("cdf_squared_0_1", cdf_squared_0_1, seed);
-    test_and_time_sampler_float("cdf_normal_0_1", cdf_normal_0_1, seed);
-
-    // Get some normal samples using a previous approach
-    printf("\nGetting some samples from unit_normal\n");
-
-    clock_t begin_2 = clock();
-
-    for (int i = 0; i < NUM_SAMPLES; i++) {
-        float normal_sample = unit_normal(seed);
-        // printf("%f\n", normal_sample);
-    }
-
-    clock_t end_2 = clock();
-    float time_spent_2 = (float)(end_2 - begin_2) / CLOCKS_PER_SEC;
-    printf("Time spent: %f\n", time_spent_2);
-
-    // Test box sampler
     test_and_time_sampler_box("cdf_beta", cdf_beta, seed);
     // Ok, this is slower than python!!
     // Partly this is because I am using a more general algorithm,
     //   which applies to any cdf
-    // But I am also using really anal convergence conditions.
+    // But I am also using absurdly precise convergence conditions.
     // This could be optimized.
 
     free(seed);
diff --git a/examples/05_sample_from_cdf_beta/makefile b/examples/05_sample_from_cdf_beta/makefile
new file mode 100644
index 0000000..69adb6d
--- /dev/null
+++ b/examples/05_sample_from_cdf_beta/makefile
@@ -0,0 +1,57 @@
+# Interface: 
+#   make
+#   make build
+#   make format
+#   make run
+
+# Compiler
+CC=gcc # required for nested functions
+# CC=tcc # <= faster compilation
+
+# Main file
+SRC=example.c ../../squiggle.c
+OUTPUT=./example
+
+## Dependencies
+MATH=-lm
+DEPENDENCIES=$(MATH)
+# OPENMP=-fopenmp
+
+## Flags
+DEBUG= #'-g'
+STANDARD=-std=gnu99 ## allows for nested functions.
+EXTENSIONS= #-fnested-functions
+WARNINGS=-Wall
+OPTIMIZED=-O3#-Ofast
+CFLAGS=$(DEBUG) $(STANDARD) $(EXTENSIONS) $(WARNINGS) $(OPTIMIZED)
+
+## Formatter
+STYLE_BLUEPRINT=webkit
+FORMATTER=clang-format -i -style=$(STYLE_BLUEPRINT)
+
+## make build
+build: $(SRC)
+	# gcc -std=gnu99 example.c -lm -o example
+	$(CC) $(CFLAGS) $(SRC) $(DEPENDENCIES) -o $(OUTPUT)
+
+format: $(SRC)
+	$(FORMATTER) $(SRC)
+
+run: $(SRC) $(OUTPUT)
+	./$(OUTPUT) && echo
+
+time-linux: 
+	@echo "Requires /bin/time, found on GNU/Linux systems" && echo
+	
+	@echo "Running 100x and taking avg time $(OUTPUT)"
+	@t=$$(/usr/bin/time -f "%e" -p bash -c 'for i in {1..100}; do $(OUTPUT); done' 2>&1 >/dev/null | grep real | awk '{print $$2}' ); echo "scale=2; 1000 * $$t / 100" | bc | sed "s|^|Time using 1 thread: |" | sed 's|$$|ms|' && echo
+
+## Profiling
+
+profile-linux: 
+	echo "Requires perf, which depends on the kernel version, and might be in linux-tools package or similar"
+	echo "Must be run as sudo"
+	$(CC) $(SRC) $(MATH) -o $(OUTPUT)
+	sudo perf record $(OUTPUT)
+	sudo perf report
+	rm perf.data