diff options
Diffstat (limited to 'src/RK4.c')
-rw-r--r-- | src/RK4.c | 8 |
1 files changed, 8 insertions, 0 deletions
@@ -156,6 +156,7 @@ CCTK_WARN(0, "not implemented"); beta); #endif +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateVar[index] = OldVar[index] + @@ -176,6 +177,7 @@ CCTK_WARN(0, "not implemented"); /* scratch storage */ if ((*MoL_Intermediate_Step) == MoL_Intermediate_Steps) { +#pragma omp parallel for for (index = 0; index < totalsize; index++) { ScratchVar[index] = 0; @@ -184,6 +186,7 @@ CCTK_WARN(0, "not implemented"); if ((*MoL_Intermediate_Step)>1) { +#pragma omp parallel for for (index = 0; index < totalsize; index++) { ScratchVar[index] += alpha * UpdateVar[index]; @@ -191,6 +194,7 @@ CCTK_WARN(0, "not implemented"); } else { +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateVar[index] += ScratchVar[index] - 4.0 / 3.0 * OldVar[index]; @@ -240,6 +244,7 @@ CCTK_WARN(0, "not done"); ScratchVar = &ArrayScratchSpace[arrayscratchlocation]; +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] = OldVar[index] + @@ -248,6 +253,7 @@ CCTK_WARN(0, "not done"); if ((*MoL_Intermediate_Step) == MoL_Intermediate_Steps) { +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { ScratchVar[index] = 0; @@ -256,6 +262,7 @@ CCTK_WARN(0, "not done"); if ((*MoL_Intermediate_Step)>1) { +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { ScratchVar[index] += alpha * UpdateVar[index]; @@ -263,6 +270,7 @@ CCTK_WARN(0, "not done"); } else { +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] += ScratchVar[index] - 4.0 / 3.0 * OldVar[index]; |