declare loop variables inside openmp loops (to make them private)

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/ADMBase/trunk@68 d576a68a-b34a-40ae-82fc-004fa1a9d16f

1 files changed, 17 insertions, 39 deletions

diff --git a/src/Initialisation.c b/src/Initialisation.c
index ea389a1..c35260a 100644
--- a/src/Initialisation.c
+++ b/src/Initialisation.c
@@ -72,12 +72,10 @@ void ADMBase_SetDtShiftStateOff(CCTK_ARGUMENTS);
  @@*/
 void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS)
 {
-  int i;
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
 
-#pragma omp parallel for
-  for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+  for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
   {
     gxx[i] = 1.0;
     gyy[i] = 1.0;
@@ -98,8 +96,7 @@ void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::metric") > 1)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       gxx_p[i] = 1.0;
       gyy_p[i] = 1.0;
@@ -113,8 +110,7 @@ void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::metric") > 2)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       gxx_p_p[i] = 1.0;
       gyy_p_p[i] = 1.0;
@@ -129,8 +125,7 @@ void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS)
     
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::curv") > 1)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       kxx_p[i] = 0.0;
       kyy_p[i] = 0.0;
@@ -144,8 +139,7 @@ void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::curv") > 2)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       
       kxx_p_p[i] = 0.0;
@@ -179,20 +173,17 @@ void ADMBase_CartesianMinkowski(CCTK_ARGUMENTS)
  @@*/
 void ADMBase_LapseOne(CCTK_ARGUMENTS)
 {
-  int i;
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
 
-#pragma omp parallel for
-  for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+  for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
   {
     alp[i] = 1.0;
   }
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::lapse") > 1)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       alp_p[i] = 1.0;
     }
@@ -200,8 +191,7 @@ void ADMBase_LapseOne(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::lapse") > 2)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       alp_p_p[i] = 1.0;
     }
@@ -225,12 +215,10 @@ void ADMBase_LapseOne(CCTK_ARGUMENTS)
  @@*/
 void ADMBase_ShiftZero(CCTK_ARGUMENTS)
 {
-  int i;
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
 
-#pragma omp parallel for
-  for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+  for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
   {
     betax[i] = 0.0;
     betay[i] = 0.0;
@@ -239,8 +227,7 @@ void ADMBase_ShiftZero(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::shift") > 1)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       betax_p[i] = 0.0;
       betay_p[i] = 0.0;
@@ -250,8 +237,7 @@ void ADMBase_ShiftZero(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::shift") > 2)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       betax_p_p[i] = 0.0;
       betay_p_p[i] = 0.0;
@@ -277,20 +263,17 @@ void ADMBase_ShiftZero(CCTK_ARGUMENTS)
  @@*/
 void ADMBase_DtLapseZero(CCTK_ARGUMENTS)
 {
-  int i;
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
 
-#pragma omp parallel for
-  for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+  for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
   {
     dtalp[i] = 0.0;
   }
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtlapse") > 1)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       dtalp_p[i] = 0.0;
     }
@@ -298,8 +281,7 @@ void ADMBase_DtLapseZero(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtlapse") > 2)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       dtalp_p_p[i] = 0.0;
     }
@@ -323,12 +305,10 @@ void ADMBase_DtLapseZero(CCTK_ARGUMENTS)
  @@*/
 void ADMBase_DtShiftZero(CCTK_ARGUMENTS)
 {
-  int i;
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
 
-#pragma omp parallel for
-  for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+  for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
   {
     dtbetax[i] = 0.0;
     dtbetay[i] = 0.0;
@@ -337,8 +317,7 @@ void ADMBase_DtShiftZero(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtshift") > 1)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       dtbetax_p[i] = 0.0;
       dtbetay_p[i] = 0.0;
@@ -348,8 +327,7 @@ void ADMBase_DtShiftZero(CCTK_ARGUMENTS)
   
   if (CCTK_ActiveTimeLevels(cctkGH, "ADMBase::dtshift") > 2)
   {
-#pragma omp parallel for
-    for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
+    for(int i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
     {
       dtbetax_p_p[i] = 0.0;
       dtbetay_p_p[i] = 0.0;