From 931ade1a9bc9c96ba40668f0a5defc7abf93ace1 Mon Sep 17 00:00:00 2001 From: goodale Date: Wed, 29 Mar 2000 07:42:11 +0000 Subject: 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 --- src/InitialData.cc | 12 ++++++------ 1 file 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); } } } -- cgit v1.2.3