Response 5

ex4/code/task_c.py

@@ -58,4 +58,4 @@ print("u = ", u)
 
 plt.tight_layout()
 plt.savefig("../task_c.png", dpi=300)
-plt.show()
+plt.show()

ex5/code/gh_response.txt

@@ -0,0 +1,16 @@
+Good summary of results.
+
+1: OK
+
+2: OK
+
+3: OK
+task3.cpp:32
+non initialzed memory for thread 0 (and the other threads)
+
+4: OK
+
+5:
+CRS_Matrix::GetDiag  non parallel
+CRS_Matrix::CalculateLaplace: 
+  "omp critical" is slower than an "omp atomic" for the +=  [worth to test that]