From 54468c2db733877021ad997eb2fa30b2bd7c3779 Mon Sep 17 00:00:00 2001
From: Yorukot <ben.hrc57@gmail.com>
Date: Wed, 19 Feb 2025 18:00:28 +0800
Subject: [PATCH] feat: update hcp testing

---
 matrix multiplication/CUDA/matrix_mul.cu |  91 +++++++++++++++++++++++
 matrix multiplication/OpenMP/main.c      |  50 +++++++++++++
 matrix multiplication/OpenMP/makefile    |  12 +++
 matrix multiplication/OpenMP/matrix_mul  | Bin 0 -> 17784 bytes
 matrix multiplication/generate_matrix    | Bin 0 -> 17520 bytes
 matrix multiplication/generate_matrix.c  |  37 +++++++++
 matrix multiplication/makefile           |  12 +++
 7 files changed, 202 insertions(+)
 create mode 100644 matrix multiplication/CUDA/matrix_mul.cu
 create mode 100644 matrix multiplication/OpenMP/main.c
 create mode 100644 matrix multiplication/OpenMP/makefile
 create mode 100755 matrix multiplication/OpenMP/matrix_mul
 create mode 100755 matrix multiplication/generate_matrix
 create mode 100644 matrix multiplication/generate_matrix.c
 create mode 100644 matrix multiplication/makefile

diff --git a/matrix multiplication/CUDA/matrix_mul.cu b/matrix multiplication/CUDA/matrix_mul.cu
new file mode 100644
index 0000000..6422da2
--- /dev/null
+++ b/matrix multiplication/CUDA/matrix_mul.cu	
@@ -0,0 +1,91 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <cuda_runtime.h>
+
+#define N 4096
+#define FILE_A "matrix_A.bin"
+#define FILE_B "matrix_B.bin"
+
+// 主機端讀取矩陣檔案
+void load_matrix(const char *filename, double *matrix) {
+    FILE *file = fopen(filename, "rb");
+    if (!file) {
+        perror("無法讀取檔案");
+        exit(EXIT_FAILURE);
+    }
+    fread(matrix, sizeof(double), N * N, file);
+    fclose(file);
+}
+
+// CUDA 核函數：每個執行緒計算 C 中一個元素
+__global__ void matrixMultiply(const double *A, const double *B, double *C, int n) {
+    int row = blockIdx.y * blockDim.y + threadIdx.y; // 計算矩陣行索引
+    int col = blockIdx.x * blockDim.x + threadIdx.x; // 計算矩陣列索引
+
+    if (row < n && col < n) {
+        double sum = 0;
+        for (int k = 0; k < n; k++) {
+            sum += A[row * n + k] * B[k * n + col];
+        }
+        C[row * n + col] = sum;
+    }
+}
+
+int main() {
+    double *h_A = (double *)malloc(N * N * sizeof(double));
+    double *h_B = (double *)malloc(N * N * sizeof(double));
+    double *h_C = (double *)malloc(N * N * sizeof(double));
+
+    if (!h_A || !h_B || !h_C) {
+        perror("記憶體配置失敗");
+        exit(EXIT_FAILURE);
+    }
+
+    // 從檔案載入矩陣
+    load_matrix(FILE_A, h_A);
+    load_matrix(FILE_B, h_B);
+
+    // 配置裝置記憶體
+    double *d_A, *d_B, *d_C;
+    cudaMalloc((void**)&d_A, N * N * sizeof(double));
+    cudaMalloc((void**)&d_B, N * N * sizeof(double));
+    cudaMalloc((void**)&d_C, N * N * sizeof(double));
+
+    // 將主機資料複製到裝置
+    cudaMemcpy(d_A, h_A, N * N * sizeof(double), cudaMemcpyHostToDevice);
+    cudaMemcpy(d_B, h_B, N * N * sizeof(double), cudaMemcpyHostToDevice);
+
+    // 定義區塊與網格大小 (這裡選擇 16x16 的區塊)
+    dim3 threadsPerBlock(16, 16);
+    dim3 blocksPerGrid((N + threadsPerBlock.x - 1) / threadsPerBlock.x,
+                       (N + threadsPerBlock.y - 1) / threadsPerBlock.y);
+
+    // 使用 cudaEvent 計時
+    cudaEvent_t start, stop;
+    cudaEventCreate(&start);
+    cudaEventCreate(&stop);
+    cudaEventRecord(start);
+
+    // 執行 CUDA 核函數
+    matrixMultiply<<<blocksPerGrid, threadsPerBlock>>>(d_A, d_B, d_C, N);
+    cudaDeviceSynchronize();
+
+    cudaEventRecord(stop);
+    cudaEventSynchronize(stop);
+    float milliseconds = 0;
+    cudaEventElapsedTime(&milliseconds, start, stop);
+    printf("矩陣乘法完成，花費時間: %f 秒\n", milliseconds / 1000.0);
+
+    // 將結果從裝置複製回主機
+    cudaMemcpy(h_C, d_C, N * N * sizeof(double), cudaMemcpyDeviceToHost);
+
+    // 釋放記憶體
+    cudaFree(d_A);
+    cudaFree(d_B);
+    cudaFree(d_C);
+    free(h_A);
+    free(h_B);
+    free(h_C);
+
+    return 0;
+}
\ No newline at end of file
diff --git a/matrix multiplication/OpenMP/main.c b/matrix multiplication/OpenMP/main.c
new file mode 100644
index 0000000..4d1f517
--- /dev/null
+++ b/matrix multiplication/OpenMP/main.c	
@@ -0,0 +1,50 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <omp.h>
+
+#define N 4096
+#define FILE_A "matrix_A.bin"
+#define FILE_B "matrix_B.bin"
+
+void load_matrix(const char *filename, double *matrix) {
+    FILE *file = fopen(filename, "rb");
+    if (!file) {
+        perror("無法讀取檔案");
+        exit(EXIT_FAILURE);
+    }
+    fread(matrix, sizeof(double), N * N, file);
+    fclose(file);
+}
+
+int main() {
+    double *A = (double *)malloc(N * N * sizeof(double));
+    double *B = (double *)malloc(N * N * sizeof(double));
+    double *C = (double *)calloc(N * N, sizeof(double));
+
+    if (!A || !B || !C) {
+        perror("記憶體配置失敗");
+        exit(EXIT_FAILURE);
+    }
+
+    load_matrix(FILE_A, A);
+    load_matrix(FILE_B, B);
+
+    double start = omp_get_wtime();
+
+    #pragma omp parallel for
+    for (int i = 0; i < N; i++) {
+        for (int j = 0; j < N; j++) {
+            for (int k = 0; k < N; k++) {
+                C[i * N + j] += A[i * N + k] * B[k * N + j];
+            }
+        }
+    }
+
+    double end = omp_get_wtime();
+    printf("矩陣乘法完成，花費時間: %f 秒\n", end - start);
+
+    free(A);
+    free(B);
+    free(C);
+    return 0;
+}
diff --git a/matrix multiplication/OpenMP/makefile b/matrix multiplication/OpenMP/makefile
new file mode 100644
index 0000000..d16dec4
--- /dev/null
+++ b/matrix multiplication/OpenMP/makefile	
@@ -0,0 +1,12 @@
+CC = gcc
+CFLAGS = -fopenmp -O2
+TARGET = matrix_mul
+SRC = main.c
+
+all: $(TARGET)
+
+$(TARGET): $(SRC)
+	$(CC) $(CFLAGS) -o $(TARGET) $(SRC)
+
+clean:
+	rm -f $(TARGET)
\ No newline at end of file
diff --git a/matrix multiplication/OpenMP/matrix_mul b/matrix multiplication/OpenMP/matrix_mul
new file mode 100755
index 0000000000000000000000000000000000000000..602ea59fa062782c2ada234691b6e79dcb627019
GIT binary patch
literal 17784
zcmeHP3ve6NxjvE;A|Oc0qs1j1Sg4bR&b1;rj&Tfu^?Thqi5mwB7t)Pv$+j3-awQoY
z=FYvQb}}f$(-df+O!}bh^hHYx^mfw6B{reqF`YJ{(>9$>i+e9gq}yIdQ(pH`bie=X
z`Q_C{$}qh%ccvZvp56ca-~Tx0{Cjry>}rqysk67s?Q#hwx42%A<X^5!$l41ru}O*w
zYZjLXU1;J8aRHD;*mS8OZk3pya9axI9%WYpPJXpikX$d}@0F4e%(;cc$!~#lRkXOS
zDqt@BmYVz)Q&Gvnuag8>m+VRvlsVfK>(iu}5SnGr{591u=5{~IO|IVZ%Wbm!m`|&`
z+@m=Akw|`TDZjUrAM<%tA9KzpHJ`R-mCxoH9TLeFxAqIi$>}$t{3et=^9~9EVor4s
z^`Fc{{*wCol;5-_!-#2>Pv*TAESPhBABG?0<((2A=54CJO2?sA<)67quUHq44Ky{b
zi%06?vBda9{X}z9eN$r~oeVU{CQ!SA*l0|3Z@E#pC&esb?C%mTDPIVT$gO17lkG}u
zly>s}#o)xl4{!bU-Wz{%)sFGgYYt4U>Ham@P(R3rY>3F7hYZQ2<BjBZN`DW-<z<ZT
ztbm&p@ck8V%2K)bdIkK!3V5~xPI)O8pK_sLbJ7@qa`v}Z!22uU(F%B71zecs@Mtn&
zrZeGG#xzB8bj%!%X3WI+sF@i_MZ=M_=-#@y&m0S<!tr=CUbMI~6B~^dCCirF#8@<y
zN~Xk6DjKC-I3l7Gv5Xi>jztq<6hV@Ms1@UPP}MpbjwQsPr5GBFC(}_emWm}ZLjp#_
zXb}zT5D9}+YO)&C7zhe#SV-{(@b2E8wsx~2u)eWaYzS;X*D)q*u(@#!S}gfuj7pN}
zkIq{c^q&Q%<ECRf_P#Y3-WP~}W~`3n#o{TJSChO*9ED6o(H;dSf1>Yv@&Llth@YTL
z*!X2JoN8|ID^h|^+@MQNb=%LEc^iI_v`{-;mxU$R(->m%*l=};lq$asM?;HLYs2OJ
zMrkw~PHkdZYs0B-rl1YyYc=VbZ8%?B2-j`+#RMoaZ1_?J!TW6ZB{uwK8}7B?J8ihn
zhF6Q#3{*2v%|JB+)eKZK@Yk1tKlv{IFJtNt9^>A{$JYsA?4HfI3&)J9XFNw_X9^o$
z09;u8W9)rb=qQmqLc08X0kryQl4;76KPKhxlk6sWR?1J1OjD=)x261bl4;76|GbpH
zKr&60^82OyFv&DU%1=u9K9Xr_lpmAwUXp1_l;0`kkCRMOp?sf|?;@F|c6nXODUxYQ
zln+Yzc9Ll-l-D31UZLqCca%!hhun{;w1%)Jir68Pd{<0EEkx_~!`ORn>_^gqxdhpe
zY-E3w3qyE{8b*lA$R5e917hs{kDlx^z6WNF-M{ud@U-WcZ}&mKn|HtLo4%j)*;fu<
zNA<NKXGZq7hkHoa03kQ~Hq<kjI)CW-(2Iw~3SIPOBXydQy{)d+$j0mZXMM{$>O8r%
z2$Q=X#%F!EK|XvlA{>&v`7A}qzMAWj{H5IW;Jw+?xh9gQo(YC#Z{5*)+mWFFcW*z9
zhY(tP{G&$pUyZ4k&-88Y49$j~H=cbqxD1E?U0=kDiPsh!$lmUMXiV3Iv5EGX6^D=r
zIixN;Oh93NX#bI+1z#FZPf7c`Uy&?(#+ZBM24n868pCzeIR18KSuw&A6~XCG^-pek
zgx;2nuej03p3L2aY#d$u2_UW=N6x8_?0>&^p$#~Y+5e^!Aa~!N3k4j-TfssvX6Fvk
zW5t_)-gw-%?=uMQ`_^-=7iO0FH@@JT*+ad$_d_0I@Ws~bT<^X8{@#06^|+23*`FD6
z|F*Pe?srR#nXkdf(>t^1H{dv~`B33QGmY-et`lc{ey=feH_*3wvrj{EpY`1$jKyLS
zc2iG#&iZbC4hPn2Y@8ea-PE%#W9q2Kn40t4U7!|?&l-0hp+e8bUyXmB&eL0W9GRMR
zg<j0xf{yP;3SNK8j{!gf)RX(o|KLFHp7l+yM7J;!jO^U^QL{^qQyQYzkpZF@*_UqJ
zk-r|25R@Z%O#p6P1>jPpY6hwqsAiy=focY-8K`EUnt^Htsu}pJ&j9_dVyD!AI6eKy
zsV6^`du;OLCk~wY`U9sPxku!_`p~I8-_1YzKz{ek>3@0b<d+^l^{EH(TT><#n=o4g
z1NePO32jmk;`HafmjCn@U;f@h2zBzYeW&i(|LPC-<@P?Fd*TPDKK=3hf%{vu)kE6p
zuiS?xb=RtzYYjAlo}CZAS}5!T|A#Y$f(OsA`@#FbzX7g+@B3q+5P|-O;1gtbwov#E
zc<-MIg?VuLO_j)X*EZpr@VizmUE<jfKl+V`z9rD$w}uQlmiP$^V!I9)|5j5+Dr_FK
zV;)?1{axPL5Be^?(=#b<xa`^wtgpL@;I-J|*k<ATQvzDO{ypwiZ|zh~w^wt2?;>w4
zkXEm!^<qIWAHo(vTic4ZJ!ISA)o%CJc1TOv_F3rTs9W~Kt4-DHaf`nLb^_Zp@SlT=
zHm`r58_n8N(+ZPC9o|0ooL6g?RknFOwJ@RFH5cRbO)c6}v#<0>Q=lqUGf>SyH3QWQ
zR5MV`Ks5u^3{*2v&A{Jc2F@(kWw&M($K!z{`X-NvW%|aBXsN~U;H1kEe<QYXr6s>d
z>G^x+RZ6C3IHJoeCeCcqWr@GzX8r4L7m}oyt+Dj{{IzqjCG$6O^GeU(>YXGf5P#!F
z&pJf!vRHA?-|FRV&?UytKFzA5{G7<&++C_XJxeTm#!p;e$=q&wUL@jg@3?;apep+_
zqU`ZooRqm<Jibb~OX;mB|9^Y_e%x7OlRCa#ir=L8km7eJzDMx`iho)0?<#&&@n0x@
zO7S-or+T~F+gr4?_>%c*Ez}qY1wvXwu%R({O)#Xb-4=~#MmVDs_19eOux?=MjiH7T
zYbfg1UoBd1YH1t`HEg&>Y`du?HJ(Vs62sc?;GnZ6(RNcyIFU%=l|@-~_eWZS6W2tV
zf*YC|v_N28dRICU9ZRpn_uuQn1F`x<G7}99ivGTq^`gJMB_!IrTAI&e*wxz-Jde7&
z6(Rb1S{g-vXG>7T`mrl%LM#!>pgviPD$9_nTw;rFj5D#p`cxEO09$PgY7R`sabPl@
z1C#G!{B+Q8zI3<Nxz>k6#Wf*K>Ka1lgwv)|>6n+Z;cLpf)z~O!!`GSjs_RWTn?6<7
z8|P^c>zJ2|vs~4|{4cfP{;WDrzqmhFo~K{j4(4<iDiK}I7KuONSNl?;KDeMcsIH%y
zQZENDXHVC2x4;*__{xOnr*dz-=aM)+{vhJ(=0=nWU+KWxftQQFMcMN)gD9yFdx|ru
zj-i6Zpu+in#rQ51T;e?<uImy~mvgaO;wy{y8>z(Z8`vm5-~Sjtgn~<~5Pbh){Ff+{
zYyZ5oUoOsEsiRa%<?LU9i3@in=XLi@>AI}A?u8>=gpRnxrGoEQjBCK@m^mp_0sjMV
zP1K4Jl^^!oQo(*Fa6jxhpKL!|!TzoaIK5*n*Um=?cl!jN&m8~p3ii(cU*=vW?$LBM
zmHd|#?Ef7&wcpv_H!Ik$zzw-v{PzRbN@}!(CTZ^%9!GoH2}jMm4s4O~-3EL)#@jJ<
z9<?d_qY^*2f{#(0GX1>@xF6%3zKgQcUW()P3SFHqJa0N9`|B-UqvbnpIDSFJaUS2t
zC=S|5A9WGQ_sZZ;N?bxEyB}4+X<c8r<DLh;%w2m<`(Ky%xyScHOxVih;qnT26Yz5N
zZYJD~*{feKX}=guWkP`@7RKN%Zf3&6_`xEUz+F-@la}QSC9uLJo(xA!>$i-+U@{(2
zW>i9if#@)tk|GjKp_I-<Q)Xt=#FCmsG>y=aq&XZ<4us=oB$G^~&G7hy7)0oJG!u;=
zcv%Z`2={U`oJxgvnbAZhwF`@S!lO|$GCn%G3obUvq`Q8JDcTfGL}b6ZwzY2VG&{F+
zV3m>Cam$w0%{}dK$8^i=G)zltbZir5ckkA=)?RaKS66@Mc5{1cTW_Z+R}<;Hyoc8G
zoU?+csKw$Rx$LM2Nb_hUoCynSc~OamS3Jq}OggP7qLoC9F05;kDq3tL6QWiql|(Yl
z<j{~jKzD?ZU36TnYI0anORS|&#cFKJtrqGjh8r4B3|jFVtF>^z(n)h9oQTAu!tB|K
zsF7I098X6h_K|`?LaVME&EUqy5)qSXz^EIbFp;<^r)c3!lv;@?n>iK(ZVsf=s@ZbI
zmR!(P64P4wRRYn{uM)sCw{LFeVNmQL62{dRF?*YeAUD=ZN6E@GT}repOpZe#0_k0&
z=o5G*W$_U%VRczFH6{Xbx;&g14~(U-JS~;kWrGIBV_04pizukIt*4&O%K({=gwrD;
z5ZRSL6&BB=EX^I!R63SSln5qtNNJcHRB0@p5dk_~L?9ELz@E-B=#p}n2BIVCyc&tX
zmnF-?I>9UtE+cSwG&YD@lZb^~KvrM{!*LV1Sd5P1^c4a1bRoMAoO&aB8VG0bbLqf%
z2A7KeJJIttF8P?-SX5i1%L1>b=I3ptmbP3b{yuE<jLZJb%AeN{>(JSW=96}6f?)Tb
zQ2xA*xXx~jGT$E=?s53@JoO;1-$a_+FYKKDAA#x(JX3S~XI1^P%Afns?-88s55u3{
zA+Z0v@}F1!lP;@n>LV4|pP%<*z{r&S`Mt$XbwfT0ot;#;BJ?qo=p6z3^Lv3V)nH9o
za{Jki`9DHW?;Kd>_YU&QdFaNJKOZA*zk<aK3Y1g!=l2{ZRE0rh$n~>7_y0lV-=^&N
zy@`h^1o3+o^1l?D)Bh1*loO6$ykBEVMa1uuEa}|wa}Wx<zozc@n)07#N5vIPGd4+I
z!QSpK)I122I!VjbIItukG%hIqO>nzE&u>?%2PA%9#%|8>`=s*c{PR4LE1F_Q#W{Zr
z|G?qT>*uF6%c7SJ6?fMEJaB41*MCgYrR6bAGEqlL8#(L$8FaKhl+(n|QM}IHY0onA
zm*DHO`SW|iGxY#LY-dReGULC8g2oQl&(B%>e-q|a`_*tP#^dAv8hpuw{dv9ciFE)$
zY-dU5*8iVS*!}srEvNh^*->$BH*+_f?f(4yG7_{LoX3l7?L^;qP%Lgg&+p^!(Ek!u
zzcMd|X8+rYh*jYFnZFwas@v)RoC3D8p^C#I^VKNuu_2SOow~HFZIDdtV?#?k=l(%z
z9=pm!vahbUBwbyHFNa9)hB@zi{i64`oIA_I+UJ@r-tSO4nj<e04io8k?RA#_wGIb|
GDE=2G!{ZA8

literal 0
HcmV?d00001

diff --git a/matrix multiplication/generate_matrix b/matrix multiplication/generate_matrix
new file mode 100755
index 0000000000000000000000000000000000000000..7691057b208dac69297b448ae480ad08837d2b91
GIT binary patch
literal 17520
zcmeHPeQ;FO6~DV75UOkz@FPk`o)&Cynq@aYhynF&Ajzha2ogZss*g?fZL;HLH{E@~
z;B;C#0VY9kqD73YGa{{)YIUqyq_#p3>sQ+uv@@MdrysE;A3KVGwGJ{<({t`SXY=;4
zOC6`vfBMeMd*`0tJ-_>L-hCf;-@E&rx~6)!%OzOc;zog%f3_+ib1y{KN@)^SEoKQ-
zC}O^t4#+gvRGEX?OtBtzn;EQ2$gc=A^2#tFwVI>9C@mpaa|=-;Z%R~3ntB);1j1T)
zO&fW|m?T?xRkXlXr8-0HWX*or`V{FVgkt)$K12;;ZRKH8CMrzc#0@5ob%g5cAk929
zk$0AOXNkvpgz96>@x++mR!#A&E>gil+tlWM>NweXJ;duFf7Yw95D?Z_2TIGR5&1mz
zttH-|B8w4&6i?PX%vf`Mw*Ze0tx5b3dp^wTsJ=qSp^W0sn!+m<$78KkOBcr@74cZ2
zyQiY3x~igTsV|-MRmz&NU4GbbOf}rRLAbYyA;g&P5-w@a1PsMZHs_=7C9q-Gk^k~w
z`}kcy`1#CJ=REIT@yDIl4&OTC4EkU{kb^!b=+8q2ZOr2h?awNggX5e$&Pxj5bp`OI
z0ywrwfz2*o0U%$!R~Nvs|M~p!x{wcVDS&^u04_AGy)&86(ndICXqxEKQ>kQ1L{qvh
z^q!a@qRB2jAyVN)M0AGZ@njn?U~IGzqca>!h-h0pnbt*DDwZ&!B6@c!X6OQ(+LN7K
zIAj87v|*}BG&D8VuF@)fOO|GnmA+-rIyW3!H;xmz#}r%bkjvxwh5PryX;37-Nm%!l
zF&JWQ@hz}XhHp}(#&g9wPDU)aybek)*r}>a&44`~GnNtyPL~Fm<+b2AK3U2vxb-|$
zEI8J|Qf|TdT7-Fi3(n^b;?)-X3kZ;?7W^Ux9qmNGiGULUCjw3coCr7(a3WAJ0{<wz
z>RhPz{gTkV#cwYbBD8(Ta8Di&^}bqiOmZfdy#?^(!r#KabiN8Hv^y{_GcpOX@D;T2
z(J^ye+ApE)Mtexw&!LTvkeO$t{UqA>=$P3j?fq!uBWY$p+TTYTA0abarTqZf_~@AF
zlJ-ut@sTmJ3G5?t6;(_e#1yt+;#Ziu7P0+V>;R-n=Wl~-=pgj<$V$w=2K<jSVyX;M
zFJbC1W(0<g+=gf+umVFe5>X%12QXfNla2jv-x=!vTd4QM>9y<Y0z-i}L$AMn+6&>i
zwAUxZM^gm!KQ?6!99pZYB6I`pjk$+piy-7Pk03BP5*;`eopQ!_AC~^#{XnYz)1jju
z+!Q)`wkYI!E%f%M#_VhfGpGbRAM4+GLl_z*y60{P^^Z>61O0lfcnu(~JC2=KANs$y
zbE*x&jfcGujo^mh`VB1;H-0)f2}7iG+jdAC!5Kqi|F0*?z-a7$Rc363^vgnlK<^&J
z{O0~KOOK)h6Dz<WT6$`zvG=u-P}}kTqoLS4^S8fu%h6BfhmM|Gxah6B-ia(a{$c5L
zLm!sjS$uMq5FZ9hXV;%?y#6oUuQc^%Uc$l;<MD6mAD{U1zo0?ACtckSHT9pGNP_hZ
zfbR(29=s!XXJ!r9;D~*aV>~;m0d&NPfD-{H0!{>+2sjaNBH%>8iGULUCj$Rx5x}=i
zR=K-XjBk5v?D+@B`oEXid+*qjd&VB?6B7pp$95ddJhdmYefRjw&y4=)XJZdO1n<L)
zRIEn}`dZ-~7#X$F5Mun%pJX2Xp%PR^AK0bTDx)txH`X^Wy7STTJ^RP^JbL2L6QfTL
zj&6TkS(ug{3)9w@I%0w9npd;}-vQ$f0ez<?C%1wwKRr2l0`zU5r@@?k5cF+ezX3W9
zy6?>7<N?qq_-+-&b<aBC>hZefT|A>?066#_4c{GB1Ls<FSUtmwm>;(FfbqMr1(*e!
zfV%gg_IT?(WnV45=<brO;-)KCT(@MwH3(P3hIJ1E_b39j9`Akb22a`UqFRr#Yg*7#
zzN0wk@%K)5Kkg}C?NMqyW#Cflff+2e1KYb2>L^Dpy6q_5HEnm%eRx!`4{yNM1o`-e
zP^|WNyO8LJ69FdzP6V6?I1z9n;6%WQfD-{H0!{=j5CQKTRpyS+cik`*lJa%ZY~#0v
zC>NV5`+YRLhm&%G%DIAU{+@Xr**rVBlI+v;{W8vVP&ogepH3#x9x75%ncrs1go=AP
z??XL8`8+-BLlA_&^(!+?aVcq5{LSJ7efP+-++lLy`5}Lkh~L9XnXXF3`0--X=62&e
z5rw}^<oe;QrtHr?@`snK(&o6r^ibM;ly8>uzuWWo<MtY>;H5Q6J?V|4qog;J-a&dV
z=_g1ZB#mV>tXfs0l*5<Gi<H1pU%(eoD*cs9{mcCUrF@+pQ9@yZGAfoYvUyjs_w|9w
z9B;^|Sh7e2H`Xle3REszF4k?VNp&X@u|&Jl-qvQXNz`tv2`3Us_(>vfcEhbT{+{KL
zD*v*oO2y|}oZga-=w0c>@O}5<aBHk0ku-E)n`l{EvqZG4stJfy^)=NOaI9~t@n0aj
zAqYjRZLC=;TIy>2BGv-CoF>E)F$3z8wNP3XN$DJ4_(s`?wN<2a_zKx<qhGOMvK$*G
z%d=tfJoMs0g8|wQT;N&~4rJGSY@}-hj(NL}PI;{J`3M>lmyyor!`G4b>3WjSr;F;k
z>;gWV$2wm*?-HN&=k?+K%)3CpxIbqu&}ME2Ydj2v1g{^{#2NgIAxE$eF35R=u8#`k
z&H<gzAAiqs3y&D0wwIHCZ2|vg^5<hf?qUt$`{>v(-by&%Ul{L!giBl@np9Ond_CMI
z@k<2XA7m!%o`MbA$@eS9Uxb89%oTjUVf+t}$k*S$O8+_Hg{xI|&FXyqC#C<T!hW6n
zSh8kk*QszJGod3cak=383*++v$75g@e*yd%z!gy@LKGk7HA{bwP-q+-rHm-xUhwC5
zvVW)aXKPAtMciE~_<ZA>9R>U!27I>r{QmwFaIDwf&R-Pp{~h28S7y<FlKx&%qNo-W
zKITFkuQr5FqujHA=Zoh>?)>e11N^-(-tc=Vt1JNjIc|^GMB|o!=Y39=@5x^K<Zr@2
zTnYYIjy=vHEC<FXe(Z%Jf0qT_EO80t?6d-SssP?w0N)GvY<Jmt{n{__^N-iT0{CwV
z;C}~vjvJo(s9)TP4B~89bz8~^_>!;~2JV`g5pEaldO}Zy4P7(e1i%d{i@@R;EgjxG
zm6eP~1QgJ!x5pAtNJLLTDsAW~&FIu%F-<~GL&=e()*er`hT~eqNT$+SxVuNR!IGM|
zZs-x8Kfi|-g<Gc<PNl+IG(BOYw!m7RaHp<Cx;s0!0L5Zyc$3X>g*NGli0oGVx?ppi
zR(JDiSQMnK{@TsK=EhY(ho>8@E~J^cq1Eez*3h)3HrS-CsjqLTTd%DT);85?a{Z9X
zOMP&C&v}c6vbnGjs8FDUrbog?nAR6%8L~X=CapNi&Crt3XgnE?Xq&?jGw2S+2}&^+
zEP=uK(e6YWJLIf^0{?VU>wu#c*M-)&2AGjpLhDZJyauX3jWIaRSSpN%R-{JantWKn
zQDHuSXk9VDwbpc+nq#h?!ZlVoWtz*dav)rdl>?aTta2btTi?8j4_daP5D0Ej7PB@p
z3v&A@DmXmk%CDR{gfG3N6FLamNSV5WQ?UL^Pjv~Oe2Q#Obo;te(6yAY#R9c<$6!@q
zEJ9GQwy^?FA0N7Rgwq|u7uk}4Dokyp%$&`7DjiEEatsafAiiNFP^v3#2%l`X&(M1y
zEr*FOm6YSir+3i!?}z}GEt6rM2quHmP;j_2)&{jEp)A<>WCdnna1ezL&c{wT&4iC0
zx@Fe^$KJ@E`oadhC2Z|B;FSI^2+lL%lAV2PSyrUVQh6OQ&od<mJM^{6m%zZe9bDF=
zR1=@q&8md6!sjfjRDogT_Yj}g&*F2IRro%}SR*8?e13j=geI&CJCe5ZZw35gm`if|
zhp7G`;&cD;y(x;l{b9gx?#%oV;*W6sE@YW)#3b{1o*M%UU70U<K1(;$KFG5Q6_<r@
z-6GD7na}g#KH_sd+<x|B{cYgmT$}BERR2EW`>7_3v9<jKi`|gGm@=Q|_s4B~uHWAO
zhlpRx1x{DZ6cL6AM-zWju=5`S4BN%!XYaT0=!3$)qnMW2R!*aS1U4&Qq5HK${C((T
z5rXl_O3II6Z{-gWe~28e=D?6Y^T-rW+xYxkd4e7iX0an_`}lpH_}OX{)$~+MHE62r
zp5w>BFCl^LXFjhR=k??(*^#up{x@y>VXA*vF$p84%I@v;{~mHYu+iy~osaPPa=Sm<
ztZ`maYT@(z=X3>tAna#aQ#Rx8K?aT;uAk>C{C@!>7#t8B|0vAk<M<I^=)!znFFd>$
zz44hhSD4TG6UesmdHyp&{7t#^Z2io41KG;w`9#QX679zeeXUXg0G7q==jZdL&&Z!e
z^^<$HXy)IZMa%@(&-yY*VBL290Rq;rBbCD@>xGcuV}mYhH>uL|Y(Ofpiyci<UaKf6
z6=w6;x8i+EOzT++;8kGGg^lCR*H!%8l4EBw%zbsWsn6L8w>1amX!-e-D@^_h8^LCZ
Fe*<7=09OD2

literal 0
HcmV?d00001

diff --git a/matrix multiplication/generate_matrix.c b/matrix multiplication/generate_matrix.c
new file mode 100644
index 0000000..e6a3544
--- /dev/null
+++ b/matrix multiplication/generate_matrix.c	
@@ -0,0 +1,37 @@
+#include <stdio.h>
+#include <stdlib.h>
+
+#define N 8192
+#define FILE_A "matrix_A.bin"
+#define FILE_B "matrix_B.bin"
+
+void generate_and_save_matrix(const char *filename) {
+    FILE *file = fopen(filename, "wb");
+    if (!file) {
+        perror("無法打開檔案");
+        exit(EXIT_FAILURE);
+    }
+
+    double *row = (double *)malloc(N * sizeof(double));
+    if (!row) {
+        perror("記憶體配置失敗");
+        exit(EXIT_FAILURE);
+    }
+
+    for (size_t i = 0; i < N; i++) {
+        for (size_t j = 0; j < N; j++) {
+            row[j] = (double)(rand() % 100) / 10.0;
+        }
+        fwrite(row, sizeof(double), N, file);
+    }
+
+    fclose(file);
+    free(row);
+}
+
+int main() {
+    generate_and_save_matrix(FILE_A);
+    generate_and_save_matrix(FILE_B);
+    printf("矩陣 A 和 B 已成功產生並存入 %s 和 %s\n", FILE_A, FILE_B);
+    return 0;
+}
\ No newline at end of file
diff --git a/matrix multiplication/makefile b/matrix multiplication/makefile
new file mode 100644
index 0000000..ed22bd3
--- /dev/null
+++ b/matrix multiplication/makefile	
@@ -0,0 +1,12 @@
+CC = gcc
+CFLAGS = -fopenmp -O2
+TARGET = generate_matrix
+SRC = generate_matrix.c
+
+all: $(TARGET)
+
+$(TARGET): $(SRC)
+	$(CC) $(CFLAGS) -o $(TARGET) $(SRC)
+
+clean:
+	rm -f $(TARGET)
\ No newline at end of file