diff options
Diffstat (limited to 'src/RK4-MR-2_1.c')
-rw-r--r-- | src/RK4-MR-2_1.c | 30 |
1 files changed, 12 insertions, 18 deletions
diff --git a/src/RK4-MR-2_1.c b/src/RK4-MR-2_1.c index 0a48cb0..914fc6f 100644 --- a/src/RK4-MR-2_1.c +++ b/src/RK4-MR-2_1.c @@ -336,10 +336,9 @@ CCTK_WARN(0, "not implemented"); beta[MoL_Intermediate_Steps - (*MoL_Intermediate_Step)]); #endif - ScratchVar = CCTK_VarDataPtrI(cctkGH, 0, + ScratchVar = CCTK_VarDataPtrI(cctkGH, var, scratchspace_firstindex - + var - + MoLNumEvolvedVariables * (MoL_Intermediate_Steps - (*MoL_Intermediate_Step))); + + (MoL_Intermediate_Steps - (*MoL_Intermediate_Step))); #pragma omp parallel for for (index = 0; index < totalsize; index++) @@ -371,10 +370,9 @@ CCTK_WARN(0, "not implemented"); //printf("Scratch Step %d, alpha %g \n", scratchstep, alpha[scratchstep]); - ScratchVar = CCTK_VarDataPtrI(cctkGH, 0, + ScratchVar = CCTK_VarDataPtrI(cctkGH, var, scratchspace_firstindex - + var - + MoLNumEvolvedVariables * scratchstep); + + scratchstep); #pragma omp parallel for for (index = 0; index < totalsize; index++) @@ -392,10 +390,9 @@ CCTK_WARN(0, "not implemented"); //printf("Scratch Step %d, beta %g \n", scratchstep, beta[scratchstep]); - ScratchVar = CCTK_VarDataPtrI(cctkGH, 0, + ScratchVar = CCTK_VarDataPtrI(cctkGH, var, scratchspace_firstindex - + var - + MoLNumEvolvedVariables * scratchstep); + + scratchstep); #pragma omp parallel for for (index = 0; index < totalsize; index++) @@ -427,10 +424,9 @@ CCTK_WARN(0, "not implemented"); beta[MoL_Intermediate_Steps - (*MoL_Intermediate_Step)]); #endif - ScratchVar = CCTK_VarDataPtrI(cctkGH, 0, + ScratchVar = CCTK_VarDataPtrI(cctkGH, var, scratchspace_firstindex_slow - + var - + MoLNumEvolvedVariablesSlow * (MoL_Intermediate_Steps - (*MoL_Intermediate_Step))); + + (MoL_Intermediate_Steps - (*MoL_Intermediate_Step))); #pragma omp parallel for for (index = 0; index < totalsize; index++) @@ -462,10 +458,9 @@ CCTK_WARN(0, "not implemented"); //printf("Scratch Step %d, alpha %g \n", scratchstep, alpha[scratchstep]); - ScratchVar = CCTK_VarDataPtrI(cctkGH, 0, + ScratchVar = CCTK_VarDataPtrI(cctkGH, var, scratchspace_firstindex_slow - + var - + MoLNumEvolvedVariablesSlow * scratchstep); + + scratchstep); #pragma omp parallel for for (index = 0; index < totalsize; index++) @@ -483,10 +478,9 @@ CCTK_WARN(0, "not implemented"); //printf("Scratch Step %d, beta %g \n", scratchstep, beta[scratchstep]); - ScratchVar = CCTK_VarDataPtrI(cctkGH, 0, + ScratchVar = CCTK_VarDataPtrI(cctkGH, var, scratchspace_firstindex_slow - + var - + MoLNumEvolvedVariablesSlow * scratchstep); + + scratchstep); #pragma omp parallel for for (index = 0; index < totalsize; index++) |