diff options
Diffstat (limited to 'src/RK2.c')
-rw-r--r-- | src/RK2.c | 6 |
1 files changed, 6 insertions, 0 deletions
@@ -118,6 +118,7 @@ void MoL_RK2Add(CCTK_ARGUMENTS) RHSVar = (CCTK_REAL const*)CCTK_VarDataPtrI(cctkGH, 0, RHSVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateVar[index] += CCTK_DELTA_TIME * RHSVar[index]; @@ -146,6 +147,7 @@ void MoL_RK2Add(CCTK_ARGUMENTS) arraytotalsize *= arraydata.lsh[arraydim]; } +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] += CCTK_DELTA_TIME * RHSVar[index]; @@ -163,6 +165,7 @@ void MoL_RK2Add(CCTK_ARGUMENTS) RHSComplexVar = (CCTK_COMPLEX const*)CCTK_VarDataPtrI(cctkGH, 0, RHSComplexVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateComplexVar[index] = CCTK_CmplxAdd(UpdateComplexVar[index], @@ -186,6 +189,7 @@ void MoL_RK2Add(CCTK_ARGUMENTS) RHSVar = (CCTK_REAL const*)CCTK_VarDataPtrI(cctkGH, 0, RHSVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateVar[index] = 0.5 * (OldVar[index] + UpdateVar[index]) + @@ -218,6 +222,7 @@ void MoL_RK2Add(CCTK_ARGUMENTS) arraytotalsize *= arraydata.lsh[arraydim]; } +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] = 0.5 * (OldVar[index] + UpdateVar[index]) + @@ -238,6 +243,7 @@ void MoL_RK2Add(CCTK_ARGUMENTS) RHSComplexVar = (CCTK_COMPLEX const*)CCTK_VarDataPtrI(cctkGH, 0, RHSComplexVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateComplexVar[index] = |