Changed to use timelevels - thanks for putting the rotation in PUGH Ed.

Note that the no-boundary condition tests assumed that phi_new had zero
at its boundaries.  This is, or course, no longer true with timelevel
rotation.

Tom


git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWaveCXX/trunk@19 7ec00dc9-1e7a-42ad-bf73-a299ee168e72

1 files changed, 6 insertions, 6 deletions

diff --git a/src/InitialData.cc b/src/InitialData.cc
index 82e1bf8..8b00476 100644
--- a/src/InitialData.cc
+++ b/src/InitialData.cc
@@ -62,7 +62,7 @@ extern "C" void IDScalarWaveCXX_InitialData(CCTK_ARGUMENTS)
           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_old[index] = amplitude*cos(kx*x[index]+ky*y[index]+kz*z[index]+omega*(cctk_time-dt));
+          phi_p[index] = amplitude*cos(kx*x[index]+ky*y[index]+kz*z[index]+omega*(cctk_time-dt));
         }
       }
     }
@@ -81,7 +81,7 @@ extern "C" void IDScalarWaveCXX_InitialData(CCTK_ARGUMENTS)
           CCTK_REAL R = sqrt(X*X + Y*Y + Z*Z);
 
           phi[index] = amplitude*exp( - sqr( (R - radius) / sigma ) );
-          phi_old[index] = amplitude*exp( - sqr( (R - radius - dt) / sigma ) );
+          phi_p[index] = amplitude*exp( - sqr( (R - radius - dt) / sigma ) );
 
         } 
       }
@@ -105,10 +105,10 @@ extern "C" void IDScalarWaveCXX_InitialData(CCTK_ARGUMENTS)
                                  sin(kz*(z[index]-0.5)*pi)*
                                  cos(omega*cctk_time*pi);
 
-          phi_old[index] = amplitude*sin(kx*(x[index]-0.5)*pi)*
-                                         sin(ky*(y[index]-0.5)*pi)*
-                                         sin(kz*(z[index]-0.5)*pi)*
-                                         cos(omega*(cctk_time-dt)*pi);
+          phi_p[index] = amplitude*sin(kx*(x[index]-0.5)*pi)*
+                                       sin(ky*(y[index]-0.5)*pi)*
+                                       sin(kz*(z[index]-0.5)*pi)*
+                                       cos(omega*(cctk_time-dt)*pi);
         }
       }
     }