From 03ca3e3b0cd46834e859c7367901b511ce80ce81 Mon Sep 17 00:00:00 2001
From: NunoSempere <nuno.sempere@protonmail.com>
Date: Wed, 29 Nov 2023 21:45:42 +0000
Subject: [PATCH] prepare to incorporate quickselect into squiggle_more

---
 scratchpad/quickselect/quickselect   | Bin 16984 -> 16984 bytes
 scratchpad/quickselect/quickselect.c |  42 +++++++--------------------
 2 files changed, 11 insertions(+), 31 deletions(-)

diff --git a/scratchpad/quickselect/quickselect b/scratchpad/quickselect/quickselect
index ae25e1c8cacc3a5009ac00b5c55a7d365c86ea13..6e47f82e41a4993cf32f38738fc4f0a74620aeac 100755
GIT binary patch
delta 2757
zcmZ`*eQXp(6rZ_|(krcZ+xFO2k=rXNEg!w>75V}6dUwT@B|<1l)kM)lOKDZ0NReWI
z7K+GuGyz_s^)HVEv?LI(Y1G0Wcu0g|2*Mxg`IUfycxchK)C3T^zS*5^pox>r{^s}I
zn>X{`d~I?!P41@q^JwBg{uj1^G@_LO5NOVU!KoQXth(Tc6|$)(6gKJ9K;qCjMPOZH
zB|4A!AI28muMKumQHXP^IJb(!QTnaBMbi6B#0FhV8sJ0Ia`*p@V9Y2CHuKB5HSfOk
z`}eVv6<2dhp1)Z!XYDs@Zq~v!bFtHg&Wy<tvQu+ygA?W=D7NI$wQ$wEn9c;jBGQ2N
z?L>)453;Yd6s-%O*OUu~Ebg?;t*<t%eiIUwJ+HIiOvH>;57UH0uu6>dklg0xEP@s*
zxNP&I%{UVic>y^qM!wQ@N@9ZE=;0YM<YDNw&ZZ{#)ha?}T4q)%pxd$BPlBv85vtQN
z3vWLlWdR$dpUWqqfPx>b-#{M&4@fPr4yA5iPG2m|XFKR9NZ;ipWJE3OgjtRSE1S&N
zbL2*^s;c{M*{9t_#dg%&T~83gdC3G`IDtPtfjgnwG40vUErcjtSpg0FML3M=9Ptn`
zx?kz?=VKS8%bLw#{_R$$6184Lo$^ki)uvvtyE<5WXE>t$4W{fW|IGuL!ejIZ6bsX4
z9zhw|JiA`$^zUL#C33wrQ;GWbAW7XAQ`J-h26fiD9X1FBa8Y<~(wA&)QVXFmdv@&y
zxOn|+JyYdq+9}_oQTg6N`P#N~I9xj_Q}^#<O3t5QNvKEPO6OozWcx3X%}FIve<`x;
z0$j_kb=6;rE=xw`(a6Qfxzs1<DJuUGkq5Q!(4RGMF-HuhkJugZ?MA!z^jo45Wj%(X
z@}Sn^Rc-7++9?m=_(Kcj<hGN#3A#vXCc!Z{pZl`u!WhhQHnxvV#?`V5avXhH?{L!3
z$S5ZXM*ikxfRR5qxx~nIPLhmV;p7S<KXZ~gC&$?4JSQJ!dtd32<EI!`Fv8_qFe%ac
zn9{W@9+eL(+6<3k5RHm`M?#6nJvjRw)F5B~Rpcz@at^xH6Mg1iRei|YtoyKRidY>%
z`FUfpS-L~dE$#AeDkdcK?2^d$q536}e?lumk$d$@S4jhkpIJKdJ}yUk>eIV@A5B?p
z!<}764@&7X3LEE0WEEN2v_@)fU9(PVUN5z5*d$FaY>=8aEszQuGcXPUr=%xeVpkB(
zjKu#|)ovyoQq=^412KtMhuDEQe^^zA5c?4wxL;<BsA?(V7-AjbN4He91u;H~y9H7D
z2m4ca%{_w-zyY^tb>IMXaNIq;y$W0Waz%%P<YH~bT8i=@qEV(|V;|NKPNAI%<wqQQ
zOpCJ8OmAgQ!j)k9V_0LTU;mJP&A9#o%Kze_@cvJugKx(T2GQVSjOoxr4!XznnHWnE
z#+ZRPJ^?eO%mA+jG6&^(aJpbodMmbqm<8>Qk*4p>fLL0u&A9IHdMoh?-A*<PA^AFx
z7z~%VXb%V;FYSkr$3>69YL8cl@r#SUEhG$wJT9S$KVj^WW2FZ!dzRBSaC=>}0Ycuj
zG!CJ%97vcfknrZwaNlLrpyT0E5e?ivr`bl(+$Z9Yj7=s@u)?>Pw!ltgcG3a8s7o+r
z5lFWae)D;J#-pW=sejG5Xi4coi6M7lOrbInf@M?8B8kJsGLgn0QnrN-0i7%6OZ@sb
z`ns`U!p9^GzPVnU)@rmda%R|U#RetJ^MSdGO)nO}Pjfv`HP2!DT2ruL-1Eec@zNOG
zu4;`OGBJQrtydl~l5E)L!`68&VLx9Y{&zqgg`MT1)r*TO;RBgNe}vQXMEVqy`c|Pd
z8z=vMw(f%+7*^JD`T884AJTzr0(ZG+dz)*ILa1CE-xAz}>}8(mVeQ+meSd&noWV8h
zrO)6NTrC$1Fl$7>(TSw@#P>xjc9gR@NlqxMn34Z0A1uL79KJMKalyt4amj3MY-eW>
zi?+)VOS!z52c{h_lvlGsP5HzTMDjpU;Btj%eFe4JUqm0R%%`>%{XZfYgj<$ksHpm^
z?@m>c_9YjDDQ<@O)n1{WPujp^a2dDnzUrw$15dw!xp1P|Yir}uHW;q1wDoY-13rJ{
z_~9@PNf@^J+i4w8OAZ|OPX@6@fTkJ;RMnK)I(RZVVAx*--8GJEBRvLwXMvM7o($u>
z;!zl>dD7HEU_rnm7#qsK2G|<#3dVUba4(z)cx}d}#FHMVgjApuo)0>p7;9Uwbo`vh
W__BoiP6Pv#?tnX?O8ok2i~a-oI<7eY

delta 3307
zcmZ`+eM}t36`#4|;0}Yi1<rS8JFtgCz%|_QWn-{`JHjSw$(9=^DQ%Phl@FVi;3Ga_
z)7mAL-HOd|4fc!FBuZ7>S`n671&c}r|51^+Hn_DC_$NKLQG}yZTFWQ_Vz=OrZr|+A
z5*xK6?f&My_nY_LyqVd#OPr>O({#u~htGKb=om^SdK)@(BanKR3wM(tllns*t6?*d
zI1QO1Fy=_HA&=BNj)yFU9D$S6EsS&4IA@K+1^Tg%_zcck?BeP!c-yky_x~ca@CYs7
zu(oDT-h21P;EAm{*|$f(_(8?}zy9NI@*m*1wa{fpXJ)d5oYY<0VZvHKTcIb(4bG$z
zn6r9lCET>S=}(xBzOKIu;apM`e4cdrM}%026(1uadsFN)ue=&0yR$Q$AZHKP^=|BR
zuu)BVJS-ERILn5$*8~Gc*XV;VYuiRCJhZvJHiM$S)jC69B~OOpWH%Hgr{ym{qBT7;
z-XzzjVBZ39WDh2%Ic7}lnTxq&ldd+$S!bU8^#CDDT7C?^<*cs2P8{Wo+y+9v&d9lk
zvb097StD<Uac9o+PhnniJpFx4S&r9SNFqdYWg_#=5Q<gUk>j>9235;ZmmITw>Bs)h
zCZi5*E+h9fmhxqe$am=c?Un!(37hB@kcFJBS8xF7-PSChs`(A;%E~{ZX>zRQJrvZH
z|7e<eK-09VwqL;;Lf&>+Iibo*@0=Wy)TIn3S;`QkWfbO^U?>HKOqQ0>N`_v%3>Sp6
z>sj`=kmD>RwE*ffw>9-KnXR?D_)ybMNe^gDdaz6S;`kJj+tH{_pnmAuWmd6fOiVsC
zA5>0!s`MsgrFll#cRkj89fQp%pUBFT+Kzd``<eCq(sGBK(LISyZl;+L8qt#+Yo3Yi
zOT?r*%16lW#mU8_PZeoie=poxhq=V01?4J^oM0pMTs>@Lyk-m)>ih#XWL{_dewFbr
zru`<9<SF$XbPdj|JbDP<?tb(-FTW+nv#&6ZEMu9@$T54Yc>ygQ8HLj<;SY4TyO>*c
z5Y!%YYoyFrddhF=b(*J`Pd&@W|9IpnbqNV(f}gG9QwXaC$U7x1V0{;MNr{ZXTbLjE
za4;XY@0($;4hcPs-<a|oI#92p`H{&@sY8#1IBteeuiS@&S*?~|JcQ}2j{YLlv#l|T
z+rU=9g>fe2abc1PPjF$H2_7!YFd>Tz2_~%P0^3wJ2N%>SX_%eNDd{|`R;+ngj_(_d
zN#|wdqkdf7OE^WW>Q5(SMH-Rg(&(QUo9z>qQP&aVva)}um;I$wbEkyvX_~(Ba%|u7
z0l7RazYqsTZK?+i*DiOgp5f4U_r`0s6%Z17evk4=sCkcaH}sv5@}OB3;>DX$-?Vq@
zL)@}<wRL4|zpLGGR<!qv*j^y{9YrE}kwiLQJsfR)wY{r1TrF<N??@$0(Qt%y52O-N
ztPgh_eg!Sq6~%+mYLTgV6UC;Z?UByzgJQTV{BpQE$_62$Lu~IpuvW*Mg*u|)SE@IO
zJOmpp9_@@C61&=uM8(c%_+?lx=0J(q(BFkIaAy$lD@_|lY(tzzj3CY-4kHdADhryn
zjCd8%gJ(xt)U=(5mk?VJXTH|7V~9O>H0>?K$X%=~;v6FX;*)F~&l>%4BcaEfbW>Ui
z+!oWHaAFMnS`pHaqTd8hVz&aHV_#}o5L~`3jYHToF9}9uo}%~+puU7Dc-xFPh>wCE
z^*cdY@BF#t>Gb4iS_<0|%>Fxk-unuE=-c?rb?n-3LcNH@q*}UkE~C!l^L?~$V7tz=
zThlF!!0c%_!EZ0YUwj+$i?BPpq-iUtUZq+VbXRx}kiBRh1xMb~wl}et7<(@7*QpNk
z{tCi{;((xTU}FjK8bC(iR&g%HLH+<8hAMw9JrCdW2ZRxR1F@u7li_WDt|0SY9PGAY
zZ2+$M_tO^e26AZwR0R&waX1$Uj3ojditZLm+^F`Ix~vY;g6O6Ly0i}6Jy-fHZG)3&
z%ODZBh_+!SB?&~`0H2roONq~j@#x!%ixT;tdKb<lROZH+gvwl2HyMW>ltw@)J4$DO
zmb<-#-{Yn~=Ej1*3W88t9>8(!$T4eWa{91C3G;oX{8>x(7WjL)-)sI5FiGqiObH9=
zt8jz8!dZG#H)WBvQ@EouJ;W@o=^+#PDsqK6zBG6<>un=^P~rCenc}%&%N68m4m)76
z!rd>UH4`UsDbrBvn+Us&TtCA$F}8CCke5*JL~qU@S8xb({#9P2ey-1EdJ9{V$N5Ke
zK7oH%2=on>3H>d5dq3p<*Jk@Mn@i(oXca78Vh=v<?Y8dCNfaKP|I#8WjF0mpqQS=>
zD|)lzu4J>0Y=E*IPkE>L@L_&>^<$1i5A^JC?|G7wdtAm~{qSISE7wcBbo%*3{dUx|
zB-fTG1hYHbwinQv2Nx~bW3<vsA+U2c{Rn>fS|f)$yYYK?WK7$6gTkaGWh_zMMri`7
zb_axEp16tUVYVy_26t~38u((Hm<bcR1CAE1ZGpwzm5u=}4M1s4<?0bMkD?5HHT`g`
zRs>J2)6vFD&<3Tol~7pgw3}ZpI9KC_U)BnlW|d5olHpRVKh<1eJbg>G&p^s9C&+c8
zu)2kPM)kU&uPz{%XVAnEn5YXtVXy=|K_?u-XJ@d)@y%a|o?zwb7iEMmPs`Xuu#VEB
Nuo9}ocWHgWe*r-^QU(A3

diff --git a/scratchpad/quickselect/quickselect.c b/scratchpad/quickselect/quickselect.c
index 4b72f6f..7bdb960 100644
--- a/scratchpad/quickselect/quickselect.c
+++ b/scratchpad/quickselect/quickselect.c
@@ -2,20 +2,14 @@
 #include <stdio.h>
 #include <stdlib.h>
 
-void swap_pointers(double *x, double* y){
-  double tmp = *x;
-  *x = *y;
-  *y = tmp;
-}
-
-void swap_in_array(int i, int j, double xs[]){
+void swp(int i, int j, double xs[]){
   double tmp = xs[i];
   xs[i] = xs[j];
   xs[j] = tmp;
 }
 
 void array_print(double xs[], int n){
-  printf("    [");
+  printf("[");
   for(int i=0; i<n;i++){
     printf("%f, ", xs[i]);
   }
@@ -23,32 +17,20 @@ void array_print(double xs[], int n){
 }
 
 int partition(int low, int high, double xs[], int length){
-  // srand(2);
-  array_print(xs, length);
-  int pivot = low + floor((high-low)/2); // low + floor(rand() % (high - low + 1));
+  // To understand this function: 
+  // - see the note after gt variable definition
+  // - go to commit 578bfa27 and the scratchpad/ folder in it, which has printfs sprinkled throughout
+  int pivot = low + floor((high-low)/2); 
   double pivot_value = xs[pivot];
-  printf("pivot_value: #%f\n", pivot_value);
-  swap_in_array(pivot, high, xs);
+  swp(pivot, high, xs);
   int gt = low; /* This pointer will iterate until finding an element which is greater than the pivot. Then it will move elements that are smaller before it--more specifically, it will move elements to its position and then increment. As a result all elements between gt and i will be greater than the pivot. */
   for(int i=low; i<high; i++){
-    printf("Step: #%d\n", i);
-    printf("  Before: \n");
-    printf("    gt: %d\n", gt);
-    array_print(xs, length);
     if(xs[i] < pivot_value){
-      printf("  Swaping element #%d: %f and element #%d: %f\n", gt, xs[gt], i, xs[i]);
-      swap_in_array(gt, i, xs);
+      swp(gt, i, xs);
       gt++;
     }
-    printf("  After: \n");
-    array_print(xs, length);
-    printf("    gt: %d\n", gt);
   }
-  printf("Step: Swap with last item\n");
-  swap_in_array(high, gt, xs);
-  printf("  After: \n");
-  array_print(xs, length);
-  printf("    gt: %d\n", gt);
+  swp(high, gt, xs);
   return gt;
 }
 
@@ -74,11 +56,9 @@ int main(){
   double xs[] = {2.1, 1.0, 6.0, 4.0, 7.0, -1.0, 2.0, 10.0};
   int length = 8;
   int k = 2;
-
-  // partition(0, length-1, xs, length);
-
-  double result = quickselect(k, xs, length);
   array_print(xs, 8);
+  double result = quickselect(k, xs, length);
   printf("The item in pos #%d is: %f\n", k, result);
+  array_print(xs, 8);
   return 0;
 }