From b8e501769f2a72e0ccfad385e69d202bd989d316 Mon Sep 17 00:00:00 2001
From: heyyuhao <yuhaoheee@gmail.com>
Date: Mon, 22 May 2023 15:10:35 +0800
Subject: [PATCH 1/2] added implementation for paper
 Approximate_Weighted_C_R_Coded_Matrix_Multiplication

---
 test/CMakeLists.txt                |   2 +-
 test/SystemTest/WeightedCRTest.cpp |  56 ++++++++++++++
 test/scripts/config_WeightedCR.csv |   6 ++
 test/torchscripts/WeightedCR.py    | 113 +++++++++++++++++++++++++++++
 test/torchscripts/weightedCR.pt    | Bin 0 -> 10592 bytes
 5 files changed, 176 insertions(+), 1 deletion(-)
 create mode 100644 test/SystemTest/WeightedCRTest.cpp
 create mode 100644 test/scripts/config_WeightedCR.csv
 create mode 100644 test/torchscripts/WeightedCR.py
 create mode 100644 test/torchscripts/weightedCR.pt

diff --git a/test/CMakeLists.txt b/test/CMakeLists.txt
index f1fa57db..0d709fd6 100644
--- a/test/CMakeLists.txt
+++ b/test/CMakeLists.txt
@@ -21,6 +21,6 @@ endmacro()
 add_catch_test(cpp_test SystemTest/SimpleTest.cpp IntelliStream)
 add_catch_test(sketch_test SystemTest/SketchTest.cpp IntelliStream)
 add_catch_test(crs_test SystemTest/CRSTest.cpp IntelliStream)
-
+add_catch_test(weighted_cr_test SystemTest/WeightedCRTest.cpp IntelliStream)
 
 
diff --git a/test/SystemTest/WeightedCRTest.cpp b/test/SystemTest/WeightedCRTest.cpp
new file mode 100644
index 00000000..80e03641
--- /dev/null
+++ b/test/SystemTest/WeightedCRTest.cpp
@@ -0,0 +1,56 @@
+#include <vector>
+
+#define CATCH_CONFIG_MAIN
+#include "catch.hpp"
+#include <AMMBench.h>
+using namespace std;
+using namespace INTELLI;
+using namespace torch;
+void runSingleThreadTest(std::string configName) {
+  ConfigMapPtr cfg = newConfigMap();
+  cfg->fromFile(configName);
+  AMMBench::MatrixLoaderTable mLoaderTable;
+  uint64_t sketchDimension;
+  sketchDimension = cfg->tryU64("sketchDimension", 50, true);
+  uint64_t coreBind = cfg->tryU64("coreBind", 0, true);
+  UtilityFunctions::bind2Core((int) coreBind);
+  torch::set_num_threads(1);
+  std::string ptFile = cfg->tryString("ptFile", "torchscripts/WeightedCR.pt", true);
+
+  //uint64_t customResultName = cfg->tryU64("customResultName", 0, true);
+  INTELLI_INFO("Place me at core" + to_string(coreBind));
+  INTELLI_INFO(
+      "with sketch" + to_string(sketchDimension));
+  torch::jit::script::Module module;
+  INTELLI_INFO("Try pt file " + ptFile);
+  module = torch::jit::load(ptFile);
+  std::string matrixLoaderTag = cfg->tryString("matrixLoaderTag", "random", true);
+  auto matLoaderPtr = mLoaderTable.findMatrixLoader(matrixLoaderTag);
+  assert(matLoaderPtr);
+  matLoaderPtr->setConfig(cfg);
+  auto A = matLoaderPtr->getA();
+  auto B = matLoaderPtr->getB();
+  /*torch::manual_seed(114514);
+//555
+auto A = torch::rand({(long) aRow, (long) aCol});
+auto B = torch::rand({(long) aCol, (long) bCol});*/
+  INTELLI_INFO("Generation done, conducting...");
+  ThreadPerf pef((int) coreBind);
+  pef.setPerfList();
+  pef.start();
+  auto C =module.forward({A, B, (long) sketchDimension}).toTensor();
+  pef.end();
+  std::string ruName = "default";
+
+  auto resultCsv = pef.resultToConfigMap();
+  resultCsv->toFile(ruName + ".csv");
+  INTELLI_INFO("Done. here is result");
+  std::cout << resultCsv->toString() << endl;
+}
+TEST_CASE("Test the COLUMN ROW SAMPLINGS", "[short]")
+{
+  int a = 0;
+  runSingleThreadTest("scripts/config_WeightedCR.csv");
+  // place your test here
+  REQUIRE(a == 0);
+}
\ No newline at end of file
diff --git a/test/scripts/config_WeightedCR.csv b/test/scripts/config_WeightedCR.csv
new file mode 100644
index 00000000..b6afe1ef
--- /dev/null
+++ b/test/scripts/config_WeightedCR.csv
@@ -0,0 +1,6 @@
+key,value,type
+aRow,100,U64
+aCol,1000,U64
+bCol,500,U64
+sketchDimension,25,U64
+ptFile,torchscripts/CRS.pt,String
\ No newline at end of file
diff --git a/test/torchscripts/WeightedCR.py b/test/torchscripts/WeightedCR.py
new file mode 100644
index 00000000..e0a869b6
--- /dev/null
+++ b/test/torchscripts/WeightedCR.py
@@ -0,0 +1,113 @@
+import torch
+import time
+
+@torch.jit.script
+def CR(A: torch.Tensor, B: torch.Tensor, c: int) -> torch.Tensor:
+    """CR algorithm https://www.stat.berkeley.edu/~mmahoney/pubs/matrix1_SICOMP.pdf
+
+    Args:
+        A (torch.Tensor): matrix A
+        B (torch.Tensor): matrix B
+        c (int): number of sampling
+
+    Returns:
+        torch.Tensor: CR approximation matrix
+    """
+    torch.manual_seed(0)
+
+    _, n = A.shape
+    
+    # probability distribution
+    probability_distribution = torch.zeros((n))
+    for i in range(n):
+        probability_distribution[i] = torch.norm(A.T[i], p='fro')*torch.norm(B[i], p='fro')
+    probability_distribution /= probability_distribution.sum()
+    
+    # S
+    S = torch.zeros((n, c))
+    sample_indices = torch.multinomial(probability_distribution, c, replacement=True)
+    
+    for trial, index in enumerate(sample_indices):
+        S[int(index.item())][trial]=1
+        
+    # D
+    D = torch.diag(1/torch.sqrt(c*probability_distribution[sample_indices]))
+
+    # ASD(SD)^TB
+    SS = torch.matmul(S, D)
+    C = torch.matmul(A, SS)
+    R = torch.matmul(SS.T, B)
+    CR = torch.matmul(C, R)
+    
+    return CR
+
+
+@torch.jit.script
+def weighted_CR(A: torch.Tensor, B: torch.Tensor, c: int):
+    """weighted CR algorithm https://arxiv.org/abs/2011.09709
+
+    Args:
+        A (torch.Tensor): matrix A
+        B (torch.Tensor): matrix B
+        c (int): number of sampling
+
+    Returns:
+        torch.Tensor: weighted CR approximation matrix
+    """
+    torch.manual_seed(0)
+
+    _, n = A.shape
+    
+    # probability distribution
+    probability_distribution = torch.zeros((n))
+    for i in range(n):
+        probability_distribution[i] = torch.norm(A.T[i], p='fro')*torch.norm(B[i], p='fro')
+    probability_distribution /= probability_distribution.sum()
+
+    # S
+    sample_indices = torch.multinomial(probability_distribution, c, replacement=True)
+    unique_indices, occurences = torch.unique(sample_indices, return_counts=True)
+
+    S = torch.zeros((n, len(unique_indices)))
+
+    for trial, index in enumerate(unique_indices):
+        S[int(index.item())][trial]=1
+
+    # D
+    D = torch.diag(torch.sqrt(occurences)/torch.sqrt(c*probability_distribution[unique_indices]))
+
+    # ASD(SD)^TB
+    SS = torch.matmul(S, D)
+    C = torch.matmul(A, SS)
+    R = torch.matmul(SS.T, B)
+    weighted_CR = torch.matmul(C, R)
+    
+    return weighted_CR
+
+
+def main():
+    A = torch.rand(10000, 1000)
+    B = torch.rand(1000, 5000)
+    c = 100
+    
+    t = time.time()
+    AB = torch.matmul(A, B)
+    print("AB time: ", time.time() - t)
+    print("AB fro: ", torch.norm(AB, p='fro'))
+    
+    t = time.time()
+    
+    CR_result = CR(A, B, c)
+    print("CR time: ", time.time() - t)
+    print("CR error: ", torch.norm(AB-CR_result, p='fro'))
+
+    t = time.time()
+    weighted_CR_result = weighted_CR(A, B, c)
+    print("weighted_CR time: ", time.time() - t)
+    print("weighted_CR error: ", torch.norm(AB-weighted_CR_result, p='fro'))
+    
+    # CR.save("CR.pt")
+    weighted_CR.save("weighted_CR.pt")
+
+if __name__ == '__main__':
+    main()
diff --git a/test/torchscripts/weightedCR.pt b/test/torchscripts/weightedCR.pt
new file mode 100644
index 0000000000000000000000000000000000000000..0c897ec3783db6774b22dcf703f5c077aae5c487
GIT binary patch
literal 10592
zcmbt)1#}%b(yp1=zGh}-j2UBQW@fs^m>FYs%p5Z_Gseu!6f?!l`7^Woc4uby-?#6)
zx<{wGucS|<(^XQd>Qj^fhX4VAfdTm!!2m%5aW}KFuyipqH4s&1GBt8BVsx~zl@Ngg
z(f)_{gP9l@xHveOSQ;3>DcTyDm{~g5nwmMuIheZIn!y1<XaanwHGhM5aWb;9x3aeY
zf`+>q1Hl*-WuRfWL5#kC_*M5?QT(L|#xG6&RTTmd-alz#;$UjV^p8;)9X%vOxc;Eh
zKXm(VZjw|-?AMr4IuGf>GpYSK_8FWI3WP9k=MkmsLxrbjsY!(5;lzvNy1Gnev{mrV
z2d^DBz3Lr9OVQNCeQD$S>XmSo9Vl<$EFZ;A*+jURLI9B5HGsqYNZY2n<H^IP?O`pF
z#(4!(IkzEdz7b-o*4pmi2}?*?1F{L63;}0>CzM2zfVfi*4+uRRyzz{@w5o`CM$#Ms
zqL`&SmROcB`L=g?r)m|dXiV!q25tuqh3KWKG8#}7nDdZbfRMx<wWS=+48k4Q*Mv)<
z>Y(xlb1FAwsfQspEgm{|U!pw7?jR1p#*tl*$+4Jm-vHbjgz85UohX*mTVdf7OXdcq
zQLgv!AJ`E%@k&8CIaLaHI*I}&`+Us@nemqr;A95DWw|VU4s<P#(7nQeTwzh~WJ7ff
zbytzsJpDD5H^%D$CKCj%q*flKz4}VN0d-?@59rOAvR&8c@kp>wyVR6$@w7@`;9BO>
zWb;mzjte+LFiuLq)Du~+E&1^e>>9fxQ;{gZzn(iVKfm2+cRdoW7nC;Vp8hOpE)zy3
zqO31f3MdkdIAU&ib?Tl<sbBHV9^ZW$v~HD5_)aS?ue1jI#u{=WvH4}8prI}mLs`Wa
zrxDLz_K}}304TRyrd{3nYQ<N?7yM&-0!j)RvT`uLBwy|uyXtd>cXM);uO4+5WSSS<
zP8mZ1XT0E`-0(I4{rt(g-OVli=0}gb5{4*8^~Be)^T7|R8z$h()$!A-+&lihgHo%n
z)cxu9eRwS!taUi|Rjp%I&E~Eu@^Tas)|yg&M$CJVr^U-_$29f^t1iz~-_R6}Wc!{(
zxk5u6y$Cl~^{zHz^p+5Q^tub)BN9mkvpAvW?)VUWb--zjoqYf02GB5cdbRcOU?3nn
zQ2)s-s6cT4vs*Blni;!V82s+DnE#7|ge2?P>PP{;*%?sDIn1)ZNv1bgu{~KjXMJ#+
zb;(%VFx?|RH%~*dBw!+-)^fH}IB&Nxc=7$jr3ES?V*KelVy`mhfHei22p{n_s6Hl%
z0C^nBnW$EtI3F=#Bm|naRU_Rw&dX35(2PRb_UGkg=FjoWG@i*6G8c9XxQ%z!(1QBb
z_9?Swr?r)F=Nxa&6uIiH*6+axaNPUB$5NkkZfxxHGCz%6qwDuJo^)M^9m#=+yAgGc
zRw%Dn46aKINfs;8v5(nS9}LZS9tMfPX*Cha9LGp#^n&(I&Ckzw!LAMy*DIH`CL2_J
zt-D<6RAselYhR>`NY{~0_h9lcGXk>II+}RwgkFeQVavJ{hn)mytulCQZ^{%0K~IDY
z#}%W><dq$_FEUD>PNGscb?WT5UEI##(A4=y9Z1J*lpLgewtk4WP(KBXF&<AUdD6p7
zFn8~)8cToBOP>9%q%5iB*Lkt%#(=$4w$gjmHmC4{G+M|D?m{2OkY8rSH7-4AQejCy
z;Btt~54h6BU$KNNMAu)bVB5ec;0K+_+H)a&<cl_(?nrB<P3{hhOoq-sz%zzwV-v#)
zq^6;y`sVLvgImqge1H>55JZ!iYZ2J+!hxk#b2-d$js_05z%3<6rIw%ZgWsn3+-_;T
zEMDGpM+<*J!~C+nab-DkU=i?@xA=%HY=XCSHi@D`VgV#&?z~SS-);k^C<n>Zxmjzn
z?z4J;$GMd{9$>4^Ms2JDd!u8Aa3Fr{1W?~oHh!;BXCKt7Nyd(z%4Iud)iAMFqE5}O
zYd!|vD1zXM+oz7-qMNr^<KwXPh#&_sw5D*ThOd@tR{}mX02MP!VFhk<ZBBvIu4&qC
z%!HHab&xJ_DjPCEke^a8vC##MHE2mmNuMLEX2ScB?#4^^xy@Of_NwS2Gyh4*zNuKd
zVO-gfiT~SObddEx!t7q3?viuCB=94>LuR&RY{up$dlQ!llle}1PVwGL;|FR?UD0Jt
z6W$n3k~<Mz(|k5kf-rWNv|Oxv-X3(ywFvqK2^f)EV-#}NnaQhY?opvL)0Yp8ZkqBl
zSv2X8>GVrw^sP(O8I9y$b}7`aI>|9|)M&qAk5^F61&4RLmgw`OVO6Wm`LwDHsyja<
zxPu=?%YN-Tl%2Z~E)UyPLqsoBoW`GdEjv1>*<llA2r5m3MTJf5TArM6qlBtdaA%w`
z;=n>u_hCu_Jk_7SR#rVh5h&g0Ft9aKAQ(^Bp+GHMu%+e)5g=)8w0g&__?$hWr|Gn@
zZM|lrq)jBRrFla)iL8ec3@fmA8x(p^c}7QgP-Ud-S^)#kXUpn^qTfb(B2YpNQ9eE~
z&S@8XGh5?J7cIKbwaQ2x+})v+$)ZdN8BzPt?uWa)L#QxAvS8vQo~mn;!F=eL_+36R
zsc7Wx-l;?4-7y94w&D!U+o2{QKDJRk+4==gse$y+E=84d>#Wrq_*uh$6S--kOz@nS
z4hji?^d2_4P}w&IC_2PbFN9$j!AenaSC1wLdcwi)Be~$4^G5If$pr7;&mG8tKgi#L
zMbztabxgeo6eRRl_;y7s5EFI<>y4HJ=7Wp8hx$nn+%NwVg8tO5^VNr#sXh2lU}WJr
z10*cV7}<xDxE=zs*G4h@nYUqMHU%h@;&G_*P{Dk;K|4-f=KH!W9F5f3Lu`+vUBMQ$
z0Vbc242FPwvFMc~32{dSBbuqQ{t-hFBpCw(Al(RwTVS7^5|UQ*l6)kF1;a+HgZBvK
zPefwGhPj9J<-0H$i+{jBipRt{ULz*?sOZf%t3HRuy`i!q)Ws%-U1RuG*9pQ#keWj6
zE6Klqops+I|48elVg6HbMbPY?Dd|neJXY=AfW(5#6t!S-By99R41<$#KMAH4sxBw-
zlNil#hcwp0dSQ6o(J8H*c!83JS`%fx4y~8Exu%(z`BeE?h-5j`E(U_*F5DBh?6VNz
z-u}x;-=^v@G8l~$G|y93Tr=tZXX4oYokgK^1|A5tX^b}V{UDepmSObFGwV^K*FfoR
za#Kk%aNKWHJFlg^$iMWbth|>?>MS%DP(6;_AI~W{D`xc~*OVq}_Rreg7i4|OG0d5U
z0d>R06;d?@SE6uGG|vw_x7t7L3FT>}ulMLpmGA|%4NhERT~l6Z&L558S3K4c>8k3L
z<`L*0S+l;XR(r-HA+p?WC;qr|0r$io{fe@WP35O@+-`lb1^&bqXR*Ur3hq-dyP2(i
z2;cZ}_n|}hh)15%llb<q5;Dp;&~p^}te&3%uXnOP?>K(+TCI0;X4tET(eb)%oDlXt
zK%(d?zJep`{2CkpO&A_<T`RE}HZCP!b$`+P9kL_5dYi%<AQ)1;W9Ch;)kEX0Vm*e}
zlfsi7xK^6Jt{%3Gc8!etL^o|EFTk{~{be2W1NbS*wf>hcLi)i6bo=E>*MPW`?YMjC
z=qC+YuXYXX??jL<@?ZEMILAi3&~CEivLCU3+r9BaOez(iK|mT&{?qQw2m<&IyEoJC
zpVUm|uJ$G_Ru1+?w!eSnqW<e9KA|BKyT*yqby-d9t%awE7~zZO=x|(W89;nw<}Pf>
z;3+CrgpmN3j`?(%MnV(e`N?i;H^|5;^ycv(V~2Z>iWJ^88%|5zTtNU>3J>5%dHmw_
zHTr|Bd|7G>R~q-_$W|e0i0n<VU<jvzNP0?2zd^{p{CD`nLNzaj$k(fn={8K%)85XD
zPq4`52VNuKi0rn5Y}ad-0#*E$b*cEn1J;SEooW`0h^V-A%wig5TGah)q~rlL2CW!L
zO|bG;mQH(0S_ID7W%tB`K_TPai4InbtOMf1T;3%uP&%1>f`Q$K9}`9qe1ifGzlCQ;
z!V&)bX=F#=BS8^@A>kkW^Ajf|9orXrXtWeYI0}DQW0jJ<L!MW+5RnM(r{*ZwiDisn
z)LQ>@hFn{KQSa_S!%~pQCw@9;JtT45OkiP9OH9tw5BRBNmId?G`Sp`0=H<~G{Rbbd
ziq7V$eAn~749)fjGvnH$M~GBX<;jW;sdT@rWGOLv9-lr<!&R~AFFB2`F}Zw-XY*Ig
z@x@p|eK&6EI;-M_3`0y&M~P+J??qRQHMDlL;AXM-Iud(n<dd*#{*##5sOQsSSS5;d
zHM*YbocavHp0##|wqYc3&16#RMiF)=sMGBsF?hOhG5KuU>99Gj&|Kf_qJ{uD`n0}$
z%08>r@0x9(Kdm{i#S61jq%#kRRPvs}{9@OO*vvpC_8ZY!QI)5D_8z?Fs^(F;-Wy$^
z*q7a<+FXr2PiPi@mNeKl%y*~7n0Jm2k}SNARscdlxYev?gdSRHIi+dE!p>UopK>?s
zi>z3}H4c3rt<9x*)@n*5uSq?tY9MyTKjXm9UN^27v1eDl;y^3uu2RIs{`y0YWO)zT
zg&8JvXu?#Bw58k-(-zvBa2mZTK+1n&0>UTk>j)t79(X!^EAe$x-3t1yEMV0jogf!l
z+$}5_b4~54iKv0$laonr*32D4JbR=3kcoDfp2WFq;n?SCQqS9i`&Q)T^Y&xSDXbqG
z%Pq}ISB367?d@c_Q05;9SM4DV<fePLN2u*Gq8W$2X_PPPM)!wSHlENyh)p0i)ZDOW
zA$mS3gFn}Tr4^fx67MBU_9ee+N|8-=7Tww>@2GY5%zSMr{W73|jAOey&R<V<7A+Bb
zzQ1}=;@jS|y}b0E8MoCEhz~X9{_R3>edE+C!(Hv*we%olYweu0tLNQxoV0IwoeVtK
z$?X8aZM7O;v8#Kap0nDKn9j<WQiA0@kGOb&gepuPARlA>SzOGTM!SF^VNQEw!~$uN
zIL-iT9;s!oT*pE+YKrh`uHt~h)N$d7oH$Juzct3je>G#)Wd3!qM*S#S4XzV0bD^te
z;qKE1_zD#Ts1vEm>;^w>-Pc=2rwy!oFJKh&RbcB%zrt@zktC356!+Kp-UazTEk)cQ
zLjR9T(H~D1eE+&A{Wk~izqhIU02_zD;Mx|O3?W(5j%HC0X62TawLtg*TeJ5>%-ct*
zLS4DPZv1}o_KZFKI9tG6qu&mRD=Zk<);6fs7Rebq2oxMBD>>-4B1H>K$|@^RE^VPL
zRYgn5Qq@f=Kx-epeB0W%@tZo6PR6maIJkTH`R4xfVJqha)6(W}BBS|m8-{>zf?$oy
zQ}p8MXt)d9%Q?TF_x{rDppac?%K5<N+jhPohT51L9nD<X;%kF8!3lkvn}s$*YEm9}
zlUc#SUZ|3|Wvu=^NP~@UIQ8RQS#OuSnVKvAH1M{f%|gc@XP1834a60oZ*O677{0Tl
zJEroMkf-jl#i39muC}ID9NmFTgil}A_Q6aZFZl#7*`gh1P^)ADS|qYwQilL-+;kpi
zqfyp|HbX6Qa8+){*&DRfo2EZJEnAzptSWxVCqMT@9aO!lJTdeH7vP0Un1)%QUIA(m
z>RBQX$p)&y=`!W-nqsM%+dc(t?bfBvLT~UOw`L4M6`dep2lhZ>HaEJ)SOMcH)^JRQ
zo?Mx1XEW*uuqbtl5K@_XA+~Nss+Y^P)|LQW)GLaHGl}^Wn$ECT;q4aqS)G%ilW&kb
zQ?frW+u<mLowja)xoRCBuabML5ghC|4XuR7L8Ip^S9%7CaIXu%Z_g_qi>in(KI<2h
zZD`JvPgvvQn3}+EB8z9A{uE##uk}I27j-d&ikOx$V)YKgB7~EE8j|1B%(-kos144%
z*0RMACbdAP(sEE9<XN|rh9t|#eH1)Qt1vn|c;<Xzp4yqqac?F8RREqpE<^vp70{_y
zInPHVXO~)2^%+UE$rVErI61=}z}!{FiH&I&WKDxdg&h&q_VH8S;&qDzY4gPCckjya
z^`?V7M$l63L+dcK{X3~C#wBmjti>@7@ceJny<xo}Gjahpq-0}*!fiY0mU^GzbGZ+t
z-b1S^E?(^20>vL7SONS2!#Re@CAc2-mbj3UO*Wqz>$;WZXQQ+j=_akuB?H;1kt}Mc
zf0Q+;2_2InhMpKhC`*{r_Op~}MPVBz5W`2-OSDr9^L~k$I0oy_sXp!tFN4ED+D*t~
zag4NT4C{0V_HITMsL&=9IHP^e(TCCylDBIHQE&2y882y%`ifL%z7PS+X$SW81k8@O
z7xr0khB4(F=S$>ARYE>MfkP@LU#(m_Ba3r*=hvyDq~O7n-3FTsG*y4Adm3|XZ&;El
zZi3aJ0d`7M?+f6$R6nbYK6icbN8lo7tb|W5mABBS%buF!Z5w4T(mZWqu4GTy>{z*)
z^7c?Ey57b{4`)a5Gi&Fg9FTDysY$V-#ekg#;8bSgmD+vMgRO&R3!Q;Az*+=w{QS9~
zgHS1nzm5g8cx<xRBq3msn7@X9pc!yzZ09fwkvx;z=N_a*x_*E|Pg%;HVo3mN%|fyL
zQ`ycBETcwap%K2ZqnArZp;`>I?1Zqz&b7qPIMB}LJ8R7j4vI{Q)C*Ww6uL+L$qfUk
zG_-jLWOr4voMZs!0p?N<j1pbbSiv20m0>~}6HQFn9R)asx(q8k`x#C{iggR)+9i`i
z-81M+|MY&7T2*sXI}nUhP3W1e$@Wg3fld)pH5KU<cOL-*ir7#PEm|~MH<$ag5f9Cl
z0I5Fr8-;yE7;V~E(6mWaoSR%#mSSpXzGI#yp#rv%F<5z31gz24HqFB*rAku_GV?n=
z{g4i{0st>P)r@B~*l=fX!IVGUAj_Y1L}>3BWo@DC{COQYpCKPLIkH~qG;Lz0)8!0Z
zmUUtB(B_DcMyjWVT}>siNiR8p>zCTOT4>cz(TE!=dFE_lA++)-O>Dj4HdleNW~ERM
zvozC?P#w0P^OiBN_af&i?d&=;>N@LLv{$7yjiPl%)K=Y$nl`9q(;q05$m4kjrMM`z
zR&cVTc?XYY?ju0^xxU+WRe0&EqF;mMMKlKPqz%jD;|l3mkbCAu?#P0Gn~+3~jEc|?
zd|JJo-eC6k%s3h`gC8J=a9eQ9Bd!SMz;Lu~sqn;_k3|dW$Y^lJ{RUst5@7LJfM)5!
z{4gOwFN3yL3|Ac_K4+@y)X~Vh8A?sJ`^%<b7$eRV6c7T`3PWzpGz>m9Tk;G2;qihs
z4tk8kw@E&5v&>xtN%3`hEyw#Hh<q^C-Yz_uLebV}*T-CM%btK8nkvAd74=|cN5DrS
z78U2A{M&r)434|oSK7(sW$Jc?Bh8IGr6%js;MTQf2aO51DK_eqUE8@PimMMToLJiK
zuG|rYRWXxQJ0|*UW&9_(ySkOGltl`CF8tsbmQvNYrSn4C1V@3_<mGy)*vs7R1uoc*
z(G9WAL4{4`O-S}D5X&})=>o9ep$+r+xBb$bji}bk(ZiE&=dT9~P;CCK_KYHla=}$g
zlE$9|77ZlbY-P`7%HMLsf81&!dx!iC+S{|1X{4R3aeF#L*!L(0RjU&{wZ=56U}9V5
zg0Z=?c!qB35y@c7*ycdmbI%6d5WsGh!Nr;(0Fk-1gB)Ue3wqOu+qAd)btEnm^+EjV
zMa~r!Za-46z<@*Eq_L_Jr%7@gx~6?Xci-^N{t7P5WvMYasUG0`<s>MN*)6B33sS?H
z+{VD5sN0WCWixt}!&#AF?TdQ-m0t%PgW%=oY&xWwHb}TqN0u%O`8~%La@J-89ZY%u
z<hb4GWhNsPhGYF0hA*EzL0Z!b+g&$P`ryv|O$oZZdf3YqInmfuBAsiF>0@4t7iuOJ
zlGOKnHIi&l-f=6dt`fQ0ZSlV$4sHr3hypEjG<^|();N#ZU}ti*HWgv?M<#Z|?s)28
zbd!?0GU>J%J(Y;fBu!*%uN8#--lu>ZRS%633k4yTu*R{0P$<=E!ut@H;eHyF$s0Y@
zEz7w890~Ou)6%+7qA2dK8>!eN)F0ek6pFCzq<CsghS#ghrQCDF@)Ow7#fMyE4q~jD
z$&)96oiS+@`6j!)fzXJw9X1^^_Uk;$&Ef?7EvJT<npr?O>fD{&pb+m!g_7J;EA1+?
zr{4X9=N#Q!T{`!pKD8Nqtp!uaCG2E;e%l8&xV=bc=ZfAt5}7jywd=RCl=Jhfz=?X!
zv;YZlhsH@=&UbeHYV%K>;NUvVcljQl4`BJ=Hqjo+7Q^gh^uGzABRYwL#u7pade!Eq
zg+OB+kX0!$oO%RZF7fJaFdCjixK8^jA~aZh+m&<j;m$^FjzSlp?l=epV*rOs2}fds
zCWA7-I8<LWPUhFV8V1trj1rDuz}tRZaY%G3v2=R`@j2Li%SnP&r06e*DB1&ZM*xey
zop(7t7R(dEhSh0WvlfNfiP>iL&OUMLn!}d_4t!)yMJSC5X9o<{<(7YSYfzg%ro>N<
zH8ybIA`Dc{@~*m&uO)O5Bd_51;A*E%iptu)S5fZg*(aqcz&Fe#KW~E-hXoFRAzP-e
z4o<L#R6vqDIDDYbQKGaWaEWx*)k!>4a9m5yKq|86md`vj>W$n_k>KHEWN)8tQCFF|
zt+1nPnwmClQpK4n+bL@*AFMIe-)o>)$!I*aJZF#$!F9`MGqY%pDbcI`F5nEsZgGQb
z!21e|hM`qi%gegbL&h-{CJM&|tW6{U=7Kt4GPBfVVBCVla=>}W0MKoh^+VU;da@SU
zXBOuYGB`dSM5#Vb5F}CUMy<~poU9f)7vwum6Ub%kuj74Ru}g5b{x0YZXeH*|NN~3*
zR4^VHW|D7S-_SlsPre#o>I&Ovy-hxgoYARHa%pe2?MLmz@+3AVTwQw>sy34I>g@k=
zpFAN&FC^UrDehGbU8tg<Hbt~loQt+&CZ(!$`dzm1V!b>%jBAL8QI-+;;>;_ansCfo
zG=KU?GHmxi7Ty{#70wu2=rPlOREFce@Ur_ZNp7(Yw#P5<<IQl{12@Ri7=Q9|Z9(bm
zFh-p1k;?MBit-Xgo!A6xR(gEUy|5N0$%s#&!w}2B`r8kZpys`v1gj+E0OGfx5kVAp
z7?IF?wMXVb55+b7m3*tp3S63i>NGQ2@TOE(jKjPl8wg<zRFfYC>WebOfWx`jXuDT1
zqq@4nKCmbp*?m4}ut0bSZN0~1kZ*p*jGP&XFq7BJDs)v$@3sLqsCreNuVCarjSBAA
zfs_~^9_mF(O)fH7#2Y<HS5N499%=J;5m%uSPl>KjK@`Jgj*2j(&xUW+FOfGTszy(<
zaKyFI;_$(IE{C-THJ%Wj)}R9x(Jzk0U`<4GQIt*Yi@N)$Zq1HwSs&dHO{$}9a(~dl
zgz~zOM;L$1aTK+=bp;$2aMIK;Mnvm<!|A+WVY?}!xqLgpPiHTZh4*^|mB^me<9_n0
zC82bj6w|d0P%LifsF^kr+LOKd%=+#mdbjP1NdV{G_kqwOVZ*%h4A(P{5Srm9Pnb>5
z6QVkGjRd>ybbr7*uyn~;+w1y;)Mu3weushJ5b@J{lxh@}UuK*O_ux~-eAf5O>opW?
zU-~xiT0A(pHkM5YPmiE2uevq+!Kci5=7FUGr2$?6LTzrem}NdZ(5ajh<p}&Qm7k2>
zaf;X(QP;@=a(0DtvMd0H5k-Carl|(w87j-YNO!&GO63wqk4_|Cu;?fa@8D-0s0(C)
zl$oJq@k;kH@bbFAu9C(^bF5@<Wi*O5kEjdjU0zXL#Pt2UILda!c!wtusg`!S$&`Y3
zg02Cb<0~wpcwAN`Xj(jUwM-_P@?WpF^?^>=v1ot2q-eK18#7j@V29e1*%Z&cd`1rv
z0(7&imS~J`;npUZu(3VkMb*7I{p8@b{rgMRkuT?b2y_x+W)|rdh@D`i-ddpTQB>>n
z9NH}&Y0XQ0brW}Xu)9p1Lw$y37y45LVm~A8k0x;AGDT({r9qFPK^$7(Q(@7vIREL3
zvF(o{DmQ~%9E7jpSekhtkvF#R>A=@aMon5zz#D@I>w<nh_Y{X-sjyK{!bB=zqSCtQ
zV<COEPv7J8j{}P`TlV|D!wg^Spr&mW_GVap%iaiEF!}GY3|EbqxAG<d(>v2C9@#hJ
zyqY&<@J%(<G8*L3HyHQBE=J8ZM!iFC#O*Vo1j#Em;+4&X{2}b<Z=+*&mX8!P;Wc4I
zbp#vK!Gnq39^eeG0Y~?ge5Sz~flVAC79%Nz<poe$gC62-4QMc|*35T!tA=Bx%h!p<
z%p+fSJIdr>RM-^5YdCz77Y`Yce~>snrq6d7FNLvMMBcQ`b9`Rh9Ke*&Qkv#!=(B1y
zI|Iu%iQ1_@c+8lc5ul8DjS6DxTA3Fc#9|xec6zQXx`5H}_&gRDm$q2^C2t^49Jo3A
zti!pGzuNw!$~K_D=^)cRlr|_jkdKj>)-M8At%-5CmJl{nh<0@=|H@`$G!}Ar$*Jcz
z+($d<wW?NGn5OoA$?htDm|y_hjP;s(1mpL<@;V9dlK*6S0?qG@eUr&eIw>Lwd>?l2
z${@YD&YV9l${A6cc>PhAC%ZWA1WNhXGD15^R6L{-DkyRg?-6jre;@F3i<TsCoC<ze
zqz>jCGm*lS#+#9BEJ8%Rr#Z<Kz!R1!Gf7rHg6ls;h1KJQXBF0sRF_;4G!?Qh;)L3o
zstdNoNRfgH5MT=Hmaog!6|N%38<t2Q)ho!*6`w5#GnHgA<cb%VRe%~MUO||~d|Z=m
ziZ7Q+))lP^=_aYm(uMY0N;IkqIwNl%wyz62qu`43*_RL*_DvNY0d+zVK<?d=GY!V4
z%-@IMk&Zz{%wCtZCo$PqXO*`n<6#wcVTDz=DG)+X4t7u}5R-M{<&hsP0Gp!!WSmF~
z42Pex-R<SdnK#^IPYd`#C7HDjn#voXowW|RbwPQfSotvhy4&*<F0dxyH;WmkD&(v5
z9FF8VE&Wlb4a{@U8|35C4B08zEWz~btFkTlGybk`Q~G=ViN7vj8W}aSJ5;}elK9x)
z7bB1k^bGb$0Y!JHS60Z6(&VQ^_6rkC)q$ad&`*u*7d)I79gm9;cnC};9pW8EEa1{e
z%c&?c_;!(h4Dlx(=oc_;UE*5D_9%YFsZJf;17Dm6AsBsRLq8#fNYO1JT{FZtfo&+T
zXQaRdjs~P{Lm@$lnxWna(6@l7w>-Xwk1ztYp<NK$$^P%fE#tExX!a8#RtG~NebeQp
z>m}jsT0`H()*GTSw5Va<inwnHa31tfAKb&+Z@}8Xbe(O4{9?BGv8Gk9AWU_HetPAZ
zCjLyU{mC@-Rw>!+6(5}c8A4ZUh`&Wrd*2A>C7k<)W8mB&Ux!4!3h8JKhR`-27hyt#
zk3M1i!GiVv<$9fnlzY(RI^F-BBe1!7h8MI0vrqKkc+2dJX}`5w><R-q46MzQ1-Cm|
z0DSx2-(Hj6Nab^3Yd_3;N3>f5^nRQFX85*X<IGQxvpY}%1Nz#$v1KU_d!yRi(nMQx
z-hM>U+cD7;M$}MyP0}4s(w4lecQSuour}D<Dz^_X+mrqmnxw#nk3$zphZmx8O>_AE
z?ZnJxXYi!xei1DitH)vfYXs9b$<m^{Dbh`NnzZBVw0rlYj}GSu>*&*>Fz;nq@5m^t
z@NnW+OZXk2?9?VV*Sg@n6nsjT@Roqy?*z9ZAPxwYHL9}%Cam`bd=Ef~Iec<ZANSII
zQo-a}Ka#084T#GVL)~-aqiWQS<{$wmP26dUF>=l?W@v7MtX3Da5IM)WWH#VVAZR>+
zQ|{ACMpwQ5MP|wwrnWhJ?jXHHuG8I+XC%u_re}L&5b!?8z$_ap!EO7%330+U^VWgg
zZ6_mhRN$Nqjy-uUVc=ckt#GObcP=tR%A=3ikO1uzf5bMt&GDK)LC6_{i62FkXGGwW
z6M-IL3)waxlD(ly@-<tTJjX*H6{eGV&zFApU3%zEHNN99Wq0~U1H}Fl0Y7X@zF{~(
zxwl4V=tA7_@`ERKj|pg}+>m?)%27I<M}#|5YiyECPT<4K2)NP;YZtTQdoa>O@yh|U
zSkQjEa_Y-It-youHOdn^tl>fb)dP$N#wHON)d<ZsHyw(gsDO;XM*)RwVgq6Wa)kt@
z-Z`ek>Yf*seXDD=YXpyEzK_{dOmWrm)rnaB#+MY=tR9d8-NO22Tl!?xY)8>3&y*hM
z+iw*EFPI+iU0c$dlv|XSaqc5t`@H);)MEl5+<^j7rq4H!0XHl@9u)je@<R4VM|0)U
zL~|M%v&iaD`^o-YA=~Zv_(9A`wU;EzsW0!e!bE<tEq?k^KjA|+;nFVgt&UV+-VFh7
z$7TLe&@P*;Xl`4rh$3A##li>fzXg3^lP=m703aYMJbw%N{(E{C4hZdEg1+|7E=Kk)
z&VS>0QT~@vU4Jl|-<e!?Ot2ci0?i=5mEV|L|9(*Pzc9Jn%$%HmrMLb+_*?97zfBnv
zh7$J|Uk2p2q9_CU+vOjewm)*rFn+7}ll(h|+usoW$wm7k9f#mgy(51;{So0`m}!57
z`eSa8>Q7L4;Qtq>f5zPZ3i`(=KjWW34dDM1=)X%c{uTR=#>D+6?2}(Ph5vo-{|@`F
zGx+B_`lI#W{ArXeiocHX2Q&Sfvj2qnqZy$82}YIj&oF;CJVhBODEhxQBm`&pN0(pO
v+u!B?T<L$eC=ig@UogMPuYSK*y;J|~8T9x6Ab-)dKtP~>{eElzx3>QSVBZN$

literal 0
HcmV?d00001


From 7782771c85a2081d84312843651a34f98bb90178 Mon Sep 17 00:00:00 2001
From: heyyuhao <yuhaoheee@gmail.com>
Date: Mon, 22 May 2023 15:25:05 +0800
Subject: [PATCH 2/2] add tests in github workflow

---
 .github/workflows/cmake.yml | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/.github/workflows/cmake.yml b/.github/workflows/cmake.yml
index 3255493e..78245655 100644
--- a/.github/workflows/cmake.yml
+++ b/.github/workflows/cmake.yml
@@ -46,4 +46,5 @@ jobs:
         run: |
           ./cpp_test "--success"
           ./sketch_test "--success"
-          ./crs_test "--success"
\ No newline at end of file
+          ./crs_test "--success"
+          ./weighted_cr_test "--success"
\ No newline at end of file