aboutsummaryrefslogtreecommitdiff
path: root/src/RK45.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/RK45.c')
-rw-r--r--src/RK45.c5
1 files changed, 5 insertions, 0 deletions
diff --git a/src/RK45.c b/src/RK45.c
index 3ca2456..6173d77 100644
--- a/src/RK45.c
+++ b/src/RK45.c
@@ -181,6 +181,7 @@ void MoL_RK45Add(CCTK_ARGUMENTS)
+ MoL_Num_Evolved_Vars *
(MoL_Intermediate_Steps -
(*MoL_Intermediate_Step)));
+#pragma omp parallel for
for (index = 0; index < totalsize; index++)
{
ScratchVar[index] = tmp * RHSVar[index];
@@ -207,6 +208,7 @@ void MoL_RK45Add(CCTK_ARGUMENTS)
if (*MoL_Intermediate_Step - 1)
{
+#pragma omp parallel for
for (index = 0; index < totalsize; index++)
{
UpdateVar[index] = OldVar[index];
@@ -226,6 +228,7 @@ void MoL_RK45Add(CCTK_ARGUMENTS)
if ( (beta > MoL_Tiny)||(beta < -MoL_Tiny) )
{
+#pragma omp parallel for
for (index = 0; index < totalsize; index++)
{
UpdateVar[index] += beta * ScratchVar[index];
@@ -238,6 +241,7 @@ void MoL_RK45Add(CCTK_ARGUMENTS)
else
{
+#pragma omp parallel for
for (index = 0; index < totalsize; index++)
{
UpdateVar[index] = OldVar[index];
@@ -257,6 +261,7 @@ void MoL_RK45Add(CCTK_ARGUMENTS)
if ( (gamma > MoL_Tiny)||(gamma < -MoL_Tiny) )
{
+#pragma omp parallel for
for (index = 0; index < totalsize; index++)
{
UpdateVar[index] += gamma * ScratchVar[index];