Parallelise loops with OpenMP.

git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@129 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b

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];