diff options
-rw-r--r-- | src/WaveToy.c | 20 |
1 files changed, 11 insertions, 9 deletions
diff --git a/src/WaveToy.c b/src/WaveToy.c index fd1a07e..38e170e 100644 --- a/src/WaveToy.c +++ b/src/WaveToy.c @@ -76,7 +76,9 @@ void WaveToyC_Evolution(CCTK_ARGUMENTS) /* Do the evolution */ factor = 2*(1 - (dt2)*(dx2i + dy2i + dz2i)); - + printf("bin: %f\n",phi_p[CCTK_GFINDEX3D(cctkGH,10,10,10)]); + printf("bin: %f\n",phi_p_p[CCTK_GFINDEX3D(cctkGH,10,10,10)]); + for (k=kstart; k<kend; k++) { for (j=jstart; j<jend; j++) @@ -85,15 +87,15 @@ void WaveToyC_Evolution(CCTK_ARGUMENTS) { index = CCTK_GFINDEX3D(cctkGH,i,j,k); - phi_n[index] = factor* - phi[index] - phi_p[index] + phi[index] = factor* + phi_p[index] - phi_p_p[index] + (dt2) * - ( ( phi[CCTK_GFINDEX3D(cctkGH,i+1,j ,k )] - +phi[CCTK_GFINDEX3D(cctkGH,i-1,j ,k )] )*dx2i - +( phi[CCTK_GFINDEX3D(cctkGH,i ,j+1,k )] - +phi[CCTK_GFINDEX3D(cctkGH,i ,j-1,k )] )*dy2i - +( phi[CCTK_GFINDEX3D(cctkGH,i ,j ,k+1)] - +phi[CCTK_GFINDEX3D(cctkGH,i ,j, k-1)] )*dz2i); + ( ( phi_p[CCTK_GFINDEX3D(cctkGH,i+1,j ,k )] + +phi_p[CCTK_GFINDEX3D(cctkGH,i-1,j ,k )] )*dx2i + +( phi_p[CCTK_GFINDEX3D(cctkGH,i ,j+1,k )] + +phi_p[CCTK_GFINDEX3D(cctkGH,i ,j-1,k )] )*dy2i + +( phi_p[CCTK_GFINDEX3D(cctkGH,i ,j ,k+1)] + +phi_p[CCTK_GFINDEX3D(cctkGH,i ,j, k-1)] )*dz2i); } } } |