From 65a104e089de7960d5ef9a2d3fb0906ac8120010 Mon Sep 17 00:00:00 2001 From: NunoSempere Date: Fri, 17 Nov 2023 20:48:23 +0000 Subject: [PATCH] add 100 sorted lognormal samples --- examples/01_one_sample/example | Bin 22648 -> 22640 bytes .../02_many_samples_time_to_botec/example | Bin 22648 -> 22640 bytes examples/03_gcc_nested_function/example | Bin 22712 -> 22696 bytes examples/04_sample_from_cdf_simple/example | Bin 22792 -> 26880 bytes examples/05_sample_from_cdf_beta/example | Bin 26872 -> 26864 bytes examples/06_gamma_beta/example | Bin 22552 -> 22544 bytes examples/16_100_lognormal_samples/example | Bin 0 -> 22496 bytes examples/16_100_lognormal_samples/example.c | 17 ++++++ examples/16_100_lognormal_samples/makefile | 53 ++++++++++++++++++ .../16_100_lognormal_samples/run-sorted.sh | 1 + 10 files changed, 71 insertions(+) create mode 100755 examples/16_100_lognormal_samples/example create mode 100644 examples/16_100_lognormal_samples/example.c create mode 100644 examples/16_100_lognormal_samples/makefile create mode 100644 examples/16_100_lognormal_samples/run-sorted.sh diff --git a/examples/01_one_sample/example b/examples/01_one_sample/example index 84f385ec57f2cac79efb293bea31786b15938706..9c9a4ae62dd96b3f3c6bec7bb7816d4c12971889 100755 GIT binary patch delta 1817 zcmZWqZ%k8H6uXOQmp(BK~L9n*COiQ08YpgE=;&+5A!Y6*1;NvEDaIe*tex${8wyTd2$ ze17)8ku$aG2a`66ri(U0l%DuSHdV7uPlG5qdFfF--=EJV!EWv7rc{bS`x1l{E zTKsNfL|kf!z`f*N>V?my9;6?^Tq#pnAu|iu;J<RlEcGOnVJ{ z9*MTWQqSGIdP>MI4o%ZUr=JaRnR}hm^}yJ)oZ*K5<6W@Y+eAm~MqbqY8!(oK6c?^-f ziexAtdio>cvL7PiJMpvQ*z{`7vq4S6Sn0Z8d22v)Zi~K(hvKwHnA?MFd|_0X(BXHd z;K!H%vs=QCW7#cOV;Q6L`>MIbs)T>id%BP0Vz&{KW2)NWfVg>B>BgG1S1}AJ{Ua{L z6pEf?ATTRdTyz!^)_||Ntf%0flMvsMlI|N0UsXx>m~XMK`$?73lYI^G$};cYY@Xqc z5%p-+7GZDhAcS&V^c;-l<{QtgT48##;?X?Slb1>za~^Bjpn;Gh=)>qePc)4@L&yR2 zirO%HLtJeV9=#0@>%#X_&BoATb~T!fBo)Jt5pP1giQx;(mM)>%>=H`N7N6P3zWWyW zUgTRDub6X@GWs#R!19elTK;_F6s)R@YZ{CeCaRsLV+JTJP=vv3b!WknIOa!uW`S;l z(bvA9)v(9!6#BB^xWAbmhv^IUqX5n?zLH$GYL{dDw#N4bBGn3E$`JT{wp9Z=A*ZjoMAlMp$o# zDfqCga7;*aI9P;J9_CjX&vh9w1)T%yOC;IPb@lv7tme8(I8-7T*2dIw2iFnaRct2R zKuab0Q_g;wp?`i8Z;9DgkXI@hg2<+^uIOcG!2Ru5R{E>Dr&JIM8e5v%m$$ant!Y`& z+|s(XVNKn-hSr9)@71-n)HSXoF}iUj1D+PCd&}-q*!4!6YFqR(RmaOmCJEo!;nw0= olLz?);;+P7b&5iau}-4H|;LfThu_V$;ZEtikcb!0%G^*!Tq4_g7=^B=?^C`+m;( zo%4I#g%Ma8fitNvHsJi#pi5wDro$-6G<0pyMqRQW>_Uw0v6{u+#-B>Izihnlnh{4o zdZp+O`Ngm!Cc~c+WBGlNJ(8?MuK@grD=`ij#zadC?8JUcx?r{8POD@LKYVD(WcDWO zmrRwqRh+YofeYV?eH##;5YvPfNsYpw_y81RpLi4Y;*q#CAv5k@*^F1>W}3gvH>u@moN|-{VW*Lv$yM!Zmq0VK5wa%cr;82I$A0KojuuKU0K3^qC&;12%a90tF9d}}Z0E5825WDhci z@1%L452w>o9uG;3-7+Pvddds?a^)n(Jp25ou6oM+rTYA^7}GAAbD@0^KZ52 ze9@fTlW|PIpkeaWtOXsSNBNgIdy6a-DLPS$icja9Q?+moW?;53Wwk!lV; zTaHTUz}0sHPElk#MVcA?8pT><=-q?h6QzElQG&$u>*$;4B6@*WX|-!vkMQJAG2jzj zmpKhOHF<>7t6V>$ZZ#RlY)0+Md7D#nx<=Q|fIdxO+IoXruf=mdJ3Ner{7lKq+4NQ| z=c49V6_(_S(m_tmydJG}&ARXy(QZzw*L8CoKPLJ*rz3oY-8}bwoFO{QDIL(7UE_2V z7is<95se}=k)_G>XRX;(D2uNobfUkhS<+vePV@7Q@1tLUan_rX?>7uE3qMF~8 zMjn%a=lml0)JUn4$F#1i(uPws!wZ}iCu*sAe1;PUlsm-f0`DQTWACAh=og%x<+&w1 z_sUv_Kg-<(f*_9<&VjAi&jL;Lwa>N(syZ4D)OA)hx3nJ&G_j5Y%{cXx1xNSnj}7sf zKvQj1Q%k)TA{&dQL0%}H4ilEG_MqK@bp;=9VcC^-aaRy_%?4|RVry|94CW#rL>4lPtET*s7@Dla z#2JIhKX9yU(h)z@*!aQK=Ok)0AdMy3M2=&jlY%+5jWi(!!ov_j(ubekcX&p#OP8xvvqz3vU+)T>iPFi5EMKs2}(K5lC z{0{vD|B|i?o}Ax9FT&XTLo^dQgiLOh$W$;S%%KtZSs0^tVK6z9`$qT=J_=TG)O>9| zmh%ox@N?Bv3B6#l1Zf%6Sh8pi?6G9h*>Ki!)UY*)5G<(M3Y$`jI3*EIr=%HgB@z

I`IQHamE&+lq}q$s|gENIQxPe&nfVGkwn-NG}oTl6t) zGjzHH+6O-%*qh|L4W~eATxhzJ35lSM1P74j_!p zV~FGxB!@!0tAB!DGdRJ&ns8+Xn_lI4Hm9#JmfF5Hy*$9%c1B;teLn30=5`^QFcMYH zZS^~naWf{s?54=ESauWESi>mwzA7)VDv|fqp6(!7=FnquOjR`!;b@`Ae0>|&fD?{Yrt1gygLtUg!on$wqHy2l^3>8`pSIm56h+9i+3VkTI~LZ z%`=iCKgwFdb$NQA&ts>(FzU(G_rA8?u&U(2b8@mbm0Az{ttjI`go2IGcl2HRum@D{xEcF zxLr{3W%%Hj<;HNKnLVW0s)jA&K*@S>S5RUMfGZ$NOz;iF^+ zA9Aa3SYs14d_+5-8_oo5-X$6u)W#IOUfi)l!@ISyuG6x-0{sDjZh_%Ij<`e9crsP; zoTfPh3nW22rm<_adU)2TeE}uNp4QlKtQ(Cz5APy-O=F2&h2PZJQ8wXz5&;n8X&Xjw913c3uo6bhnW(}nT>mDwiMky1EPDCp{9YOz(* z1+}@yX3`F{NDx2J*zpW?`%idF%svHgk)R7Bo6fqTm!Sr~Z^W|FyK-j{$K};F)VHi_ zYOdbgu)e;bsjg;obz@CaP2KC&%?;JH8%d0=-N=B)MRHg19SVD1Y?du6Z&GK!|ghJ6RkSr{%qoeJ}07nHAGUy135n&)x$#gJW zW*|d>L!xeO>ih@ykKG;#6_SN)(<~+mEQ`sOsUXRegcyFv=BC@*dA!?YVdsyY_r9NV z?m73od%GBb#{uZigy2Bd-AJE>t(nh)agb~1kI-KGWe-?|1f4gk_fS{nhhr^A^9EBU zuZ~WQY?qIOT{jwjlNfXNMRiJYS^Q&wchGIhfHG`0wZL7xW6BZ6%y`l)Me7b7TsFoO zr|XxEyTU)lZHYmc#)ZUBVH@^}*@970qi{-m4i@pYxBx47Jt(r_aACh?1P;rt)FBD_P)cVrB~ zL%g4HH2HiaV+Z4?#EYNN({5wlit%{+--_Z4MquZrA(+R6%|X!Pjm=k{eNP{XJ%R(7 zBF*Pg=4}2;z}UP+Q52)KSRGE8sQA4Z@CQf@{e1(L0|pN zw?=ChWB4%J2}Nkn*;E*k7+W!>{^_hL?JKHJW6ZhJ({sz|^YpAbE1f-0eeOP6CCP7; zX8olWeW1uUZCeHGGE5EHgZd;}apjW@1_uW1X)umI*$d>P{1mXbpD4;gIAg=)bAW#?aL(!&tEMp<=vg%yI}4CjNm`XEv_z-XN9&my`#tGr{}W52 zV84)ddKJ9|+x6*GZGtc^IpXCm!=-R^6uAYrP3|q)sZ;$uXS4tx;nFMTK*E7ytKhKV zSH;cHjvHL3C;-P?B@R(7%2**yR<%FK$5Hh(wK1|#pMs$2fpkrG(7S3TY<*RfGC1v` zqorB~U#(E3;92`^#A%9bqsSOWH=|sO4E=T@l!(#^qOpR+^y}y+m{0TzUS-0nWxc?Y zN6}v*=3nDk+4QIa>JIH=svjn6}?3t~cO?5-Z(EgF8F7(S}L$WoiUJ z!WJ{wYxa>%&9W@PWEtQ!k^PcB< z-jDO1bA|?K-yrRA&{%Ko*t`STBt8ygLppWq2IeZ;0jYvoGz)$jU$H#*vy)qLs`{srGJjM78!T!u^YvJq|@4f^Dr zJEl0k)TW){Hl!>E$QYu#pi($We}vxzr{;7fQ$b66X=JBONn5U4~YeE12dEwLg+&-;fK&79XL#(4I3y2O!_^`D2T< z39St#9Rh6!UG5FfW=hDfW?7awODXH&Mq8`cHKA>6V0g;^@h-^9YotEdpVv=&r6Ol1 zh4rqbG!=SWw!Le^BDd!pn?}Upm=%r4IgG}~qPJu@%-y~?&DC89bCD73T%}DL3m}Ze zqlx5IBuB%XV|a?&`0W(8DYY|+H81cL>r(`6vFmI8r4i1uKmI7Dx$FlR+ktFqPh6R{ zx73!2&KLlr^Rcao=sf1w$SCEy0+(2pSfx_aEhN|3v=|&i)s2U_J>z0G=B&Gl=BPM) z(<#G(LyUd3zLOj?b%nh=+;y}cwTaNOw=6vN0#rT5E7V?WquAJ*Ku{= zDn69uRw=aPcyb*+9r9r^A$|CKYFQhidTYI2OH62nL16%nEMiqqN6(=<5Y&-ZwX$@Wl7YjZ=;<9=#Mu_v)>2}O1>z?%TV4LjMs zNSeE9nwHc@BRd+}+Jnu}Es68fpu*Y!9|agY``$LDx4iV75s5yYe!H@fWv3 z^o2Q6)9Q;h^Z`IU3nVMD-+PrWyg F^DkU~=BWSx delta 1577 zcmZWoeN0tyW1gHyNjxRT4!)!88 zV<16+Lo%1l!PJ?W7y+M))g=>MoFpd8W)}TN-Iz-RgJT1kM263Md9TLc`Qx5@e&=`Z zx%YR^^-n-x0tS*`c`*5A=mjHLvkOLygf#8wMmZc6Jz(OtBqocwLTSU5u=#>-IvxzY z-*@A$p6!1C{EX?A9@vQ=TZZ8* z7VbC^R~kx4M+OF_k}I$gZ|)d|)7YLEfB@Q(KHHS0339D?KZ%DD z)Fw|m9s@#VjglnkOnFShXA|1ougqy0au}}uKVF4#DGlJkS5wB|w74s^2Q-ymUDlaS z5`x867q@x`2d$RfL5z?^-PXTc6-E8oD<(o*2RyysyUIPiORoK{-iPJxe%l;{_Z20t zvZ^y9x5f0d>wuHmYeV)aO{^_%|H>oU9GtS7;Udo3GsPrFJQ&|vk)$~_AI7O%FP_UhrN&84k6GDn&TkWYva3|g-*&Zda0birMxYw+=b5-s8%E~0 zz~|VVKSl}YEGWbe3Uq>-BP5d6(_2`1LUL%#>=}$G;Dtj<_7?3LW)}Je%wA?$+Q@AD z2pxsIF+iDmLYux+W(3EtDQ?CIE~@JWokeg0(X3-kGn}4**yT0!ifeX{;_PH*vmEb| z9a#V`gC?GT)ve$K$di@iLBo7DhlKmbv-lJqmNuDpS^!r_=J_4mKIDMR1#- z6C2&>LZ6ZwNR@LJWk=5Agqs&GDQ6mplHqhEE8?`9_gzxfmaeP!or2exWqPNw>L-fp ziGk?XDK)F`%~ULAjmm4>Ww&x@sBlcjDMd&6xz=nG#hH%dO7&pxnehb=5015QFP`)8 zLW#oE^YY3=$_`a{oidLxrqHKw%!)?A-QGd)^1_Eome4Q6EMImcfDSLOW*xna^#a^M zpO=Tlwb~&#=(RaAc9SE%=C-EVYDd~m=cCSG)a+~NWI$FBsH;87;tX-#y9lEBiL1bu zdWf&7zV^ApzUm_ljdktSEv<)}eNFg#aZ^kXKj@=dH?`KwDY2r2{b_E4m{fKvgnKOx qL(9`b{mPx__sA{30$HN=fCO+7wUwFLNw0hr9zB(r@B?;NCjJ9%%7(fC diff --git a/examples/04_sample_from_cdf_simple/example b/examples/04_sample_from_cdf_simple/example index a2f9d8571ea28631192550c6023542c58f525b11..40d44cf06d4369687b700d5dbc7c7aa20ee3bc6c 100755 GIT binary patch delta 4990 zcmZu#3v^V~6}|6G1}Blsy!@CXB$-JvPA0D(h(MHw-Gz zWL58(|3eImsvJ>~)(VLUnyFVa_0Fe8hD$5 z-(ldj1Xi?toF9vi&ptt+)M~GxLz02-H}J`!FJc|Ts85V}RvPY@%9#JBeVwA*?_cw= zWw+8tuFL^FD38jW8Y+_Fh2Rv%Y>rySF4?lKf{u@p%~A4oWWLe<+bH>blzci${ya(! z9*QCkM9F`Nl0S@+--(i6kCJ~ECBOJXdD5uojf8Kee`945GSk1YdSxkNz8N#8``4Va zM1HPtm8q&~DO8cwDpjXa?WNR}m?4zXi-|d?`w~;)xlu?p+E79#6K9PmVXUz*GU&It z`*!eSF=+ddd2b0#Ng6x)ml=G}c!vVj+rO7(Ic0n!W0!j{m3B1z8_f$j+s ztC4mjWmptN<>*E_mXskh(055$g$;20NYPCHn(BkR>F++>Io$tD^FB)0oG zH?2r1NOnUvtG@et-%Rf5G8C)a^hV11;ZR_xu^=9ChiXS831W$xTI??I7B>a$%fw$5 z(gk~wxT26=%1E=ttb7y`sC8KeZ7t?|3aK{rPhxigT~B>MJeyBj(%j;nd~QhLwBh3R zeEKHstoVmKI-2el*X9i=mgG_C$Z*onMq%vo{y|%U>M;Jii?TDOipisBMMjZ$!G#WP zai5F!Wo)+R;$%G8g*6V^=EaY8ag@xXU9>tgU&y6bGPA{xa%oRymhB~2Ms{7|JKAF2 zoJ*H8UE=CoO3xZCelhQk*mPy`DaX?(ocE-VMvL;FLwzNGby-syV~-|cw8|-HAy+Qh)^;LU zeLIQ~VkGeMsf?8$Z0tURcx_3g-xfR;YgfF*|2maM748u1bi8m!T`5S_Pcd5XGiY18 z4SxI^=eaugiCfYMe>R7gy$@vXUFP1qPn7pHZj5gpfn&7MTr+}>x<4upK=V~5Qd-Hx z8~qn^-h4B^hV~UL5vQe4W$_#0%gLeB#ZL${tn_x_C*&=41*f<9ZIAcyYPWx9#EHhT zKFmjooH!}V{Cr%p&Dq;#TecfJ+32!F_JUS6-vZ@9DEr%N1v`-ff5+jDUynbQh|w+{ z$8T_(Bax-KyKRoAuqO4+h3^7r<2Nam3Bg)dDmudf4{zkC?j4bq=5dvpiOu`?D3K4; zIh}%XhRXuya5X2~W_vvB@5Y#O4x-ub-+9vMkL_1hvb(KRADr(kKl8rXjT@^U#?8b} z-{B13%=+$ceKUOB*Kie<+zWnsW8rlkk%)t=WrLxMvC^=i%Ji0Joj^V1sltc!YI#-g z!~5@xy=U6hal%N`owW!+K33`dhb*r`yXl{@ycHM#^4l-aFU$P6^bE){4`SCrS?&jJ z1)c(ifmeXizhqg0T`#bL<0rD*2Mhy00ZNB(rvbf(;?Kc|Aj2y3v`~6 zb75Mu0)56>Gg9orquGI>E40ht%3Z7x(^uld6;M)#w zxDS-Vxvlu#23M&tjPd?&;|qhsrHgE}N(-zO=0tNpzANB21K9+tv_YI`B}+w$FhE%q z4Z;Dszhc}p1TvfJQMmse%pnFrEL0ekdj*^`isQ*}7r2F*9c~b|AKa#cbhTn+3c~(? zL30fn%O~Wn9GQtYkA;MJdc+fPR!jY@%n6YvNz-9kSeYw)N>5jo3IC>dD;>g7I$Y^! zLChCx-uN7k!FDgW3Ilf$++}c&X&j$V4DupyOqTED0Rq|Oayq^pDF!U?9o^Qt_tECpryTJR;V9AtS zkwhJ1i!D!q=>S9j@h0IIx5tpD+Mih9;T;*novswmd|F@a5a!cs)j^>((4#>nXQDNO&hWYn722&VlC znrq~liP@R{bIhnPu~g>Op2$DsvDSc*QQIb#&ZGdnKmPY=nwCeLmnM9IGCcjbri|}R z%>?lSKMmHD>A%+$J5>C%2>)k7hZKXGR~AwDn|J8Hil3Voz;7*)|3Kp@j6z#$CDG}j zw`!k`JA>?X@<0Y?d|j&$plx*%!~zd}T{qFI+$D@DS%Pk*FinJyKD<rug zIw*$>{9N))u&42>Qt6e)`ykq|=4BN<1noMlZx8;g^D;`?zAaoU{86d*(ZB?I;&J5Yw0?pzBn$YBLM)BY~n{c3IrskeE1RN-5Fo z`?FP9zXa~k`COGR#>l+Xt9&90<4Ma)uL`V%c22aD)1O+tQtKMfwp2ObuzFVGd$kZO z)_QKEY`;D2T}|78hZisUhVP+iey0>rxS{F zFZvanr130jdifeJX|Z@r?Z!rfmy{Es5gOdB*?V!?d2wq|u2FbK6YCRjKz#LfVHvN} zy1?^@MXzN84b}dnq7`;wq2v(;8`pWV}eNeQy&PuHubm=dV1PP6E5@inPqYMWr&luZDv)hSDS>F zPD1jjQrJg1v#O*TRno7$P1HWi6RZ7a$1|C>&8pI+v6JcK4aqy1a%NY>YL}UJDlLQ* z&<1LvZL_QN?Undlg<_$B@kw-T_9%-+6e+}A%|iqRZ?1|xGDo#*SQs(gqKNYW|BbOS6p{)~>d>G&p>xWF=Y zw~ptF)mo2^=aXw~i;mZ##C82TK0!0)^FOXrBwkQ3woAt+>-ZORe2R`A)bS=AKcwT; ziSDz|V%59%+fK|tj==ij-uQd0H-Zrm%|OXbr* zk4w{YtHVX+WFhn`#!PNMW2a5oXF(^&$c7mCJThKy|3i#?GDbccBYzbmhdzrTK8lh5 z5hK49BmX5vekn%&V~o7>Qu(^+kLd|tW1z3ximVLuRc*><%-7JkG;se3WAqW3E6-3= z>!6BmR;U_ELlP}?X*n8eb{14=0;(~Xq$_USeB)`7^^Lajt5MI!$bT?jGI;( z_POW=Yvt^BGx?;+^)ghw@UtXIwrl4zHqnoz+=5P^&#;QaHD|hQfvtw@`7~tBGRlhb z+4JduHA|RBKUuR2=E3ouMU8>`t3Ku3K=0A+q`=my{Scx*O-fSqX83;B@X_* z!9sP}2gKDDDspEg;N)|Q4Hl|*mkFtKpWBW5fcwU5S(93uXU3lFiC%x6yE4n7*>%!P zsh$O5WhyQ6ypd(%K2nqNqO%e7K<}vbC}zVarf_l2n{nc^N#Wyp4T6wJ-uxEacjnJc zK(T{E+c$|6$^VUzK&xjwiu;)vn@gHdAn&q3S-uQC$sDT-%r2gkoRYXIO9reOdT+**s{DxscLwCSDxx<-FmS{7U+u zXsy_6pz>KSi!a89kIuSJFgCrj*W~VDjD9Zm3Ad5=%KS==Qc9z5d2L_S0$vf2qx*f0 zi+g|YHTZfbQ7?+;fL>Z#a9%DHqja!j{8jH3At$5Q-*7G`!Heaijq$U{+pxkryl>Ee zuOcL5#?A9%sWYJR9+afJ(GMPyq+P%^Ag>c`5lP~S)qhx$+93BFk)#peF5uU|2=ENh z{9j2j!)_2bhvV0hGz^RY{|z)B#pwrnfsIEylW%hNYX4|@LNgp0{f$q)B+qi zh6Mq;j!V*$T>f5?J_LG3CFumP2bhG9iZP%ESUQG@frG#R=TAw}%|OqDBn<*1z!O~l z5%&C5es+N*^`kR*R+4rD+s;YS=Rgm#H3Qq+11tsh1O31;U<;5T<-34gz#+<6Fw=n$ zm4IjdDe{{%58XMqh7ZuH1v5flAuwF4HvY-fjcXTphQMi!+ib)20=&6Zl5@Lo?FLt_ za(w)5ToG`EoC~Cw+f$6pgYJj8&Vb(mWD8Tw+r&jFWc1sF!<6m6UigY`^;a*!M&@um z8u#CiHN;~R3uH#&o&)EJ;dnC42e(SK%K&}|?(rjZ)}Lv^w*OkEnS{ppH5Jxmy0GDu zkg(2D@laxlaq(2hH&pT9ae2oWSyARwP z9XAGU0^B_+$JY~&yl6WiN&H;0IG2#(X+~r)H7*Zy51gbM7G5um(vF4maNoc1D&Z6* zF4`bW&{KB0C??#jc*A<=J+pT zn*n~C&ejjpS>U#S>!2LJm3jiRjNM=c!O*FoH9VGZPdt^@1uRW+Rx#$_{?SV7=;1oI zfDiq;kWe2kyumI&ctuR7rutWeU^p?jUvOxHlp)GJj^U%34ll8VtW5hexKy9UW?rRd zv~e~}dO1FwY5xLxZ5*>R?SJG_K7)DrcVe-gML&)|TQG){c_%Xm>Gd0*cc@yvqr7oZ zV0?PIcWDWa#nOf1pDSs6X^HmNBHQ7{PlIro6@4Zf{PM}7GXKj6os<2jPaxYUx0T2! zk^iRfoIs;p4Q4S|La#PFnD`>H(<46|UEA0y1nHT^g0itZt`#C@{=Lw=c6h zeg|h3cpA=^eK9W_@HUqQpuG+f*Y?{LcDfplo48h3FW>t}SZ*)hsaZ}_Eajts?(^Wy zmx=an@E+ulHMAGKRMEk-pEqSZ3NLs*cblB%jA_Lw!tHrtl`57v6mlC??MHGHSv$16 z>T^FMxe7lEGH+9rC%3RjkwOeAz%uF7a=VaEOPcHj*Qnv?x>8ZfM-trzm5-==2MJAX z2Y+wysk-p}#+w%EUb?TzW9BKOYXkI3lil&AI)th5533e0(cvb0=zUc?oUYJXvHVOM z6|<~kDsP^q*lR^<0=${=+}16CcM&|kV6XafTd3}@LFO54YECV4qtD~LMg6tsg^jmq z>ex1Q4O&&Ipxc}6LV({LEBILKZS1_A^hvYba3{2K4hWBM>#_mW<4BCh*R+|`o$XK- z!1hwv3cKTDRm;@ie6H^AA!>*AoTgRBE>XuO(spQ_sDJ$Vht(XDUu>`_2=7^uFNUYz zWESng#jJhZnvPW)+gtBgdrN0`>-r5FZ);yib6VF?aOD)!9jn{db+oSAutwp-Pqjvb za8G+APOwpXM@gbq7;y@BbW|j$S<6c43rK2ABW2|5tT3ykP%DXBsJ*i^L9K(ORkWkC zLX#@1=!=Vzw~Cx=DiV05L-SVC3P{>kduhi-$?T;s)>I2#a;~j5dS!DqNbPGY62`7m Rq$(PKRIHA6QqF4k{{c~<6wLqt diff --git a/examples/05_sample_from_cdf_beta/example b/examples/05_sample_from_cdf_beta/example index 5f5422beb6caac14f96ebed5908683dda4af1171..4d138455e4fc3a665906ad66cb6f575450bfffe7 100755 GIT binary patch delta 1758 zcmZWoZBSHI7(VB^Yj;ui?gDjpU2rdYS3!|o&}10XSVg@!T7DqxqPZ{(E2aX23{D@4 zic>=~r1)B+Wu^XTvS~`aGmT~xh*PG`1Wp+>lOpwF%`FfslPPcS*?YDhb$*=pz0dQ! zALl*i{@G7Q`ssme8tBZqk>XQHV);}^rFk~RW_~JqDfJ5J^FCf$w6SX}XMcLvF* zsyHfztruWU##!OM3cq9=rR(9P%vnOA%4E=&Igz^Hhs-$jz~w2kgnLu}g$E$hdsBFj zbx^6I`=DCMr4e{pnMLQrN6I_46AB?%zFKkZ=!7@ZlCq=FnJ3d`IGJ}dzlsuaE7LFx$y>svvq5dt+Tvne zEyJz<$J^k%zmE2R8u)-Z^{)eaDg2RNL^~i*pti1!YSON5w)JW~F$JCaZp`M0-MDTT zQEB*#yHe$)s8ruODmhoEV(b#a*f+XJUPZDmDrNVKN*gYZN*gDQ+1Z>z7uZM^qOY}G zbgb`{y!#T@eOH1-`gT^NiBFuNmmJDJ^qH8wEHnYh%@6;>torPienT;f-c$xwnc?v*VCJ9YgUgv zt`{#-X=dYyVQAwC8AMx$w(5>y>_*#%_L9D5@vtSFjz^DM2pPC#7!_;`4R)o*E|Ls% zEmlI5pD~o-OYDv|p~9{R%k7S!U1V`8@{PzhFd&r>T9_SLJzby~U@Lozz6~eKynC|-)y}APGn!@>~4>V`tVV*brO!3 z%k&J>(B6b{D5;!oD&})#9B74Xy^pK{r-d{sLRLsF^^7CLDY%JPWvY{VHp8)7QfwkP z>BK5ToSmWDXMWG+cnNe?$nLeA?em-PDbChFPlXIh$d$i`D^~FtB(4rsM>&3s&#ab5 zzQ?+x-B1+ryZgB&kZ+Rbxh4*KLb7|5pZfqd%@%KgjKS%U+-k)Gu)arlY_hJKk4IB^4bZcVYnb`M3e^_?E z&SOupw~K3wX~?@kD`k3-Dc#q&W)+`%az2AlT`9YBI6IPSh7RB?*os3ZvZ9rrT{Y3w zrXe!)VE!5Uca;l;1vL%zEwxR}k*y7z>KmH2Rd0J~3wpd_Diqcn_@-4#iQm8PrruC( z*WT1OYUZ&lYf1nI{+9JVwBU%jKzP*7WVm2X1}iR^*TIJKn+pUb=RRIWvt>m{$a%}! z42LmnErfP_!&(4U_=WYj?j$`1ZJ;@cUD;1S1%94wNf9YqJ$pk@SWhGVlgUZs}he_*+pt zQ;9fOcmz`A^rC)v;FGkN>?n>Gi|{0#C|+@G146Eu6-Cios@Zh**!#RMuWLKpjBEbS z_u_O(2V6kA^At$(Mdu;V2G4$_w+s=2*Igbsh|bd7tuaK%E&b-ZoKQM z?>!B^;om*=p5Z?myrV_c6z{Id`-?SQRphVBN|Vu$v~b3~3fUP&)%EwXb#4;1$ToE) zGb)$d(+HInvb)j_W@|)IsuKtqqU#>Ij?uM`u4F?|rsS8ZZp8Zx)LV}yad?``#35ONy2-Z3g98YyWELct{~-8F;~4r5&I^~7-`2zItgZO z>NE3Vs$Z0XJUV1jse^tOW?i`G7h!;D;MJ(|XiS-_D%h80Hi52ysF|iwE1nk8V4kUc zOVm#lB%*DOitazw7^mkUZWYmiMePsD*(qlns{b(O2{;=NrA={abfHAWlRWCcjerPF z3}#eE`=c40ag2CPp=NSOJ5$xIlV(enpo1}`EAfVUIOT1C4^c_=LK`-vhCsxNsXpO%2fpgC+u~kmjfzr{ z#}XA+#h=1?Q~Lq%`qX!!9(T#v!Xi8O!mvCU-1xJ69cJLEv}{3>|AR--X}^?glPOz% zI4UL?p%TYYas*%x);WBz4);2;VF?~}^jP~+8KZ>peYhdL3<~hmbeBygEYdfwX}pGH z{59PRw^i#j3E+ElIA_2$T2D#M2HcF8A(-47}#u zmL5)EtTmZ3Jn|h`K{5B1_)5z?P4oBrux;jku;Y|j{V*TDopt2F!5ZKLkk|v zyp+Ed82insX`1LN;w@b5*{pP2x2$jAc>MqI4m_5%9){5~`y;Tc=VrG9PGm2DO_=BN zw56iteZT+y_kcw_jMSaUKcHVskd$~~4yKfzQynM9r?G8_~$ z2czQZQ&DkEVnQ5mUiAsyr-?~aT*p_^Q$wQbC4DAti5Yh&b|$fjX}V`pTahP?4pRV% zE=5kqqDz!xHK&Z{sy@xLL=G7>T_Lj4W1-*_s%A7OZXZ=TDQC@DGQ-N?Ww&At>o*w+ z=ERbVuEmTO;IAt0D%|a2%>QI*=lN;=m8G5I{z`x6-IYq$oUMefD9`$b_cM~K-t)~B zUdtQ9zw+F09PRl9mgCQ?vo5c=vj_wEcb#kh*0jy>jD<;elMdX~v>ws}q{r08xr^dz zrqIP7r>h>I;i=+H2-)2AHVd0ZW+mYf!q;=W*d}!dRW`R!W|RCj3qRLF{QL<`Yv#OS z%OlECMCKkIZvvx_K4CdW=1FqFn~PJ_0qcc$oL{I2X*p_p;hH$EFZ*o*RN~Dg@6-3P zugE1xIoMb9GJJsxir*v$JX>6nIfo};OwU8wB`}uEB-ZQGVW!td=0Z$bD%;CrtjAy% z)1FH$`iTvh><;>3xOizMxwLZ@P+$A-1M>WwJSAqQt{EP&br{bO`y;VQf}L6JW3M8V z$g6IfSqA92xIAnq1~hWK@`33F;3GdSl8J}2jD4@6nA^lhdAmun@JNX(f1!y6%=_!_ zM(ipRKVbH&K|jM-BMz3xunjMie!{G~JH>!=P4+Xa;H-Ie z%sd;*dmeWZdx=4U;jyOa`aE=Ayq;wm(J#qQB)-Bh>~oV;MeM zD#O#<)81P6M;Rm&}Oco03 zLrtv>%`LSXL+hGC&70~r)^4tAuG{onZA+-OzLCY~`bG}i@1)k2j{q)Rwm}_UHUjD! wfyv@xpb@{ZcfmOBcJxCL{^B@c zxt`A0u}rG)=r!6k8cOdNPPhLpN&e9^>?rDoLd@Pc3O@XL?`JBCw^2s z=iLH~-L^`SWUvL~7Cvh~D!jU&ZSu*u=KuJ99CtRsG}=opL816tqan zbNS&pEOQraPC|^`G34L%*VXh^?kix-|8%J5nm-ciS@GBUdsZVtucwCKJvGkza?yE7 zob;4)IH8+*vushD=Lyu_E7bXLv8;gpI>TFmVehKdcu$gI8pe*2K1I5d^Z@Ba(zD{> z@@=WrMw*((QyDu)UDnZFmTdaMguP2xoeZ?-3l8eFO!jRe{OAApQYYAP!e&>Ht9U~D zIc>FnN)j&lvcw79Z4GX&6ga^nc2(|AmGw=(o`V1`J^4O3urFZad>*_U=!Q;wtm@zz2NlI` zs`#wjSrV0Izsl!yvJ~Y{54aUJ(Q?{aRD4Rs8}NM4MlVArl$ydSSLRY;#udeyahJfG zPO0_cSqfaLvM=K+GONZes12{zN2gDkgnHkq~qT)q-orDiyZ(;;4ShS>Eu#AU%nglZnH$Qj{F(Ob=PVt36X@Zk`P#tt<+-xjTZ zvFUJQd;O7?wilu?+_NK=o5T-BV-5AOmO~09e!62C#CUi%g+{lvwlrPvC^cznFV=nm f*TvR|1WuSek~P0fW+14g zvEkQ*ov?Jm`}P7=VNXalzW_=_bV);9zs20<#w?DyOslQ)wY zX>C9KXStuuoO{nb=iGD8IrqH#?#sOI6qn9TPfJr|F)80vsMUN0$GAj?txQw^T*`DM z7vDE3Hz*^}&R{s+@8STxyu6(r53>Z`1W0^&R7nD_=V(PjsUbn)OI6xN@*W8lD{m8@ znX00l2+iX-orJPKHeLbmR}x?^8308>BOi6SX0mR@wICk$CfG&t5t{hQ1YeoplTZT*{ubbE_#~DJAd>0%R0{lqDe(0v@S9TLrvY!nC$X#mBAGsADe%7lete=3 zY9k&SDEy`=o|5tfo*G|^Z&h7uz}Hf~U`~Bgqp#e%vfih7JgXX-8a=H6Z%e@A;iO3k zeQnK3ebXwX^?{avpM6u)uQWHUQCjO7X}D@JoTs^^t})>ERM+0G_}c0M7zFaG?8oc%OP1S^NG|)cK+SEeC)oosnzpl|+U-uv>NRowu4PK1oZ}It*d8H+@=XfSN zCp)JO)$bUp&WKf~#A>%Yr;uUN*_tt7#hLLL!CIO4PFLt#!FL8=8vF1y>jE}-Z5^%@ znaW|*X&rC+-alh`jZpeUJH*-<%2~9j5Pvfu(NY;xSC>nxpDvw*92zUT40uZnizIKrs|I|J0YBP+?=|3)d25>Cs-TV|P|J(;+%_GSGw33AY2NH!4J5 zOFe^;BhhGd8*9^sU@&qZ)~1cXV5BqFrj5#An!658Vu+G7*ixe0A^Li-CG2N?A^lh8hy(Ec!?eJG** zJ7asvD`e2cFW2CkVnN>oOK7huI{c=E_ONj&KE_e=lwUm+_G@Zmk4`tguaNI|(U-P$ z>I`60|965 z$kDdU{5=}_J6}YjdJRoEp7nA>MdrWWebI*;aRgxSL@VA%UYjVE9?zrhQdaDwc!kzF$S>q zYN5gLU`5k}reW4;E1j{TqOux`~# z)JK=vvk(yKW;%B;YTl(kdLbI6Df}ad8MnGLtUBISojs@3Mem(fM_1W1)~trOB6}9m zwSZ3l4>0*vM-PcQtu8IHXP~PSk(}4f;{~d993?yK(+)(mkq1bpeXzs6NHsY*=F^n+ zig-Y$=VF}x4ui=kUSKPXE&PEXNMWzcw(-W9ntzN%G^kd=8kL>Uqzf+p5`XpYkL45i z$AK|n#Ui^E3Nv%&@c>LQ_9cDw{X}17GwqnKa_r1i-y7zt!=H0keS|!a!c{qI9)hbf z?BuB)km={a6X&TMKL5%hdmgMnGwRswfGzkMOYOIN$#2_I`-}3Lz8s@P9zspZWB&SE z*zfhcSgfbC-@pBl+wVEh{t5Qm3^M&?@Lbw{k=-EU@suj=|0Mg3gP39VTRey*sfErX z4sCoX>mQz@nOL!^F^_{-fxRr$)zM?%`K6uu9`tKlOduI!M=!X$Fp1YFX5sW6ECWRG z3iE$ojoe+>%OaZ=ckWY!NRshD0jnCaZn6|~rMmg4#IT+Mfe_47kd?*tzYL4(Og5WJ zvkKxAez{)V^d0cAG-RdV?!OAoxdDuE&bdAo+?Y;9h8?kOTNb$mWb*o#GB4eIiTJ*U z;``g5!9V=+#n%hP_hM3k`=_p)^WEo`W5fMOcLb7W0K9I{3wRIbQ< zgbB|7qxxMs^hLDA;^QkN#Vgg&)i7}^0pS8FFC3Jv!;y6Keo6i@asO>G!K4SE?ZV~ZhnHnDy+Ty z18Y`*3W~yMHhtO|vOQfcSiolK6BxoeZUN`0qIFAGPe(Thz<~*Zw{qM#F1g6=00x)R ze27ptZpSW>2?CZ}gs3i9hLl~z7AxOx;5sY!8(3_(&fWogfXVzWjwzu_Dd|$tso#a& z8PKe09rr7AvOGb-gpBY`oY9hOe8{@5CtO(RtX9AIQ<6KR9Rwxc2=(vR?*MINBU6Mn zCi)u2OW7>)Ic+S4wg>F~X^9s40=WPJc@^eg;?NN?a_stK@NRZ7@CD7iPXdT{}lk(Pnf7{}3AJ=)2hQhcs-Yp9C`F zRklCDS35$Jt^aP?6`u<^wsjIT_LMdD^=19mmu8?M6peI`RgSL`K?+QPXIk zlSY-9eW7m!>09&X>+3qZH!YAU6ozux;@$ctY}0Abdmkws!&9R$OuH*KFVY_C-UDWj_Q3(yA>lkc!;5cl(TiB4r{Uz);W3RfN_D@1HXtalBxRxY32Jm@b zNF8%FCY9;q7yj|9Rd4#$O4_o%<5$NrNY*-C{uahQ`)mWNCem%x|cGFZD)$Yh#^wUiNI*=bILFwh6uw{!X z!kaosD%L;Q&wa1^TP2}4-AhVBU%1QNq2Df9tc6aKfDH%L$B*OY6Oo~={}>fASs`60fBVUc#P5p*{p=(TqNI{VYh}7VTb8eF2&wkCk+OL`gKmUK)z(jRQ=}D`6O1 zdqaBwL4rBMLZhuETSmmvi@0%%zYlYFTu58@ZTHr4`$cx&tg#qoeu1Q%^AOYnGb6i? zF8>0z;;$m0KMx)&35~$UqU~W!Z76@gjHeQXXF>gXe-7-gb?i-Z&-irREf7?62KUPu zXAmA(R{5Q)KPU}-=8l$|mHyG&AS>`Z=z*FOEuVQ#rWpO^Q*3H%OOdy1xx31}+`Yo> zk#QbNcf;~<6J3{L_dCh$HjNAas&V7%w9sMw!pUgV{c47th)P~X<^)tqLREHa#Zqkxmsq3S&W8x*Yme-vK}$n@%iLAn z%*jQe&eG6pI30N!!u2~pjz$?;zXvt_5X2j@ouB0^(*J^taTgpq?krdV$69bwh2Vf1 zhG`Yqrx%9h{Rpib3Lm=UB7Ch;tW#?}#V8^bSsA=oUHLlM(4{`HC;Aq<2*>*Q2FMD2 zb)yxNcCA0DuFpgT3`Zd(cNtW9iHv~I`v7@iKk~ur>JP$dao75=y5V2Z1Iz*S2VH6r z!_ZxOA@U#X_8^3|5bWZN{h0%Gj6^zCVNUSv(50JE=buwIy~`x#BM)({jF8gI7S0cVc?1e-Lp>o&-DdBzx*P=w|6h+X5rT+^BB&8EBA) zT;xYAh3=vF#%KZ3H2;&;RTv$ym_N21qp!XKJyX}A2Fz)7QvmVk-gO&A(bn9b<3O_T zU}@+hM8V6b>)*oOp#M>dSIpuFK3MVBSSjQIgNG|9^vr!Y=M|w>iOY$zXDIB( zGE!%6%-tJIYzPW6u1hIf5l|L-A=aQ*Dbove$YywGXWM`)A+z>|7K4seEE?&SA^*I( ziQkXpA67S>Wx6_g$seR;WC^*5>;_7G^@HR^)N)rqa~>g8r~xymZh8hT)3#RGn<)>~ z1~3}M1ItD%1tI%TPReaV%0X(X24Go9uZu*Z6-!wT8kU**ZyAc8ef{hG+mc2+Z7q?z z$j^qL|2_r%BfBYFMBKZVV~g(S!>5PMZVgMa?V!UPoQKs%{e;+XAHD@0S%PT$B6l%b zM8BNjeVAcnoV3pn-bZ-zaC`^hPeO{mODChQr(yYifbt`(X4AHeG6XIHx)28eA8DZ; zt+Nm3W}U32;~LCseK<)(ZB*NQWEyt1iDXo1egKab^xd>OMixNyHs#n$LZ~4%;7kCx zyqJ-KUqe=~!8()^n4d93#}RQJs%s;A#EEG|xc^yQf!=IA{&n^3#N!eDw>U`oYwKI~ zkXWBX@MnVMHQ!%!h)QF zQSv99Kwd7FE%5)$0%=ORGQp1>AuD5CuG?&NjkXD`HphgTiHg&ikH6Ljf_1A_)%%>) z*f}@iUrtwh>uuFdjRE}oW?kbdTcFlw!#|v^@&#<@3Hoe}O#z#~DcD%Ug!zMw)q%RE zMhIT*Yiad)s%!jWyn%0sN;q#)QTh3su5s~p->hH#{i>1oEE*-+1GAQ}@;m$!0+sQu zhCLUbe`s*KXb-#sWM4G;n^y&#$N4{>m3E!!8}z^AuQMNuZZVX~t-1t0xTG&(?di*>FgZ@!wn zCaYbUb?w)0z1=>Z(8OPb&(8miMvebpO85XiJAm&4l!fFzhEFqmR8Poai*;i1cY02K|R?HgsunwQOF*nq2t|L-Sa6r!I3Tn}&*E1)k(KMeYkY}QBdbpr8`*k=Q2 zJ^HMOl_!Z0^Xo{z%Ti#vJ#7@5S!wodK$GnsLmZXE-;%~fg$LR&u#p9;_yA>MBTROVW&|Byc0lJ~f&*kaq7QEYKM zW-hW6Y#iaXOz+5?XDKsPzzEabmIAlMG23DT>1>O2wk1o;o(}eBA;W>S@dD~3>#_7A zii>%c?HMa9J2K{5c4o8yy3t-b%92;aSYe<7koST2C*XAnJs9ssliOnJ$e3rj&-5?` zwn3bmG2N9-u{Q>`--VMt_S+?lx0Q5%%v5BtZNwms4)c7=cC%@Fnq`N%*s{}H4~VhN z>B~o193XN_?cLewC8Gq<{On}a;%pnnTaWP^r=!uGsPp;FLr|@Oql;M_63q_Ky!a(F zcY}}mj7GU!E?eNT1uk3QvIQ<%;BR69`8zK8 zyDi!ZQAvc&9IfOC$Sz>0%ikQGqs<=*{-uo-`FpYJM4R?iRIcS9_Pa5<4pWi+XD&vY zsJ$I~0V?v^xIxGlzr6FiS}2TdB`dVQqk{i|i7SFlDHFeIqHPrw+V@hS>pm6w zmu)Kek0(|fLP3eI6&0y>SkU2koE16#4WcdSa0|ochv66s=fkY6i;3YJnYE=}*v7GT zSoGr-iM8cCg!sQ2i2Y7GL1F=ItWjJ7mI-*DfXxE73%E_d?E>x;Ff3rNfPDfE2&jza zjrx-lWLipBtsqZ#y4ON^g$Ec|83HMaFkLK2wqFAReEk$oPoIk5pt_#p6dQGJfOn z*@}$wcsy3U^pBy#gw=0Linl1%Bz`n0s$xruhY6Mktf ze~)kMHPNb`Y{W?Vv9aH^8hcI3Rk8dOkH1<;wmwX_bC&reo*vn-%!uK26SAN+DPCsF zAs$Lg8OkToAq=lko=b+uKPMtOt;+TUds03kd~X8&YXUFxKbow#@ttPSQ=~^WowOwO zSp#~KFRv$D3jVW&$r*0vXMj(pf2w)?l*t*R?4Kqj55Z(|*gYh&zx;i&q~6cu4EN*j z7{PEob%q}vSLd0W;eNjx>pY#Vk_c&I84Uwo0px(VPRi@rO$=|O;M;P9!w2KQWJ9~n>et}nnJzK;WqcC5x^EL`>lv)myi2{$% zWrdw@p>~I)znSA*Xj7Tb@VsvPUMlGC6B`9~<_3a>6!hy-;GdFmZsgM085F9|0Z;Z! z^w*2P+c;_5cQ^2()2~tF^-p&865fy}an!)dVc^L=+r`FYiol!}s=F0nC_ zex!qQEAmxdLY%J^^wwLsV4fEg<@*BPJD%g&*&M1r74wzv56b7=$}@s~Tf%&QD&+|I za$cPR-zV(6R>+|RC8~1ujFKQ=~?LCJ3YFbpa3F3)9%tf_J=lAluu}>IpP?s_SCi)pdkOB$9GCr($^Y zCv^qK(72OlNQNN!Q*eRAZ=ktD8k@h~hTa}y=!G?T#Iy$BS`ST}zh(xF_?O~%YFe8- zwcf@WdSOn*l?)U^eIEk4oY_Fmc4rnDiXYwP@hX)qeS zMJGX*`c=M_End7lr>Q1bZLlIOj9AD}^orN##IuXQh?!Nv(^b!m0uP+xudDGjR{Idm z{86leA*y|?tsY-XOH)fst=lInUZDdfBTU*VX4v5KHbRe;;NtPlUd7`vyy|C29g1AM z5@-lPflDs}ipO|7coUIWYC?PX(vQgblGhNC;o`*z9}by+r4JEoj32Un{9A?MB}mhG zo9n&RK4{u(=#P&xsi5L)UE2`wt^^FUa9AsA%*@S-6A!h0&Q*;;C(PZ9b+OifS{bav zo5Jd91j;?TWKzJpN@3)+-qu>hS+lkgLvR>q;XU*;9M8q$8Xoku`0BmHAZpF^0maDz z#fhqORTG+2Y4uer&Va8CP0E*`Zepw5>8lkPvbF}&WQ)t-sgiTZI)?H#)Ik_$M^GqE zWReD?UYYg3GNd~)Y&Ef`=S7|eNGQ**L^q?Ab%HnJLq|oDe>$$)R3x;q*3cI_6dGzO zbh~5ZuMqqa+6}Z+#bqz9l%Qqgm*)@?+Ty5(^V5C?kEoiY{_Y8YLZ``lSnAGwTjeV z(n+`u{B%Yq+wwd`!geu|w6WA*pp~ttkWVDPJa>^WEC}WJl3&jM3Bf;G(8==|3FWyZ z@e@xX|I@&cO{9EmK6w$4;@emflafAh{Pz<0W9ML#IY(~-L*n=!f{flDCelvqyiC52GLc@k zCHw@us)1jg$7nNwM3Hp7HPn{)FVI1;E611TZ}R=3(lE6Cs7OBP$NvP5>=)x-NL0AJ zp2Z=>m$68G2}gk0$RE2ecX2K`j+7z!WIG4Ugp=dT`!B1DQzX*USjIxc_e>m7k!vFvP literal 0 HcmV?d00001 diff --git a/examples/16_100_lognormal_samples/example.c b/examples/16_100_lognormal_samples/example.c new file mode 100644 index 0000000..84156d2 --- /dev/null +++ b/examples/16_100_lognormal_samples/example.c @@ -0,0 +1,17 @@ +#include "../../squiggle.h" +#include +#include +#include + +// Estimate functions +int main(){ + // set randomness seed + uint64_t* seed = malloc(sizeof(uint64_t)); + *seed = 1000; // xorshift can't start with 0 + + for(int i=0; i<100; i++){ + double sample = sample_lognormal(0,10, seed); + printf("%f\n", sample); + } + free(seed); +} diff --git a/examples/16_100_lognormal_samples/makefile b/examples/16_100_lognormal_samples/makefile new file mode 100644 index 0000000..ef385f7 --- /dev/null +++ b/examples/16_100_lognormal_samples/makefile @@ -0,0 +1,53 @@ +# Interface: +# make +# make build +# make format +# make run + +# Compiler +CC=gcc +# CC=tcc # <= faster compilation + +# Main file +SRC=example.c ../../squiggle.c +OUTPUT=example + +## Dependencies +MATH=-lm + +## Flags +DEBUG= #'-g' +STANDARD=-std=c99 +WARNINGS=-Wall +OPTIMIZED=-O3 #-Ofast +# OPENMP=-fopenmp + +## Formatter +STYLE_BLUEPRINT=webkit +FORMATTER=clang-format -i -style=$(STYLE_BLUEPRINT) + +## make build +build: $(SRC) + $(CC) $(OPTIMIZED) $(DEBUG) $(SRC) $(MATH) -o $(OUTPUT) + +format: $(SRC) + $(FORMATTER) $(SRC) + +run: $(SRC) $(OUTPUT) + OMP_NUM_THREADS=1 ./$(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 diff --git a/examples/16_100_lognormal_samples/run-sorted.sh b/examples/16_100_lognormal_samples/run-sorted.sh new file mode 100644 index 0000000..b7c3750 --- /dev/null +++ b/examples/16_100_lognormal_samples/run-sorted.sh @@ -0,0 +1 @@ +./example | sort -h