From a9aef810e4cdb1e9456ddb7bcd37eaccfc3735aa 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/IDScalarWaveC/trunk@6 071ba4cb-2270-47cd-a45a-3e542e115664 --- src/InitialData.c | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/src/InitialData.c b/src/InitialData.c index c03b04b..0d9b384 100644 --- a/src/InitialData.c +++ b/src/InitialData.c @@ -72,7 +72,7 @@ void IDScalarWaveC_InitialData(CCTK_ARGUMENTS) 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)); } } } @@ -94,7 +94,7 @@ void IDScalarWaveC_InitialData(CCTK_ARGUMENTS) 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 ) ); } } @@ -118,10 +118,10 @@ void IDScalarWaveC_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