diff options
-rw-r--r-- | src/InitialData.cc | 14 |
1 files changed, 2 insertions, 12 deletions
diff --git a/src/InitialData.cc b/src/InitialData.cc index 35a4733..1f659c6 100644 --- a/src/InitialData.cc +++ b/src/InitialData.cc @@ -59,14 +59,14 @@ extern "C" void IDScalarWaveCXX_InitialData(CCTK_ARGUMENTS) for(int j=0; j<cctk_lsh[1]; j++) { for(int i=0; i<cctk_lsh[0]; i++) - { + { int index = CCTK_GFINDEX3D(cctkGH,i,j,k); phi[index] = amplitude*cos(kx*x[index]+ky*y[index]+kz*z[index]+omega*cctk_time); phi_p[index] = amplitude*cos(kx*x[index]+ky*y[index]+kz*z[index]+omega*(cctk_time-dt)); } - } + } } } else if(CCTK_Equals(initial_data, "gaussian")) @@ -116,16 +116,6 @@ extern "C" void IDScalarWaveCXX_InitialData(CCTK_ARGUMENTS) } } - /* zero out corners of phi next */ - phi_n[CCTK_GFINDEX3D(cctkGH, 0, 0, 0)] = 0; - phi_n[CCTK_GFINDEX3D(cctkGH, 0, 0, cctk_lsh[2]-1)] = 0; - phi_n[CCTK_GFINDEX3D(cctkGH, 0, cctk_lsh[1]-1, 0)] = 0; - phi_n[CCTK_GFINDEX3D(cctkGH, 0, cctk_lsh[1]-1, cctk_lsh[2]-1)] = 0; - phi_n[CCTK_GFINDEX3D(cctkGH, cctk_lsh[0]-1, 0, 0)] = 0; - phi_n[CCTK_GFINDEX3D(cctkGH, cctk_lsh[0]-1, 0, cctk_lsh[2]-1)] = 0; - phi_n[CCTK_GFINDEX3D(cctkGH, cctk_lsh[0]-1, cctk_lsh[1]-1, 0)] = 0; - phi_n[CCTK_GFINDEX3D(cctkGH, cctk_lsh[0]-1, cctk_lsh[1]-1, cctk_lsh[2]-1)] = 0; - return; } |