diff options
Diffstat (limited to 'src/ICN.c')
-rw-r--r-- | src/ICN.c | 6 |
1 files changed, 6 insertions, 0 deletions
@@ -120,6 +120,7 @@ void MoL_ICNAdd(CCTK_ARGUMENTS) RHSVar = (CCTK_REAL*)CCTK_VarDataPtrI(cctkGH, 0, RHSVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateVar[index] = OldVar[index] + CCTK_DELTA_TIME * RHSVar[index]; @@ -155,6 +156,7 @@ void MoL_ICNAdd(CCTK_ARGUMENTS) /* "This proc array total size is %d.", */ /* arraytotalsize); */ +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] = OldVar[index] + CCTK_DELTA_TIME * RHSVar[index]; @@ -174,6 +176,7 @@ void MoL_ICNAdd(CCTK_ARGUMENTS) RHSComplexVar = (CCTK_COMPLEX*)CCTK_VarDataPtrI(cctkGH, 0, RHSComplexVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateComplexVar[index] = CCTK_CmplxAdd(OldComplexVar[index], @@ -269,6 +272,7 @@ void MoL_ICNAverage(CCTK_ARGUMENTS) RHSVar = (CCTK_REAL*)CCTK_VarDataPtrI(cctkGH, 0, RHSVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { /* UpdateVar[index] = 0.5 * (UpdateVar[index] + OldVar[index]); */ @@ -302,6 +306,7 @@ void MoL_ICNAverage(CCTK_ARGUMENTS) arraytotalsize *= arraydata.lsh[arraydim]; } +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { /* UpdateVar[index] = 0.5 * (UpdateVar[index] + OldVar[index]); */ @@ -323,6 +328,7 @@ void MoL_ICNAverage(CCTK_ARGUMENTS) RHSComplexVar = (CCTK_COMPLEX*)CCTK_VarDataPtrI(cctkGH, 0, RHSComplexVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateComplexVar[index] = CCTK_CmplxMul(Complex_Half, CCTK_CmplxAdd(UpdateComplexVar[index], OldComplexVar[index])); |