From fb110a35f3f6e50dc655e7502c00c360b90b6211 Mon Sep 17 00:00:00 2001 From: NunoSempere Date: Wed, 29 Nov 2023 22:51:58 +0000 Subject: [PATCH] refactor & recompile for function definitions --- examples/more/00_example_template/example | Bin 27296 -> 27360 bytes examples/more/01_sample_from_cdf/example | Bin 27584 -> 27648 bytes examples/more/02_sample_from_cdf_beta/example | Bin 27568 -> 27632 bytes examples/more/03_ci_beta/example | Bin 27336 -> 27408 bytes examples/more/03_ci_beta/example.c | 2 +- examples/more/04_nuclear_war/example | Bin 27416 -> 27488 bytes examples/more/04_nuclear_war/example.c | 2 +- .../more/04_nuclear_war/scratchpad/example | Bin 22280 -> 22271 bytes examples/more/05_burn_10kg_fat/example | Bin 27376 -> 27440 bytes examples/more/05_burn_10kg_fat/example.c | 2 +- examples/more/06_nuclear_recovery/example | Bin 27664 -> 27736 bytes examples/more/06_nuclear_recovery/example.c | 6 +-- examples/more/07_algebra/example | Bin 27296 -> 27360 bytes .../more/08_algebra_and_conversion/example | Bin 27296 -> 27360 bytes examples/more/09_ergonomic_algebra/example | Bin 27296 -> 27360 bytes .../more/10_twitter_thread_example/example | Bin 27448 -> 27520 bytes .../11_billion_lognormals_paralell/example | Bin 27336 -> 27400 bytes .../11_billion_lognormals_paralell/example.c | 6 +-- .../more/12_time_to_botec_parallel/example | Bin 27496 -> 27568 bytes .../more/12_time_to_botec_parallel/example.c | 2 +- examples/more/13_parallelize_min/example | Bin 23152 -> 23143 bytes examples/more/13_parallelize_min/example.c | 4 +- scratchpad/plotting/src/example | Bin 22640 -> 22631 bytes squiggle_more.c | 37 +++++++++++------- squiggle_more.h | 4 +- 25 files changed, 36 insertions(+), 29 deletions(-) diff --git a/examples/more/00_example_template/example b/examples/more/00_example_template/example index 65013ada734451affdf5bf5a1d2768b9d82c57a9..b854df73ad441cd666c8d88d40bd7a11fce471e9 100755 GIT binary patch delta 4284 zcmai13s6+&6~6cKa6#C+3z!98tY(#k@bqotV=pYY=-QCTA`0>lK}}6%@KG~kT^;gp zJK040cUqH48rw+|(@CeMn>M7GrXV3vYa`-gG%>moA6Ik9 znXTCS6~sLur)Y|LGlktL*lOfXiPd9B(bzVBvzhm5X;&_;(}8U?t`nn zaOJ(ac0Bol`zw1OKaL}_JesRR>)(lONqB*)3?a`2qNBvyvO`h z&C7mhz}fgvjh=fK?A6>WRF9s0syBH0GMxKCSaieFsk?Ij?ah00e;Moha%aCwOY+W_ z$|R3eTI@Td-f-RRdBwQswsR4@bCH-#J%vyFodO9R>ZxQ%|nBM0|pbNE?wyjD#dH;|wuathU}?;?$D@S zIu(#gW4T4GorhDn#76geH-Y!_7q7Xxhd-}cJogaj_FZ@SdgPKDL|_jU5Atz>ffa_}AAg47XE^L#>47h=4PoKDOx27`3p`aJ;lA|f6|*SkH$=X<1{;vF(&p%g8DWM`6BU`k8KhjPhicDG2v}9OU_)M zYU%1G5G(r%!xw^=(CwTKeDyZ~kF3$K`cV2m2?kGPCx~1wu(T8<{mCcB@Hs z8}Bo-!mJO4R*~J#YI$NLq?*Iwn89I4eDU{i=<-+E8bOa?*$f0oCH!{^RLuoxbNi#30cv`3Oprz8{B>FeEYFWn6h%yRQcz(}IV6vIP(GLBmvE@O?j*aOlOtG8 zv9#Rjf8z$m0yY2S% zx%^qSc(x(WDrBEv+w+E%_QVm}*izy0PY+_TNbBzN%dj>4>}gHz_D}!V&T#|1{vZ0~ zyybz*_aL6bN|U; ziqZ-U0&mg!vZCZhgMq7X2y6uI0|tTpz><)n3?GCwupC&@2RpVYKP^SxuP8ahuyYWq zm*Z^sOu)y%+Mq25v#OX>80bZAH9jSbFR-L6!(KQW7-abtJ|!%#AaTTdxP>ubnR6Mf z!k$>8Vb)-dEDqt*%4%UT9QQCqv&a=r!lHOEhu=aj5;6zdS}=@X!v0k7s1Q8Nx(c#+ z6XTt!yn~Ix_eLf;6Vs0$1K-#{gIpiSM#Gni@~kT1RA;@At!6)W4&%46SDjDr>)3hc z9Nx_)7S80Su{DJSi`0~HJ2gR+2ko3~G`%FajL9!V%tGkx#^geysmN#`yGq#EajPM+ zD;lty1G|5R?Oaej2Rjb8zKC6NCVQR6+!vx}8K*Q4nqw@9$mESvBx9~*v=n;IU8U{PvS z^ReJo8m)CWh7=Q^XoTW^xf}8p$d{8nU0aROK=w-@H=$NldjueWjVX|?hrIj@`&-eZ z=oUogr5YsJD#(|=VO{&o!m;^v`IeEIJp9&dgWSXoWj|YBPI*hW*{+`zek9*j<8RlsxpvLNBK@|y zRP8sqcIS9ayIBM73Tb(swJo&N4nen=g%r+Iqbqe7rOF?yp@}~rpA9*2;-DX3jGt)0 zEnr(0TKH9T7XLhbr@WzC21jTn^}6Ld?CL_xh=8u$sFO;+ku=&NAFbBVN%mxgg}_eMc5oy?@S~+gcri9D%`u0&FpN&t2I-D`Ew3;(2KZ$P&u?e^)01O@6`FD^ldy_j MtgVK!wk4_m148CS#sB~S delta 3475 zcmZ`*3s6+o8NT=OaK(7f^3W9xVbN6=R9F;WjW%3Z^rFzLad`>KA}(r5l;Wc{&4319 z+jcGKNjrHpTF17w(@sg(OiiXWAhyw%2k|xd%nDW)H4#~>q(M)=bMM8>nDort`~Bzp z|MUO<`OkmuZt%3Y=d{>tmx9Y*O1l)}9V=+VJC??Y=`s6F`gFf3us}=~Q(THZe-yJs zozJ4LN}nWz^TtKZe*MO)Rkxn;HSRuZ`IV<~SKx(+1!U;jY^49 ztQ0BkBGq?#vFiIut$vX2yRWugi%^?xj;bMNB5k1qHyF$J{VQ+V=M7fHUj+I60{@B* zlj^&b-}JtTo$OOc?EL=1+1EDQtRc6K)sQujwlFb|`0;2#08XjqsF1Xeh`lRJ%G!v% zleo{3^cx!+Cx^-8RxCOXySdr0BSpb5N$Vu`aF`tKjVAUZ#C#(-M`E2eQ%D6<;Jb$zwm#-MV|~ONX@hSr6xd8^(~(hqWJ!MwY3K(Y zD_4DoAs=0;RH$uTR%NNOOdZ;xRMx0uO)E{YCp4nKsOft`-1CCs&iCCZT8!G7LoWDM zeg9N_?FIf8t|ZAqf$y5y)NWGPX~fU%&!2s1!%fxq#k<)#M0n5DE(oF5zG6(D_BP0Z zk@Gf?I1N4L(w?t#K7pLm^`y;6wgZWF{Qso&l)KUjrGNC$%t_0nr5wh6XE78Ft^C!l zm7niS+rsV;W4!8Nvg9ZZr5y<%SZF6`y;!t&QbOb2kpdDeJ={miqIN+>1gm6OR6{bt zWXJjkj420i!Vj+9bkLb6_$m9RpNyyDGRmZuE3_daDdu*Z=6Cq*c>3Fnze(0Oy3BF# z>Fr>3hg%`RVISGC7YMrWde0$C~SoS0_#z?Z0CpGBF&iSZs^7A5i#$rB(WDNH@KQ zS$~kun3g8J(M6X}n=vHeGR%J-q`Rl3jY~et8&#pHyaaNC4P;1j(*0I{JKq)k^=wSF_Q#ADbDwu;kW8Pid*pw6v5+T zpcfaU-seRL-MSc9oAuJLfNPjX`fGN)n6yw3Hb14c8(!FAd#onyp<^d0BZ;wyA$-a@gK7Uex9R?uJOuq3lp)dyZ_KCy&aJ&3Q8aMW2IbGuG_Gi3WJId9rzz zlp~MwM=X%7KZ|tA_U%!M>}ZO1%Gpw#>?oA&r~?dOm?N8bAzQKDhE6Eq>(|QGYIx>+ zXAu;&k0>68q8I!D@T>TOJlWhNm7uf>WV50>G^1rA(NgtXugVUnOtT=`!Pj=vwD~il zS3vCTrZw}^qhEnI@NfFs{AZ(0Xr0;+jZr5=XYnrL;@EU6mU=L~gKn#W%tLO6Plu$$ zXfc&xS{=5+CcB4jSEtI{%8BO)sOBDeSj`#w7~1V1`ZOw_4Eug0H5HahrBegT3mYPV zv};ki=%TUDSj3%y)MqA2(ms{`aq$G&Dn-&4+$*5F?4BqJbZAks9Eh%0cO=HU}jkT~d3!?>x z$1pHp61RhK&krpQ85S;M*EAR*P|*ZYj5j=d>8t>09L%=|L3zpuEdmS5;qCTkEZ@ zTvzp6U9Go1@b`+HB3-%cOIla?v3;d?edTONW%b%MHLF*6*Hn8eSFc&`t>0KxE7aAm zt@f^4SLv;>1%J*A(3C@fTubbA-$bUb5^Fulup%?#GQS9j9o>J$*fCCICN=5-eAcbthjnH*#^Ei~H zPA6RJR!6h5QRnQgySXjx^vB)pY$)2Xw$_pSgMSJz)6uYSu$4{Y?(e?$#?-d;o-^-v z@9+EW_wM(*-@OmIjtFhN!e*=Bcs9K+#$zQ)^H}M4!4_i}rw$F{{3U`-m|UN)j&HLI-H;kQNA;i@l_sSH1R>ISQ7@tD_ zp4cbuPW8V$AurzBBz<#Au8@5fLsG_RkIVPHYu%u3SX!G{mqbWNGUTQzShc4~vi3Gf z6~~$++wkhyA(xcdg}KK;p9?u9LsZ`cLYgE`5S9!$S}%vB-6PTH>}`^zijk01)7K;| z><&p)10X*}sE-iJdtvnq9t+{W`4IdGymI;rYN&%e)INkciZB(Qg`~=V!UMAxo(xHu zBTnB*m+$M~mMe^*3+4)}!Qe)WUF-JfFxXd%$6mRnTS|RWk!o1F^*gyjN!fSN!!!13 z8RNB4^Gt*w^SJwgd>PXeR5A}>0SB<$bE{{CGw*@RoWVq%wDl+w@8zq7pEZOP5Anpi zG3gTWK8zLbM&46TD|uNDUJW)QRIlaU1Ais=ZAwJX$4WFrIu_3TQ+S+2q=9R4|IOv= zbKivV4Y||o+Y`Ktx~O7PWFj3JK-iw@6N3#3h2%xh}5TqD=YOXNCK7D@(1 z1%-=!hh-1fb?GRV`U=}MOlDfph^ zeOB1wT97za?RHvlRBA9CC&duFtcFHcDG+{&urkEMt!CX7JfWIvEINd>g2KwI-MA6d zn~45-`n-e}@zTFu#3lX@{&4yFf|&zUv=vX74MXs~0mk0vH+--My1NyG;x>uKXHhX)qT939WbKp)Sv?-++;`AxlCp-S4XZ{OlXngy4!!|pptF1U z5pK`@N>3Nqt3=XA->^@LIS{ABjT0n#$o{Tqj-#8hci#OXST*4D=pkT8+VC0zy4q_j zEzn1|d@feVYxr+`rAji$BfHDp9%sM|cYD<3SS`Qc))VxFJKnx41DDAa3CoT_fOX#sYjbg%jk;hTCFS`kR9>Ak}K%`;z$v30}^|| zOIvyS!@MS?eLaub^UYr5%9|G+N)u8d?&RO}OpLldVXOmv^M1zm0qwwkU|WE(?IIza z2N>IkaYqkhBfu_TEO_!MV`iWom*kj=x02P796+XHL?4gvdt2|A3wWvm?758MJI^iKs-1T#HQU=ot(*@9_8IW;;> z$u1~g4n@Na2aY_0Uk8<)@wOe%KUXv1&LnT)*9$I#b1uE1L9ZibbU(&#C_tZZrU}jT zx6V=_Krc9FB?hr=_rb&){;nc?3!QXMVM<~pHtIozS1mhXd5f+nOc8piuP{Zta)|z} zFim)bepr|#Y^VPzOobu2C^bC^8`Y>yvd>>?%! zdUF+4i&Ki<33%A)vqg6c20BnwDxT`2w)t~~gLK9Gf?5k4_+F3jiwJ#EqBlGvy7Wow zBj)SPVz)l2P;V&G>$o4Ufi8aealfLgex>mHB<$ycsRDjmP$RYU;ra12Fh9!c)Z5ob z&)4Uz8&j$;kH{4Cd9vOv>#aq4bG|+)UvJ1)R?8=MVDb(W)&iahUh-=Z9-ATx_2%|S zm)`nJl#9obc`PSlrTtVazJr5&!T3BcS%d+3NN-+>*Kk85ic*KFy^+5Kd^z|BxW5Z| z)$4WKKL`Bwus;?cLXkZ1Ztz2==^L(@(HqfSb%vgE*`nV=x9u#wwfJYznW*BSujr=Y z+oE~l=M18oW;0{zlNdgN`fos=K;!jd3F68zdV{Rj$>RtQza2h%25C&m9WhP}PhpsJ zo|cs4kIh8&J8)W}{BMpM7U+n(Qgr3{?=9OD38W|Qu?h9mc&|wa_|xy5CW^78^p6W} zqg`SoT~o0Pwt7MZdyX;wS_xa3VK zF5sKkJN}k{U7?8se!L~nXEBRkYeq4ig#T|+ZG_AM&u4Yu?|a^+YHKdxvxLRU^H+4& zB2)6yVb3(xQ#m>4K2I;e^HuIi%PS+Zd9*I&3ix4d|N4oyp9F8FZ;J841$s}V$vdjO zgYYi}TB-^}kKBw{Cvm6*hDw{H<=j769f!|(Mwp+X@|l_zvo(H5ySA3AW>T#1gr2E1 zX4i!0+EinBSAMDSfmtfwsLkC#2P#d%HY!({l6Po6o3-`7r1^YCAFeWa_h?pwwnKi+ z`jVElfu9^SZQDst_>Ji*A5?+#X?)uaD;(B%v$jKzYW}0(4P+u6Q*BI+!BOKo)DJOv zOVjpSq7Bui`twl+Lkca@}}x=RyVS|3{4xN>D(V{7X} zt>gy}%rQ2U3jP)M?HBxxW&IIi<5Zfq+!o{F_Z=GUoa*I8xS?KMo@ET5?l3mff#q4_ o0|lyU7#(B$X%CDFf{RXmaMsw3HL7(RZNO;CQgt+u4m76z7hWiQ0{{R3 delta 3684 zcmZ`*Yj9N69Y6PG6P76V-X)Jscr=?3WP`Zb03o4*Wr6@2C*`m(Xt)$g^3e;;+po+PXHM&J-N+yW{ok?)X5N;1JTw3-$dE(JSbD zFMUz`C?T3Y@ot-9$(sqZ^=@O+HgcPXVsK5EJgH-#B1^&KpJ z@~du}%U=QdUAMn!#H{+hF712Q%s!cL$(GfVYnPvY_)?4f^?NOHn<+CQD@hlvgaB_l z#?d7&{7|-^iO4mBvVBzcd@egKcJ(ZZ$cut7ItIO!`=BRRgd+07QQ01e$Op#bWc#q} z`6eQJPs)`;>eHr7bueth0Z^(rE!#hkJ^zp$7rMBcVScu3Kh95V(>xu4r$cb_S)A-S z5s|Ckm+hzJN`$az((Sumn!4{)KJM%=#hYTuCA#XiOPt0qgApM3a|>btzZ zpim}nmk$vVe(77r^zC@JM2wi91QXrRV-BDDf6RvwvqKMBL9inTbmXUl)}`I>!izzL-#kkPo%7?nX7Xf|q~9daapyMiTMO;Xp<;`LnioHw zodRLkOo(%nLAA))5qbMVHyCs6y97VDg>PC2$vMos?9W|ep@GFUV*P0RGQGTH-`&^csKcwE;BZ3T{wf@*es5+ws%CD@Lxgf#0ol(Z)sjITIYGC7 zMmsFr{@4uwEx+);Ptt2k|F~ivG%w6$jLB(5JgBHAyhh|g?2q{+oc|RL4evNfzs$`R zA2~~{%T|gh`fl1;x*faMPSAtP7K*#h(1~TsrVXEf`kg0ec3xKU@NpEN?>nJ8*X1w1 z;r2yT|7}w@)3!XDc*hCapEuic;AUQKm~2q}EB9f8yTtF&g;tn;e@177*g9 zEcG4b*Xo+f58C+=WBY&wz!6~o2xBjZgakik>{09o4kJC_2yhChe8QLl zc>&N43<2F7kH8*S{wZT$0DF%zMq&sV2POiYqpXKHL9~Ac1E6w@v39;c$=GgS5Eum7 zPB9h&_5-g1M}TG%9ETX21MEM|SUxa#ma!ec5GM8rP{I3^7>k4P*3<*X@s#fY4qRaD z5>Of2z?dlx2adrJFa%rytiQ;Z2WY#5itznq#)g3fR~TD54Gw_yK-&akJAwVc0ltqg zHqQ6JeEQz11wsM6vdSvt)9FPEgm#)!WKC^FmCI4J=)&Zao%jW5Ls5d`1JpX#5aBr> zZ{ary9`6rPC8a?!6C216@KZjfuNGwpJLx+`YlLxXDb7nWqa}?H@s)z#4jj}(Hx|!M zdK68$U*qc%-rytDUp!lwq`xbkEzUVg4;E(#FVG9c*}@^3T9O9A@{+Wyd^DxmFyyiF zNq8Q-rQhT20WV0ml_bnvgqgjaXBRV9kZd)G7S|BJ;}vLr05mC)lxM_3Da9XCCfRnHiX{6}vGCOgUj<^jBa) zUgdt+{|)w4+8#2w?&-P8k}T(KJODzpT!Ftut7j601Q%&+qSXlGL7G!# z%{gL7oyod%LI=`~%5bZ8ht4Q~9$#0etJ0by;MQBj1QPuhsLKf9&j`2Ljek6R&F7_H&OmPoxjM)H_oU-{}P=atzI&Hu*RB` zWk~Cd+G=wN#iD5SjJG3tVLjX`wbhXJ&eDD68&&!xeWk|Q^Ex6?$SmTl*QL?#vxZf{ ztvW}$Tp;%t`~joyQT`)FMhd?;SlK#|{f0jILCpWu`82QAlJgm89)nVE7#VQRunnQG zd}SDsk7;?Ca9McwN3wKjqajV8S8A=PnT9lElzfM-M4VhI6!LFrC3JXCe|T=A#{coW zb^QeYX?o8FYwSyq+UQ;}L3m3$;~@Vwp9EK%Y-0kNjX1&vKWNh##RmVX!5?tyj8S7U z_*=mJ)bob{>j5L%B=9q|Q*)@rYmF`7dfL96y1drZI_xc|t*u@&BW4{)CwiQWH{TG16xwiq-t=m(F5OOF!R}O}zMDnUThjgu D3@J*w diff --git a/examples/more/02_sample_from_cdf_beta/example b/examples/more/02_sample_from_cdf_beta/example index 8bd4d50dcb4d70169e14688126db404b632e5578..b09bda41894a9325f5c48dcc5c82e1e2e9026d02 100755 GIT binary patch delta 4341 zcmai03vg7`8NT=C!REDhlc1aBVHVinhK6Jl0xTLN8+MajSjs{k1ju94kQmx1OvJQ~ zl#N9(TUkwdrX3%nMQ2*k+A+9Au+u3C0VmL?fII>k^YC>^BZLjsh`{N0?!AkG({g9_ ze&>AWKmYmv@BHW7?Ks1ao#ETef@MSdoKIpM8Jzq$GDs4i8Pgi4bX$WizKD-$ka|RMNnS(2n8zKeAH1dwV zU3hD1(4C^?g;%DK3d0JacsluuVL(`&8r(U>n(V6=FMeh#vH35D#5rG-I{nw2+ebCy z;kOy_h-PyeQ-pGdeN%BhVuuyQqu08}{{Lsf3X^ zFf@M`5)V#9qj9+2CYDTu#EOA>v8)@ehCueh)ls8Zr)u&tt7 zOnqEda9q6WimgOW*?+l|W$aTk#PWb3SJZ-ml2xJF6DZ_cF?x&I(ADaZo2GeR~NvnApDlP{3RD z`@Q1oCjESJ@!TY!+y5WCe?;E3lLqg26QoU9 z&uZ$Rt-@AmtFl$2vQRQ8Dk$vopSG1UU2Cs1sjo79ON*JpJ}G?&8{8Pq{lm&7!SyPY}la&cQH(9O`^X*jKDXOaaZd;irWYzw! zuuBwgHwpV=U2!6D7^N!{HwqI<<}hI;P)`3vMcDNNw-J^MXRQ#f4UlY0mGB~iSw6K^ z;S1t)hQxJU(y33vZa^Iw7uWSDLy7Sy=yeCBy~=pJ3G@jJv*X!?6L~ZwRwY=p?3kCc zlNog~0XHXh7yH@5ZVI-iSkH{H1P+NKYPa2lqf&wCI4L?QPRRgGxj@({LS*oVTg}`h zjVdcwm^5%}28GDX-MA5yn~3}|d&%V42;iN|xZD59UQYjjl-Ex}*{(Z}AR4K+pV9;D z;+KZel2-JSRXJ&~Rq-qs-5sYBNK4LRLPR`?m`8*Mcl zedLJ+d-LYN^89p4shDbKL4`eKF+S?Y{0hc@g+(QEA8B5gCFqC9s|yzk`_IU|R|kpS zx>#_YA!XJ%Lf#~Ub9;Lg`?s%Wl9$+VM z5U4vwX|KR>USKQwp5v5G0NYMcIt{!9;aPzlzzU!TS+_9kL-xQ<;8oz*SxQqP;O88r z8Ngj1VFx2PA2;|hr5ixBpjcahgSe&k11#O+~<^f zfR-90uKHg;|xcUATa*|rIVxYML}@2*40x6GFH>rn1ZBQ~LtqQA7dQr-0D5u9Ow*tbEC+hOg&j$^-_0K*tL%Dy zH`$2xAZfSj4Sl@`9%_akBN{vL^N=q4th*XNr1T{=BQn{^Iz)SWxlBfRzo0i6wb83yk5jZT0 zlP%r=cZ>|dA_LoJRV*yR04!W_TnZ+949-Im3#am1$lStN!q_>ox-gSBlE(^Dc@KFR z?>5p^n4He;oLo#`4_o|i*!vV*Q;*0D+Q1%gkCSU~wwoA=R`RW6W6^THg6u6?f`q;< zT3Ts>72E91*j$OtT&C6SL`Eq)B8m`tw>AYi7HT!jjunCR3*5oXuJDH4a@f5WwsT^r z0(LqSNhNupD8*;jT6RPiX|3BQE!TP?GI*`krnT6#=0dHpK$}va)fF)LI5u%NCYn%H z#Y`70@QeV9rz0F%V<6J0HSdhV_UwqZF<*B0n)$zc;Xf>QnFWi=tXhK(w@GWP#WHLl z9wo_2_Ez2mz6bn9X3w_GtJN_3Jn%-8wrtPp18P{omw;~=Aq!j!qTA3M9U)Dw%;-aC zPK=U2yM7*>ixO_QNbX!VJDOGf8>~4q)`+P)S&?5P4=$S$%RF_Un>R*&w=5^tj_xSB zZDZubvVxdgRQZ0Kj;LT)+@n01aF+|#yx{5;Pe%gD@Tx*SKbXFHiNHTjHkRz;pCU6# zQw2*NDJZQHI`1K`l-|d`NQO#R3Z8pNPFXXX^m5t8`{X6uSWsU!*pPUBUq!|xQvBOxwdVbGPIyU}229!pdDN>(@@WKX%? zuv8r#T%?RPD8Lnyp>ln8eRyt;LWd=NOyxW8Rroq}Ry&);|B2Yv>a#yX+&XRwXYs0o zD!i@E<KfBjzQd#t z4we5})rDKA5PfPTo5|36y}_=|;#6;|)LF$Ou|jXC#b_e5GkO(TqmKTP*eXK0hiV?! z+*sFCy}72Trm?ZEF}SpHFF*5EmJ1X2*KMg@WT~#*82qLxf){I=nra@YzCn_2tB0p- zZrW5^_s~Ptbxlp1nz-*RFejp>oDVkstdkE~HVj4xrRgN8Av4yRr$mHnVp+ouzL6Yn q$V$q$D`Xg*lY&Vbuk*ZrMyOlyTQ`V?7(UUIF{n?jD)rT!PP(QbVJ delta 3583 zcmZ`)dvH|M89(>t!6qU1ZUVYlFlK=THw=Ug5=eMRHs)~?77{|9u#gQgX{ZrGKoE;0 zOa$G^I$@4~SW6wCoes5j8eCMg9Z3jQRMbTUgTb&wDT@XO*#K!pPJid#yS7uuJG1Ba zo!|GJ?|hGQ&JKMcbbcXp*u<+FQpRoyOc^3g|rx3oHn)j3WWk8#wBP|7yXU+ z$@u=Vnpa+Ls!VU_s`mHYJkavHgC8{cwnuy*$LyX)h~kZxr}>J$zD|~aoX*DSACz6Kcm7@*vj#1$Is}#?1$R}4jTuS$V#j)D4 zMv3t|Y8n)IeHT42YkVsTjG7)JiBA$oS)pg7dv>?O!2&{czT@PPOg+wh|_af zY40&O*ct5q@MvLH|Ko#-=i1xZIWjp`(nCn_rSBQjrb9a z8Z{?t{B8vE)(CFbJA6<)7lQjP(9$`D_ohM-oDpE`7{7YK97w^xT%d>NOpSd<;vpLE znQNqhIXlHViPmLM@#AAt<*bLY?tPNDcC=ko?T%yNo%G4#E%a<|5GPpYkqy zGy3A`@8?#F?yu0R{P6ha&og`TtgTHSP%IXrjD;c`3Z70%QjcvXTM&ibw|NQm_L zojKN4{<=;FQSuRE9gsUeX6z`i3+M+L{fzxxB*YDT z0`rp7NDbI=hOx2W?PnRY0XuyAv^sA|vJ)*^##MQe zaknTck^$|)bUQX^rrmktl24&YQ~7p^WD(zo!-6=~;#F_~dIJ^<(XLF*!Y*EbMPVFC z!h-YQ+;md@MB!<=D1VAb`smvHG{H{y=BEg5dKlw=dOCkx>K(Lg4i@l_Z+-;5wq0Os zo5~Q?*-3CO(D;Ihs6}c)vG5|@Rv*m~S8+#i@?+b;tDybktcT_GbA?mp&f}qyx+I1-mb`bg?K{ z>AJEvg(^C!JVo?p(VX%%qC1lwEPqINl3pw?7Ki51*%e#(qE{;Nt1an#FGAE@fydD1 zL>z)nyGj$aMiBdT-kz-S>Uo36p9OC3o~`lOnDdX(jTI(Arh6+*=E=D6+@@oqrc+Nd z#0=V9X%ZfyH!v3S(Y!uvIZd-%z)v;4C+MUqldwy*Jji_tZ_u4pCR$PzO>b2uw)SFO zf@mfNw`%_5aEbHSn`CPg8moRI`ab@t+B|&t1sL#1{SRpT3_TnER&l;lKd9T3=#^8CBg~iSygN%+47~R@l;c-x3Ma&xW zl&%=j`2fnxKlK#@{!7hk9IgVdSl3KdN!IF}oJ{wwG?^dNwSM(l5%R2NL>67-RAVIy6RaTg!WFNJ;bON{SP4alG9PT2IQtIqJT_Ns2_o~oJD>z225zSF7^c{=ac zbJ?#S1pcSTzX)79o<+$FZT^^^#dwYHjd)jrnTRNvy;wEA_yx1>fC=xX;x+b`-j*JRmi>Nc!zSocW%`nvj> zb?Z0Rw`{3xBF!xu>gqRbs;O^j+0f#vs;v@!>g=y;y9L@MM*5l`2}FoDr_s{JwAfz# z0OIEVx^Xui?Syr+679KKSYaHE_LZ)?D+oz+`J?kTiEkGmz zDJz7z|AhEL*eMxfW+Zc5v*cBv$$+8l`>=FiG!~sh&9dYf4NElx%~E++SgIO=@CaOe z3Rk|X>+_f|xW953@-MMucK<}xRWEbZ4L2v@#?v2`D*pxx%q~A2ma;~j{?ji1W##1? zl%fmC4QLf*JI2V_fI)XHK8NHQ+M?}#+@y6jO?_Fwa|jD32>ST^PvWPr>Q-vqKS z>Qz)T??C{^5blNbvm==whsnHIYNWQ4NW70l3p-mF5f8J(yD;hB$omUKd;oc;O;stf z_1OVC~Wq`uw;t1;^Vd&l$`UFmij*_+4AkCbeF$N4wTTcK-(vV z=9~zHZvQu~&0Pk10yBz7OA7{`{+96xGNP@Fyc#ePmZglZ8H18WkMAV1X9e^ljq0Ye z0#a!lyiZN~-V{D=W&3(Nfp_z_OLKJxe?GN%?jX?Z|IX>Zpe(-52)w|G2l+U`z{8NZ z{n$2;qCEEvrO3l$HQ<|I1+^2Gv6?y{*T}VUom`K~Ldl@0z;LPm3)#zfZM?~-zQOqQ zE@ce+l#_6TQ(Cg+?4~r8SErS_18s*4IF%fhW0cCbae7JTG^P*fFj@^L3q~s=I{QO4 zN87zv-7!ID)pyUzd{*=L|A?KhdAozypYnx z!nQKu?Q=vKh*aIJsZ0qd=;{N?)0#fP3h{BM z+40mnt#tEB}NJr*#!ib3^Wf7syMD}L2 zvo_-9(mVR|%mZ^84&5@ySNj zZ(@OooXCD(*dUReIr|>{Oo})>f*u@(rER}qZ%UxX+5&zIs~$##tcE{ek*dic4@bxy zNOWO_I}noxQ1eqZ4v<~Bf0~;C%kxw4KuIZ17F5JjE|bUnR9?p9-(pck`HF0{=L+T_ z^1gkZut6KWIY?$W<_U(cNV6kDm^(;bay&T3`W4Lg^pm5G^!sM@gMiIH@SE~+o85H# zM_hsI-UK4fvF6Z_1`o>vnX-7s5*m$d_fTJ3nB@B5NPdE3)b;SS8pv~k?6!~QS! zbl&j7<^2#Zb*HP(r$SE!_srzSzE`-E$Eo4$#jn1?+FSFR)pmJ4X^|~HXw-(BS-Umn zCr3TReh|PT*Zvu$?dXS&Qo0w|0}KJVW0YP3wgXoITaHuO2kZbE1?WGg)C#osQd$V? z1bTpjCn#+JTK++ObO$=2K1$zV0^~XfYyookS`_&4U-R0x&hcS zL}@p$<20q?qv7ui76RIF06Kvdob+ynxB-WO?FyxHV-U#sMoPD%Q+NS}z)s*kpz$K5 zJ-`;A0<>R3Cb1Y_rnC~c_$nL%L%?f5dl-R^!F-?-XdeM%Yz*ushl^(LMsl^tEQ}r@ z@x>2d;3zgH8$hd_(Wp06Keq#a9i*W+sUQI)w=qqLXxu84iWXswXe<^DOvEe?;%`F_ zr7M^K4X+&6hreF(cJbs13by1IOpL}f2CU;G{l$~{RpdYT?jRGK(|H3aaOU!J$ZBT_ zxEY zPO_wA5&slcDb3Jka>MN1==VG34_FGF#+d@baIkgsO;F66abG%)+;A)kZlRqYXg07m9O zz8muPugTPf561SP8-0#cFU*NeU?r=NoeS5*wxQd8f!tj9KJj*r<3*z^8sxDYb3Y0j&qeaAYi7I?!&xZpy%)&`uA;a$bbBy6CMb{nbrgv! zs}NlI!4+kjqJiY}k{oQc=`k}u987;~hQPOxC7u@`8uHBG^GKXGMJUWC8QwafCy)Hf z`!HWg4tr+{%kzj(UdJ4)C|^*uJdH)haZj*~KzkZhG!C_H$)M`dZB`-1tl?cRH#@$=$Rs3jk^hM_9mmC?I9LpwY*WzE z^w{nv4HafSn|+hr5t|gvMm-VemXjA@)5vVda7DC(`3cR?Y}ATIw*zzk*Sf5tE#zGx zi4T&iSR=&j5D8bP&BxxW?I@X6X?A?6M=)s7r0o2$c2SoP&(y}NG~m8xF*!syVKEJ6 zYT^iw>AJjUmL_k~XBCilE6vGey0+b{iRb8y>>)#yW?!AIHR>z}^||b}R`(LWDp@mtt8`8D0LWvXV`qyhIPxmsfxW(QHOZ{IE;Y2aHf*f_QD&__I;o}g$;PI2>*||YTc2#@eza&@ z9h|&;5g&Y};TTUsLUhp8I2a}THkD+o%!zm9Yca9wVdVl`CU33GHAN0+1l{C3q{Bs; mX#_ih>8mKu?<47JX2*>#)0DHx>NT_Db604iNjzE8IQ74-32q7i delta 3571 zcmZ`*4Nz3q6@KsW<3j3tZ&`3bC9I0NAU}&j)T)&S3qIn~U|jxHKm|#PifByRaX^Dn zx5=iY*UZ>X#-W)^nsk~eb|+JtNezfKL}NlUT1-e>{J}?Jlz=8AQE$Kd-XqDF^v>*l z=kK0#?m6e)-L`Xl&so0B$y=UHx)6D5I;XE&(@7Mc79mF%tQ;ud)A?y_ioxECi+Mx7 znEX+AFZx;Qf!AWsJn>?3+1Jlqx%uj{juoERz>bjlQoQP4QUX);Elvvo_>QZ-|7hD{ zK>D+?$l0iQQ+Tl#?x5n^O)rQYE)YFiB!_9J{6j% zc4-a-fKmB?wV^ z@qOv`eX3Pe+@=%+M7KvIt2@!6sgQVlUxCBXpBkccM3O@j;hPbjL<_1KJ0ew#)k)1` zz_%+^-%-##tCdx1*F}f2T3Mq;^eR<#sNrL)Xj0)pT4EkV27JU z&dm8q2+<8A=N|3*oKGX?G$Uyl$@U`A-v6Jpk#f_VLLqgCbj?{KoYDX@e!dcl=1%t2 zt(TvjOh1BzP^{{sQbiA%(v1WVOdGqZX8@AyOo?u(n8YE`vZEtXd00MZHkXy%=1^il z?*9hPoX}2m*_kEUfC{hb(!I&rBI7)sM4Uig#Z;_MaCov#a&~D`jT|mR>D9XfPZdEN zG>FE4JT0?Bv(}(3-!+BAq%Mq3Lo)xmi3@dveW|_+TEie&n_4jGeOxK6e~{8Q*~QdG zKxrEXNmuI3$Rp7#MLA|Lh770f5>`f&t!YF!Vk4IH9qDCIZjXXGs{~Y?ls+o$-G-Z} z9KHrWxOP`;9G7vLb=jZQY9oR4a=~$#h_3j^-$dztC;r(+vRwZZ`mE%->+nMlKx#P5 z1Ohi?gi48ImP>p z$jS|~N<47GU~re%JsMg7*$?MLDfRf5JeACGN=bq59J^Mx6tCjDT@05yAYSQ7y2GlX zRec%Wxl0QUA# z+RpfXN_&Aj`)Lau1`&in#?xN(F{O#XHsAuF{0SZ zG#12JPwwPS;rEd^_Z0qbB+EU6cam~ z+mIEJFWk2HH=tR>P7o5$i;i+VIS<$r)+;hC?;*aEY{^^B?C)_PAXq~V^EDD*#iwg=xmnu4QL`S|DnXuf zpuqvqT)rqjBjkv&{*VVn$J3#1(YYrKzFOfcTWAzrg`!gx9Wa1lj%dNN!v(S5gTW|e z`!|b@b?{ul6k|~Azeh0@av77-4|zS?fNr!3rRc_kqO2GWUC^gs!kU;~6 z{g=oemOT)@8)WQd^6s*<@WUVz2Z^9Q8t%j-?bpabb$X-|_Fk}kjGYy}6Xet(azdRP z$&UOY*!m%IU7d^D&v9;iE<<#9h$I)}M5M5R!=#10(MNkn0(pv?`eEo=@2l8YpF*v|I`k{+Hd@INCfyuXKNz&o4gNJMdh5S&Yri`NK^8RU1xzu>dT+r{$*E`taq zYnX$zCHYA%ydG@nM+CN_!BH5)&LmyT>7@br0^v*TaUG_%agL$Q2G3SEQA+F?XTVtn z8#n$h&(Lw)pt-A^^#W@q8!W7^*)BekJPXfBXc3Fsj=yPa>B9${N4iSwe70`08a6-W ze+!vbW+%g?p)pn=jd*1XPyx-7_kNqoT7qsk2VbR-MhIahaJuGh4 zN2a#gtcgBfF#OV#H_kP7)f&KsV-TXqA;HE^Cf0I$O9Sj=ZZhYpHVsUe!|b5RW;24N z&pj~Tly{iD8J90K<2h`)IhSq#=QZuQd-AI{OsNGhz&{ef=^QO8-*4J)YS6lVSrfs!N+iM!Q z#-`2dYPW8!s%>i8+!Sc6Ue5p6*XuQ>`9RycV2F_CAWa+6A~kagak2lpAs^4c7aP*8 w<5zJU7n6+obZc*pk>)sN2e#Bx9zPUY9-SZ2Q)wt?kVD`;t1-CAB%?Ore?5IZ`v3p{ diff --git a/examples/more/03_ci_beta/example.c b/examples/more/03_ci_beta/example.c index 1ba8207..6bb75ef 100644 --- a/examples/more/03_ci_beta/example.c +++ b/examples/more/03_ci_beta/example.c @@ -15,7 +15,7 @@ int main() uint64_t* seed = malloc(sizeof(uint64_t)); *seed = 1000; // xorshift can't start with 0 - ci beta_1_2_ci_90 = get_90_confidence_interval(beta_1_2_sampler, seed); + ci beta_1_2_ci_90 = sampler_get_90_ci(beta_1_2_sampler, 1000000, seed); printf("90%% confidence interval of beta(1,2) is [%f, %f]\n", beta_1_2_ci_90.low, beta_1_2_ci_90.high); free(seed); diff --git a/examples/more/04_nuclear_war/example b/examples/more/04_nuclear_war/example index 6060499712c5334cbb8f6ea50729fe630bc7a199..6e51fd6fbeb0b56287d51e64233b474ed1556f76 100755 GIT binary patch delta 4489 zcmai14Ny~87JlzxFjC1&1PzGRsG$Z_2nveZstE*qS>vdcC<;O-B~q7GO4ZR0C0MOB zyJE3-I%_*l#m;U^ckK3u-Bx#ZT#(w)R#$4RKgbrawmu#4Cs4~;_3n4yd$Im(*E_>I z_kQQxbMN`ix%Y-`$N8?~e7%it*qYQC<+X82^V-M|K4nm2lsar|a`UPDXh)VhyjK?T zDqcul7mmbkG`($^8$a^hD}^&Ub6+gF6nVt(W>a;Z@bw?s>_%;0hYP9C=HIVO2LzgPHK@YP}xa;%)^6B&`H6| z@?tV7HZgMjO-jX4WNxf^(DAX%L3HuhrWLVeLj!*C@>yq|GuRgrlfU;$!JE?hi~3%% zX6Vc~jthy#j33QRF8M=O56rcUh!d7wA z1|<`{;=|XSc}mE^J}-+nphb*jb!J40$fIW|5gR8`Y(tB@4IB6p+r6-QYB=(vm@<9y zP!_c96aqiVb_+XO2s^&Q0&m5n9})L)?06p>CGJ(?vK_naY;wq_#eNj_O6-p+9=(SY zZ}8L}j{PfGoPnqAJ7RzHr@Ld{g7H1EleF=%fkn=RPOmf19XzIFxb)TnHxk`)E&|Dp z3yG!yR(w2Gy`1trrS9NoIWP6pBuT+mXI&1RQ&;_&-k%bt&=Wi-t#8%SuQ4OL-#x8! zRS)A6WDKw}@@`=NA*Y!99ZTfMs5QO7Yx9peK0E6}%wNyUqx_z={X)xWK>+Af6!722$j8XDLM-c~(8XgRP)8u7K6l zK4+nGsk6vgjLJgEps2vGJ9ymbWxST&VpQK?{CXEKh6C~`IFf?ehf6Mr$J^qm>aou|T2A_8?Zz4$)W*Jg_o{s~-O!vAs#1BfM|@pFXtz9r5HMZz2=#s;)ng)flP z91@qb$giCVy8(2tS6p&HRhb5(pqK2ESE>5JR=`K0X2-JyC-U=FdUA9( zDM?9;`dDO}F%Es$MCw!42}C3@sYJMCAdRWjsr68<9s+g77WDjLYOnb6O1uNk-Mxqj zmwqDPw4>|@*JZRC$O8Lfq4^u~XZy&gXQLIr(WeaLYx^fcTr}B|wtMa~Kvjno=)qw~ zeEBW*u+$Y=E5VPxdUqLAG$!j~vYm@%^STZcn zkEHB9asvu$gl`+wiA0sLf)_Nn#BZJp} z30HbDLZvJ0{e($_`lS%#xfxdE`p>7EC9 z1-IUINDE$-4#@g**pLo>W$2kaAHdGh>AFDar;v_|lo})92+{`0^-$`7tc1*mY=f+X ztiOav*!T*iU63F2!r?$H)JJJGWIbdfB=-ZQhu9cWkHm7~Xx0bZ496^AA&MXN+oX1mtG?8Hvd?3ifIEI>>xif-d%7l)guvb&caUkk|2@ zM_ODlBX=PqiEO@Pum}gE4Ek9QaYhU)fHjk=uo!`yWU3Z+p%WJFhdBJ;;1Xyf|$T3*z21o|mW}U=*IhijP ze2w`fUo5kKfiD}%u4oYk+zNxG6w5HhAt)-LxEDVQd^hl?m_19I&!A`aH-H;aw~9SB zAb^cb;ISyz{BtB%ni9!Mx#JvpQA&y2fUbO=d@L=G>|+KO$s_JZA}uJ-l#YfvA$jVO6tlk=h92t+?8Uyw4Jo|scJ z0Fs=UpMr`oJ#OY3nvxzLC-8qD^YdN-)SWku_me^1QGz3rBzud5wkhPd-WhxqIq01# z_@)qHVG(oi#KK&Y&xR$@o@Sh*ts*-Yngg$(aE73)EE82!$q2eTHT;8A6<0_>w^_p- z8m^ECyfchzx1dt-z_F=h`vdP;llq^X%^gA#m2}rm{mkM zXXs8LC5z1bEamGbTNjz*D`9Emu(U5cSCua?o3uAHKh5llMN`Te@=w73*SdSjRroxt zXyZA(kInhHwUg2L=J@lP&#uv`CVW7zk@?W-P)dx|D;cV`T$?_WbmyDhHcjimZ)VnZ zsDg89cvp%FggGzN@V2QcUZ&YsFmLg-nznYh3Wu5BLrja!fvuX>sPW&V&0VeSe;%uR zG%cOnj__JZDlV&n+ZA5rp1n0|xMhrL7`~Lg0B+=llXn)I;!kKc4M;d^?`aj-q3r*f zrfnXtYL{wS9VuC2jvt1z#h7$xv2Rk1xSJ$xiJ2e6Z1{BcbxW``3GtlKri6 zuJT{NUmz_@%(~6cTF8e&4F9&W#%|@?Mn)Hy)7qh9duq0*OXO%PeXrrh$trMCtMzMS zU4c1%2#z?LWz;HWzBX$#`4HMPXieO3&f-&Hy{63}-38wopDTH~vZAc2xU!_Gq@tp% zqUl)SCjLH@6ci;K+GaZ5I1Oy={SJF zudQy2bMuF>$<*mMZ946=m~9%8$wXGtL=%&sF$Uw;f>GCRWHn)odi%TYJ(5m7?hNmD z?(dv??z!ijd)eSge$Pq1!^W5IPCFlaH&DsKBUw;(W-d+RGktCS)j@y58@wa!eo(OWMvF!El|uf|<$J2vfK8n39~$WAL$YLga=2*j<>#;0 zOSg{HOHGm1uw)$fG0p>=@h^t5H0=|~HV~GU2P9jcRPvRSer-eR^sqF&8;efEZvHOV zNu|NCG_6mvg~QUJ!Dz|$zEtvUSgJZFmGvk)BdzSH*Z@Y$`z6~^spKOm{mKUB=6$xe zWIMxl%u_v`fTzQ7b1GUYIUAPBKbCC$QW-*+KJ4;-J$mv0D>~A@_#Ob|Lx)ixJhPvTf@%~rw_PBhVOesy`T;9t{dyhe; zUts_IVdva)&tF%(H{W6XaBy)CD!=zTO10@=i^%CQZv}~JXfYqY_(ROc5p%j0v=PDf zBGBGn4q6MjIlV}jaft-dtA$T40#bgxl(k~F+K+1b*~!!+NZ5*HJuFomMp1eY05%H^ zf2~FUk+fwbwBC#75NP?KA!&KkCCJF2v9d>RN;Jr$_u$M1ZCq~xvqTw?;i`b@%@}%I zJCDnX6Bw(Qiq#1YKdzT-fzU)PhOeP4QL6(_6+s*{h{8ZTZe)jMr9oM~FA5oFpPTS0 zf_eK6?$vwjqj=AU8ZVGKyEE}QC_-m0P^#Xp&=921nhPXgPmN6+%R)4c+ZIm-?VE+L ziLA>cLbZt)W^T$DHgVjBu~27y1F2rh2us`7<2uS+*VW4#w-t2eaaLttR*H%EXD$~k z{Y0FV9D7Vs{l>{Ak~QlsVWUKD&gxop60GL%1SB{dlD7Yu-EUu&)q|`V8nO|&G%1&R z2dEte=g91K`Qj9;aQUKU0o3z_PMsrnXTLcw1(sKmDSH+i>@67aWMh2PPx(T|-@v9J zik`Lxg6J8x4Lv!Yg+Hw}=^bEr|{N#lk*Y*T^ z1V;wW**YRKB!)aoD?cIParSE4Kn#9DX&YqwXOzAU%msD>J-w99#%pakPH6+?#?LA3 z19k$3fo&%!O@f~Lg3=v85AapszLS)80gYc$dK9?04=Dn6oT4-`0{fh%t<(mg;~b^s zz#woVF!wyAdw`w5e*g!8-9S%}(i4pDr*se)M3?49Vq;v0T|hhDm!rUaml04DHoijX zQ@~!}exUs-rSAeefG2?5H5a9KAhcbhbWAihzK%u%w%w$36R`F+rPqPIz@$;o58>?u zb^;sNJd8}SIdE_k*_%Iw?<9eIGv7c?sZ&RAt%CNl~r@Lj}NxQO3FUMXD2H&R<`$*m4Cr+rkhxHuw_7jue4yCT~1L`$KVlqVVsL<8*S;V%g#+mBNE*=~iR@%KWW znBL>uw|0~qFs20Wi!Io4b8054de;P?o>0s3j6e z28vl-O-sytn?G&IRDqvOic9_gs=s6^Ka0eaP7;DyWLjyp;F(Q+SNa=%5jj$tEpW4m zP*%+x)Rh&s2FI~|Iqt9qTnkD*77byap%H36Aik#a-S}a~hMGsmcUZpv;Cn0{M z^PY5#&&8aLvnsF2_>x6sF_|~OnS^+5yzH=5ljp z540?rV2Z|0*H3g@=R0O;{77a7*xp5Vb#1Pm=RavioCxR4Tdkf?wUMDq(|J#pX6V-K z=Q3W6wTQ^e&8^$u*T_A>*=uzJ9jf*9xsPdtdch#RsPmnAT}Sw(dVGu7?ZL>>mpoZ!w;oA^DhY(VM1e>U4|Xh(*?I z-6E*JQ9Cq4?jt?k&P0t^uk#N^R`L?rU18=^NT9;(%)k}=xytAN0$&cERi&3br5Mj@ zz^x{eE6x07)iYt=@qgDXJ(D#Pw{H0+*<5MP4CvZUyGHlvy>d+FYi%08MnCreIZ$aH z|4`Q&b)Ut$&p6T#Z89DM6S8fo)o6pR%_OF(3r5e%=9P_&?#7yRE7y1$-Oc{U6+3vc zxOQA*)ALPh$>HjIFdTdBL8nBhuITR`|{vt#DhX|qSjX7{B30S7fn+a delta 41 xcmeyrma$_U?zNJmoEdm}1>vqGm!}9OGk? zkZMobG3}(a9ouN9e_%IFOedo#l3=Pujao=h5MLYAh>J-XHJ<*?y?2dbQqK(cch2{H z=R4o`JKs5H892)~pW}C!c*oA@p5QH0Ii+oxN=EbX!Rf*3Fg@5kpO5FG5;N5CaJ2KP z+)myVI>HNHnJ}vK&!v|Pi`KnR`_en&+tJQqcZ1*a^4R`^<_w9mHO*!dJlpgf`l-7} zdsxb|u@5MJKrHNMZY#^%bRoG9 zt;e$+W2CObpuH5IHcJU|e_f$=Pb8;i?BqbFk5eulZrKBfL_mzY!n1QRjJZjnW*~r&(KWZj83MS}jdkbYWA&r@ESW zABA?;HQVM^9X*Q~nL{}Kzh~I0GVy>`_T(8pO?0SMmWz$&!DKt99 zDz)=;3Rf7}zCKFeaQ>#Kt{&shsT9v+1lnESSzY~}MRys2{j7M9j};8efyC~@wt*DS z;;WP*k6fb;-%&WF%#Eaj*ylnjar4ClDcS@IaKRSy`|_Zh!=OBusX&ptS^ zx!SU}T;4cE;nixS_NIn5oix$RVL3*rW6M1IyNXRYgbjFqfpSxj(DgU>_zY_N6&M!l- zQ^4;4N`9W43#2XaAG2a$`SmlDQang3D4(YUCJ%U0zL3dpVo}e5^W?*%1VMk1TuYiQ z*woScy`&&{w$Ry4ev=#{%<3galG8@q>4y25ZbFlz$NHmh+5A1fF%`Gj1G{U`))Y5u z6q%D^5PsfGDpJPh3Zod^Fe}EJ`WYT;?aBOl?^i6!N1l!Xd(-TDc-$Wz_HWtKc|Q*> zZ$P}%9<4l|o{lbeN6M+7znib|cuSO)nEBgxS=%?aMrq&7{Z)p+iF8?WYDOGC-tT$# zs+%AfI6-M6`iUKs?gg5GEkMUfN(X=grzouj_MWD+6PWk~rRf64?FQO`C(l3+G@qrk z0oZts(p|uo|4{k?u%e4PX(tmPb)6r_$>%9G0UbaauoqYYG+dz60o)CIoymJB-3^rc z5C||4^*PHQ4pH+Bz>`-fZ3njC(i>oiM?9e=!WvF<+;&7uXCO z0?Icil>^}jXa~x-zyxsMHXH&gyp&ED0X@(LG!If*$Hu_DYiJso zuF^l8sgK+MRujRNuiWP=k2l9*rX9Rb zVIFLvK{ysH>gJ8)aIUxskX4XP*07nMU`Ntbi$tLZ7C9q1Hk<8iFWaSIIstw`y+nT^Plc=nTQN3N1E{2xN66LMJJkhXTlr5r8S?@k(7)WJSSaidel+IvQXe<&V zHgzv#uPPEwS({ZxIC(E?JfBWJ&U%55AVJv+gc}z~e)b%`oBTREv(yOl2J9zy5vPYu zoiEC72{ties~=VtTEz%NkuB<&-5%Jr!_Ld>vLD(lfZYk7oeiqvu-lC~DJ3D<5l*X^ zxHWK|n7(<$0&$UFEH9>8#6*i|$`%b7Vnl{0XCOBai@yjI&Ic&0rHm^q@Jj&}Z~JA5 zh9-ZTXnHHa#(Y_rFDraa`YB&*e^nDj0nPcI1mAD z%gAnl1T8EOHl(?iFWl%4Bz=qH`9u=BM9(+6qnAt*gcU?ay(mF^~p+2cwi$I>FN@KO$u4{6G zc1~-xS@&2>Ng|xEnB-wRMr(2|O*OC9X3Zq~3-yr;HEpwUGjW`9t-#_B|Cbj*Rr+vE$kclo@TjwqH0pDS-wwhFViQv zHLW>G6|>iswN_0oYqdLH6?6B<;$nTI9w(cvZq|G%7d|Yel1nmR^ z|3PLf?-=K(UsF}Rc5Qieef^qx?gxviRqpW>3wZbX%5QnMqpH_WI6IjvSs5R^L90z% z9dE46#zk{#WkRU$@cPhBLTeL38#C3|eR#mVxR&z#0kZhT=YoccRb?63iBbM?bu^BI HR!{sN`Ga+# delta 3655 zcmZ`+4^&iD8o%%H#{uQNH-HXQ%8;mo;V>u|tY|zK@C^$0Hj3 zaki#|-A>2dV;yTdZS8T-vOA~EJZBTc3j0f@u9eAv8lzbfv}T)j_jliWv>j{jIWxcS z{=V;i_q*Ts-8;~8LKr(PbUOssuJobUIcY>6bJA$M;EWxJF{T4C!BW8~q}mIO`SqAD z7<@kcz4$@mJ9QaF57w_(AGdi~=*|91ir8}S-Qd=UC5rVcEwHo#w&8u#3DE*M&(QZ0 z7R}E3fica0Dg3}SVAo+KtGD(pwNAZDUDMy}lv5spZ#pY6Ddt29UldoL61*+(jf76M zTCGwmRbQp%KfY4)f3L0k!Q;QC^^HboT^Ey@}sGA@>|}( z`auh}^K`vsVIPgFvUT~Dvb@pzE;h@Tjx@^~Bhw~iW%^^uga8-Z&ry}VutJ8JdvipOhUF6Y`<4XxVXC_I*1c*PWKD`nAU* z)3m;@9UH)C)v)Y%SN6RvJI{A;H;4J&vf~urF;@3<44wwy=HqDD_t}J8^N#Eoma7m# z#ucysvfFht)&sa%QND&cPU z*8Fd1{(f(uhb!e2;`NVeUHul79moDfS3G&A@4Kk^zdBG*D3gQb{V4h$zGuvsj%|~O z8S}#+Q4J&J-J}1D`6yy`8bK=vb`XIM{{Ntjke|;e6R(ZZuQS$&&v2OitCdh}>*Zhl zdiwdv^dk&?&KS@7guHqHMd?QX*epB%+DI1dS&-Oy0DU?x*Bly`S4F+Zdlv0k+i&zF zZ;)C4finlRaosbxCCY#d*Y)Y%l<-vJJl-^%K&j&@UMDzwr&)INg)bU0oQ2Y-R|lR} zhj7pk3Ip+|!Vk?$gR=Zs6k6=eOMDu^y!H*Q((C+7^ACmd&(SBHo*8*ig!i3e>>$6B z;X{zZOV81-opWMuPvjvg(@UmPd*%)?m_YxSNyS+SG}pB?bAJLM9q~}-`~y<6oH-#s zx&;?g-G32&aP7)KXC38L4&;0hPtUkkiDM&lS=P+hn0Vc9oRC19vR)N)t#sDH{dYWV zopksVBsd(GAAO$R^gvx&8?t6;xfK!0DFx&Jk5oSl&XLvY4a6x};SEH6ubAfxzZ{}- z7QK=`8`d|gGjw_OQt_oBeKzF`eHF9Z zA$mD`zPM+E-kh^|if;(!ADpK3Iq5UjpGE?@t_j<6)IjkMUjKv^Sg=Q;&*s?0_fON_ zoEv>>C7y1W)#7QJ7e|Qsa-Fw;`)Qx3)!x9;+H^wHDv$pJzgGWNeX9RzC0u$zuk@y0 z<5>x3o(g6z{xBw3KD9yod%>uXyn>MJw;@9W=j-G+Ii4)EcM3N~Ea9*3Tu#2vSR3Sm z4;b4GEC==gw+}KF0+OSQ)d0tU1Hgih7;}h(^Z@gLp<~bk%a1eG3hX|?*iPW!zZiP~ z*!D5j@%bsnt^n?gy6Rx^U;$q9r+~e{eZcJ(mY`OgxX9Q|(S+>1#Mt9N zWeg61+kt0+V?fIk%)dte0mp#d^ns#z!e07Zkxkf1U(@B`B3ttA_ZV~WY?%X>hCGgs zom$0=b5TOB{lQXDZ>Vy_`c8uGiT?w?LOhzjV(=%+}&eVUViDDXHWW#wwsC#gg)f zxLmT|FDa^I(f7NG6$NxhaYFhYSY>l73?4}usy6`c0i7ZG-nrni=|78a6l&nVoJaw`)p;R30%ykaEi;mnJz1C3}gK zR46IM2!oJ(d}&*=f}p6GHk+qCm^zGrUqR*2S;hmOn zR|&Cn7d)q( z_)L4(b@8*0p#sTi>OnRl4lGOuhR8#T-EuX_l|;_@=GzTftHm=5v z(7%5}Q)_LF4??GR@+>2qVl%2IOg`i^h?8ayKWp;ca}EAp^PGP*dByDBQq#{U-BD{x zj=){zyAL|Kdo#3FWZmAtb$GuFAdBg6?HQ%5vAwaiwW+m!bK|{jtxfI0 zV|6606ru}NmzpF?~v;u5H mkrC@8wgh*zFabZ1yY9I)W}wzkmeLO~i(hNZZlD!SDgOn#Y&uN< diff --git a/examples/more/05_burn_10kg_fat/example.c b/examples/more/05_burn_10kg_fat/example.c index 0914434..6a71ea2 100644 --- a/examples/more/05_burn_10kg_fat/example.c +++ b/examples/more/05_burn_10kg_fat/example.c @@ -41,7 +41,7 @@ int main() } printf("... ]\n"); - ci ci_90 = get_90_confidence_interval(sample_minutes_per_day_jumping_rope_needed_to_burn_10kg, seed); + ci ci_90 = sampler_get_90_ci(sample_minutes_per_day_jumping_rope_needed_to_burn_10kg, 1000000, seed); printf("90%% confidence interval: [%f, %f]\n", ci_90.low, ci_90.high); free(seed); diff --git a/examples/more/06_nuclear_recovery/example b/examples/more/06_nuclear_recovery/example index e56830d123fa25065ed7d6a3fc46344d80f296d7..3cef207dc6eb81ffbe676d693ddb220b9bbe5371 100755 GIT binary patch delta 6019 zcmai24OA3YmaZyjY!JG-`D>aVX|TZt6&e+PCZaSn*hGw|pr{}i18PE6f|M}lN=^yGn+}`&yd|UIVTB*=xD}FkT4k{W(9x7(nTS_B*ZAS->a&|XeN_7 zr~AJ5-TUsl_j~ufdfknkd}k-$7-e!qaLSXd3U2@NDHnotcD1LoQz)Mn*c7C8o0<#w zbUwk6qjv9?HN1+~&;jAo@GoEJU%C6}Z*MI-yzSCoJ>QN!xb}^q=Iwsdbi9#{>dkyx z^J)EM3xJ@C*0oF!iFY> zi1FE+d#~>4$uRf!Rf~Ns?wzgrYPa|=r=*$?jvEk7xxtJ%i$k zGtxfTsqxgIbss}&Jl{@K;<=Zi%VQZu%XOFOZ1@Dj>oAexl_~qAwGOqlf*a zQ$ZtET3Q>VM!%LGii{2{20S1db;a~VWDK884UyvmV>K+8-idsO|Ah4BB>s1lX5PTh zpqI_IprBi_%<>OSp@p%QF#n%#fd`jO&#(d1zLZJB=C%AwR2em&PoX!XUgG!CgO>3@ z6~pSlwG-GtYP19gxF^7VFTG-UoWDWF=&i#1WSVbD42u6oHJy+=YN}6$>3q5vojhg} zv?yZw)RN7=MMW_a_$qobCZ5lu{V~a*X-S|YqYYK81{HFJ+G2|Mw<#gk7P$NGs^jPA zq1cH5kJ;e!5P4#g`M=U@u?2hx^~c)GN5*rU_=|c}@PJqr(vBYEsGTOpP2lTjNn9dd zL0jVD4RfK9rVObbe?|x6l7eFuuHT239XKTWsez7)PopssPzfuAO*5W3ii&n73};-Fa4U@ z?Wue%-LcQ{WI*?S2c&HhHsIcn*6AGgt`qh7K6W5#KI;!=#CZcsW&cW_p#@Ka;<#I= z%*)PVr*~*TO#NGl+k3~ovsXVT?uhVU-y+@28H!x~QMG9AtQL#URf}nZTV@To#Yrs~ z+XeXgfXjnL8iP5i#S#gsR9M=-8xTJn4nX5XwNoq}#(G|;7MHf7R(k-Rfv=C@>#HqO znJ>8SI}P|erkwmsw~{GH-`TR=5(W{RZ{Xy6c&1~ zvX_MX7mUaqE|_&;^9{x)*zmV8_GTl;0jHRHg;Ch2WMlVyf$UiY-Q`iqoK=v*V@tS_ z^-ojyh>fl5KO=Cn>gpZ>JxcZ5!(gHJYnS)BRD=^1k`)i~VQ0r2AcbBm8%UAPUS|8* zSj>)CW3Ua>My+Hu^`W!OS?*lrtUzU{MFoZnyq`Nu7_Vza7}d8Kza>i;!yYN+3VSo2 z$b0tkjwFRwmrX8g*nUFKZUmT)QR@62r@b1dx?@O((JDbQS{c!kf34(b|0z~$!ZcPr zKUtY8RFD6R*m%|3J;eTyS8PdCq;v}}9HlZ9UZ*f(B|x}gsKh>X=_iP#Urk&f+_^+O z6ITf@F)`kw)hc{}oR$G`#c}D%MV}i$r*XG9rK(K9DCiX*N+(o(unllK)NFem$Bz7X zKwK4a%*6J2DchO4O=z>#^}YVt(TMr%(1q*?Wp<`QmknE`48z$gK^m)OpuP@@nM!@> zmEu&hAD4R7nagZ?xU~a9WcF5^2`-YzsGJ`iV4{zFaU~pobiBLEkF3Ph*Hhq|J{15$qw9mT_$1 zA)u>wbtYVE7G8D$P7k$y3t>a}O<>nkAg-QhLAZraz)r z+C>4=Qia*wG-cZKF`K$z-gk~ROiPJ;<{Svvd57PTiPLPP&^zdEnEZAS9hznpZl0sI zX)(pya50n4z-%R6m9^-UtK0jU9hmRjV-?|CR@jjB6fQZVoBi+X-FdqNE^7cUDNIq` zpVGq4=DC@R^nb{{#IxTij@%ayu;=x;)yi{o?yqvJ9w?LtDYIv^&&T*Zzmdsvk3cp- zTK^(@iuDlZBhz(1SjdMkli>4FGChx^~kbA2gkUJlro7MOg&`Nby+?G z$>AOp6#z$&4#*ZrH>BeRGKXx0d=j$grYyhB`u$jP$l5{p8-s};pMfkILZFZ>kR7a# zYp79=z7CC2J#EWP#vSUnxn|x?@8nM4E2$^fyu1{=JHV&!!LH->;nRN>zbG^iM<&;V z&uw7I409Wd+Y>(j+-Y(1d6{b84JA;pv7P-PKlb}gb72IsdA$;_ImgQaaudZ1B zQF_reE4!%FG+j@2BN? z>HI13=EV!c=jiu&<5Ff~TU2QySwy8+@E3qh`yN)qih;h$3yYg{UX~wV9QeiY2J3QU zh7FC+Z(vbG3Hb|o2R)HLA9vjM^5+Ue-87P)Q*MI+tA{9@e-hf4cVFw<#pI1ToEgEhwxm18jGDABUPl^A>m#kYiIf=QkeVDU@8JcG5t-)*q( z)47>1C-dckuf!ko#a8=G_+kaTY8g7*27`4CreTVYps?fWbU%I`crEbt%-)T>Dh+yO zKLq@g&mI8?U}PrnU=+UbD$Q|E4{)RT!d2SlP7BzMX45q~BPNEu7;UkMc~nTP*Po?r5FOG}b!Vq)un=PEwTvYQf$)A`@fqT=V^xu;-2}@U5hG!^#M!)8uEJKy4w=PZh)W)ca z{n~VUG_BFD;*V?i0S#yWq_9UR6Ua9-yfR&ZIb~(h{8+>L(^OnpSv0TtaK-?~C2}*h z3BJ)3^$tbBWuwoYAbNh8xgZ)_z$AomzW=0D9mh=u&XOrIt4NaOEnCBnXbYr~%57xc z_-EL+<)CkDqB>UD!)U%m$6@m>vmtA#e+Ki<{+wWsvIQZp{zyFE0k_goAHqX}siUcvGf|uTurLSj9uRY%TR< zO+j?D)SU5}&z-g|zV$k;;q_XzRB5BmvQeHZn${7m!oEfO8tEFjXwFfonQP)}&Avse zie?ROhCNHt_phX;=KQf;8D5{F8cLeINyS%ij{Cw*(gw6?Y;DyBus0-M8I8 zr{`Cgt1d$OkE-|UAHm1slu@cq^X7Xy=4*Jpx_LO?+pPq+Qg!?WN?2*;pJv|$2Q+J| z#!;0P_d%a`s|w_4_+brinW6%n8b=?{(UoRPmrtvWTBMD-&gM$E1+59ZTPsyquW1dq zNw7@1$x_y3+_?Jj+M24mirUq6t7~eiYMTAa-{ceTM>!v`wrW$wEJwwf_05M@P3Fbb zb#<$^R@?=(P?AJ- z>(YXJM=KZar)J|bd&FB@yQy>Sf3u+l&faqliZl*y7eQx5EV>yk4+COEmPHr Sbo8+qf!o%oy=cm)j{kp+j#u&k delta 5202 zcmZ`-3s@9amaZxuHV9qayqf0OLIaH=G)7PpF*E{JBd#b90a2PLL^noJQDfE*I>BV@ z7-r0zb;cxSmT%ViHoKF}II^=b?#6(aL=%&chb|^Q5H*yKL)iXFbI(2Z+;h&YVsMb(gTGj_Et*qXp+>lh7cwt~?oHE1d($X_&kk{hYC~s_o6q5= z*ow5_k7)z1;SKbn@NtwU{;i-VTvK*@T_4##p8xRF^2c|s?AZ}8Gl=Ktj|L0h+jG(| zDKz36sjFai0LRGzu^cCL<&{sT^-)<5eKIE2zc0rsit;uW@gyhTyQe5_-#e18zr3rP zA$z3VLB8pj;=9!%-s+V)`wcCUxc-E^nwCVyJh;Nhac?tw4zmaJx*#jR0WDZ>Fbn;7 z*7vdHKM^o2x2R%!v)J+-XR-VTauY9W5*pIeBmRj||+-jJkDKHLzng8&tA^YT;RILf#ppTKM5K zeImAf5=*onYZ2Rv?O}(p3r=-K-*p?z2h<|2A4cz#{30}I8|gsrfBCV6IWh#x-S0@F z>d7mUndNms<8&!1!}uoFm2*^+ZpqsyCfZqDYGW3Kf1J!nK`nSn^e@XeZhQ}hZ=xf0 z1s!5!&O}=@mRZmG&-Av1C_`ehmwB&Wr;zul@_A3w%h3sZ2_243%3J}6+cK?M>MTfM zUCDQ~Elk>7kdA@8{jZ85uf~$^2IY|woyi}d%$U>&?^MkT@}r@e_1kD=OkBunz@wrm zXf-_*6VF%BD>11dSvr|cJ52CT@5M55Aap{wGf-5E$(nV__{}y@TGx#mED?XX8 zrdQ+BBjzDXrQ87omu6J>w=@)A!5^pG1bfJ@ztt4~1FcM$7W{X1=sZflNJ!^L=-(6E zeBv(??UvqDjuW5uB5+2TKY;pWwdT!dD7lgoD^knJ*!T1f?>^^=W{fm7uPxQULK{#_0Pyy{tD>2NtsVKg8 z&tkNz;227Yd>2f$0LhopgQhVt>3i2QSBcZaBO9)h)rr4SlLMll-z@8_U)PJ6KwI zsr{fubo9wPG#jpfS*30in(BwJxgnekx}TcZrrFWpSbjJZdM9U2R1%DNa~w%>lszQh zMfucaig3DP`Xb@!xmWxNvlKn~!^^bMnHKtbgj$rM&qPwU^C=-Rf>LuyXgAa0+#NX* zm|G)YNnR;>En?1?__H?L)ULxf)l3XOZRWUK+_@uo0oZJ&RWs@Y?#vtT(YGIxaKy1E(cBe%D_ zD}0~Gaphe>ognJ|=HGitGf0%xopFio!d5AZjr;GWSG2`dpfX=&od!DPS zd2-HGQrq zcR&k4qrmroIzYXkb3qS+?g6!(R+Iyv6=xKs7ql1jb5L$j=}>N?Q+QTU!U8yM=nF-0 zf)-v-6c6YgP%r3h&~8xAMMe2H&~DHZptd1JImhrziZTJZ=Q3g>5DO!6{~mM*N%A>p z1!|5h2n(a8uLJes;&=^o=vsrK97m_{IuwB(1ib_5yn!7BV`0z=Q0F&_vJ=#V`)(Z6 z3u>DLK4>|p2RCUe8-pHX<1rXCiC&tYj;pF~zJ*^yXXj_2&lg+jp7;>EgxmZLAox;ya8;c{YbHPbtFxp)F#_@CD zV7JmM#R+^p9Vt$Vj)lx4OaKj+6eq!&Hu}0aK6)39_AwRLIPAEuP-02EU^q<=mc;Q{ zR9%wFze?Lna`-{|=aOW>@Fl%pl9VZZsVL3*OlDCn&h-GWnfGDqfDMqjG%|56Li{1- zgMdWdXsg4{khS>jhlGb*rHlDOdZBb7UM0?#J}d-YqND{y;SvOnI*N?8Qe$k9(Nv;pXTC5Y7RnJLJ?L|#M$^xQ zB4cbG8)*4RABxo%wNXlZ!p?cL32H0 ztN|kmL3ckR0k{`<1MoF$LW$ATDJ(PEgfgSar77$N|26P8F}`GUf@PSE?qR0yDjkyY zgU8X`ca^@8vV$Fn%mX8o;{My%D=P}wkFEzX9@OKSd_^ErWgu4Ldda+Kv9K$z=aEJ2fuJXT2l23XuedH!SADx>SV#2M^md&dTC>|MKH~z=c?!N&(q23Swi;%B-E_5?#j~? z`k(UGYkfRvaZQN(J78uZlJkF-`?;ZrT&ANCFKMi$>ADkjyc^$@Dy|W(g8G(O z_$KzWdi`@#G^Uz5=srh(hRm;38CDqouXRt*l3EMdY6CmWf|WCPboDy9!e+rPu22Km z`Li~rujrkkjs0oNK|u9)s)lD_x2%ae?tDNS`#qZM$747(06DdqK?Ye~XUW~HvpM}B z)C2)avrPUAQ=+aY*C$1CtJ2-f`g z9(}9+f5GHvZ0@Em-;FPB*8H|{$S~%DUfng@H8wXHxQUCSN%dx{MVIM@XIa~-yD5Qk z>n+xNo!g49SJocUL!m_F*63WX?(lPRA>p8ZfB4>?IxNSfN9ETzA&*>C$qG^*xH#TqDv|&@vnx;SSKh+kq>L}kc zxOyl+xSB#o)?|nJYZxl)=$ca0zd37j%>HuWrwz1yZBB%DzGktX9t`d|vi2S?$dl>F WWBDPD290f`zdSaJ&!^mV8UF+Boh7CK diff --git a/examples/more/06_nuclear_recovery/example.c b/examples/more/06_nuclear_recovery/example.c index 4caf57d..572120e 100644 --- a/examples/more/06_nuclear_recovery/example.c +++ b/examples/more/06_nuclear_recovery/example.c @@ -50,7 +50,7 @@ int main() // Before a first nuclear collapse printf("## Before the first nuclear collapse\n"); - ci ci_90_2023 = get_90_confidence_interval(yearly_probability_nuclear_collapse_2023, seed); + ci ci_90_2023 = sampler_get_90_ci(yearly_probability_nuclear_collapse_2023, 1000000, seed); printf("90%% confidence interval: [%f, %f]\n", ci_90_2023.low, ci_90_2023.high); double* yearly_probability_nuclear_collapse_2023_samples = malloc(sizeof(double) * num_samples); @@ -61,7 +61,7 @@ int main() // After the first nuclear collapse printf("\n## After the first nuclear collapse\n"); - ci ci_90_2070 = get_90_confidence_interval(yearly_probability_nuclear_collapse_after_recovery_example, seed); + ci ci_90_2070 = sampler_get_90_ci(yearly_probability_nuclear_collapse_after_recovery_example, 1000000, seed); printf("90%% confidence interval: [%f, %f]\n", ci_90_2070.low, ci_90_2070.high); double* yearly_probability_nuclear_collapse_after_recovery_samples = malloc(sizeof(double) * num_samples); @@ -72,7 +72,7 @@ int main() // After the first nuclear collapse (antiinductive) printf("\n## After the first nuclear collapse (antiinductive)\n"); - ci ci_90_antiinductive = get_90_confidence_interval(yearly_probability_nuclear_collapse_after_recovery_antiinductive, seed); + ci ci_90_antiinductive = sampler_get_90_ci(yearly_probability_nuclear_collapse_after_recovery_antiinductive, 1000000, seed); printf("90%% confidence interval: [%f, %f]\n", ci_90_antiinductive.low, ci_90_antiinductive.high); double* yearly_probability_nuclear_collapse_after_recovery_antiinductive_samples = malloc(sizeof(double) * num_samples); diff --git a/examples/more/07_algebra/example b/examples/more/07_algebra/example index 7ff0169d2c3d4e4674dde3370c500d4beca4df66..d4cb9044c0ff0b0b4fa06abd75cc51fe76f9b26c 100755 GIT binary patch delta 4325 zcmai1dsGzH8J`&*E(B+G1zph)7G07Bd8#Njz6KUtlC{CGASjPT1Wl@yVq)|h>tc-2 z?Qs*+n?LesQ__0!IH!M*Y;&S{v>-9XB*xX46s;(t=mcM|Sj$o4?eEUal1P$z&zbq| z{e9p4?)Uxfckk@Be9iZK&F{7embw>GPeu17b872LCgb>wsJ19=*ybF(1pj2@N;1_W98C-b=H_TwZ;Qs1C;3u(U`g@wIW6V`2+#H7+S9ee<{~Q*WUDg|v z(uVEcUWfOJ^3p9z(FJk~T1DB3v18jO1`Il@@H!w@c1VfOs8kI}_uY_PYRcYgr7UBQ zo-vkA1<16`pL_TPmmFAW#pwH$Dk~w_Av%8 zsL9wo$B{iNpzmo^H=PxbO5@%MYSMS7aHW;)>kkPW&EHzh)gAnK)Z)2=K&SV*-P^Ae z-(Upxv*JNMb}%p>h|`O011ZYjc}kIo%WA+o#tLdDE@w5hU#^s^BNO`{s?U}L-o3bvoylC1J-w^FBX^8tgDxRAqgj8gdroL+-kOVB}P+MS`Bb(1BA$I9k>y+n}{r! zwjk~eEV}O+F7dy!m&1En*?gA#ep=ClA}EwgXDR&$yWtf9St!3fO9rPUN6#9^f-#OS zizSw{t-|Lf@|QFs%rcSO^v#)_CXU-M4(gmw(QA-0hol`(+@@6CI|M(t^s4}8A7w|_ zmvhrd&ZbugS!c-njD+YxN%cF{WFqw$?+e8exi53?BX0oJ9PU964ujH;Kf|HJS7~ho zKZe!w5h1JL4_KsXGRPyl&FLF!!3?J_;vNV!Kjqs~BsuH7xv8-HHi1$qCE8g~Ay1hM zkN7dZnBkYPs8W84Ovug@%xB1=?0JIHt&ZO8A@5_fv721aP8H_%kSWt=jd`&f=2ADQ znVvFX&q)xl?K^%;4sNsCPVbPzm;Ton^4fHZAa;|s>62YeF^q1QEqD0lt;V2AJKo_PTThbPYRC8ZnDcORy-1+o!xFJw!A(g0-O2&Hbw zwxg5|KxTbKX}AFUW0YDT1CTQzjo(mO3|agyO4mRRbWysSjZab!Z9}I8Irczupo|N{ zLXdHgt&ldztW%WEW%y}I%OP9NQ2G{RTQ8;8AO}#UGsCek3ceY#7+2<7Oyb51K(^wG za&H6z>8ErxWFurVWT3x<(wER#c#+ap$QH;>NcSaV5{dB@N^>Eb2jLKM0P+~5JBUEX zzyYKc(mezw*ckFHl9c~I{AzHN$r!bdF*0GrGZRk>sfBj$Axc%4snBq`xhL?nlGpRi zd?(qLpE$k~9#$|V8op;-8y+jUoIiPdF8KZxMyxT-0G<~~f_*Z-m(0QYBq_B|<;%$? zdnUh=ylqbee$bwnGJKrUdNv7-th=!d5)=W$Vu?d%;c@QHmEP)p-6`}!M zC!RQLl#$2>N<-g55@JdX>IvgvTysJsY`5 z%-uHTVX-(YjTdueF-sP0g`y=-oR}vX^VGoE6Rv=n@;SSNZ zE5gBi$;_7>zNY+~FSc1$R8=qK|X_ebtV zH|`?3Q!*@=`VofZDCAo7IdJ@}@DTjU(V|fn4e~u4bKe7-Jr~Ibj_J{M4DUrv4__oV z9C=YI(Y2ubBmD9`C1J#{xLoMU@vm6CB^;9UF3E@=`3u1O>n#k)JNzk+J|OTPktMF@ z1ul>Dy6)#^kf_o`A&^5-ORELB8AK?nW#;lzua+x{% zIz}b|pSwH_4ZWjKFD!;Ec738H>#^NSYL}Y%boP#ShHNaV38x+ybj!#Ku&HAQjef9Nkh~0 z+%D>P^9&7NtqJZXiz!=#6BbjW&R2-X6diZxXy*0$tU2VJa&!D*U29~w4qC3x$e+lW za!{QrOM?!YO2=P7=hSr}OS?(~Zh8z*z>*f^O)@hRbj;o7v&ZR1nbemizm&zq_af#wQwLo#+QEJ}~FK$q3bN(Z5 zjz11MyKGYyJ2PgU^L+rhHjswoT z$xxQ2AC=65eyO5IGWScxXQh;z>)Mk2QgSyIor7IgAMB)(0lzf8Uo!js(g9znWIiYr z-}Ot?7p2l(<>$ePN{`os4PaDpSu%ei6(5#ThSxDS2ie||`8?aPRrPcdp8DYCOsG_R z!7r5`k<6E+QiPB^YIlBPb>8sSRNki)0;u~VQlvF8%&UO0JMUtKqu=vVIwBc@6XBZ) z_QWuy^I+VFR5{upwT=Q0RVmH`koy+O)k@ExNnR)~Qo=oQO@ksecM~}&rVR;3PWN+B zZ*X$4&3UI}0di~bT81fAoS!MqUc0M{DWwF??!2ya_J+yy6!y;@wavY__Lk!O*C*LI z61Tswm*c#@yho`v^|gwe9`i8vobSYZ95JV8K^qXP2Z4HiIA|?oA-O>Km6sGH zFA}mCO#aaVC|bMOt5;1wJDGZfw(FF#tox(%WkuQY!Dg1wXiUe~;$4wBxMGuFR@vl543Dwp}4REE6M&qEx?Q zS4ET0EPoNYB1vo7-uaeIl3HU0L^%I8H9L zIZv}|^&h!dcHS?6%Uhrq*pnWxtay)|cfXUlD#Bg(RITub?7#7`uW(%Zi{B`gq~%LW8zDQ6QMwJ-4%`jw@=)3h9QcaTLZIh3r31k16O^*nah#;o zAV3dH1e#7$nhmV|C#BUu?hK{vY<`~7E?^hX12kNqZPW*${Q{-IK^)f&j00xl;K{&t zBt-#sU8eK{pl6WMAs{z|je@c9RZ1Pe0d(l6z@4}Z_kkR~0Wl%i5SR?~0PVm+d~=op z4fv$*Wb<3;l%9lOxsC5&C>#JCK!cCc?Z9?mH=F;PQf>_9K!wP;)8aQHC#l$pvyL$^ zNW_zYr<-&@+kX@pVYEJQ9k?}kJftVr$a~3|+=N&!JS^7ralmdoRx;K)A+{X(eeF9M zALKuhEb9cmpH$#|pRBb`=GT)Ct*QKOa>$whe%P9j6nB==dVL++-Gkb!0mtHEgKI zKS{FkbH^&MorhwK4&WRCt2SA*YznoCb2`TqiG@MqcyUgFXi-FSj%dmk<8nkpp1RRI zEZ&U8&!D1i$Xvc?*ec|Raju}3MbnmGt7zU70$-i*l`S|#OMz%sL=y~Pm?MUsN}VCGcFy6a!EUJ*IdD{3!V3E9h~yAWt-O3T2{6cv&>an!`A>&IagPm|hVr zP=(EbXaPTTg&cqB`OqeajYC9KQbKn^Tyu@&D9b|)XrUo5xu`rBVTS!YOwGfXP7Q5@ z*fdNcY~v%?k+)#FcbKHuX5pHn+<3kK6?TIxx8;P#pwld9%@DWzRBsRoD<~J%r@I#y zv<3soU{OkZ;7<@+92h}d3}<#H{dAhZ+euOJJAzLpgT+&Mj)a#a2#$0zy=0NV%_2W5 z`7xhK4wcLnwxtoFbP@BpxHLbhd%EUHJvGEm)ISQnQmFH4w!rTqe3>zx-Cw5VlC@Rr zp0m+QZ3hO(JZ`mq7g1igQNvIC7=L3 z$Y-lI3t{sl|EtKPawGAT1;<7T@myG~<{|KQk)m>A{C3^H1AhV82;7Q~Ne8rhb*(8) z(*{o6$M(%Ur)z!q%Ce!h9Pk^R=h9W0tJi=FMQ=orJwh};oic=3HD3J!K`hqow>=hLqsEhuDk9@oLuN4;YPGe2 zcN_Vt!Wi-vvt$uZ+@brbou*kV)@${M{hYrC|D3wV5zH;HYbzkL<*Sy-iroLreT@&YMS-GTs_3E1X zmX?(*?(*7K`0xALRkx2P-9oUtcF91HaLPn-mZe0jPuGIPoqc0jKJOzP4Mx(xEG{x| s(*l@5_{P*obFOCo2;Fic<#`8@m(LFGs?wBEqyw`9b=qt^;p-Fr2Vb)!kN^Mx diff --git a/examples/more/08_algebra_and_conversion/example b/examples/more/08_algebra_and_conversion/example index 1346ed0c297772638a15cd10b3f628dce8242f98..e34b096548118c2016ae0740f61c3c04433ed8c0 100755 GIT binary patch delta 4351 zcmai13s6+o8NT=Oa3Q#N5p_XhS#{Ng@K#Y9iCkH5iOYz}A`0>lK}|`e_^8fU7Kf;8 z$4#)OogwxS(lLonGab`yr>4`SJVL70DB`1PP*!~1fKgnmI%+)q&b@biU^1SWd;fF3 z|3Clv|L^?g+}-9=y!RB}Xc0{1yJ9bdPMXXqZPH{if=>vML)4+{$>tOJ@us=zaHzEM zDsCn33dcs)&-n9=3kiFgc8~hq`&;bJw5yX3^PVjMxuURt3^_XDXTn!eM2g7BYKY=E z+X0SiPKTW96T_`$U8PO@X$!hxW_bEo;>ZmV&Pws-WhqZ`-S&yGR~ZSY^(o5hrz++cy&lc?PBy>3hRU7kM2u(O3JR=4wkV~ z%NWZBC4xm{K75JFX1r06O6H(5?>NHk+c?9Yc^9UnT{p6%8Pkq4ELujp4-pTr#M?3H zF7iH&h})3&q|-`XmXnX2P4$&(xp%={$^EyAN5>Jx8$6x!=Y9efz3_DLsoekY;Mv^Y z#&}5XgpV2NTq5O54k<6&byC^k@<)ravC*y90Knjgbqm|luF@-Bp_Tn=H4&^ULb@dE?PNjIBA<*u+V|De(x%U}? zeXMwp4-S?58W6h++XhnPPp(pmJo1!t6=nvtk%g?L+N2_Bxl}Bbpt4XhC@L_V?K&+v z7_a4z7}XCMzmCO>VW%8&4IAU?$lP|hc9O!Y)lBW~Egd>BI-SFEj8f?-POod61~nob zM(Y7cMk^yaew~t|i1*b11QQM=IW)A@+G*aUyXT<@@&@7bciY zVZ_RSY_6Lsv4vM&KrCsC$q~+9A>pyb!Y`PZ=+tTzzCaExzGba)$7R18K!*m!Wd~Ii z?nLO@**@_c+*e4bbZ;W6UTmd_oL0iygyNms-V>bocQ>tU8O+SDMDoMhg;3kDj!kRMP?n`S^yz3OFM1^^(G?DSLTe| zjzvGcjqkzd>}7NH$g?kzu2-^0$3P)}e2LP3v5$q^hMt^ui6qBO3LQOy1!EYM8bM0q zwhDg=Czs-gFgl!Uir;4sBRNgjdtMoRae`?FhAB!-b;xY zy|Nnw)PKjS0rOvb z>LItg(+a0>+~c$UKlbH(=zvE%;KlYB_Z$huBSJ0S;-Q(6k??V|J^ z)1RO;Oo09*r54Bm$k~udrztIfto|>h>mdziDQ#loZtA25(P=@BGNcKN=>j;e5%Oio zPRJz4+zU9!4DZ2$kS!M}-4E&QrSu`B0d+n*5DVi#H$Yb7#@x>&F1$04o%o=bgAhm` zrPYv)kPVOneT$J9I;H){46+5X6SD3)G6}}`2Bqnc&9~qX(%_@C8?p{kHw+FS;~?t> zDP7IRko$?0F(s-293?Uit-}mVnDI=*Q%AN!+w>Ks3d~e!INaQNJdLC^!^rO?XEUZG zd*NXPQ=;K>#vQ~n8;nl^)0|Cj5)Osxb@N{0n4imda!Il^V%*2riw~#fN3wRy*G;==RO@gxK@JC2yR>pFcF59sIFLsvAd|hwYDcJOp+XLny z^mct@rrwaH*D<>)*d4%*`j}nTW4rmV`^s-;gX-6?Ys8K(C!b|TI<5Mo?ZNZ(GwX-V z*XIVr@%ovPK1tGBvh=38`pCI@!(6=%lh{NrCho!YwV3gO1>O;0@peF_-sBFn=`A~h zY|NL$d|BaZ!VmdkyA&6QU{P+xOEBQp>P_WXhAFI2aHzkbcoA?L@HNbyZCj~c$L!Yw zZ}!_G00C^w1illve1&YEml}K!-KwjkZC*m~Lv%OE6HMY%P9~d zQ#~tkY6Brj@4^J$M8X#t`PH77MNc7S%>6Tyz^T_xsk}1KdPD5h8dG zxrH@am>nYF($)FcTcsT%;|q+*-C6|u$E#?ZmRrAuw@g*>)vDkgu$YqdaKd8RAE&~8 z9%D4TdAf>MXtQRKPYR4tIhrnryxmaLy7He9swg(=?hU1E~NDQow(Uf}* zvdX`f^6@|cDq4=)>|bR_Y;Fx7h*k~#S4hEh$biSY{ld<{+i0sM8~jQ?i* zhGp$1s^FHAUCWI8I<^+SjlH9OuUV#RSyyP5ACX(jjN?3-)|;flXSCgK2R=cmp&qiZ z(8%9VJU${@3ys2vV6wZg*RZba&FZR(nv&|WnzE{@iYm{}q6U7#P?Ylnt18x)%t|UL zUqd2_?~7$MHD#Mio?u}m$Hr9GtSzrtx2~k3re&oc=19&lEuK)l5 delta 3543 zcmZ`*3s6+o8NT=Oa6!FydAZ;tth&a<;Iat#s4;S3!Apb&;_?#Yr65{Vf@oc?|Ot_Xfz_@dk1tq%dGt zBuqFbOOdXa&}|_Na;}^sXUq0%#nHP;aeSk!eU#>SsC3@(QJROsE2Q(jCb#`ArD=|T zrtJK@-Y7e>L4Ta?T-T#h9QV?iKi1K+BeE2=NsH6|O9KmtG z`6n65(!5V4^Yu|Fzf&?_knEQv%fGibM2||*9oTdkddm(%Ps;8emF8WL%%h{ykzs$y zd{nZ3H7XTelXAM0U-+7oPPY*bpp-v=vq|<7l4Wo^Gjo*rmdsa}$0F6zIaun3nT!6C z{pzTccTzGBNI3{0dMw@XrOk2IUA%gnQV?K`YiV^W;G{dg!U{+K&`l|F>WiG-V??@1 zr!*f8ACXp%RY>(?z~cpq;|S>PHFBZSdEF?lk&Bdo9=TYp4{3@DZ9syN)5Bcw%baXa zbKK8fh1}}h!?3M5jwz0=bY~kAN|9W;JDk~t zGXC}(O10&1y@;x_R&ehG6V=dS9=r2h%x4g@MGIPwV9^%!JwF_@7IIQ_hOov>_C^;8 z5AOhG|7aBi^&RZfrKX>Co_d7kaK*AdDy{BDQMwQSoVm?f1RzP~{LqGAwLSBWj7a%@ zkHXM?%}4bGT~g!@^(u@?PkjwzW=JRZtC~0xGF;fHTGP9K;anODg2ILt%J2Wc|%JOYd$gTNHLT@9Of8587Jj6bVqtCtX7Fl3P zo8o}LJ?S>3huL*?&jRIs`xe=2nG?7@RK06qZ-tQ#OOr4$lqAOxVS5O<5c5im4#Mrh zh$O)uRD~2XD!o?wh*G(I2zGGo-VWio*fXrk&V;2Qq}iG;WZodxtTO^*gH^kevqFep z?B9i!AaW?S{e{~>UV}G5!QhDW+5s3;oQ0-pWKCDO2oXw=iQES)QuQ!6M{GyBb8-oG zq&xjq0@U-l0|v;k1%G{ZHZ%ukP)emp8w<*7$;|kEPx4b4|8qEX*Y%UX#l;9ey+Q88 zEfaLwYQhbY8^28Wbsu>vex9)JIyn-*aKih2P><;&G(Kv|(Q8OR^FzKqL3XA*N_UJZ z&iO+kc`m^ySo%mw!ZiDPB1<>66e!MR2eH6iV)tmU5?lYjPvqQm=dur@I8M$c6sGc8q7kjQ7Nx2|DvMr7ge`;6Y$Z52YPI?hK`wz+s@9$s8u7$3V1RrF4=H$Mpk^z!Drh9@vVcWCD8z zC~XA}U#D~kXuJtWzHoet(mEiAOK=p}0rd5Q9B<2XU?VUd$le+|uzHZvDxeil`a!lH zilg)jh|GI<_55J~SO>HYQ~C#BE3libzo1k%0c&6;DNUXmaS%C)foIM(LC?U18IRdY zeh+E;r)UBAv zF*`Rgt4WqEWX31R{bEhSD)d1kZKAq66ue>5WSb_0_LBX!S$qrWw8iipAp+{xP$kmh?4K{Ic^Mouh(`O1VF3u3GifB#}jj3XIlBiD+brZ8PpBmz1adZ6Ff|LC&%mv9%h@9q6=D z(9Gp7*t80Apz2eesOI5K==2PdoU{eF?}!ld99YJ?q(*ViFuxq9o*USVk_V$9TIQ*2PR$p{1BA~tM6`LNiP4%g z5o9GZYk?+ zj5yMnYv2=Aoi)(;KKn0`nRx~>oa-AFBt&q!5_RLZA1BGmGeqq3*l)uhKQ_Fr|3}gU zX}d>S601qQr|xEc6E1tCR!>gUri0>_9)4n+wyRtN&L75uMo(?L{3RZLde-R}>bnE9d-QuecD(UbYJBoeeo(|# z&@2Xh3ETPIXds{E8~k=KO&0Nt{T^Gbb2XhZPaXdJt;QjbM%w~SW4FiH39p9V6)%-G z42sp{aj)Yl|H)3#qxZkB!gEOHYJ(qp!i{QPe@seO8~76STSsOV7&2akOzqpJv~W^9 zQN8ZrH(E3X`@LafzlU#qT3c`OobxZ>&r*GyBuxbd{*vnBE_t`WAUyITCkn3Xt4nK2 ztE$SYino?-s;(-paeY*{n@nYa`J~srk1o$(quYQ)15T7)doZpP1A}*4OkzbLavGei?xmtPk-m$yC!N<&&=-c zoZt7I@0{;(lFnXyp2!XDQn{g z(ZOrHgS;gi8}o~0!Ql|zjW|cC_0DhWU1(uUn)OeCzZ7MrDZ)JzeK1n5XwEYemaYV@B>G|pYcn^ zv|26H2^Q)I!gL`_(OJJ#`VV+ucF6_5lzPkIx#09%S9Xn1iY}NDSQTYE#?H-$3>X}# zz}I29yiH2@g(~E*H0icnq^9h-QOq)S>ltI&kgu^UGLIakWuEpa#@H^)yd4`jj_qDq zKQoYd2d2!~Z(>PXyO4M%+b#U8e(ZRVCEkWf|3%)Xu;YUWl(a$3%X;Eou<3r6o_h!U z)!biDBNl(DMnk02f!x1>#|1>{`Cjh3@7$mJYZyO}J889!aW9oiZ25S6j==Tj<7)5 z*fLgA2jz0PLavmnP+2G$6crRM@SKv1nXWZ=nbadp-{K`qVYgE74GzY0IB(0P4N0o5 z4ht>xHXSxd2{{~=W0K0>)3jgLG$b48Fj<*GOjagz+5|O6>!Va1GgepC{m9BZt-AhS zVLuAj>b6JY_}0W-Uvt+HMZNEi8r?=eF0DTNA=VFJfhX2ZTswRUxGFl3~k-rx?uFxAo0BU~9 z<#Qy-{_*S-cwUaDlu8K>wy8j*bjF7~%r9d6H7u%>o+INj(gjlw$sv%7_!VW&53Mg-ZNX%@t@ zq$P7g(e2SpZn%{>y>kXJsL-zWWhc16yQkgqLhqb+rgGf9%l;qs(>YR%kUKyxDNI&> zK9$v{ed*I*3jb41E{}`F*>m51i?t8t)v4|JJmQkgZdlZYoLRd!<|hO_!G8G>k=yw% zN?%65{TQXqz;56^pz%1R?ZBK4N?kzvSCrlYHUlFB_@AKE3d}i4>1^P}QuJ{CO-0G+v-|FVK!M9Ra#fq_cyu zFbaP=umw-&UWR!2I)UcPlv+Ztkt>uo0K2c0!4aLDtH=zv7q|~->_=w6c3?y(#sie* z00)N<5a{w#dIs1IG>pRhVM8tX9_z;q3$`(g3kpHkOfW3)D*+y_1?7omZ?IFe?g(+RSTc*{K&;6>7KK-*I*fDaKNw_h6nj_@VmgTV*czrT_StmxwYU21OC{605!6~?*O0OM{X{7JoE^< z`};^{zAbbF-6IOA%YQz!9;KRnm7L0-61ofBCf1!CZpPF@7`9v`Hs|S{qMtnD z%nWy6Xkn$?Pu_6O535IaJE}dzCy)MH5HTz&6*{wg&o0^!3?vtxw(&bj)H5c&)0h0r zR6!^zARjK7LRy4i(p2;uY&}I&`ROFAI6-L7A}Pg{f@=o(W$|3p!O`NG!r*iwlvJ`v z&z3A$HkgF9(OzWYqpd^jM5D6yrf6gJP|)3`^Eqi6uj+{Idph5qs`09h=pH7=OH6TJ z=~kCbv*uvTS|5>1Oc^&ZiV|YDz<)WKjpM>l5Nr>2_Ew^)cMUj^R4+C0Y3w_$FyLcm z9vt_CzW);PI(%xF4;fqvJ>bpcBO!+Oks++n!u+rkE=QY> zLsr`{l2~fWIHT{NAyH#f_1v!N{Kn}TzfuG44%<_Ph#+iF4Z2a%8% zk4jB(i*&1zJtAo8+XDOz=_xh2D|KtLel|YUhvUliowVcc3f9!O4_F^qi~XRWDH?FE z1Xg)qZ(g0>o1i%c-kM$Djof(hvt?0ndvu>>$i*7_^=9n@b^mR;b#b!BR_N9a)p}94 zw(AeucFll=UP0li->DF8s;N+6x+R%gF1? zP5c_R7Qcmk$NgIO?yeSK}cZ%6rde)5AX7a&&DuC1D7uc}!^Vk&P-)s2nS zuT*`dc|rWNT?fm;|r2WbV{IbY9D$BjOICFS{M!aqMVth-DzP9ZWzq0efg@uc;+ Gg#Q6UgkCBD delta 3527 zcmZ`*4{#LK8GpOEgd@oAB?nvrgd7s&f+4w(M4*zIO)g;rmwJ)Qe;^5v5N(7rAf=9z z3kH;{o{jTpM?03GjIFe`Q}m{7bUKm{RM1vVAQUk`fCvYT5_42T`SbdFyLYyt-I?9r z_rC9Y-~adS2K$BG{X&;b99;dv+%xgpXA$k$K8q#_neic`z6==yONC5fW{0A$@5Ed| z=jYOw#3PfkuG80N26JljcLdM9DT`e%rFmU}Cu7QG(QBbc5>vzm3;jHCM&jTULPq6; zyVCG+q!~lIEp&!CUu?9{m1f`gr&9@0{S)QXz*Y0(X0KAER4NsUyF&H#FIRmx)HSzD zeYe!!OEGHaa8k2;!jKl)d4;i3-@C=zKWn!t{tD3VIsI#cM%8zxLICrNDIrR$Ni4+jmjoN8fbq949JMut?3Wo+eebBgKBs>l7s{!`>AR$M z_8Aq{5C27DrHjryKCJpaf2*KSCi`7|goJ)^gE4(M(k}5n^A>Lh6W!2bUid#TA4SZW zde9~W8$_VN{~xp-vY1gO{xU?%GggWD9AqrA)h zybFo6$X+9MU7)+{De*$0W_O}srtjGQBH9yZdv^DOKTn7n+zAQ>N9AYsz@X}1k>)|w zj4cZgp`2Pk-r|vJ$smvH9;bhz5j&jz*sv**U#RpP-9GZ@@Vs3kXFJb9K@=FSygK2M*@y>HyLvryk~ zmhR6@pJD*5blwu$^AvyaZKp4+`seMG=b?gQQg4jyMLNrb$gu?%3>3A73r{0C!efu4^U+rifZjCBL8rx`m5>^RN5 zYz)K*7ilzx` z=|smYVUT7!vSyl2F;)yIDV9u6i%TTy4<(Z#8Oe08UeR`Jb)*U%^j8jZ`T%r!%>3a# zhaoq9!q|^B4$-_2v}Y4N?U*LKNUu8{6t>XP;$?!9{;c=`VKMD5E}E!9+kvJuoMUV? zEh|p(I;4f2tO_L^3@u7D`F} zn5B~SIfFy8Jsqn^_RctmR3LgJdzoZIfuI1zLdnSM8N_}A2BM1Z-ym7nz%u8DDj;y) zA>io~Ovp~~TR4sKbc$7|`BJSue4`ceA;{afT$St)85e?F2)=WOKCuDeNlononhdR)?L? zv0kC?s`GILppG|?Z^fKWg zi_jfNe`t;<4yp8a%kQB*qJci?UJKbd_Z)%H@fE3}Cx_0hSS6DA^d}YH7xL)giUnd< zHWe#ZaZ9T!OVWGh@&E|Y3IOax?-MZ@t_WW%9Qb>vP-RKp7m<<-z6Z02HGlf3QfbLP zhE;-SCec54x{Q#EkzIED*6>H0ec*MzONgg0z;Zeo$o<-_T!WDNF_-+wp}kcWp+M7F z37zk<{}`QFZJ{GohRF$HGBGx28^8YeX?e9J`Q?ax5B?PLCwc~cOSeMW9g$kIb!jx~ zL)>rPClTog?p^+9(+B)A!Y4Tzn>6V_;xHzO^mWlJOr;4m7H=!`Ok^stH$)U7a4*+` zBxuta+Nl6p65$uz`%V6$h`*p&=WBNnNOwfP>yG$ZbUuBxCLwtvXdZ*9LEjtwdg(_s zme_4vlSiDgD`Kl*j;^sf(uY5Mt?_n5!!uvk*by=HepDlNAWD^vfZ{#5)9X!*|LsXJ z;ShXRhiB8?T1zay)~#Aye?S{+Ekc9#-9l&9S<0S*Oq<&odN{?AsGg7T3o>;E{|(~L zt_Xi*w!VHelII`6zgzS10o__>5l(47uF%)(EaL50`hMMclc%w*v9-0SwSHscqn_5L zw!mLk>=5X>RoAJf{zKckrcL#W?DcCltZ#mJZPWTSP4y42-_+E$xv`aa+BU3d+PJa4 zsjY28TcEmOnec5}dm7&s_)!ittQm|E`>nKaU1t2296d;!*(cVO;Dq222n)WOs9ag_k}MG;Ye11l5MtCsipIu_{0z6vzb@h@mZ3#s2vK}-2S@uqPB@a_X^^g- zk;`Q7wUCr~&aHUAQCcpW2Bh_gy~%`xB+KMqYgq9~gJeI}AeD7DNZA8xriT$+$v@=VQ9Y*l6DV9qqC<$mdXZ0Qe|(0w4fs-RrG^A09X6r%JcP_ ziQE_5U)=}32U||KV$)n5=C1a_&2hLXI}wt~{{aiEUT`WTWeqyLrxfp3>gJn_p$q0F zW>s|~#!Aa4CJgp2#@8OXvO`MQpoude{rI|EriJXi=H?N5jEJ#q)VqiPk;k0}@?-o) zLoM=79N-Ha_u`tV;mEsSGQ+8z)P5X+ALFNmojrsT2YKKfSo9yneH16&jkw2EX%g8F z-*z@JRAa>61$!;_fd{nM4`|-tsXH9|Ay}M(r~W%)f8)m8vA=@xJ+afT?ZZ7!$O~k* zTvqHotlepnF# zWQQemdCw~?9VT`JD+&gSr}sX0h3g41BGb5d+n`8DmNL(B1?TGWc#b1_o*x~seN^JvL)mw+zszn$KTYUy<@;6w~uw6H4ZBVO3Cy_o^=J^0W65 zmcE@)B3?O5h0H3^$6>6;$W{1)IDwF~aF=@IOxO*meFM_Mz50|j774v@xB8AgADafc z3)8$jccCKphoq{wqy%2(<-9VZBGKo!i~DTmF9p9->}SVZoT}MbUCuO=N+p(~QY`8X z-8H&eg77NB$>0yKnthkrqfb|+nc&tA3MaF7;6>11BD(UyS;Nh6{N6P@+aK{S#oMd8 zPSf2F7C*2fo)GoJ(~N!0-*mMLJ+5)4Ls?XO!$Ox% zSeN6#?w!yQQzfDIDZA3s`it#C^tv6 z_0yl_&Jf#AXrnj#=y{Bk6Er<9Q+&9OKAtya$dVH<@9m+l=4Cvvx(5QZd?T!%jMwaz z%R8X>CJc|IpXa5DCwpixZ)92bFfKRD7Ad|Nofs_E%lkqSxbVBDSGmhKWBGVOZWsH% z_@B<3Zn&%jJ>Qj~{rOa9b^EI)$}w-}&lgZw#8I$y3vYj()u6SnX3<7D)q@$WA!gog zkA;yD5Ana|4`R0YbH+BKA2`TZ0N4pU1|(lFmM;Rk7;DBj2n=%m5Mz!m?7%)~P!4oJ(&vG#z|BB%lCigd z`M>}$06fh3J{$z-IE4d7Vq+xlCSWrP;s7v+%1?@d**Omkcn448QegfC#&!bRE;2@< z;ph?!fi*w}uoI{N9hVua0S1BF`1mSgoj^+n{)S*9;0j<6xEWY8z}Sa;3>@U+Yp^rH zKRr8p?1<+(8OyL1kEMLKaNZEgWJmEQ*H_{N^AVz-^~Lb&eDcQ=fCJP)c_= zbA)Ev?@ST4Q=u?rW)Rn14@=(offS^ItD7;d1#g_I1ou3+wHiY-&yC>DbknNB5rT=X zEu1T=C+PmdM}!Y)V$m$&E$S{RSj;uuiInI>YN&keBWBAcQ86dI6fp;pbeNM0&6Xmw ziQDzSZgW3gb#7O5+wM`=eHOM;Fcld|NF!2gG5xGC$>TITUW%S$p4KwtQFCcTmSCPH zn;o**USv+4ZBCkPw#+v38|Mp4v2X)&az2*{8@wvQ;+2R(bE+>=G21sqDcqOLeL3N4 z%n$kEnXnDMHY2aAp1^>tG^f^K8$PiU6MdMtm*)jwM*%#|?fEU%m`&XN4ESd7nmrC6 zLXBzQZ-U=@mi}tal;~+F4c9r^H77fI1-j)bjaHUL_n^D?B3+}TN1IRsv6s*t6O)Ri zlQ1m3L_b$X#BfjL7;d{nLrPwZ6T{ar3|^vn#j}U@pqqqRit@|DMn+ImNx5j6;(xMa zeI$^cnwKrKQS0M2(Ny5id^}wguF-jAuY&6@8!tRUhq_ZlGKFTktHjnx^cU`#!b-Z& zJylGdM8ySF+{2R#iYijmu_xLyTy(VOkTvnhslH5otX&4WFB^R01fAFPK=)4uPYhnu z5Bwhb#R6OIA!D|cONkaPz@9Llt@GM@iSAV@V->@!RZJk^|D$vpoC_Jq&sM`0S{|9py*1l)8!fNVeG_ZHZ}$(t zr_vo_g0P>CF1ICr$?b3+lCQ5vLTK$AomXzl9WYMNHd<%-OV67KGrsHgiTZe{4kVuT zW7Rlgwk=Di!z%OVny+e|tn=abpp>slUSiDV^9Kb@duM=uqt*(W=Q(4x)3_nOW<$uT z@Cor-hd1r51-69u8n4kr)q%8zce&?me=zvI6x}fV2K^Phg^ZwYR#=liFl^c&7jHLo zAOY?C-NtNlhE6XwZtDUS7TSbJTD#Den}~Aa3J^nqjfRKO2H$36Sok^59l{5O(>!NU3A*As`i=274=QkjkQg+D^}F6 z@L#UnD2%z6C1K*J`qkCb9o2PD`&U(21*x{FsdjDkZIZTB9~{xxw6d;#)vD_Hrlyrm zkKJQ2dwk`Fk@wTl{9rLARngPIZUWM0<{gUk-0j3Y2^$ zPjTd_uFFeR*U#$eCt0q?YWr=O~*1P5orCF|jWW0T^-lDkkAfL!~ujw?Y zuKQVyCrtF>kRpXHe3&!q_Qtyv(*08vQni2Vuq2OKH-_VY)0Q$+q$yuZ*1ln>s9myN zlN{ejwtsDEm^v&?ZNs8(VK>(fJ1MVcSekN8vJMYR#|8r>>v75P)38*0L(1<^-}8@E z+r1V9fKgGeWIZD}{wdi8HZeEHS!~JL#UkE{WS%a;(@D6w9w<4!8$h2RaZy0yOk-WXfE4zTW#zxDfBYpFL;&C{rs_Wlhrlv{U(cBJ>^M3R* zrTTQRUgV6N-v)`(&~r9Dlk)}SY}1pLk!&Xt?fm7W^^~Vi%@Njn$+oF0g%Ac4XD@}K zzKy*)wDPl^XN@UivSbS)mPaj%JoP22{AXUGpY;f9vyx9mvH@#m)&1Pgh9Rc5kts!%Zk1 z+Uda4vL0--2L}W3g3LC}jt0l_^ias|UaG6xd*FLA!Im{T6bkS7 zTa+GUx5_&mlK0c^$u`@h;IrW@MS0ZZQKZ$jO(5YUF`fud!pO<^E%AGx+!O|N@&}MA zr1)X!ow_HKDo5_Z53b!L&?5Hq+a2PSIWF5Ba1o&G z&pW1<9G-r7&O}%aL{Unm=yaBp&y$t$0Wb2IjK7DV-liV%$HaKy*Zt&b;#?t8pN;4z z3zOyw`@6}Gq$$EDedJ)$j1kAXVLrc`+)s)dedz`Y(D<0IPgdL+PqJOZs(V_UNNmX# zVO}@MO&;f%EwXZB$#T^_*NzG961ztOo3ZrSK2Zv?-E&XJah#H$<+{SI)sKopaXrd| z%dbE$&5nD_s^a~9m*;4*670zxSt=Y%UCKwj#&KKcYyE~7_SoQ)8InrH(!h`NdzF1e z;f8BIOKB72&U2Ku0(*cbf%17u?E>%urB2Mbi`kA0yIIN3LFnC1#Sj*0(Szr z%ara1I)PsSTdz=h4LEp}HqZwUny*tD=?9_;0Qt~K1w%%J-|o6 z7JU0t0x`c!X%(>TKBZki`yd>SfF76xlz*VKkj;U6*n9{HjUd8;$z%6)QW}pGxbzWd zMUIQblf;CtN!o(R4Uwp|#xn{d0#6+t8N8j<*ci-aC!R>0f$@{bF|Gwq1h^0&8;m=| z5FCC_w86NOc*Z}&rDNX3aRPUZ3@lj4_mG72@TfU3n8bprqD43rESeUM)Eu(L$R%a4 zIDn&1(JkykH!QM7a`>^~wu9>=yVA$;56GwK6L}{&pB~S*5I!TC?<6r9(evbQC@qI2 z8!RCe=UoL(tA)W8gWCh{H_Qlwd3J*H|CZEejODH5!;Hm(=?1xzv54;?$(i%{Gh}1t zf)EvUhj9Lb=z~?{UPe?yx|q^9Vu_gRH<}kyazwi-TGK>JrWlzf${C^w_H*!X=#2>U zjEBjYqWr#)CPupb7KxUf{^_E1M}Q*Q8w1nDRKY3Qb3`l30s|PPi6(YXsaU@o>kHWW z8qu;Eo*7@52SwL&it$i5!N;OGE7^h!QEn6pQ0YaYtmqD#p|69!mg!Z|4wWed;(q3@ zp9HHj0&hYry-6}vTVNzQwCWalOMNr231XL*+)*b79)f5YfH=Ws!O|W~*PxxtAhVD= z;L|QB!CFjZq*e!iv~=qLsmPj+uMawvJ)pQdWM5WV&`yX}T!a9RGV(n?V#+BL_9uE) z<<$EFN#7D1-%LVZHw)f0&y?3E3;cGn#PK_Dy^fc8GYQIz7FrU?l)ROKeL8tJZyujb zPUX!M+7gJ6zmj=am7f_mIE4k~IBny=Mx1&W`eduYYXt+pkMITNgjPe_%!-ZSyAWR( zynU+9r(({=8KUHy6S^=96~Z~+e@9)%arcc?rTE!nL%RV$f$SRH`j}`uYReR%jmMG5ojS zUn(2EyU|8IgZ7-E%}vm?zM6EixCujs)?!pl>q!`fpkG+d_C$SExei!kc^5LvK`zxyp1wXH{4z8DG0QgLu;HiqWnI0>sQ+<;Z#DYC_u=c-Rt*~3mg)Mc zCd2n5A8)MkeL@3B86uk5JL?tP9(G>OXEjsfU+ksVet|C~G0V*XHPBkfXF@pd(qeq9 zy**^xa&yjKq0>74B|W_iBg<0;@3iRz`@v-6y1`o}>+|(SE$)Dypv4Fzy@)vqw>M0i zPa>hkW_}ilEAE$_Wwm8hRpnJB8_L!@tIBIV<5x8Co)ar&p4=`uYh7DjS2D|9vbv_a z;*B-s)vL=(-l(oCuiaEu#W`zhR+n$sP*Pr7TT|;fS6air=xpw)qx``&J$^z!4Cz~I z3*Mio=ZMQbwlWi+hO)}|P~SE5p@ZxN>rK;5ec0mZt)#rrIg<3fIU`6e*0nh#X5Gx7 N6|3~wShBY~`hRxELYe>o diff --git a/examples/more/11_billion_lognormals_paralell/example b/examples/more/11_billion_lognormals_paralell/example index 2850178cbf85dbb001fca02070b1df5eb8218005..adcf8bd219a5d46cee2dab6b5536d443d9b01b40 100755 GIT binary patch delta 4278 zcmai13vg7`8NT<1#3dwqcL}ly53|8UHzAK*AQFl|Zn7a8S;j>2eut3o5X%}VBOS^{ z5X{c77q3l!iZvqAs9-h}@g_utjn#nu`&8Kur+IT#= zd5w3or-XN>960?zNH`Vy#*CIaa_AtA;XdZN85K{@yJ*VI=$_ZSUwgO>j&D!%AT-TGYIk+LcN7h z-iw>&QY?i3@=frkuw+)e7U~3rI*KqKB24+|uvq;zJTSZROjyhubNbG>d|$|WFDnYV zU@pTd%Uu||w!dz`;AjJWdZqdvG5JwdoKf+%OH#R-vhS-3%Gj%CjAf&}SQeSb9{|!e zR;#FH-iHmmhwTn;njgvhI9wL3N~EN1A0qL7+AaKSVeI$-m@m5Kd zQiW7r;ya<9aO;&7B{=BApGD!Y6T_l;LMndTsRnuRkfN0M-YDAf`POuouSfD1Da-v^ z-Y~RhM(`EY=>DpZhct2GpCQV(n4QWR62Y*B+5~X!f{N z(8pn>>vsF^gl7koXC!bYuv3huAj9UAN5p1nz%0&0ce)Pj#M#0NcOkx8~wC zyW;kZy8KzE;@B03RrvjBW}P#=e19Ct4YxX%f6-nH8no*@B%h*PUiay!a{Cu8v~k?k z%l;qwbY8AN$V|{H-0AA)Qx-oCJU+)f>5uu#c-$J!UikD=)LvN9uC_0hux82XMXuDm z@K75+J!&C+#Sw|?{->gJqTlzfqC5`_0fRu(dx|mw>;$d>dX6hfKd>8U5HRt5MVSG# zpHP%M;69)m82SMIK6M-E2fga!(un(9R#c}!2BNNyO%m;Fiw4C?> zMezVVgV+$T`;4N*O+dhZV9$*M)A3-KD7=NKCrNGjQ2nY-UM}YP)HhK%@1MddfM-}B^8Uusum4Z8Y z6B{hB@CNoR`c{@wXh}+hs&=RqX()Pb7yfm#vciO!GobNmGTm!5?G%J6qXE0Z@KyZl zWRDakEI-w!C}}j`WwZ)!#u^R9lT?opP#16__XDZnv&V7pbK_LyS%`K-*D%zM~__&vz>I1|$2K7wcz3ypR{??fKB%_^go>dPjDG+yP{nbOQKj# z#G3VEvG6E!^H5jScpMw6FyJ;At*uyw3{;5wU>MJ@0xyGKOa63*%|-+HZvmf=!d3mT z0Rd{{fqw#g*QYF|_^y~fbmPvkg5vC$38>CQnQbaw7qb!Fu3>hr_>P!;=pLi)tXM0i z9>dUdo-KAwi={{xFswh%R=VcIIx);d0Utci_PYuuZbY{qH6I<2CjBmo8J1TIdAWg= z%ePJdvNL7bd^a=SXW>JE^!t_w!m1Lsu6#GRZ_DpQ#iUdu3-(-AP|+v^=d%43i}^M- zSTSEHoy(?GHd2J2S1$FIreRIA2S{|ZZK#>as3=p0HdYS<-6wQ@U#7;ZdY~K7dArW5 z`hhh0$yrHx)c>tDDXQO}%Qs zbk3mQR)~3%|O0)P{Hdt+OyoAUmZW?EA){Q#cug~Q)iMfEzuVqq=#j!aO`3{Yib?Z)@ z_l(=g)9k?-OVV?&(psiw-Cc)6>^C(Q{x7tvB)SKb%R!rFY|u+;K($`bt#&=5E>_mk6Kv{U zXKLFh(svWodCsm8^!=mO$I9v~Nh1-fe$r){k^7D!Cr!e2HA5q7vql?qYce}jm&Gq& z19eAC54AktS>M*tS--w*ec(`i7e8w}%SDLQZ5x{w*qd6{vb@GiVoOIy%jTwQ zoTKU8X`YS^t!)oI)YR6|v7v*z=Fz$~U~68+2Oez+@&Q+CC`x#87Ry+j9h;Y{Z438* z4Cu>Pz delta 3520 zcmZ`*dsI}%8J~N3xI(?N2)d#YHVNvYJQhX81{*Fcc!|pi%0pfe5Wz=L5g$oB2Jp}% z?1`IVr#UB${y~zq=Oj6&#crCWo-_e54QbkBW7J^kWAV8p7>T5nH0ty>_ug%K@|bgW zelx%Cn{U4PzHer4*C}c6q}1(@EC-U$M_s;$=1=hhcEz!g5~6vhq$m2)W1op_ zSJ$ZJYMEMErg=`5Xr8aNhC79x+gktCFs*ZRMiV_1o)padjPpXzU-KWo)@E0|Wf;HX z_OAC^G|%KGrVN5>+m;}|WyIZkWN(eeT9 zsqiGNKWN7WP^uiJj=$5=ztOaj-NMWf*w!&FI9QFUO(}v6I78J?LN7&&7Yvr>+z`Z{% zfuOBVbOU<$#mV#|Bw{NO^>JD?h@=d_0c;ku7#?8AT2kWLtELk;T5faNnsJG%i@;6nLmGti}R-FG|&0q-YaZ%YT@+L(-7&= zE1bV2ZdFje%6agyE9`LUoT$&^gp1Z`SEjK*>OQ$7jy0q)`OR21k+wIz0>a%>AzoaL zQ4>ucr%&#@!@2tUDD2?cU5h0ooj`L6Q zPz*cgeEq>CU=4$r7{TBT`s52RsCjFWS`anM#-;F(CgzYM!c+Y)I7ddG+dHKt0=uwi$)sNfpiToUF*AAkug@a@owud7#i}>|j@vV?c~i3~ zTbpT@l>qC=oK-4Q5pJkeYu@D(3L&^l;vS7?;Qqf()HQDJ@`LjUQOgTGr^L1Tms+ZN zZk560e$XZErZziuoC=2X>s{?L3G<0E8uskm9qP&jxk_oj~h3JWqm0AnJe}!*B@L zbAj_)K>tN-6poEAbKV8)Ma`Z7wxD9;A|MB*0G+rbE?^IEJ zQ6Ra_c}*k?0K0$z;H$vmG0um?95`bV=3l@eb~<;S(gKym?X%Q}b-`S%{(& zQ3*9G2K({rW!v*&6S5K04MOywVtqoE*C_ad6B=*e*TbIAi!Hk4=X{P>uPJu<-6+Mf zB3joFcgTu}{}d?N4ADWa8C;g3@08C&zc5;t%?8)aB-bqIRhHno7oMuFbg7PQcO^>Q z>;QT{JM4;0PWgoMOkoR6cp@Hy3UH6=jIoaeNjk=^!B&t>&VNuEWR>|Vr9Ws5pa8{O#ezqM3}JmK z%3+OIzfG|>z_Q@S${;xOEx~LEEU1r5;5Ui|`HHm@Ik3wQC|1=l=zu&Eb=4~5n&O1W zl7(Ro_(K=jo)t?X_hM+j#9mvG7TJqo$`$tMip`OB)X$+HtJ3C1x0v>t~;&}kWAr?iE500?p6cLr472%B1%J2?rp7leffpBgom~=hT-o`-if~t z(e&B_{>9|$7U(p|!CW*yQ_F1`0nB3LSQ7ehGi3NMGFNqrq@n4R4~%37*4QMQ-Y;SO zYi#j3NT`+I%cf4~ebb7BPQqqLzt4Ufcsn~R$4Xwds=^lkoY2F$NRD9v`P17=Y+r>f z<26Xdnf-QyO)<~#w#j!m4L%f$lVb0TpeZdjBc<0O6vJR#SgTuTtZFoX#Il%5Tl_pz z+Jb+5(L!;2l&O`rj4V^?bQn_oeqdN)@;zp@LVN++j#c%hbj(~8@?@9ETk-iTTIgxYTg;Qb!S=!CkSQgm%~ht&tNM;NA&nt3iM`ICEv7V@ z4Ogd1S&Y^kwYJo^);BjdHrMW`-_p|D*y5|y~wO_GkYTtKkXxv%5%vsy8ZEMq} z^^IE_8f!Og-Pzc>yS|yUv~Fu?+_9s!v9)zutM9S8mC|=@y;I*Wu|7H6*VPaRlNa0B z<_&34LuTdSiodX-01v{*hV+=w6%1iEOWv3s)0b;_9KtN$rj1;Zy4a@8OD7Lj8^T2P U(&nX6SM`UCkeL#Dx#6z=0#3LhLI3~& diff --git a/examples/more/11_billion_lognormals_paralell/example.c b/examples/more/11_billion_lognormals_paralell/example.c index 62fc0d0..b9fa6c6 100644 --- a/examples/more/11_billion_lognormals_paralell/example.c +++ b/examples/more/11_billion_lognormals_paralell/example.c @@ -9,7 +9,7 @@ int main() // set randomness seed // uint64_t* seed = malloc(sizeof(uint64_t)); // *seed = 1000; // xorshift can't start with 0 - // ^ not necessary, because parallel_sampler takes care of the seed. + // ^ not necessary, because sampler_parallel takes care of the seed. int n_samples = 1000 * 1000 * 1000; int n_threads = 16; @@ -18,12 +18,12 @@ int main() } double* results = malloc(n_samples * sizeof(double)); - parallel_sampler(sampler, results, n_threads, n_samples); + sampler_parallel(sampler, results, n_threads, n_samples); double avg = array_sum(results, n_samples)/n_samples; printf("Average of 1B lognormal(0,10): %f", avg); free(results); // free(seed); - // ^ not necessary, because parallel_sampler takes care of the seed. + // ^ not necessary, because sampler_parallel takes care of the seed. } diff --git a/examples/more/12_time_to_botec_parallel/example b/examples/more/12_time_to_botec_parallel/example index 750df7e800edd53ef5521603e9afb701289c7d52..b61f3d2fdc2013ea894b9e000bda7fefe8ca4938 100755 GIT binary patch delta 4508 zcmai14OCQB9)Iug(TMQgAn1UmPB`j7Fbp3eo+5cNz@X!b2?)r?5J6N}EPk|gJPaMn zoYpZt*K^jk-IndN+iCY~QP^4A(-w$!(6#)COpW9t(#h^_FnNz3FnHP-XENsFB(RrmZ>kqa-3f>F34B0Y=6CEJy|amovxQ`qZ=3bol;^a z=6()(*zb^xAxBIcS1%PSuq43I`n_M;Hx`P{fqGdg8uLpP{q<5ww_hq71bG;t4k481 z+l^@~7Q$aX2)++XCOxZ#>SdvN5atxZ6n){BO8*59%q|)7ONnC+?~v2`jk5bXrRajW z4y&SU#n`!dzX5}uN<0q872VS8XH;=UrCHbHA~j|2m135$N6#3`hCCF*amYO4H$b+H z6)LKkAH@caVY|Z{7X~sv4wtOHX)I~$DI|W9?G}DkKXyFA67R;O?~wNi?06sYo+~0C z_}=kgX@0k!`*HAU?!Q(e79UolA=2qU?w`P82$2SF%Khc#uIWUE90X*tK|~8 zSS~8?_Npgbdt+4r4!Y}n5Dq)WFGU7h@NiiSil|Up;QcW7*-I_)PH(sD%caYGn?5wu z+X520yyu;pyAAYH%*YulSk%Aa64Mi61Y4MRE}WfTmJ+^V3Q|oT&naZj3g|}})k9|m zq{_H;s#-g@m+(OgJJ*{9oTy)x7U~xLJZkaWBB0CrZ-;kSDZIu69A?FXd>l~VJ}@pX zjt!zHHRmWr9z_-dp5Ycq8?%Ph)IPaFu9U0fYE%|V21Nyh3%n=fVy0{D4JP$ab9kZ_@#TG&`Haa!gYBCQUEunuhcs9VY7-kjct~PWr2w zqxE*Grbp|ldTv{pk7yD9pRnhn)M&Q|`%_(UBXJp(Yh&Yu2_;jRunH)r_p&DJy;TF)x zVP@B}3petRU#g1enZmAlDZ7~=d(da6i}P${Zwhv%SkDSu#BL;xUAr6>T$Kt;$4xOR zFKMCCRSSgOBJ2$D@K&>SDVsENg~fodR#4cPwHt2)?M+0E+`TyZeJuL!6}-0JXAh^h zU%5C$;t~s{)xe-sd`amC?47POqo-UNB2OpI4GWEAn=wv(Yzk>ld|LR@NYavs@Czg9 zN!pY=VC1-sQ(#_j2EBSIc~p9C{S8Xxy`zYOm;MIODaY9r_AQ7ql6AH+p>2SCY?~R@ zC8}{F;*I30?e9W^NZv@^yW)MYTEO?vLqNau+$#v^^i^1zppT*9KJ1Xy@ITm2)nt%I zN~g;gk%<{DU&tJQnx9g3hP-2c=iUT(UYJQKm1aBGrUH?Y86WafemUbWVNs>$3-afb zWWh8@K1<0GO105z10+8+OHlgAmed5{-U0Gh=qXC=KszLL16zTu zzyV+zu<#6}dw@HDJwRhWrG1Paz=nXWLs&Q%3!{j40gvKh^a8Wbd9cF}1U?T3hWPX} z0e4)W^gZCnFr|i21o|3&z;<9J(0q~7LSQSf0hoP>(pTB|GNmVhZa?DPfrWtUf!U+j zD6k#4pN+3j${8>YLC0wz*OtuV+sIuxCccL(%9)=s@DZiBCgTk+QEokcC&AmHaG)`t zZTOkV&YY;MJC9Ksp_`ax4}LM=@W+$d7}tZJWdhd+&Q1n%qWEp(YEIfz1MCl|GcX^R zi4`9tS&kV~e}$8N^d?IKEWKo{V}>yH33=WzoxhX($uWoDO%6Jek;;D@vxTu!WP0wb z_{`Ik*6B0ZK3E;>0+)Uh$7-pQ)Z|9ZOhnP&!}1K8$&2PP5mIq3{IKK5Z zpZVpD`z?dtM*%-4Ou?v5fw$Mg&+Y~a(`ABFVy2h``cIUiMb|s2fH_m3p)N~AY<7U*L&q3oY z!cxSQ!$hMj8sy0w3ttMKvtN_2g48evhPyD#yhxT7ESc1Yt{e9#q(hz@6hsWmONEAv zjt7>v1Ov&?N*mu!BJVfx=8pLL{~!qAOUa?4?cmCa=kw2yr;BF`?hNuuag{KVPK1*C z_}7VCvQTJECr_1Bu?VM2@;t5c*w#31EjtOc4^c-`I3t&Rw>DN!2;Gx9e>6$sRZ-FX zPUqb^uZjwO8p?>paho)Mc2iih==_LH9fon)2*BDx(zV)@k`HsF5XA-lebQ_kSFTU% zWG`&Axq7&@Y!=_5K7UWXU2Tee4W3307u&5mjz9LO@sH_o2CSMvt#Wh+!T+yyuacFe zCJ(FiNFkOpxW~h(ok+$bO)gvGF;trw9hN+2*2aOWoD|^WG(J(!{a&3n>W&SXnRBT; zN6wc8Mz-r{Cb4Isn9DiITe7T9g#OjSIu{ z8Hm@+mHKJUAv?=WvFWR{5XS&z7IzUn1X@n?n8q ztC_qdMDhF8HBP9{Z^+6urj#+*)VAh7X!6z?Q>;B;ovpDe zHQ*MLL0DJ78VQZeZjCnR)>;x-k(8Ta=hoG%ZECEmPO+t>r=$f&Yik~Rk`apnM19>8 zY@DHwo2#4a)@^894=y!>Y_8a6)WS4ZH`O%PG&a^Xb`)1`<>Mx@8rXX)igKPJX!rC_i6Q!=}pa#*Ve6DS6$cK{CG3>qetz-9f#e^_zvp>ojmCjf;*Jj zfgs_)9CGL(TUbMemLA?XLPH+jPX!Ihk%0>yz((Q)YtGSJ131!guz~VIOgK68NXDen S3eEZine^yF{s@sDo&R3}DUGB6 delta 3694 zcmZ`*3s6+o89w*&a3RCJ%gY4`VNEt>F}N%WiitLRVIg38CdhNpQUD>1YsWzgGR2xka9Pxc|`m`@Dv`laav-=cdejWM* zo%hk-h##bUdDs2smyb2q{CV@!mdrb6z0fl9yG5Z#qqdskVy%Rz!GvFAg;?^J$z94S zrAnz(Dk{~$v6X7z8@2JOC-AL0crHrqxscW*ACI<0qR%tt3H-BU`^8SH60C&$rZ>2L z*rWz7d3xS7u_F_ToOth4+2V7XFEq)Q-fEIJMB651b81F9A;5W292I%)d$RrXq+C5H z+fT?9CuPUKwshr8$~gmA^fBy~hSLafveyK=?bvSWM;cXNR6E!$7=9e<^HItot*;pU@Qx#E*adG$N8eN3)G z2su;Uz<=C<^O3rmD~v$^b!9?MblXf36^u9V6=q1iCBoQ*Y>u|UHy6A%liG72Z9=Y@ zYLYvrfN#~Rf&GvVu2FpI;AyL}Mp>)I4J&m`s=Q%<9?DMXLV;1!*GbY=qEvVSUsSF{ zZOxHB_*MgNsDUAGu%9dCEaDBEQ+tL?3Ok1Vi>Ew`M>bzj1OI;2RV0(w%ZCVw{Pr8h z^y%eJi5M|&2Z?UzG4DS2|CkRWW``cM8Nm)C(BT^gt%v+{PMMexp<_8~#lM^b+;-PW zC^`rDYe*|U&rD0<(r1kEs!z%_2ho%v1c1#VCqV1PqJ8sHx|XF91iE_vgj^l75Hc!g zt{sxIwrMAJQoiMDII}~WV4cY=(FRo5H>i0tM>6z$yf$QD_Hh;O6CA$NB-;lgSM(Up zLRq182cBxeNHmPbKs;>bN%PjAE#EbT7C08C+=5_U`2rW`AN-{TMkC%&>CYUVnS&{) z^lP6o_By{+k;9N8WuMY9$L;Y=$vi}J>hTn6&V5WQPNp=MiZ3NmyYtap70NA1P#2U! zYLas&GN$ai06(~P$C3!iJIuQrEZCGpk2tGEa)vtcGUI0_YJL-TCDF>fzljIU zRGPnU`5dsiLmMPGoRFV*9uC!@&(?;jnVJ_NLOIJtUgeQ$$smus0dFwj88ac?U`!7{ z%P;aun5HgxdC6_C9M5Eo$ysh5)U+o%=VQLlmva6hHjUWBG&Vn1eCP~al)qHmH>%BM zoS}cjOdX}4<ZZ?FUY7kpTNrvvsyK{G>i%E62C{|Y2dCsC(0^saA`|6Axf1eaGYPOFO>=wO~z4%vGqmX$YA>H?A=M699*flM4qz^0BE+Hf8ZvGLS z3!?jd#(E(Se8AW~;4p9)NDkrKAQIv_%-F-2JC7hupbK~fSPnEp-v+b+dx3?(eqcAy zdX%wW0?Usv)(;#29srULSr-dKaGii7VE7}(tWlVM!dN-57uW{84D1EgqcDE}_5t?+ zt*H0AoDU-)VDD*c9F2|7Ft!aCK8uZj<+v>77;KC?7lLaXj)48Z z3qU8n{n4@580Z2zFEZ8vG+$=yEU+7Bo&h~@39$Zi#_IVT*vID+2xtbKRXj%+ph~et zh|ucdJMst$jC1+jzeQ{SVG@$-zMy%L50CyG`DaJ;nQ8E z$<%dkOqU`#ACGlQg*`J?Naay81*x!1a;lQONV1knX+@H`MDxA`iwRDE4<{qUFyP>yBu*WPdyczIx!xCALY5e*0190H%its84Km8lToZtm|ay;WTp zdlusEvvjTMh)u(J-1`}QTK!4vlMrp=)Zw``-VS?ya+~>dPHY(Bi{rG$Gb^5_%+K!B zIPLVzmoOu4{(w@?(^ox3aXTQk;bg~zl$$n1QB&D!(Uc!rSJoK~q^DOngg%;hzeOOS z?E9Y+#b-VA?TVe?sw?jhR?x>PvqX|lpRZgi_Rgnb)sKY7sZzB_w9cnntJZQ4C#y=U z%kShH6Qbn^Jc<@45i<$%4-JM^8{nkD*IP7NJCs}!&~a|xo1^nC*zu7HKuP z(lwsYhpH3vatxbpqo&$9g!#P&-@8Cx>(_y-qCc&+q_-H_KK>2B;GY3L{!V{iZOMDo z(AML3g%9o00Do=p!$x1H`JdB6SoMOTwO$v=D+cefa$kfzq60ZV^J*;V!-lpG<>F(B z5y>$v&P#^&@*TQ%ouR#|X=NOEt)51XPwGY_gYK%a2=lbJLO-vuc*_GK*O2b=qWTms7^6fXbHMBRhw)$J^bbsxk z?Q5PCsI%_BbaDOr_8J~fe8aHledSJbOL!-a$feju0_AL#qq^*5pqkq$;I)8in z#`e(chWmu?`)XM?D1@Br!%?DtHce`F#G8!EfKSJ==2CoUb~WcFPG9Y5^wC3LNwJZ! lfqkLG2bds6ZlXyo3*+oQUE4@kv@D8yYMnltMGyJ2{s#m-NNxZC diff --git a/examples/more/12_time_to_botec_parallel/example.c b/examples/more/12_time_to_botec_parallel/example.c index dc3e758..94958ed 100644 --- a/examples/more/12_time_to_botec_parallel/example.c +++ b/examples/more/12_time_to_botec_parallel/example.c @@ -23,7 +23,7 @@ int main() int n_samples = 1000 * 1000, n_threads = 16; double* results = malloc(n_samples * sizeof(double)); - parallel_sampler(sampler_result, results, n_threads, n_samples); + sampler_parallel(sampler_result, results, n_threads, n_samples); printf("Avg: %f\n", array_sum(results, n_samples)/n_samples); free(results); } diff --git a/examples/more/13_parallelize_min/example b/examples/more/13_parallelize_min/example index 92e31166613a83f99293b2f867c99ad93e83f6a5..37cd51e9605ade354ae8dc9b6c39b01a1d426287 100755 GIT binary patch delta 76 zcmeych4J|o#tkVx0>z2B1v#lj@db%Ri8(o`Ih#v-yqE=HqUot6@sr4eY)2>{taA*28R diff --git a/examples/more/13_parallelize_min/example.c b/examples/more/13_parallelize_min/example.c index c7761ad..13adc11 100644 --- a/examples/more/13_parallelize_min/example.c +++ b/examples/more/13_parallelize_min/example.c @@ -29,7 +29,7 @@ int main() int n_samples = 1000000, n_threads = 16; double* results = malloc(n_samples * sizeof(double)); - parallel_sampler(sampler_result, results, n_threads, n_samples); + sampler_parallel(sampler_result, results, n_threads, n_samples); printf("Mean of the distribution of (taking the min of 1000 samples of a normal(5,2)): %f\n", array_mean(results, n_samples)); free(results); @@ -48,7 +48,7 @@ int main() return sample_min_of_n(seed, quotient); } double* results_quotient = malloc(quotient * sizeof(double)); - parallel_sampler(sample_min_of_quotient, results_quotient, n_threads, quotient); + sampler_parallel(sample_min_of_quotient, results_quotient, n_threads, quotient); double min = results_quotient[0]; for(int i=1; i