1 files changed, 13 insertions, 13 deletions
diff --git a/Carpet/CarpetLib/src/interpolate_3d_3tl.cc b/Carpet/CarpetLib/src/interpolate_3d_3tl.cc
index 6fdaa854d..f605f67e7 100644
--- a/Carpet/CarpetLib/src/interpolate_3d_3tl.cc
+++ b/Carpet/CarpetLib/src/interpolate_3d_3tl.cc
@@ -6,6 +6,8 @@
 #include <cctk.h>
 #include <cctk_Parameters.h>
 
+#include <loopcontrol.h>
+
 #include "operator_prototypes_3d.hh"
 #include "typeprops.hh"
 
@@ -124,19 +126,17 @@ namespace CarpetLib {
     
     
     // Loop over region
-#pragma omp parallel for
-    for (int k=0; k<regkext; ++k) {
-      for (int j=0; j<regjext; ++j) {
-        for (int i=0; i<regiext; ++i) {
-          
-          dst [DSTIND3(i, j, k)] =
-            + s1fac * src1 [SRCIND3(i, j, k)]
-            + s2fac * src2 [SRCIND3(i, j, k)]
-            + s3fac * src3 [SRCIND3(i, j, k)];
-          
-        }
-      }
-    }
+#pragma omp parallel
+    LC_LOOP3 (interpolate_3d_3tl,
+              i,j,k, 0,0,0, regiext,regjext,regkext, regiext,regjext,regkext)
+    {
+      
+      dst [DSTIND3(i, j, k)] =
+        + s1fac * src1 [SRCIND3(i, j, k)]
+        + s2fac * src2 [SRCIND3(i, j, k)]
+        + s3fac * src3 [SRCIND3(i, j, k)];
+      
+    } LC_ENDLOOP3(interpolate_3d_3tl);
     
   }