diff options
-rw-r--r-- | src/RK4-RK2.c | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/RK4-RK2.c b/src/RK4-RK2.c index dd1071f..c310b9d 100644 --- a/src/RK4-RK2.c +++ b/src/RK4-RK2.c @@ -85,7 +85,7 @@ void MoL_RK4_RK2_Add(CCTK_ARGUMENTS) case 0: for (int i=0; i<totalsize; ++i) { CCTK_REAL const scaled_rhs = dt * RHSVar[i]; - ScratchVar[i] = OldVar[i] + 1.0/3.0 * scaled_rhs; + ScratchVar[i] = OldVar[i] + 1.0/6.0 * scaled_rhs; UpdateVar[i] = OldVar[i] + 0.5 * scaled_rhs; } break; @@ -93,7 +93,7 @@ void MoL_RK4_RK2_Add(CCTK_ARGUMENTS) case 1: for (int i=0; i<totalsize; ++i) { CCTK_REAL const scaled_rhs = dt * RHSVar[i]; - ScratchVar[i] += 1.0/6.0 * scaled_rhs; + ScratchVar[i] += 1.0/3.0 * scaled_rhs; UpdateVar[i] = OldVar[i] + 0.5 * scaled_rhs; } break; @@ -101,7 +101,7 @@ void MoL_RK4_RK2_Add(CCTK_ARGUMENTS) case 2: for (int i=0; i<totalsize; ++i) { CCTK_REAL const scaled_rhs = dt * RHSVar[i]; - ScratchVar[i] += 1.0/6.0 * scaled_rhs; + ScratchVar[i] += 1.0/3.0 * scaled_rhs; UpdateVar[i] = OldVar[i] + scaled_rhs; } break; @@ -110,7 +110,7 @@ void MoL_RK4_RK2_Add(CCTK_ARGUMENTS) for (int i=0; i<totalsize; ++i) { CCTK_REAL const scaled_rhs = dt * RHSVar[i]; /* ScratchVar contains OldVar */ - UpdateVar[i] = ScratchVar[i] + 1.0/3.0 * scaled_rhs; + UpdateVar[i] = ScratchVar[i] + 1.0/6.0 * scaled_rhs; } break; |