diff options
| author | allen <allen@90eea020-d82d-4da5-bf6e-4ee79ff7632f> | 2001-04-18 20:37:27 +0000 |
|---|---|---|
| committer | allen <allen@90eea020-d82d-4da5-bf6e-4ee79ff7632f> | 2001-04-18 20:37:27 +0000 |
| commit | ca0e8efe2bd2797b6ac769fdf5c93b9fb022dc54 (patch) | |
| tree | 8f8c10a288c4d2692467e547111bc425438d513e /src | |
| parent | 75aa8fa949a06d574d96185a5e44e77bca220660 (diff) | |
Updating timelevel names, and some formatting changes.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyC/trunk@56 90eea020-d82d-4da5-bf6e-4ee79ff7632f
Diffstat (limited to 'src')
| -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); } } } |