OpenMP parallelization

Author: arneish

compile.sh

@@ -0,0 +1 @@
+g++ -g -fopenmp -lm lab2_io.c lab2_omp_orig_parallel.c main_omp.c -o pca

eigen_py.py

@@ -0,0 +1,19 @@
+import numpy as NP
+from scipy import linalg as LA
+# A = [[108, 23, -33],[23, 121, -6],[-33, -6, 39]]
+# A = NP.array(A, dtype=float)
+# print(A)
+# e_vals, e_v = LA.eig(A)
+# print("e_vals:\n", e_vals)
+# print("e_v:\n", e_v)
+
+A = [[-1, 5, -9, 1],[9, 6, 0, 2],[2, -5, 1, 3]]
+A = NP.array(A, dtype=float)
+print(A)
+u, s, vh = NP.linalg.svd(A, full_matrices=True)
+print("u:\n", u)
+print("sigma:\n", s)
+
+temp = NP.matmul(u, s)
+print(temp)
+print("v_t:\n", vh)