From c5d918944cd5a28b42aeed0b61da1ff5a3061de5 Mon Sep 17 00:00:00 2001
From: JosephPro77vn <83563375+JosephPro77vn@users.noreply.github.com>
Date: Thu, 14 Oct 2021 23:16:38 +0800
Subject: [PATCH] Create matrix_kij

---
 .../Autumn/LY2021035-Demeke-Yoseph/matrix_kij | 131 ++++++++++++++++++
 1 file changed, 131 insertions(+)
 create mode 100644 2021/Autumn/LY2021035-Demeke-Yoseph/matrix_kij

diff --git a/2021/Autumn/LY2021035-Demeke-Yoseph/matrix_kij b/2021/Autumn/LY2021035-Demeke-Yoseph/matrix_kij
new file mode 100644
index 0000000..78c5303
--- /dev/null
+++ b/2021/Autumn/LY2021035-Demeke-Yoseph/matrix_kij
@@ -0,0 +1,131 @@
+#include <iostream>
+#include <ctime>
+#include <algorithm>
+
+using namespace std;
+const int N = 2000;
+int a[N][N];
+int b[N][N];
+long long c[N][N];
+
+
+
+// i-j-k
+void fun1(){
+    for(int i = 0 ; i < N ; i++){
+        for(int j = 0 ; j < N ; j++){
+
+            for(int k = 0; k < N; k++){
+                c[i][j] += a[i][k] * b[k][j];
+            }
+
+
+        }
+    }
+}
+
+// i-k-j
+void fun2(){
+    for(int i = 0 ; i < N ; i++){
+        for(int k = 0 ; k < N ; k++){
+            for(int j = 0; j < N; j++){
+                c[i][j] += a[i][k] * b[k][j];
+            }
+        }
+    }
+}
+
+// j-i-k
+void fun3(){
+    for(int j = 0; j < N; j++){
+        for(int i = 0 ; i < N ; i++){
+            for(int k = 0; k < N; k++){
+                c[i][j] += a[i][k] * b[k][j];
+            }
+        }
+    }
+}
+
+// j-k-i
+void fun4(){
+    for(int j = 0; j < N; j++){
+        for(int k = 0; k < N; k++){
+            for(int i = 0; i < N; i++){
+                c[i][j] += a[i][k] * b[k][j];
+            }
+        }
+    }
+}
+
+// k-i-j
+void fun5(){
+    for(int k = 0; k < N; k++){
+        for(int i = 0; i < N; i++){
+            for(int j = 0; j < N; j++){
+                c[i][j] += a[i][k] * b[k][j];
+            }
+        }
+    }
+}
+
+
+// k-j-i
+void fun6(){
+    for(int k = 0; k < N; k++){
+        for(int j = 0; j < N; j++){
+            for(int i = 0; i < N; i++){
+                c[i][j] += a[i][k] * b[k][j];
+            }
+        }
+    }
+}
+
+int main() {
+    fill(a[0],a[0] + N * N,1);
+    fill(b[0],b[0] + N * N,1);
+    fill(c[0],c[0] + N * N,0);;
+    clock_t start_time,end_time;
+
+    // i-j-k
+    start_time = clock();
+    fun1();
+    end_time = clock();
+    cout << "The i-j-k run time is: " <<
+    (double)(end_time - start_time) / CLOCKS_PER_SEC << "s" << endl;
+
+    // i-k-j
+    start_time = clock();
+    fun2();
+    end_time = clock();
+    cout << "The i-k-j run time is: " <<
+         (double)(end_time - start_time) / CLOCKS_PER_SEC << "s" << endl;
+
+    // j-i-k
+    start_time = clock();
+    fun3();
+    end_time = clock();
+    cout << "The j-i-k run time is: " <<
+         (double)(end_time - start_time) / CLOCKS_PER_SEC << "s" << endl;
+
+    // j-k-i
+    start_time = clock();
+    fun4();
+    end_time = clock();
+    cout << "The j-k-i run time is: " <<
+         (double)(end_time - start_time) / CLOCKS_PER_SEC << "s" << endl;
+
+    // k-i-j
+    start_time = clock();
+    fun5();
+    end_time = clock();
+    cout << "The k-i-j run time is: " <<
+         (double)(end_time - start_time) / CLOCKS_PER_SEC << "s" << endl;
+
+    // k-j-i
+    start_time = clock();
+    fun6();
+    end_time = clock();
+    cout << "The k-j-i run time is: " <<
+         (double)(end_time - start_time) / CLOCKS_PER_SEC << "s" << endl;
+    return 0;
+}