diff -Nur MFCG.orig/src/Grid.cpp MFCG.para/src/Grid.cpp
--- MFCG.orig/src/Grid.cpp	2025-07-31 14:05:35.124111300 +0200
+++ MFCG.para/src/Grid.cpp	2025-07-31 15:27:01.244919711 +0200
@@ -15,9 +15,10 @@
 
     //always pad with halo; to support Dirichlet
     arrayPtr = new double[ rows*columns];
+#pragma omp parallel for
     for (int i =0; i<rows*columns;i++)
     {
-        arrayPtr[i] = 0.0;
+	 arrayPtr[i] = 0.0;
     }
 }
 
@@ -28,6 +29,7 @@
     }
 
     arrayPtr = new double[rows*columns];
+#pragma omp parallel for
     for (int i =0; i<rows*columns;i++)
     {
         arrayPtr[i] = 0.0;
@@ -49,6 +51,7 @@
         int totGrids = rows * columns;
         // performing a deep-copy
         arrayPtr = new double[totGrids];
+#pragma omp parallel for
         for(int i=0; i<totGrids;++i)
         {
             arrayPtr[i] = s.arrayPtr[i];