Small optimisation: substitute division by multiplication with the inverse.

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

1 files changed, 3 insertions, 2 deletions

diff --git a/src/RK45.c b/src/RK45.c
index 41a0462..b2f6c88 100644
--- a/src/RK45.c
+++ b/src/RK45.c
@@ -122,6 +122,8 @@ void MoL_RK45Add(CCTK_ARGUMENTS)
 
   for (var = 0; var < MoLNumEvolvedVariables; var++)
   {
+    const CCTK_REAL tmp = (*Original_Delta_Time) / cctkGH->cctk_timefac;
+
     UpdateVar = (CCTK_REAL *)CCTK_VarDataPtrI(cctkGH, 0, 
                                               EvolvedVariableIndex[var]);
     RHSVar = (CCTK_REAL *)CCTK_VarDataPtrI(cctkGH, 0, 
@@ -134,8 +136,7 @@ void MoL_RK45Add(CCTK_ARGUMENTS)
                                                (*MoL_Intermediate_Step)));
     for (index = 0; index < totalsize; index++)
     {
-      ScratchVar[index] = (*Original_Delta_Time) / cctkGH->cctk_timefac *
-        RHSVar[index];
+      ScratchVar[index] = tmp * RHSVar[index];
     }
   }