diff options
Diffstat (limited to 'src/RK3.c')
-rw-r--r-- | src/RK3.c | 9 |
1 files changed, 9 insertions, 0 deletions
@@ -121,6 +121,7 @@ void MoL_RK3Add(CCTK_ARGUMENTS) RHSVar = (CCTK_REAL*)CCTK_VarDataPtrI(cctkGH, 0, RHSVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateVar[index] += CCTK_DELTA_TIME * RHSVar[index]; @@ -149,6 +150,7 @@ void MoL_RK3Add(CCTK_ARGUMENTS) arraytotalsize *= arraydata.lsh[arraydim]; } +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] += CCTK_DELTA_TIME * RHSVar[index]; @@ -166,6 +168,7 @@ void MoL_RK3Add(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(UpdateComplexVar[index], @@ -190,6 +193,7 @@ void MoL_RK3Add(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.25 * (3*OldVar[index] + @@ -222,6 +226,7 @@ void MoL_RK3Add(CCTK_ARGUMENTS) arraytotalsize *= arraydata.lsh[arraydim]; } +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] = 0.25*(3*OldVar[index] + UpdateVar[index]) + @@ -242,6 +247,7 @@ void MoL_RK3Add(CCTK_ARGUMENTS) RHSComplexVar = (CCTK_COMPLEX const*)CCTK_VarDataPtrI(cctkGH, 0, RHSComplexVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateComplexVar[index] = @@ -271,6 +277,7 @@ void MoL_RK3Add(CCTK_ARGUMENTS) RHSVar = (CCTK_REAL const*)CCTK_VarDataPtrI(cctkGH, 0, RHSVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateVar[index] = (OldVar[index] + 2*UpdateVar[index]) * one_third @@ -303,6 +310,7 @@ void MoL_RK3Add(CCTK_ARGUMENTS) arraytotalsize *= arraydata.lsh[arraydim]; } +#pragma omp parallel for for (index = 0; index < arraytotalsize; index++) { UpdateVar[index] = (OldVar[index] + 2*UpdateVar[index]) * one_third @@ -323,6 +331,7 @@ void MoL_RK3Add(CCTK_ARGUMENTS) RHSComplexVar = (CCTK_COMPLEX const*)CCTK_VarDataPtrI(cctkGH, 0, RHSComplexVariableIndex[var]); +#pragma omp parallel for for (index = 0; index < totalsize; index++) { UpdateComplexVar[index] = |