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