aboutsummaryrefslogtreecommitdiff
path: root/src/ICN.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/ICN.c')
-rw-r--r--src/ICN.c6
1 files changed, 6 insertions, 0 deletions
diff --git a/src/ICN.c b/src/ICN.c
index 7f75a35..c9a6114 100644
--- a/src/ICN.c
+++ b/src/ICN.c
@@ -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]));