From 0c8cad060c4fca72586a9a7fe4c67a97f5816367 Mon Sep 17 00:00:00 2001
From: jimlambert <james.lambert@softwareq.ca>
Date: Mon, 6 May 2024 09:56:39 -0400
Subject: [PATCH] Testing in progress

---
 VERSION => VERSION.txt                        |   0
 examples/multiidx2n_test.cpp                  |  12 ++++
 examples/single_qubit_gate_apply_test.cpp     |  50 ++++++++++++++++
 examples/two_qubit_gate_apply_test.cpp        |  54 ++++++++++++++++++
 ...ly_function_two_qubit_15:00_23-04-2024.dat |   5 ++
 ...function_single_qubit_13:32_23-04-2024.dat |   6 ++
 ...function_single_qubit_13:32_23-04-2024.png | Bin 0 -> 12198 bytes
 ...ly_function_two_qubit_15-00_23-04-2024.png | Bin 0 -> 12012 bytes
 scripts/process_apply_test_data.py            |  32 +++++++++++
 scripts/run_single_qubit_gate_apply_test.py   |  26 +++++++++
 scripts/run_two_qubit_gate_apply_test.py      |  26 +++++++++
 11 files changed, 211 insertions(+)
 rename VERSION => VERSION.txt (100%)
 create mode 100644 examples/multiidx2n_test.cpp
 create mode 100644 examples/single_qubit_gate_apply_test.cpp
 create mode 100644 examples/two_qubit_gate_apply_test.cpp
 create mode 100644 original_apply_function_two_qubit_15:00_23-04-2024.dat
 create mode 100644 performance_data/original_apply_function_single_qubit_13:32_23-04-2024.dat
 create mode 100644 performance_data/original_apply_function_single_qubit_13:32_23-04-2024.png
 create mode 100644 performance_data/original_apply_function_two_qubit_15-00_23-04-2024.png
 create mode 100755 scripts/process_apply_test_data.py
 create mode 100755 scripts/run_single_qubit_gate_apply_test.py
 create mode 100755 scripts/run_two_qubit_gate_apply_test.py

diff --git a/VERSION b/VERSION.txt
similarity index 100%
rename from VERSION
rename to VERSION.txt
diff --git a/examples/multiidx2n_test.cpp b/examples/multiidx2n_test.cpp
new file mode 100644
index 000000000..0be183ad1
--- /dev/null
+++ b/examples/multiidx2n_test.cpp
@@ -0,0 +1,12 @@
+#include <iostream>
+#include "qpp/functions.hpp"
+#include "qpp/internal/util.hpp"
+
+int main() {
+
+    constexpr int ndims = 3;
+    int midx[ndims] = {1, 0, 1};
+    int dims[ndims] = {2, 2, 2};
+
+    std::cout << qpp::internal::multiidx2n(midx, ndims, dims) << std::endl;
+}
diff --git a/examples/single_qubit_gate_apply_test.cpp b/examples/single_qubit_gate_apply_test.cpp
new file mode 100644
index 000000000..aa65b0187
--- /dev/null
+++ b/examples/single_qubit_gate_apply_test.cpp
@@ -0,0 +1,50 @@
+#include <chrono>
+#include <iostream>
+#include <random>
+#include <vector>
+
+#include "qpp/qpp.h"
+
+int main(int arg, char** argv) {
+    using namespace qpp;
+    using std::chrono::duration;
+    using std::chrono::duration_cast;
+    using std::chrono::high_resolution_clock;
+    using std::chrono::milliseconds;
+
+    if (arg != 3) {
+        std::cout
+            << "pass an arguments for number of trials and number of qubits\n";
+        exit(EXIT_FAILURE);
+    }
+
+    idx number_of_trials = std::stoi(argv[1]);
+    idx number_of_qubits = std::stoi(argv[2]);
+    std::vector<cmat> gate_set = {gt.H, gt.X, gt.Y, gt.Z, gt.S, gt.T};
+    ket state = randket();
+    for (idx j = 1; j < number_of_qubits; j++) {
+        state = kron(state, randket());
+    }
+
+    std::random_device rd;
+    std::mt19937 mt(rd());
+    std::uniform_int_distribution<idx> gate_dist{0, gate_set.size() - 1};
+    std::uniform_int_distribution<idx> qubit_dist{0, number_of_qubits - 1};
+
+    std::vector<std::pair<cmat, idx>> gates_to_apply;
+    for (idx j = 0; j < number_of_trials; j++) {
+        gates_to_apply.push_back(
+            std::make_pair(gate_set[gate_dist(mt)], qubit_dist(mt)));
+    }
+
+    duration<double, std::milli> total_time;
+    for (auto gate : gates_to_apply) {
+        auto t1 = high_resolution_clock::now();
+        qpp::apply(state, gate.first, {gate.second});
+        auto t2 = high_resolution_clock::now();
+        total_time += (t2 - t1);
+    }
+
+    std::cout << total_time.count() << std::endl;
+}
+
diff --git a/examples/two_qubit_gate_apply_test.cpp b/examples/two_qubit_gate_apply_test.cpp
new file mode 100644
index 000000000..5f398ba80
--- /dev/null
+++ b/examples/two_qubit_gate_apply_test.cpp
@@ -0,0 +1,54 @@
+#include <chrono>
+#include <iostream>
+#include <random>
+#include <vector>
+
+#include "qpp/qpp.h"
+
+int main(int arg, char** argv) {
+    using namespace qpp;
+    using std::chrono::duration;
+    using std::chrono::duration_cast;
+    using std::chrono::high_resolution_clock;
+    using std::chrono::milliseconds;
+
+    if (arg != 3) {
+        std::cout
+            << "pass an arguments for number of trials and number of qubits\n";
+        exit(EXIT_FAILURE);
+    }
+
+    idx number_of_trials = std::stoi(argv[1]);
+    idx number_of_qubits = std::stoi(argv[2]);
+    std::vector<cmat> gate_set = {gt.CNOT, gt.CZ, gt.SWAP, gt.RXX, gt.RYY};
+    ket state = randket();
+    for (idx j = 1; j < number_of_qubits; j++) {
+        state = kron(state, randket());
+    }
+
+    std::random_device rd;
+    std::mt19937 mt(rd());
+    std::uniform_int_distribution<idx> gate_dist{0, gate_set.size() - 1};
+    std::uniform_int_distribution<idx> qubit_dist{0, number_of_qubits - 1};
+
+    std::vector<std::pair<cmat, std::vector<idx>>> gates_to_apply;
+    for (idx j = 0; j < number_of_trials; j++) {
+        idx a = qubit_dist(mt);
+        idx b = qubit_dist(mt);
+        while (a == b) {
+            b = qubit_dist(mt);
+        }
+        gates_to_apply.push_back(
+            std::make_pair(gate_set[gate_dist(mt)], std::vector<idx>{a, b}));
+    }
+
+    duration<double, std::milli> total_time;
+    for (auto gate : gates_to_apply) {
+        auto t1 = high_resolution_clock::now();
+        qpp::apply(state, gate.first, gate.second);
+        auto t2 = high_resolution_clock::now();
+        total_time += (t2 - t1);
+    }
+
+    std::cout << total_time.count() << std::endl;
+}
diff --git a/original_apply_function_two_qubit_15:00_23-04-2024.dat b/original_apply_function_two_qubit_15:00_23-04-2024.dat
new file mode 100644
index 000000000..00819eb9b
--- /dev/null
+++ b/original_apply_function_two_qubit_15:00_23-04-2024.dat
@@ -0,0 +1,5 @@
+10000
+2 176.545 174.842 183.596 178.394 176.159 174.526 185.764 196.628 190.345 175.515 179.988 176.591 174.93 178.4 190.857 176.199 187.605 182.297 182.269 175.799 
+4 189.184 181.322 183.454 182.137 183.592 203.339 247.339 179.266 176.318 200.083 181.82 182.155 198.882 183.476 180.763 185.164 181.143 185.028 179.992 192.943 
+8 320.075 263.933 258.172 281.16 344.28 336.755 266.493 267.674 272.531 259.63 257.905 261.764 257.426 259.785 260.897 258.122 259.11 258.94 268.418 258.673 
+16 23787.9 25444.5 25880.3 26085.7 26617.7 26229.1 26678.7 27656.3 28167.8 28822 29575.2 29536 28568.4 28914.4 29033.7 28474.7 28631.4 29014.6 29148.5 28824.2 
diff --git a/performance_data/original_apply_function_single_qubit_13:32_23-04-2024.dat b/performance_data/original_apply_function_single_qubit_13:32_23-04-2024.dat
new file mode 100644
index 000000000..79843bdfd
--- /dev/null
+++ b/performance_data/original_apply_function_single_qubit_13:32_23-04-2024.dat
@@ -0,0 +1,6 @@
+10000
+1 177.045 183.401 179.482 177.491 177.767 176.061 179.892 176.771 178.809 178.344 178.731 173.207 177.846 182.742 185.616 194.434 175.579 182.799 186.081 177.499 
+2 174.65 177.879 177.104 178.997 195.871 185.893 178.443 174.799 178.844 176.05 181.729 178.526 179.879 176.963 177.265 177.043 186.762 188.596 176.317 188.308 
+4 184.946 185.677 178.205 185.135 184 185.485 182.079 182.476 201.217 189.12 181.365 183.386 180.097 182.053 180.133 179.742 191.665 196.478 185.7 182.02 
+8 219.992 232.964 217.757 241.844 238.279 237.996 228.896 224.703 218.453 243.654 220.595 236.702 218.217 271.035 223.113 252.001 228.766 225.088 225.37 243.18 
+16 16772.6 16456.3 16817.3 16767.9 17139.3 17423.9 17191.1 17600.5 17299.9 17661.1 17426.8 17547.7 17839.3 17510.6 18148 18134.3 18854.3 18829.3 19560.1 19014.8 
diff --git a/performance_data/original_apply_function_single_qubit_13:32_23-04-2024.png b/performance_data/original_apply_function_single_qubit_13:32_23-04-2024.png
new file mode 100644
index 0000000000000000000000000000000000000000..f4107ee5ab66b5327e6bd10de08de0bcdf046b61
GIT binary patch
literal 12198
zcmdsdc|4ST_x~jottyq0rA2jDLP})qZXptp$WDdqOZIK3aCZxLgtAtW-B_Bjjv+-F
zOBxKu+F&ro62>|+-!paJ&;2~V@Ar@A_viDxUdp)Ux~|XpoX<J$bI$vm$Je#gwrmvK
zh#<%o^(&Wj5oApef^ZG+u7@i!tvy5V>zw=LTkbcUt=+vW-K>ynmhLVN&h8GjxA%Km
zx!t#Qb~-H~Cvir6|2=nim;3UPl8*nlLBiS1MpE?fArE-T2A3;__Ys88lJk!%QzgR|
zK_vOrFa4_LokZ_F<*?5tWno4%dEK)k-=1^rZn(Yf!$k-FH=WmYgX|4zO%l%7e09Rg
z_&5pNI+ef1K>j$H{OpF)q3Gs8oa3557_mH;jwBx}X%eaDdLUb7+;c+5oqrdjGUPOW
z^ufc?JI<Z<uM#~UbR|NL{lZ3WmU8+uAAbtfXAOb`oyvb?3eq9F!x3HtvDmX_KU}or
z4T4Ld2N1YpPw0Qa7mQ=O!!^Wq7@3>%RXj_N)ylYeJoYUMkHg^x*e0Y1_k9?$qfBB`
zFg4!{m)+ajdwqP&vjXlmIQaB>MQ;j0(S&qUIn_WuSsyOI|H`THv~q4@lf6=7xoL!S
zr8td7JH9(SBhG3^uam*!OW)cUQV!m@C*7fMv%1XL#n^>+v-K^z+eeto3sJQRO(QQ*
zvkP0`p==WEmDcIcZ)jF=0faLv5t*W1*1sYXH}dok<GKgRd!0%@?*Bsf1aaM|&BE%@
zj^ysnPTrH(lQVR)lDf0xA`g+V7qe`fowezn#u2I)UuceeHShN@O()vImSQyyx7*3=
z3cn1W(}`oj13O22pMRK~bSe-pGtr9GxV&%QKD)Dt8X3fHGC7gpA4Vc%I^Ik2#fh{{
z$VDEH{oSoBshf`NAdzgceG-<!Y`we;Dk>^qr|?rVn@?Mh${}&~@RV$B@}_BJ%amwX
zEMfSm=vsKHs4G`D>zN%U!GR3=66G}xPqI!@ClM+aPU@&kHhy?W^+_S!JEN#>maaP^
z<}>VeZX!AkHrU5tF!b2PI2?V~#*kpDd~&3XVTV<mr|rg&2k{Me-%A`SfbY7nN+qP0
zVYfTv_eJ88@q{MFzA8g(jwMUk!rOQjUZx#Qvcnt+Y)VTz1Rw96d7EwZ)iIMe<{e#^
ze+NIJm2o;AHdKn+6!L&77fJO&cNewZTN|a+V%+(?v7#v~f~>MU<)_X~tBEHJYtaRS
z>i1?^J|#Xm-Qwul%~X9BrBJ*dVZA|dH{F%y-bc%f@M@clP^~50l}>r{bYGk!nX|3W
z!$lx`_14!W5<JfOgd$auN2|80crSV37FVIE(b0R~1nbV(v((hoQYdm}8N?abBRrr<
zT%=W_ZT<ZGh7ym(G;E$I5GM%we3>Uv4eRpm6N<D#{8d&L*v@4pjN`Ey1zM+HG``@u
zfs$v}qv099mpCR##OVgg2!B#lxHH4{PQB~;@frI)W^V^QTsyGxb}#DJX|iSfDIK*w
zw8ctCRGaT}2F`-aDEqamC^z_Fs}ETRwMUK?(26&!bG(>#SNg?izALPiBK@2#il~ki
zX6{l;5$eTV<3vWlh@=6#3f*qE<x#I&nUJ_wTI+Lf^HDr|l~0m|9l=i?Bv%s0)`}La
z&P!J+J?U1#8l6O6|7Gp0=q^Of!+nu3C`mh*vDd1{cP|Dd*jM7+w^ZLg=;$3iuhigr
zRDw51<>^ib-eET%i`v|FyF*?dnENm)*Ir$p2D>W3;R41HCo&eSqEhn+rQHr=LB09Q
z&Ie1N8g`w$VuYyj;dSWK>RZy1WKpOvYqF<{)C}n<N-F8DAWC?V)-{D2MG=?LL_9@=
zqyu1}<W}F|+#{cQ23cm4@13dCNdIOS=V5(CPG0`ITUk#Qcjxr#a_tk;%V{s&G83nP
zGFEcaJ4b(a_n}=I*R^2jFH=|yIx2!CJXkt^_IL9P(Z*(JJEf$O$fI41Ft4HF+>AjF
zrQh+_qKrja_6e_FedvNZNqI0->UhtjiOBHR`W`4M!!~3{5}oXcl`WwUqwvbB)Mf75
z67jP}c^zdqx4<Y{Hg+V8i|rxMHNx`E)H@DVwt=y)s+bX_T<@wyKdMIeQ<|b~>%Dho
zw|tWvrGt>kB)U{+x|s;~IVY7=DLEHWv`n@kIAZE0Rr;+k(G`VSmUc*Oszm5RN@hgE
ztQp-s)aw@}E}XUR-r~RzN`9PAhT1Om*=KfvRPQ2N$?lF&z4O|3)kCUsb_(sYh*W^u
zuq*v2dhruuYRvlTLAG$d^D~_c@w#=XgSMf3zO|0>SNAQXt`QLsTpb?g(<3I0PkUjR
z2e5!(8zU(XrZbHcoxg(-M*1hr4~brB*~D9mK4rUS&qCF-%IFh$?_{rTf=bZVg-Kyy
zVIPc<<3I+gtl1AAmB@Zan>$GE7Mm)z3_QuAF+|^#l$78!(%r0Iij?Fn;YP!q?&IM|
z=KUC}R*LrQPII=4#Y{=^peHSrdNAW0fSqgG+0&KEN4uKUJgy{R{#1_(|8H-M^4@A5
z;Ws&NXBa2_qOsakgLB?9ELNaPYZk41+_U1dbus*36NRc_#vIjtz`ieJGqeipAI6n?
zmzkmYOoV7`vtM?uXHAVXCDu$4+ob5L0ZakI721a)aT&pgthakv6{y*WHv?-JtYJYf
zOEi1cM|R?g5jCdVl~PwtAyhuTT`EDwswq<|uIyV`wGAaV^kbXJ`t@0Q$G<CT1KNjN
z_%)xFeM@|Nygi57NB8QGPq|AL^?W-2Oe4KH*F&GuC&0tQ6R(Uus-5ynY75`$Gs5t|
zR7_R>N>h<mMpuUfB>-DooTia}kWeIMYms_GA|Cvwxl-Mt2Yc3p{YIjDBuh^t&I%T^
zg%FcZI}TtZPBC1mFBTFJQMF4jOe#Socq*jc4Y(R3>h-O95$U9RwmY0{cPfd!BK9EF
zv*IMl<6%XHQ(xsV_@E$2n=k*qR4KNuWq8)d6z40msGdQrQS?vi%VqWYvgk7k3JToB
z2=bTE+L|{qnhL1?%!hZ5#|oAb?nm1$?+X9;a?9$xg1>Bk=|eZR*0i`l+38$%geW%r
zU13YAwRiFpqxj-n*?Lz!iFiPRwlVi-4R@clP+M~8xgLrQSQUr3`FDkd2gX*n(eaFT
zTvk0RwsEhl`~$bt2ETU|@<BIGewwWe=F|H+G_?>|oJ+LDV#2FhG}5DC13t$qucZ;k
zp5)%?VAfwI_Y{d^j80lf?}hLR(UuclE3UV6t7z$xUF2rTjYYX3TBpNjyREwE{?YTj
z=2MsZDs(LUk-4xdneBFVX9R1s^YmlpvTZKv7Vyu9?ffH<^=1^$rI6fQ9Z6!(sDwNz
zuBw9L_i%TA12GI^ls6D07J+&=hkqrZTkA6`DuZWHm;3c%G>I)($1ZH1QQmI<ujA~y
zvJK7jq6<cCDy^pO3Ndw4m36|oEtK;79gjbJGE*jFAHF(f-D{{-MCVJiQd;84W~E4*
zvw3R3ZJZYOVT2=$SjBw;F-49mJt7^0W>M<>Q%q~qjIYCv_LRli0nH+2>Tr{6zfhNF
z%$!A2P7^aAjp}2EF>(_+I{GTwda96Nzehx*4&Bu+DQTDB^V~{7#W`RhT$fs%`-Cru
zwN)ve(BPUhP><?C`AWwVY*8rlm6a9l)NCgA`{~He{!_^w)3aX9UXwk0J70#@YsI+R
z9h$2hyKMLTE-P(>%0{D%?j_iUg4S>5+-#@K-$)%eNf4~yW_XnsPp;*tkHTa!UzpkG
zNE9V6tqKpgC@hgXd=x%eVjT%Z;`E7$C#w}GL8p!z)EI}AT5D77SdC#lzISu82o@4D
zH&fX~m(j$#3WE71X6V$Fs^=d{69NMHDx4D*{P2byV^+@gX`55AbkB;@sLtos#>STt
zSMF#p(@7}#uc>Tsjaj+K``MEqLongh?sH96306P+LqKKnQiJ4siT#?aWA&~){HT|J
zAn-{~e0OY<^}}C1R$>WOkg$2k3tJT%3JVF{gJi4BWU_HfPydEEjdY6w0tDEk-|b4s
zzn-7>8o@f+qfqR&jAlcQW|*l1MO_ZldA#Ym>Mt}^BPnY!Nit1xkmbRJUZ;h&VDZGJ
z$cfmeX36=p+@Ic^J$v?DNlAvgyl`c!WHHZ|DVg$#;CP$fX8$-V$U79|@(Rn!_Ey$G
zE~ZI?xPL6M>3Fi%iFm@gjCK1WLdWNqI~j?O9LfBSph81!))-AU>y35sUh%6X<X9`i
z8oyU;_R-!*6Oz(AptwWi>4vjs@*(L;^TkOJiU5WlBpm7Of06(Y?k++$!h2Q-Z|>|x
z**5291T&d>g}3xFh^-?dZy5^<+Vd|ZvYTYjo*gRB$6~P$QruDG1$3~IJ=;44G~!RV
zrku*xhh~bpfVK_a_NHnTCoiyDd}zkyovugScDtmh-(4yRiQ4z=Jw7o9^_-jO9#8oN
zV?^(9Mbx%*eYtn$f#u0n<1b)3Sqhh#uCwBVjZTiDPVZ;!_6%1pE$}sc;|rJaR~rx`
zCSkTtoVru6<YgUUmu4a%>-nW`c7z&5Z^x&sp4!~9+^-&LdJkh%vn(X**_ladzAK#(
zXVZ3{h@Hs{ez{U_Qj8J|%j<tbj%+daV%Zz1YzaQ)Z$KBLS+>b5JCjA{<%=r?6sK;k
z%=5`m`r9Y)+RLP2AM#U`&D!k6+qGU*QAUz!*jKi8o<c{xar!wj1rsq7&UH<wI1hOO
z*0XP=QIv7<YgTIb_poc@&R(`I4po;iKJ3!)<+0{(tc5ni6}A4%il6`{dWtV=$H<wq
zC!F`a7@R-X;HkVSPt9{~Nt`I<7f|J?E`^jw@Z6ia(%@&*Q>bfFW<(!i(q%z5{lp7R
z>!C%5j^|vj3l-n@G1OknRfh}@gym*V+Sk!VPl*!Tdrq+gtL3Jhui-st`Ys0?i@VzT
z`uaMjP*3>r;ekC_i3(nmN%A+Zqe)g|hDjB9ckA|#j;O_aqy3s`jHu^PRtl5bkUQU6
z)6}?2FgC@8PzG%3U)ufFsU+vpLDxx{U%{%w!os+{o`N%v);#3dl)cl^sLoNeRQ#{!
z#S~j+-%G5CsTl24FNXTWO^DjKq$S-{U|sAp{(e~p%DLjM&FlB(lW1#^+D8~Vcku*y
zcNU8M+tRH%%2uvz2vRMTv<5f;sm(wCSbH8S;@UmTYx@P}zYkvxU^5qdDG{>Iu86&t
z;6V`0>RCl46w3s^+>Ul1{{2OO<=W~^i>#F?3>LfFb1ziZ>-z&1k^`AE^DtA@UwDHA
zrt9`OU-89o<>lvxnruOk=a<kVVt2QZi;(PguIV+a#HyT}9L}n@>xKv=jfyjHkA>7f
zm*=OHk!5pfKDfLwjoSie7C(1<Efnm}-;Y)V-7hogJPZ<xC$%aC!S5r|<vcD%M~@zr
zKX-0}$bv6!TMif8x<x7IMrT)-j`WMofs2B22(nFWL(6PsK{&7R!KYuuaxC@u2D-Zw
zBBP>~j(|>Xzb4U+z0ed4SRL6py;EOb|6Nhh4$kfnf%nG-Qr#<h3$DT+w|LY1ki@39
zf)Sx_3`h%(PEOKiHY3OZznQN@gO;uH3sjZyc<0J|zWnsL?$}G6RwsAmb5$bw2M}by
zEnumMZ+Whtb%RY@SxQ(bi9yqdQT^>|Q@DU*vjdAw$1Fg~?dE~A51rU$O50dsWB4dk
z8(pPjeppuJkzUt#4r>wg1?N3BHfhev*2y@`q)xJi(Y`a${PFR@95Um=LVNK=rpqd&
z6Ym!tD96$<S5h!B(YgBM=a&mqN^QS5anOOI{cYvuRWr-!T4BZ85iaEXsB&rCS8@YM
z#;G9wky-DS)f!COL4uLd&s$0PVsm5lQyO2hwvDQ=<|tl!`KF_{ah?<NpX(!+FRVQ2
z)7Y9a5^!Rve3Ger-<E>BpLSjYC(-)7)zA6F%A~PB|EHuME4``oJp2lE1;5vG@70JK
z)u7h@G(KkOa_727O&+H0G><MU<bjY{-Zj+kw1Vw6M$B+YrMBd>+o}8RTJPLjo(|cI
z+fCTy?47w|79<ZI5i&Lp;)$*dTno2!M;w5rNY~|VN!LAwA=Cdj&M|bAG;J=gir%L!
zq9P)Z{;Z*)OJLwN6f?AJR1rKj@_AIG+O=+7ZE~{k;h2U28%IaS!Vp;9WKA0PODn?+
z-<7R=eD!_n{DE2#78T`D0T=CElf4!}N=;OG&eyxDD&q?aPkmpCjok%qx^)#^`XQo~
z>kF11T%S*~yL0EED;#dM?&;+t+tjn$?VuR%?&)c(so}cAF~aotX~mXJ6~C2^&--c|
ziIG=U9&wjKkW)1?<mn$LdgS1L4a4?tF@-PRcvj@i&(9YSK$Ew4DgTUG0|J%;+ps1Q
z0-hV-%S~4PeEWfkw0$zgl~mWz@K{Sr>p?xvO!9*rd~1D?Z(ygVr{^=xBM(zC7)-%c
zjyV{;pH-BcI4D8x>r0VOZtlcG>^K*BdfEh_%qQpeKcEa#UhZ4?8qV#)oUhX181=NM
z{;NJ4ORxXZIP78Y1Jk7AN$WwtomyNAzyn&w!;(^7bv)J9#aX&pVq@(xfvnOXxJa-&
zhUPNa$pNnYfpURFZVrVi^v@XmA4#E(z)?tjrqo9bm*)fh_I2!?1G!@=aG|1r0D8W)
z@#z@X#Juk^6~(eKEpxlP<Zsq)cQaW_4gQDFjK(n%ocgQs%+La!S)ioC-H%18wUw7H
zMTscs7T#)Ot}N7R8K6n?75z-DE(d)m`8MT3aIx_D8MQfqzZ-%Md2aSAZ61!;!{0x*
zy}iAF0IthFuqgcWMzu{yS$&i!jqHq1_V$G}l1}d6`m!)Lc7r)0i*E(V53A#snw*E^
zIKs2$&Dk<T<8WKt#W!m#US0pv(Gdp;@Z9v^MPU(<@nJ%Ixz|(*wxpz1f)h=TUHT??
zzCZSa0VJG#K2+h}QrC;lRZDLtdn+<1lskYxIPaNnU#VSAB~a=fR)-SXcX}XLoRBbz
zDhmapT^s}<A<x1((#D7x09(F`)2|>Hij6hI9OVxyvIVO7rI#z@_r}J?P&(azA%Gy~
zcU~IbIn7{T^Ud-3IR@zU4^>pWr_Y0eLM{w`rIMUVctT?uhTS0lcQal1N?ixBY1h+n
zF7cfk-GtOyQj|7L&S4yZ1HL2S#j=*(tK>bs-@AwId00J5Z~U*@%hyXeKtVwFRFIoj
zocJ3sGcXvLG&IMIol4hB-Bt>Uo?D`(=0LT(N>;=()5!E}D29BP{msphlYKaa5<$@B
z?gP49)3tmmdc|GOC~_H2%U8s9@Z*~~x!$&8>L7R|WwWBe@z}O|XPPU!hMJQVmTN^>
z+Y}jo0Gh}*G2M(Hv`N`({haiLeiZ7}Eb#6(CE{Zo$u%B9F+Cod0@4UFQ#+$@+l$La
z`b9{8kPey@+i>^6v_YOe_3RL26ri(go@vlAXyD^*9bmIFSS)S!oOZmt3dsX#ex8(K
z@vFeD@-tOZ_an$}qLKaEE7lk2W{r-HURjvxS83;eenm_|3F*e>@eZ5A*@gV3hve$#
zK+YI1G6$QSF1}*Q{ahnmt2cQQd`;Muh1|;015ROnip`n=V9yAIyIfqKmTwsP982^$
z^YqH6S;fOUHh_&Be3V^BIS4JKZ($U#5BXgFOW&=2W@Wf5$(-b}diC-_Bye#sZ24Qb
z<dFo%U~SlPJ!-Z7oB_n3g35!g8%y<bY;A1R;LIW^_W}Y;p@zjemgS!yd2sj*hqJ)8
z>*eX!N<a>_AA8V)|I^G&A=GT|@OVkMh3;u(7ze)Flqd_q7hI`U;&k=pHxR}6>Zg8N
zS6yn!z2*mUsSOF9hi!=enfY?E8KjnzxFDM^A|Q*DrdS=<T)+yw1fm>p7xMD*8&kB?
zb)kg)JnZI6P4;?!g5$y)E*(Jvzc(C(UJH;|p{&Jvm6YWbVurQ{n-u^`frS0hbuhDP
zx8FX<1*HZsx;{Kq6crVvnX`t?4{l$xSB>{i_d)g;IDiyC1npb<eV$~_r#Lz9JGc`N
zsNm5-enf@6dS5AU^%l-dTRTN@?;^RTrY7q;0B}&(rI90lfCr7Um@D8dAPe}|XX-3;
zPB@-lwfJollTm`-QlSYqcZLti1f}|ae|%O}7K6u&uDww2>OR$nE4aF^S_XFhLxA<w
z8QQ7i%!PhF5fE0}J#6}QAr}|2O+^TQ$Ysj_=)8)_AO};IYg6Cc7(|n3u6IC?>L(ZY
zf*S7X0#Rq~KlAYtQN@4pdb?eF=ZCFeF*;96Hw<;;8T&5JWjXbfzH*n(Yiha%i4vix
zjqVA^TeUoXLdlPz35HEVL*MEyN7~<eWP<AJ>$CN9#@*$Ce=HKeI_7Q^Ng3_N*jvWS
z2HC^9Kdc-eIv${t){*x!*c!bu=R~KVkH7mP{&S;G0sq`r81w<B!7x72aR!|{xEl;|
z+l{AuK_D;2=o%n7nM__@rlCEU%N~GO1y{M@-RYl9x#b-0-@if1Hzi{1N$)V2^cLk6
z%@k@wQSM*2EFc<RQ(1Ev0}CpwnJ9ThMX2UC6dr`-7h(CzK<09m^ZKstZrz;rjai(_
z-r;wViK&?wvvtbqzN?Fc%5b>qzxO;30y}iw-hgl_DGu5TFGddV@QGxC0|VvxD=!>S
z_1wX|h&?7js9xu6;uc?SU>Kq6@+XLfnt-vL_rENS<fG9vJd1{(E5S2QmXa5yr>{>{
z&PIwIdwS)Ofe_r(c-d{wa|(qru;5g3Lj8w<-v^u>8IC#%oTi71%gaPrXW&0gx@n+?
zOV`Bw5VUE<TfO7yvi<%y*eiH;4o44g6xyiem@3w?y-$jl?JU?G1ORXnXtKfQct$sD
zNoQw+NHvXwCb$Zz>q707f`StPssmG8eBX8!T#Q>OaI+K?d~v4lD;g)?U2UVygPR{#
zmv!p-GCn@uLw(U?9oxcbb{t8l`cadOFD`l?nCC)*!YPZwts>PnEh)!RTL3K7p;Mus
zr%!Z7t&qza(Y;PU8<}@GFdC=5+}vKr#Ka7^qDhHYl|Zc7oelPgsIq@1(8ZZXdD%vJ
z<M31>2u3}N*F2I2bvrcBGR@FY$xZ-D4R@s%4ULTxjjwT@-2q$D>Hl5`3G)nlicU{n
zGOg)yI9pox#l|oWIA3T7Mn*=enr-riY@0dAhvJ1YRO!Ngu&LxL>#KEYD;5114&T_R
zbAnZnXZb6OZazgxb@Md!8oDlZ;=~CfL&F!P;5b3Q{C*d47g7@qTzmzM$-}sfcEfH<
zU#T$~>Fa&hkm*=)!fm1R01|Bn{RF0qZCdoQ^u*L<o%`BWrtm89OKmg{2>Ke)5Rc(h
zzVLgCOdk9b1m%#eD`4+bPl4&&&Y(XX>D;<i2!*Lj3qZmG;|Ro@cqYeWhj%?jyaCmE
zD}81gZEt_ofwXm-&r{x@VK-2Uv5nSgMF)r&h&BJ?;Fm3%P|+iv{2T|gJ;*8VkOH~?
z4YL12HT*OB{s(<AaRUNA-LqB#oWvp@`ggV^<nG?yEfo;1gUSvetPdWhU0C5Cq``eP
z<e47Cbi(laAGaR}-0sO`N#mfJ{cBroEiHF9w;deRzn=~B85V(13qWXg6S8@2>-Bf?
zf6f5-7nh|SFwkHr;MxVz4|xANmQrHVXICMa#A}WY4r`Ty044ME?H!^9GDWH(*@1G;
z-os;FS^$osQ5{XXA^qaa@*?h8M$o2^_aBJ*Iqd^pBoFY%w4Yv|l;CM)bs5YsUiT1y
z#jUs^9vA(`mo6l0*?g-AURnz!{;JFBh~Rm3wA0#|dPF>)0HNmqCyq<0%pqpTU+Kpj
z0s(AK*nh!H1!o`s-fRf){xF>YKhZ<wPSXm2*+tUIi2m_id1&;Br&ls>iGRk@f7Wn0
z=8sIzcY<F-+e1vfNc=?4ZSbMp&{+Cq`zlzFA-oGFEarVpqupIxTt;jMVq*A{lakg(
z`QiF5*nzrBHK9<u@M_x}C#Q2bn*)u0w{BnndU1YU!R>I|CRLB9K4zVlO|a1MEj*wV
z4ag)^|8#C#VnFlAp|P9J1MN8m?PKn$j5U$tTqefGM2PSD1_qKo$3Sw$ziY-W?pqxE
z%?C)ew?9$3;X}Mex~dkwsOWWUZ0sfySmm}OX<@;Ij>a&K;N9~hSt~sTB8p?uC%BTH
zSs8LNlF7+@ae{6JQ72D-Q-wIkNu=cw4c#L#8}hk0))U;C5d>0eUP3s3{r@AulNs}M
zJjg9&Ni*N&@vc`a)~Xuwdjx5|2J_6UBAT~Ce7`hr5O5XWfj%20qE@t-+EZfQlk_3l
zG7wnvKud+)+Ph7dLxze+NaB}%En)YF@S>vU`s-E38L8yu?}?%pvN!*>d(D0*ohpCQ
zZpWQJ^soL&yScUP3^M$&#Pt73yPdzp0mQ3=Q!Vvjc*D2caeIMrE!Y{jR7Z=K7cOZ0
z`6!6nr&X2j`Vc2JH|{h5)ZuQ!-P_awM5}eZ=^7VO16T9Kk~apkDB3E4@iPe5kCzUm
zE%bHZ`fcX)-KNh>88?tTrWW?rajgw}fBp++&;KCZcwN?W<RAN+bUT07b>kj5aH;<T
zr)7hY>OTY|l_vMED`*4xzEBCyfnR#@t2b4I*sB4B&PUBj8KQW?9ycUzeMZos(T`WP
z*{w^cWfgO_Ia7_w9{YGzaw1x4UAJUa)#?aOFYCVF{#Es_C{cDM8RbNSxAA640}C~W
zsBXsb&8V(?B!$PPaA-`yT_h#k^6|n7)KX5-RMp2Z3#3qu&*IUSm2J-N%qs)tC#Lnt
z^Rr?=Z>ClGdX|>HMl7^=gSNYxO+B2;UI->WZoRbfvA;hxiVsz+$%E*qB2RdIMZZNb
z!tp<u&VRke;b=X-dqh-}OTIE7-4Ap*iTPB7Sioe)_FtziR<;JtEOf^{Qf4lis>|<M
z!>osMkdO@ml3;E1Z!Z@P^xeY`PlJ|eeN(4>Pn_G+)~`{O?}u-3IVu~JULz&7eQ-g!
zhjrh0gYAIz8eh@ZQ4C`dQQQKSvo+B2Ejys2OHaD0LY5DiG-SR;rlz=Oq`N|8-f2u^
z|MH0IYUI8(1M3>?xxe{cp+qpagYF$d+&12dhT;g4Xa{nf*ttl4J0-NK`l#qv=TIYV
zYqb?F<aP66;`xN8w=hNmBR4RiRCsH1h1d`8tKuoRB@Q{5dNf>4YvnEDMEfP+&%bz>
zLgw%p`o3gGkBZ)bGr;8mlF1`jn6A=2_JiikgJ}q8N){4Rlaew?v@L*H9Lg{wh>wpC
z!~vYKGCaO>W4Z5qc3tlt9K@FjuPKRA4OoJJc~)CJL8oUy)y?!fn;1)xgV`*|89BKq
zEAEO|n;iz)*;>s0zv(b-Kpxe8)vd_+9;wZZ(F|nw#*oKv6B)jE4rQ5}+o-z(qSf6!
z0b5@Fm>jzf-MKdw;?+7h6bMGReasKZ(8+1nJ0AP4sw!VU=c5$t$ztr}uuM?+`#;Mc
z=UQzp2<Fc}E#9aW81zb!7TFJV`JiaYHRTN&NM#VUCG~1_M@Pqnd;lm-+R`i)xK7LX
z?Hhs94~myGABVh)w9mQn6*O<M(DbNNi@jd`Cx_=Zh?}pTT#+vFpDCFjeEQzh{I4Ap
zgW6VZZ!Uv;lU6i}Nh@o<<Kf|PtiMw^PLdTDu-C9^(U@U%cggyWT}x_9V<%)<e>X?A
zV*rD-RhaI!c;Cq3aJ-ed>Yq<3o2&(^;sh>{^YH6mO5Dv9B=pn&!&Y~MWJLUKYv)u#
zj(!8d;YUdg6vXaZ^f{Pta{}I4KgVOl-PJ*%{|#YJJtM}fz4vcQs74pZlRg%ADv~}b
zX<R~Xuc--!If8T85jpJO8G=APW1s$yxxVCdKT*}_urDxYHL#$SK^b-{yZYkn@wZuJ
z*MISjR#b=Cu)fyT)-}hIO!M?%5-3i#!T6U*N2(9;<ei@G76XZRp1>*R97|WndRxwF
z_YC7;_6=HKoB<{nl=JNVYe5=y%DFjPA<}Wr&9&8kxXUNsP`FJfl2CNi(=@&h^J<Wx
z-!7?@mlXO~NDsZX1EF*VLe}*1B(Q)-8))wG8E+%(6|{<<qzJTlI=~3SRIG-TLv@7x
zKCN<pRzIa9kO>ZDLuPa?Q@TTSKFS0gu>A*bpY?JRs>3$zc3i?qf>C+?svNdA_ghCG
zbgTJ96nRCwDgZ&+&2;_At_fCiAM=ET>o)|#U>Qsr9?A&+u|icTCKYXlhAwpSxoak5
zschh?#8{Euor8meA#dTJ>?b`dfVqh_GtGduI&>E+KN~-Xb;Q4Q2Qqn*n9Km+gTD{#
zrFtg0?@eFC-l^#AV6mzo0oWMw{*&#up3oTa)SdkdNJ05SI^;?$cWy@xS@*0DYz4VZ
znwSFL|H>(fOW*Gnk972kD=kSy4Jb}o$Lp{K1Y*SRn2=A!G<?Fs=w<q2-{|F!ueKaD
ztt!QuV~mteccvWO>bU1}(weQ0q5GhWff<^})tzUW-szljY=U7*{EC;?yo!K<t3!q1
z2OJZ|P6Q+BIz0_TJr$kzrhAv9Ss4Nxz|SG>Uwg>YmBTULx&oW5Z$MS}{9boyiDlb1
z<l<M9sTw{aB4X3EeY=QlJx(RSk9&+S73bghJsec=?pt3VXM{6WqsY;?*U4+BMrM*O
z*C1V2C=4At7=IO1WN3RuZd7u0sI`rgjHFX@^mAZ3lnA7{38}rgx%m@q`CDfw#Es1j
z(;j<&v_TFniV>Yk9?CBH5LH{>5x^iH+QQr1%1qwGqx+|8&P+PWe)b{kAu(WFv}Ozg
zom@Li`I)I-`_B;elz3`FW9QY&UME6PuEj(}Y02O_f<OuQQb;-EajRE=G&V^%`=sb+
z+V=djhw07ZyP&B!J}_WmF-`DSMgwu;=<J-T;HY=4*q>#Afr%yJ*yf2(uQR(GvRm(Q
zCMkf+fY!|Be1c&pHs-By9hlNww&~Zgx$66IR%5^hesz0Cc?Ul1nF}J?)~Kevg5iR-
zAy)VPR0G%qsNqJIp6;vbycY(f-&yMYDC9QR9<%sgQ1Z_`<A1A=YlO)V1bO`$h*l`P
z+wHPos0%0xn7E47I2KO`F%>~f5`HQGSWl>7A*F<=BT&S-$TU!AUfxNV5Go*W2VLdJ
zv#2f+e*>@{V&5Tk-|s_@f$q7GyAhEP1<986a@vE5IN*4oKL7ALGVmQX@=`mkbq}w9
z;Qk&mc~jwW$UPwM870y^9}@A=?v%PM5Y>XI$o|RF4FL1iQa|gMe=q-s4C&wPUOs2@
Zz5m*1c;XcO6Fh^cU)H*mf5GCx{{T}cQLO*~

literal 0
HcmV?d00001

diff --git a/performance_data/original_apply_function_two_qubit_15-00_23-04-2024.png b/performance_data/original_apply_function_two_qubit_15-00_23-04-2024.png
new file mode 100644
index 0000000000000000000000000000000000000000..8e06a68163f4acd734dada048272e43016a51ca3
GIT binary patch
literal 12012
zcmdsdc|4SD`}akK7NNRR$&!RhNZFUEJ1rzsA~Keu#nM=lZASa8M3JpRvZaa1o@I>E
zCNUz!#E=+cFqRqHn3?w&-S_j{&-?!I{+`eKdw+jC_h;^pX0Gcxuj4$9<2;V<_j{g^
z)|Lmwmu*;vAc(lxA=6_BA{>DrLOo(4@Qe1B??2&>e!%{d0mm_C1A<)r+>s-$0X|-s
z0Iv(D6)w5^`Cq_X+@-FkzFSS<d_aJYzk!B^_djk>$M~Jokl(U75C&Q5bI9HwK_pxS
z{|h}s<y=4zb$c_@y*8I0k93hSGUrlh;}R)qKkvNLDw277cdcZiPNKr{09E}Omx9D+
zNlRaTbhofitR8;-bw*;{R+C`=W5?;0`wVUw|5Uu+=%TVtLF@wAH{xSv_{X1GwSEyU
zeW~xGW@>8`HbkLAt)=%JQ#$eBknMdkXNj0(?kWa8N3Lzsm3w*3q!D>Fm>1={X;>Kk
zYN8kKfj_m9h!BFT*oQ2JpH>J-B1nzqf4~JF(?^dU-F|o4!QOsp<%9B)s;U<wBO~{5
z8foM*967M-5GKQ;<hft@clP)G`oxmTN*QzGy9f4d-@Y9ldTxt96F$Ob2fjO9knUVx
zUw;No)68dhgomTy%gFS^z(BTRyX<VqhPyXU7krEK?!X68rYO9OuT}6KHk-W*gFo6H
z<80yiC}`+=8GBn|{XX0PiC|}IYooSyA38hj$%tP9w@>wq+1C$IW(r93c^uBU**C{~
z_lDlJ*Gdm$q>qevcU6a<iglxu+S=LeTW*t6*YBrt$s)S%W9f+XP#tX&jNBG)Zmt{$
zrc%l<cx%;SQfytOTK!OndSd<Ez?5slqoW?d!S<~lC+GnV<&KtkJRU4@aB{LWNlW_}
z5?ba>8bKW@IO)|stQYN9{<xW5{JizKt*vdN&7>+rqGfb6|58v;NqKp0Lqo&VWIm%5
zXREgMgxcE8(@rR|4)bdzh7sRzmV=diW{Wq^I%cl4YTKG+h2Eza3<lcO^I`ppRmEs!
zA}(q!X{VH&@b-A~R3(kX`rz0qFCCdnX8v`h-WX3$&rE7ZXXk`gvPHXpi9<t?n&;!J
zyGNAc#$)T;zPT5BdwU0$dRw`p=>Z`gjG><-Tqq|}BgHZ_w#qzf)r>ZSU7=D3LvL<f
z6~Jlt(6)+Aw9Mh~z7xuA@kiXau0vg-Zd&Wfz8lGt(Pi$rHXf)cBa8=*oLLr3*vrx3
z8*k6W%Z72g&qenIv}!iEizK@_Zwyq6oJiWaba9PTz@-T{t&Y}SQE?9<InjKJMzTdn
z?A#?wN-5D%o#tL-PcD3spJJi&$eYx=qkhQn-faC0iqPtj+NbA4bdqqaRj6%qCbS!W
z@U3#l8*m3p3SaEn(N}iGxuvbmma6Nfl|0-@w@+vd>Y1YWC++pVXx|ojk9mvVNNqS<
z<fNYHS$07yd82^>r+(HM+jZ*HDKv$G8|n5Q(qLfjV0&F2cn|AY+8JZTiMZpKQ!5zw
zT2{46NLd!DW~w?Zq?U4l%`WujeC@elZ`<w?;2=M^Vo~i}0IgKEn9DI!j@u`y#5X>n
zx~3vke#%L<lvgw?bCnf$l)9NmZ7ruhx7mhv)qT`Z${#(JyXG=#fE1)wKhaDN2$iRk
zMt7G5=h;p|)||7qtwe`#_L0I*WgQ!bNQQE`!Y-t1Yk1UNeqM`XcG}*NW(>8t6341)
z`S7;eWMd}rTGaZ<tB9Vk*+}DqhO?ChLFv|mmAjZ>XJ+alUIc5^%P?0J5g%~0n-+ab
zoEiL792Ajkkxk?^Qb&fzi1jmz@-n9wi*wV{TH2`k=Q)dFsRwuw8{F2D#OE(f#Y6;D
zk8|7U-)AKZl87ZGtwFSJ7pxzNibwF62@ps~!~}P7@eS;duERMkDk87D*?~7ra4*M$
z0_Xu1jPgd=wr7-APK3~8i}yv|m1ry4)JWLf#uWOKONK`<2#XUzrlzKF+vMlDSjo)l
zZc;&6${am!%lUO<qP>j{B1`J_QgGG{8FJzh6HX&#0?!^8QOT$}7;-I!D#!Cr5!wk^
z9B<zLRbmrQs--bJQ(k1gJ!L>=nz9LOD|gK2rcR4f%4Roy&@YH#lu=n#qtjn@C$4w1
zu4GZlSQir&&YZ>|MtMYIjc-XiHN+x>*w}~kGXK%=VWzbq%6<y_J6$}QYmiBkWxfgt
zGEZ5Q%U2$j?U}banUX$GykSn?n~M!)$$gr3uC&><#Nj>WTDMSU(Lnd!5LlHxTNm?7
zZms`1uy>IQlR4GaZbfR56GYWLM6|KnJ0!vTd+@T$)?;IQb&`g-M_CI?qz2pDh}uTK
zHQ%Ts?%=F~Q!w74z$!i7q^bGLio>3ksj@U2=`TCX{6Y*Nj{dgnMn0jD#3-H~dcuzB
z#5A^8T?OCFmdp?+`iC&E>dPH3o_``&9K+&t7MGPfga6qv!?U=)m*J<*afw&nF%5~7
zJpYsy$FoV;@B-&Fx;I&YTSa|+-F0Ln`1Sj!;MYeKgl44)9O-n2ge+LafjP$wQbt3n
z%F5cSeOd$2+8e~Pj&0T6aCa&Z_nMn=Fs744)3c<q;pJuq)kwujwM%Fh1}?46tpbZV
zb0&k!Rn?BFz&g`zbYdT!OvpNxySrJ7h?>S!-~U=IRz1nw5a9cgKYu}YnB_V7<l^^N
zU)Yx)>mFsbxl=AtX<@XE_a?aH?vbC<c`b88&*mTfK0o5Do)~YQ^}s@zLYCdQvAA8^
zJe;^H@@+X|<eYVyiJsm%PKBItNG7zDe1ew_g%{j3?U7T&dUUk3R<C=rp`qcYe!oJs
z!?6`p(j8v5$rfh^j|g9zx=&<_26lU=9#V3%Pbe)bizn0iEApY9wo~=Is5+l}dtEaV
z^zXDx?Z}v#R2bQ@M1I$RkdS5s^c}zQ7ak>@+JXB^9O4-jBZe<^baZ}JTnsgcJYDb%
zI-;2es;s&?zdli;j5JTzFfcqt(yhXMNW<OXHa1$e`&~2e#OR4R)IMWOvt?XXN1C>;
zk=dI-FD@_d=oZ0prn;r{&bu(b-(TsUT3K%pB<E^K)|r~EXO)_fN;kHngQNS3gSvTL
z`6=Ytbf2$tVbHT;EhTDr%;*i)B1Gv%t&%7*Azjy6FM2a`X85iVnK7rSv;;@EWN?_Y
z-rex=v~eMkn0NQ?-POI`Yf2pQFtX=HYH3{q3GTVaa#fQpJpBESLl1{0I9jnh@LTnc
zL*B&R&(yOC?neTXi}buJqdLjgii6G(X~E3KWx_59Z};*Wzwd=ubD>Hq<r?3(LE1c(
zv*jLkwwO+e^ek~`Z=~=iP<1N`3kww~r6&mOf={kZ=e38sXO6VZW*b?F=j4&7a|UUQ
zurng|*>T>KLhnjr^dIi_=(LsBh7~^i*8G&0J66TD*if}-<L;iA+V}4dw1>IWxgE>3
zG0$f_cT~rJIbTe5K9@_De`XN#VAcTjw%@P(c&#(ZqzUtKU0cpIL5ZLJ8WIz-IDEks
z$+oERoYza6gKh<nl(f%XtW-GoQn{po%Y8@-hNAAbF(>@GxdBJYlXSD2?mGXF8$MK=
zkbJ#}HQCZxdHjLNH>+hrnh4JNP)^Ip=G71s`JmKPi=38D3Z<pH+o@f}A7@SPO6C}r
zxu)sN%gzG=sWfn>F=wtkE885@NVwFt;pz*R^gNYqwc~D~1)-tg$Ow~(KdF|gx6H-9
zttRCasn3KFU2KaFH14KQ3|kl$^}%z=WNPHpMHLbzr8LfXIf;L^sAG!SH0wyaf3v79
zvA*AOqj3T@x#p6)`+h(psBg1ZSy#FpZom7H8*cE&UTw~NMql;OZhT5(Gba0sRfrWu
zXT#lH^2tjJV^~MdJmUCU%~1%G0huJ8E2p$YFWUd<D3M6aKb)Lh{#h?2D^2tQhCYAh
zy+v|o>&;yrO*d*y_&;3uio=bD>o{G7`ACq^w3E$<a{zSQ_IbOuv~2k?jcPX2XENvD
zV34%Wy$vBt>`A03()X8@5w~{idL}+UqLA{~wubb`vuR4-dnjga3wlQS*0~~kg3Y$T
z6530uxF0PQ6~2#lqlJ4>m>xdjrWO3-cWyVWvG8svO=!64k^i@@l#K55q2Oyd78970
zLt#hw<=Eq1k8Pb%SjtREkfnAnib^jwX<|lO<79~~YuMeAO_kh(xgAueI=5!%Y$Jh7
zy6jk|`1YIz`#H0F(}$o%7k0d`!?$)icX;WHbq#Q5DE{BjH7%V3DY|hpx|6KN{p;QM
ztxRP}5_f)tZk6%WM0c-+XY$@T_B(8N_H_Ob)Y-|vtHnW0I%(32YmU#z1~*xLVA`Y@
zn&D_oR_WGyfpRt2YbOD1F|W^4zIm2(MjGzt22ZNhdoX{71QqABHjytdn2#U6qvlR^
zpIGOX?CR|zk@TTK;ug2`s6mRH-Z=9DW0EGy@KV98taT;H1o%u+{F~;_oD+{o)t%EB
z)s+gMRe4XP4frv3k0OJlb_0~2xXd9-Jm#KxR<riAQ?e~Xtf_*Qxd_w^OWOQtVo+Xa
zcbzu{!-+Xo)TcMnjt6K%+4QmcbkiKU5Tkb(r9Wc)!-~_38)VEIXy86&G>&wy7tStk
zFr=YbwPxq|@}-zb@nsQ20#h_iGm)O-T={q{!eytlkwcwn#ZyyNE@$3XDIkrJ(A;2k
zt~R{AclB1Q(_C(JNOZ+#O<v9f7Ai0`qj+yjAyWNhk@Ubd`NF3aY3qtHwQfnf&-#e>
zJRZYfNWqh5Ij6rv9S$==jo}V8uHxhNudG<!#K#q-+<y1eiz3J2m%a;Q5b<lpj+4S=
zlxvH3UPJT_${Ag53iJ%y7Jp1RPJg<~-_>LH2H-NCU0tnhZEZD%MF7-RN=Q_Wno>GC
zz7RvFvu(B_QO<vk(AWBMdCiPcoLQ?!i9eYXuoeoYwdxAuq{c@rE;a|!^@y(Kp>%8O
zaz|UO<Oht3&Tro|2jVtH6piPIB!3Oq($`cjn1|lkw?^d2H|he0*;2dJ<^N!qx>p|8
zLib>|SfS|6S%<+w*ThA@#m9`S$Z4o{0+ziU{mNyaKvTF!f1cH1<>273ta1f{_=uo*
zGZ|Pe!}%8afT~yo&4=E+#q5LZd0*e1KRpm+k7S<hVfb&v_q*xh8vijXj=H+~D&V$B
zJC&WQ-3KBuV9GVA{|>TCxQT|b_T)%w7>9FM*0?ZUkrOH4C=-2TFn)Q#%le1g1q51h
zy)$mwB8WF2&m!k#O!?dWkPf6}l(Bma?S6mZRqeBE9fFv?X{Q>u<Y_!d_gcmBn&$1x
zs-F!pn5j2xI_5Va!bn6D@iLtrN2pG0r=8>Zr-iYJPF2HYix2D<(DiiA6g@pv8z(p8
zZ`E8hK%SeT_p)w3JhSlj3Y3ZbMqXp##+#de4TPotIGI@M^XuYtzH-18?9pq=p^@~N
zi#_l)+rI6);;%dXHf0;i6L*-Y{QiVL(|8cGOTH&bpc3p+3Y%&Ct&8h@4D+_qYO-32
z^sJI}L`$M)-A<|92=Ygi1hxF>R|;k<ZI5ykB8DJ}{qf6payuqGT@!b<>n(vOx~vmc
zK#&Mmu?X<sO~DkAbvFgK{$DRdNkjmV1%B?D2zrI4_5snF(*>$@y~Qv@(hta9F#vZf
zE5R)hxx41X1aFCqG7;8$o{_OxR#vt|zRV#{5)Y5+$wVX3F6Yi29KT$}Eh#QubIA;B
z*?cEBm39Md`t5W<($(GjNSh-dmppVHIXF5#c8di2EFAI@evH5R(x;mfaq#8Kmrrio
zxbd<cp45}J6ylp@d=dNgZcPIuq2ApagE?$aT}@7t20_Io)mOYn-$xu?UUX}mGPXGa
z9#a-L`@cH*f8SmI(@zWtBJhY&v2E=g5CHNS+v4gz7QT2{FNn^%zC}HAa~0d-Z)f&t
zUY;B{0@(i^UT|e&hLuodK85|Q)H}Yu|K4-(P}0hbMJ`Ze6z>m)sNX;I%pvcljTDGd
zKi7vZd-df!E1Wkr_L4by@ZdoM0|O7##S(`&)#c#hr-ERAeB$z|s3_yIQ8y3j#>0!Q
z;OH|h0Z6eL<L#H~+S=MHK|)lJRe)M^*?lo6%XD+dlW{Er#2~t5dbYK-wVJG<rUo^o
zfFT&PZJjLQO&JXtuE+-!CRr4M`B^&yBYZ#e_w*+q(V#0Y9T623{VW6-JmB<s(nw=i
z@$3nF>oFrayP}St0h$Q1>ncS2c+G!a#BU<1xej?A?(SQLwP1GHZ)T8_YHPR1-r57-
zR#$DfB;*R{!!;4UX`Xm5FC@?x?~Wft4*CWJWM*V!O!e8~e><%Wm(SHo)kqBw539<>
zl#vn=65gmU1>d{vBu`eoytA$>%!J^Or!z(aFbo=vvGkpH0S?DmhQhCki!;v_<-U9O
zZt^#AFyLvc$5y!*8DlKDe`}A0@~->s=#wX#-n@BZ>Ldf>Ip{6<p8n<hbu=Vi&D-_z
zNPT}HaHB1<&V4lUatNR%VPIlI&4u0#aKl-NKHY<YUJ>>GD2zg{*n^UgH-V1%yww9}
zM}jvcgwLI>RTt!!y6J$J3%<44nb+?;^ybEL2WRJY9<SP?<l)rxbi8E_`_p#9U|ZfP
zj}iyhL{P3ZTNHg4k3r(LT@#f6PL)eG;du4<y*quBN#a)@%YAU^(xvQ^n_z<Te=O9q
zw}SghOLYqRRk9Y`gYWS2fQ~urhah{9D5{kH5+k+G#$hwosQIydJ^e6aO-Ci+C(QA_
z4<U@aGS*n{KC|Xi7~+->zQv$zl6VZ#%5no(G_#VqNK%~}44&QEdZLSjbG~)!R><r`
z)@U~opO=^SvL3p6((?PSM4Ezn-z+K5$LmtVj_Y3h-sDLjW>*F=pLcj|yMJi#-5qJy
z#BVn|Apl9nU{Lw;Z|5kbawzNzQ^j=7G!t-|SuAj1rTW4*J63|vEKP^jjJ35T<U5lV
zXs8OVI*6GHa|hk&n_vQ;ZyZVyPEcwOWV7fkQz_LFaz^*oG#KLaudhIAhF#~QU699H
z0u&a49P~C{=s_kK#did*$Ry7&<slzw(z(dYYf<@Uj_UQkk5}u4zdc>xb6ZfacHLhV
z(eGy;Z@!;d&AIq33vT|Lcc^R0V+-Z@LrQPdVahJ~`-jAO0A2y97f7d2Q&yw+cStjL
zB%ce0(10O8i)DwP)551ZA~#<tIJwX(HAUb8hr@vay4V=}yY%)~?TARz-=q7kqeUQ~
zj(2R3&;5Ra454znsJX2kDKju$&Fuam)D_29K3XyFPwxy)%{i%-Ev?~n-MZD~v(n06
z^ZlkBq2XM;Vxd2kOXn^Wqok@si)$oIX7)H^XWgH;GpQs3ia&{xr&Qoup9@sq%1S)Y
zeGo?qh&qrfYxe(G<CAXP)Yupei>6;M8f0S7roAW$j0!)r9e4!%0x~nq;LNF@az}ui
z2}&1&g6u#_0>XJFG#uviEb{9Sg@&`oXTzLALqp*)P}~#TK?e5r@=}efgY_Td9S&Vy
z@Ju@+7D6b{?H%%zDWz82aCOBch|82F_MjSaGySC6wiL_51}m=vB7l1f)fiu?S&M`a
z&CqqE&3?=DEU+_inM9(oe)KNQdFcQ2!NF(Tr1GUgNR8U<w?`C?Hl5!MT0C%i5%r?K
zwT#m@Af|mg6LTUaXRAO7sY%N@mb;qocM)_03Wc(dz-0t2_PKDOZgAg62>u33A@T)e
zV=GAI0I$6$7p)f6YDd5YmG>A+A)ve4k}XPWsXQL`$(SgtExcmj;;~!uSZhGiz>cgM
z6wM&UP}#tBnCjVqZ*4RyPIDb-%7GHKDzGM?5CvTq?)e~(djXM^(Y-j>6+|^2pZpLb
z?}@=!6`ZsH5onXm<yEp`NJK~NjO=}6hMtoqlKQi}!)qvr#6c_luJF%mdicRTV4d81
z<{#aQCw|t&<$thv1-MMp<o?o#G#sy7&Lj-P{H%j9a)I38o}M04RjAuF9i{`$uhG?E
z_K>s5utxwO3A%Y(dwWSm#q*L9?eN+D#JL6xcJ-xXaNjPC*meXvKmMWO(#XT9sj1E@
zO}-nuh@rZ}%@u#4Aj79ZV^Cwz$6xP$)UYQ=%u41jBd>TOj~cSEGrhO+X(7OY#c_GI
zcVza#xQ(hTWL?aM6Wiih{fW*Z4m=*0(@CZ3>J-4bdu*q4QMf#KGFwnOz)$C4=c7Fz
zEl^w34nv;hGlnP_Qfo_&m0(pe+6lMN(6(wm-~O&00VxT23-uW&6Wk5@kl)6Gb_6H!
z$|`k7c=JNl_$XClIeytu7t_5$z(*p=f149T8cE9kC|V;bA$LA7@I=|gE)VLCFg`E5
zt0Hjw_cw6)tj*6wF1a?xAm9cmf>n`f6#Nukk9SnRh28*TTL>&5=A$dK@pgovp;o-Z
zsZ%XEY1nP#NjH1jZ&`whmibX0$`qB)n;v2EhX_1#sI_43CR;&MtyGjp2+kzz57f06
zd3hjnU0(o9UU%G(VWS@^vZ5wtAv0s?Cj5DSv)kVv3v4b@j1?M#L~4~l%{Q@8SD9Uv
zIgaW%0#$FH0$e{JU5A`hUGq~f8YC#rd6$BdHad@5*^KtU*@^yp+u^NA==&=oVC~K2
z;O11geJI`;&P8f%?LOk{z@ykq00K_FFwQRZSfQ64UZo_=?hUb-A;fNMV_?^XM1j@1
zit7=A6Y)7z`F5w!Bm4;s{!=Ry*Cx7eaGJy?f|nEpJQARv%h`8Te`IgvtGXD;VNlaO
zuNTGFyN?VHUq4}+2&sD6DE}7<cv%mw5`y}mhm<<(0JZ6GENHM&b>gO?!oLa!Yt0zM
zpUvYi8ILg4U!Wmn#K-?;96X^kTnN?<{I0~89^-PR*@Is)@@W^T)p^wU0o~p;0(%cC
z1}}a|#m$~Xqs7x#2^!$)I^-MXVPcL|+D&T!EH(3yKexJk@MS^gI1d6<tPN<p4QKb|
zL*I5WF`W?Wsbp3|)K$onKZ8>zW~;h!eC{uRy!?KONOpUiij`u}YO%8Wi~g!&Gm~VL
z{M=*+F`vr~2YuKgS5ffd-y)Hjj&o@j&YU?js1;{mZLglVMV3$Qy>@MpppYh8Sh@)a
z>QTu^B&p$S3AF!!^&fnd3rC@#8($vV11=sLvxZGNNTJpC(t_%4`a}6%%bT=-urnW(
z)_^j<84H4_8}uRMs8l4fF1vB_9Z*4GOLdbitToNx`+7?@mcj-qR9KME-L;ah)c;jv
zU$uH1$uI&t3$%D*Pn2NiMzE1t;&9~MY2O|{s8+W>^o<a*va+_?g4Z<VcI}erQJt8b
zZ<5i6WTQSZ|I}Xy-F%20LiH5l$O<_+IB12)ZGx2~K8ZBaxAwT^_!>e6<PivXC&J-a
zXh8w<q(3-EmK6PBW#hkr$0db<TBxjw%-?(k8cKruP6-7_T2KGr#hTCDv~us>74M%6
zHiH))%R08Sd>fcs_psg*$@vEmxtS&mZC6w3KW0b#Uv~Bly)YOMm487CQz?X>2~a*L
zNOu9sD@P+UD*230r6UUjH%U+UNhYZO0S2;{)<Z@dANTmS<a*wh^FZ+S)XV+ZTLoc%
zw&2H|Jh2`Cu^T@uEglBA@MHWk$MeF43vXIsc}6+^B1!Leq>=Hqc=7q*ZSk>vdV+0f
zG}>EU$PqTDvyQF5Z2FQaRg;yqMf=f_9qQ`QlUDHF)y*CyYhpg!l_Sw6*7zvLEv*`c
zzRC#&lmHDJBDy86Zgssn&^K3CS2b+_N7voeV~|wXus)t=yF5%+Gfj{K05TmHsnTZ-
zNa9x3c<JChN;dfIhvDOZh4}ZZhP4hP{m}wPfB)|Y6#vtV`@bfe{N1Yez|it+PvqHd
zF$49S1b8IL?!9OY6n?0J?@sUZzIag(^TXh|q-S=u>Iz=Jz*X!5_J-BqvJ6z1Rj(!B
z`nwtA+X8~xVg>|wpVcDyzvV&J0a#DIS|-R@t>h`GJ(N<|<FdiG&Z$>IpsY-0-ERNy
zxgfE^H}-*f7Brm9zRHNWy3EYJj4O#6&bE2BjL<dznhdi9LA09wnG9kskPO9u;9EO8
z9XQI8&Y&Gs&cAtcb?w@<jc7phdy=A3Mc$WsyTf)uA1JE<1`S6xA(BlbE(bT*5W}~{
zVCT!`%5eM`=pm9#j&edPJ{79v<>$kaumd5~77tqJb@W=~bCk3OuKe~OCFj<hJccVx
z^xL%d^nT=O(qE+`X&BTydLwoI8n2vPU-IQ&M86g5qJL=t5k<3emFFth`%AXZn>EWn
zL^0Z=PlR*3itcS&B)J$Wj%KVtTbyW}9txWKh#Lg~kZEANLcsx|8F4=hG)MA7l*D|w
z3f+G497vEpA^U~aEnRr-FV!!CE|C2WOw487IfGtV!x9Len}eAizk@~E=L-{)qM0IC
zrU&1QEp$Bw1h{o|+oexeOn#xn1}H9`i(!VbbS8y_r4AsAk*ML?p!~s6Ok7BK0f&0#
zAF`i<r0)T~iX@=JaF7ML)BL$KHkCK68=J|WB;}8%?TJ{w7ao7yk#&b}hMFJ6-9(E-
zNcjqM&;OnyMxXJn(uj+JQx%If(z1^`K1&YaHt#0}X{Xz%zkLG)TT)n9!6kw~3e7ge
zO-jts+KzI{w~yZuOml85W&0*u5kzcJPQ(@cX?74ffN{se4oCxPH?9p7**S>NrW0_(
zav{yDf>%tJ;9B|N6a0mL+(~r$w#D_ppX#g~N$w|Vy3UT;@&TsySU5UK@J_4hIVXry
zSz4csY~v853_3hF@vo3s<4Wf{S<4n7*Kco*kouAKoFV&C|HQ4LxxUHQy-YQLG6Rc*
zG!^wu+%mh&ZGKmvE{Y^cEUx)uUB8X2<kwKU)fKBq&-TU<-v)$5zJ{+YYA{_bVHZ=M
z=PTpJKHW$f4P%umrHEbO^?{b#5Eda&Sbtre00lZ=y3CAN{i=yQF9ZJIP?U6~g(9xJ
zTCS%!e>@O(p7(v)pe89=)XTE9SIs1hcKaBH^cp8Dgd9>t4&1~&ABB8<mzgj5G1nlf
z781vTEsWfIyS3-gxK%4rifPhNp9>m+rsYPV2Y`ZL?-`o3l794RhaHcf1Ad4*8-}PQ
z*27^T*uDlD2>hqSVaL~pdMBVJDc8*)uK=Q(XDOP2lnC(R=WST<%JqYu2)ezsw>Mp_
z{#c2_iHtAj(?p;sCE0-#3{V7sBnZM#YsfUH_@e4KkeO`pt?VUWh-Z|^lCaB{RoerB
zxx(YhabMunSAQ^<PaTC8WitCSZmXRDjX3SZSa&A7&W2$J?`~WWoInKdgmYi)5F)n-
z7?^q1Q^6kHbaZrdKtRA(>~{~4kpI%5#(4s`s`!^uJPn29vzykH1Q~V+>BfyEJJi%3
zS3?WhbA5fd_!Dmmv^<!_z04`c2yQe_fb;K^8wu;({-rJEy?F8LSgxixaE+u;|6uW+
z;=j1lzv_waXFL&dEh@C}dJ~MW68BwacAFk-o*3bHj<I0d8ow=IlnrdI*cr)5R7P3*
ztp2E96*hYe#bw2*qbmTc>$J5`sU-8pKgi!D?$K85+@|-p&nvUicv^4>Z@wJM#yivB
zEZIf4*(cVc51PYX5|2$#sfd2}ZvUlKLSJlO|NV31d7g#w1Z~%@Dd;Eb7Br596a+5}
zHl`TPGkcNq_@00L)csHfn1aAA#lN5xaLtNZqlis#<gEf2-nRHpP5qATCkjqJ@TNfe
zBwHK+w3_qr53x>L0DG~Zyg<sq*}T}gpS#VYq81MNshJt*={12aBQho3!f8|Wb`POW
zP;$SVAMEMbBO@cTh>vaIv<iME4VAzV0yt7vM1CkDk*>OCnG<Z8UiG>PM@V?pw)kK;
zi58o+zfsjYzRu0aEL8~(ml2!@R@F6s)SWFtQ-h|`oU*6$j&^=xRlD0KKz759SQVTW
z>nJt$vH!<QellsdO0LkuJs$&>js<Wvcs!L+fsXN}T*OuE-rrOjn+a{Xl@SP^c=IiN
zChKkSCxG@8%jdLQ0F8My_s~<!-DcZW;c)aXfGBZGHl6O~K#ILu7H>)#`YBj*a&wcG
z%2Rp1uxK6sDXAA*$@D-+kLrq425ep#MBsw1eyBOpMCnV}feYKy>cT<wnwtKF_rnu+
z5uMFmQ$C&ZBTf8vXMa`;J-xoE?EM{2WOJkZgAiFcIm@zp>%aWM0<z96z^?UKtNV?H
zzDJIZjMVu)uD-TBeEa27W~?3bnP>Fw?mwkEs{5ebhpWY4Yiu+{;SiFIxhpN?b3i;~
zb5j1~RhPdo!<!x3lz26Mop*O}vE<!*J8xG9hJF1CiV+4QYH+%xt{D4Uao(U^Gmpn>
zx)-*gWv=a&&oaFqYV`)=^ntOlrkR!{*8-e2|DbNs8GM)@Tn2=N^-yuHO~cZ2dt`Ng
zk=kM4K_IFz4ab@ZSYHRQ6AuUh1cb#73Dw$8P;wbn37IWJLqiYX0Z65il9JVAiOOtG
zrx7?`85(-3)EkgAOEzrqlwdyeaNu~1iT?VL54^3PuzTOJu-fLQlDi3EATfNRO~MCK
z7HCj7fw7B-xho^2DI6mDxY1B5l(57XZ3hQl-eZ60=j9c4cyZ1xh-#g$N<kGIX|Vzs
z0|pt<d4dU)$;z`BJr!yC)N?G)7GJ1_56{`o9g@x{jj!XH8>|e3b8_4*W<c0jd0T%;
z!Pg^!Sea+V3dFsVk}S{m{1SqcJ;;eAp6$gYC0ppP;%I*Qf9`1PQ>lA=FrzpgRX}jJ
z2bzdBv|IfuAK%~KpHGAONpwmyr4H@(fO(a{4lajv3WyY(ft2{j!e4~P=-f~3=XXCK
z(*!HYLb+w-6;@u|pI~}He*Cw#{{D_(f1%zaK&J@u^y!!Lf`hXlKEgx|jEq2Jhcl?=
zf?Wqu<jiU~@VvZI4Ef{g0(r^n?|1L_)7OvAfMZUeS%HiM8g)cL?;1^H-pnAy7D_Sn
zsLf^}tik(&Gq*36A@2pyNEWpCO3UfB%a*_?Aj$2*5#s=@_S`y_`}Mb(`8`Qw+QOZ8
z1b1$@dw|0tD0RK^-3Ld98Tv?4;CI0WrOUc#Kx6+OcKed}5*iNd+om=ru0q{J%=TNF
K=Kt<;^?v~dgpeKp

literal 0
HcmV?d00001

diff --git a/scripts/process_apply_test_data.py b/scripts/process_apply_test_data.py
new file mode 100755
index 000000000..c56f1f10d
--- /dev/null
+++ b/scripts/process_apply_test_data.py
@@ -0,0 +1,32 @@
+#!/opt/homebrew/bin/python3
+import sys
+import matplotlib.pyplot as plt
+import numpy as np
+plt.rcParams.update({
+    "text.usetex": True,
+    "font.family": "sans-serif",
+    "font.sans-serif": "Helvetica",
+})
+data_file_name = sys.argv[1]
+
+file = open(data_file_name, "r")
+
+lines = file.readlines()
+number_of_trials = int(lines[0])
+sizes = []
+data_sets = []
+
+for line in lines[1:]:
+    curr = line.strip().split(' ')
+    curr_set = np.array([float(n)/number_of_trials for n in curr[1:]])
+    sizes.append(int(curr[0]))
+    data_sets.append(curr_set)
+
+data_means = np.array([np.mean(np.log2(data_set)) for data_set in data_sets])
+data_stds = np.array([np.std(np.log2(data_set)) for data_set in data_sets])
+plt.errorbar(np.log2(sizes), data_means, data_stds, marker='o', linestyle='')
+plt.title(r"Average apply() run time for Two Qubit Gate on $2^n$ qubits"
+          "(10000 trials)")
+plt.xlabel(r"$n$")
+plt.ylabel(r"$\log_2(t)$ (ms)")
+plt.show()
diff --git a/scripts/run_single_qubit_gate_apply_test.py b/scripts/run_single_qubit_gate_apply_test.py
new file mode 100755
index 000000000..480df2909
--- /dev/null
+++ b/scripts/run_single_qubit_gate_apply_test.py
@@ -0,0 +1,26 @@
+#!/opt/homebrew/bin/python3
+
+import sys
+import subprocess
+
+experiment_name = sys.argv[1]
+number_of_trials = int(sys.argv[2])
+runs_per_trial = int(sys.argv[3])
+max_size_exp = int(sys.argv[4])
+
+f = open(experiment_name, 'w')
+f.write(str(runs_per_trial))
+f.write("\n")
+
+sizes = [2**n for n in range(max_size_exp+1)]
+
+for size in sizes:
+    f.write(str(size))
+    result_string = " "
+    for trial in range(number_of_trials):
+        result = subprocess.run(['./build/single_qubit_gate_apply_test',
+                                 str(runs_per_trial),
+                                 str(size)], capture_output=True, text=True)
+        result_string += result.stdout.strip() + " "
+    f.write(result_string)
+    f.write("\n")
diff --git a/scripts/run_two_qubit_gate_apply_test.py b/scripts/run_two_qubit_gate_apply_test.py
new file mode 100755
index 000000000..eb9d6d8fc
--- /dev/null
+++ b/scripts/run_two_qubit_gate_apply_test.py
@@ -0,0 +1,26 @@
+#!/opt/homebrew/bin/python3
+
+import sys
+import subprocess
+
+experiment_name = sys.argv[1]
+number_of_trials = int(sys.argv[2])
+runs_per_trial = int(sys.argv[3])
+max_size_exp = int(sys.argv[4])
+
+f = open(experiment_name, 'w')
+f.write(str(runs_per_trial))
+f.write("\n")
+
+sizes = [2**n for n in range(max_size_exp+1)]
+
+for size in sizes:
+    f.write(str(size))
+    result_string = " "
+    for trial in range(number_of_trials):
+        result = subprocess.run(['./build/two_qubit_gate_apply_test',
+                                 str(runs_per_trial),
+                                 str(size)], capture_output=True, text=True)
+        result_string += result.stdout.strip() + " "
+    f.write(result_string)
+    f.write("\n")