diff options
-rw-r--r-- | src/WaveToy.cc | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/src/WaveToy.cc b/src/WaveToy.cc index b730fce..66e3be3 100644 --- a/src/WaveToy.cc +++ b/src/WaveToy.cc @@ -78,13 +78,13 @@ extern "C" void WaveToyCXX_Evolution(CCTK_ARGUMENTS) { int index = CCTK_GFINDEX3D(cctkGH,i,j,k); - phi_n[index] = - factor*phi[index] - phi_p[index] - + dt2 * - ( ( val( phi, i+1,j ,k ) + val( phi, i-1,j ,k) )*dx2i - +( val( phi, i ,j+1,k ) + val( phi, i ,j-1,k) )*dy2i - +( val( phi, i ,j ,k+1) + val( phi, i ,j, k-1) )*dz2i - ); + phi[index] = + factor*phi_p[index] - phi_p_p[index] + + dt2 * + ( ( val( phi_p, i+1,j ,k ) + val( phi_p, i-1,j ,k) )*dx2i + +( val( phi_p, i ,j+1,k ) + val( phi_p, i ,j-1,k) )*dy2i + +( val( phi_p, i ,j ,k+1) + val( phi_p, i ,j, k-1) )*dz2i + ); } } } |