7 files changed, 191 insertions, 0 deletions

diff --git a/.gitignore b/.gitignore
new file mode 100644
index 0000000..eaa8d95
--- /dev/null
+++ b/.gitignore
@@ -0,0 +1,4 @@
+.envrc
+.ccls
+.ccls-cache
+.direnv
diff --git a/flake.lock b/flake.lock
new file mode 100644
index 0000000..d24a185
--- /dev/null
+++ b/flake.lock
@@ -0,0 +1,27 @@
+{
+  "nodes": {
+    "nixpkgs": {
+      "locked": {
+        "lastModified": 1633083543,
+        "narHash": "sha256-thXKms0SvDirYb8/Hw0zqAE1TsuSFWd0y51mSzkPGLU=",
+        "owner": "NixOS",
+        "repo": "nixpkgs",
+        "rev": "92609f3d9bc3acffbdbe54fa1c591a885612aa73",
+        "type": "github"
+      },
+      "original": {
+        "owner": "NixOS",
+        "ref": "nixos-21.05",
+        "repo": "nixpkgs",
+        "type": "github"
+      }
+    },
+    "root": {
+      "inputs": {
+        "nixpkgs": "nixpkgs"
+      }
+    }
+  },
+  "root": "root",
+  "version": 7
+}
diff --git a/flake.nix b/flake.nix
new file mode 100644
index 0000000..3308ce1
--- /dev/null
+++ b/flake.nix
@@ -0,0 +1,58 @@
+{
+  description = "A very basic flake";
+
+  inputs = {
+    nixpkgs.url = "github:NixOS/nixpkgs/nixos-21.05";
+  };
+
+  outputs =
+    { self
+    , nixpkgs
+    }:
+    let
+      supportedSystems = [ "x86_64-linux" ];
+      forAllSystems = f: nixpkgs.lib.genAttrs supportedSystems (system: f system);
+
+      nixpkgsFor = forAllSystems (system:
+        import nixpkgs { inherit system; }
+      );
+
+    in
+    {
+      devShell = forAllSystems
+        (system:
+          let
+            pkgs = nixpkgsFor."${system}";
+          in
+          with pkgs;
+          mkShell {
+            nativeBuildInputs = [ gcc ];
+          });
+
+      apps = forAllSystems
+        (system:
+          let
+            pkgs = nixpkgsFor."${system}";
+            execs = with builtins; map toString [ 1 2 ];
+            mkApp = name: with pkgs; stdenv.mkDerivation {
+              name = "${name}";
+              src = ./src;
+              buildInputs = [ gcc ];
+              unpackPhase = ''
+                true
+              '';
+              buildPhase = ''
+                gcc -o main $src/${name}/main.c -fopenmp
+              '';
+              installPhase = ''
+                install -Dm755 main $out/bin/main
+              '';
+            };
+          in
+          with pkgs;
+          lib.genAttrs execs (p: {
+            type = "app";
+            program = "${mkApp p}/bin/main";
+          }));
+    };
+}
diff --git a/src/1/main.c b/src/1/main.c
new file mode 100644
index 0000000..975c495
--- /dev/null
+++ b/src/1/main.c
@@ -0,0 +1,13 @@
+#include <omp.h>
+#include <stdio.h>
+
+int main() {
+    long long int total = 1, i;
+#pragma omp parallel for reduction(*:total)
+    for(i = 1; i <= 20; i++) {
+        // printf("working in thread %d\n", omp_get_thread_num());
+        total = total * i;
+    }
+    printf("total: %lld\n", total);
+}
+
diff --git a/src/2/main b/src/2/main
new file mode 100755
index 0000000..b999f77
--- /dev/null
+++ b/src/2/main
diff --git a/src/2/main.c b/src/2/main.c
new file mode 100644
index 0000000..ee7f2c8
--- /dev/null
+++ b/src/2/main.c
@@ -0,0 +1,82 @@
+#include <stdio.h>
+#include <stdlib.h>
+
+#include <omp.h>
+
+int main() {
+    int **a, **b, **c, l, m1, m2, n; 
+
+    printf("Order of matrix A: ");
+    int _ = scanf("%d %d", &l, &m1);
+
+    printf("Order of matrix B: ");
+    _ = scanf("%d %d", &m2, &n);
+
+    if (m1 != m2) {
+        printf("Invalid order\n");
+        exit(1);
+    }
+
+    int m = m1;
+
+    a = (int **) calloc(l, sizeof(int *));
+    for(int i = 0; i < m; i++) a[i] = (int *) calloc(m, sizeof(int));
+
+    b = (int **) calloc(m, sizeof(int *));
+    for(int i = 0; i < n; i++) b[i] = (int *) calloc(n, sizeof(int));
+
+    c = (int **) calloc(l, sizeof(int *));
+    for(int i = 0; i < n; i++) c[i] = (int *) calloc(n, sizeof(int));
+
+
+    printf("L: %d, M: %d, N: %d\n", l, m, n);
+
+    // input for A
+    // printf("Enter A:\n");
+    for(int i = 0; i < l; i++) {
+        for(int j = 0; j < m; j++) {
+            // _ = scanf("%d", &a[i][j]);
+            a[i][j] = rand() % 500;
+        }
+    }
+
+    // input for B
+    // printf("Enter B:\n");
+    for(int i = 0; i < m; i++) {
+        for(int j = 0; j < n; j++) {
+            // _ = scanf("%d", &b[i][j]);
+            b[i][j] = rand() % 1000;
+        }
+    }
+
+    // printf("\nA: \n");
+    for(int i = 0; i < l; i++) {
+        for(int j = 0; j < m; j++) {
+            // printf("%d\t", a[i][j]);
+        }
+        // printf("\n");
+    }
+
+    // printf("\nB: \n");
+    for(int i = 0; i < m; i++) {
+        for(int j = 0; j < n; j++) {
+            // printf("%d\t", b[i][j]);
+        }
+        // printf("\n");
+    }
+
+    double start_time = omp_get_wtime();
+    int i, j, k;
+#pragma omp parallel for private(j, k) num_threads(1)
+    for(i = 0; i < l; i++) {
+        for(j = 0; j < m; j++) {
+            c[i][j] = 0;
+            for(k = 0; k < n; k++) {
+                c[i][j] += a[i][k] * b[k][j];
+            }
+        }
+    }
+    double end_time = omp_get_wtime();
+
+    printf(">>= Done in %f secs", end_time - start_time);
+}
diff --git a/src/2/table b/src/2/table
new file mode 100644
index 0000000..8e89bcc
--- /dev/null
+++ b/src/2/table
@@ -0,0 +1,7 @@
+     | 1         | 2         | 4         | 8
+-----|-----------|-----------|-----------|----------
+ 100 | 0.000862  | 0.000813  | 0.000492  | 0.018057
+ 250 | 0.016979  | 0.007088  | 0.008280  | 0.023153
+ 500 | 0.187274  | 0.095052  | 0.083152  | 0.081629
+1000 | 2.522172  | 1.221258  | 1.139272  | 1.361240
+2000 | 47.200504 | 26.900284 | 19.395751 | 15.864646