diff options
-rw-r--r-- | src/InitialData.c | 12 |
1 files 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); } } } |