diff options
| -rw-r--r-- | src/InitialData.c | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/src/InitialData.c b/src/InitialData.c index 7eaa6a7..2261b97 100644 --- a/src/InitialData.c +++ b/src/InitialData.c @@ -4,7 +4,7 @@ @author Werner Benger @desc Initial data for the 3D Wave Equation - Derived from Tom Goodale + Derived from Tom Goodale @enddesc @version $Header$ @@*/ @@ -98,17 +98,17 @@ void IDScalarWaveC_InitialData(CCTK_ARGUMENTS) phi[index] = amplitude*exp( - sqr( (R - radius) / sigma ) ); - if (R == 0.0) - { - phi_p[index] = amplitude*(1.0 - 2.0*dt*dt/sigma)*exp(-dt*dt/sigma); - } - else - { - phi_p[index] = amplitude/2.0*(R-dt)/R* - exp( - sqr( (R - radius - dt)/ sigma ) ) - + amplitude/2.0*(R+dt)/R* - exp( - sqr( (R - radius + dt)/ sigma ) ); - } + if (R == 0.0) + { + phi_p[index] = amplitude*(1.0 - 2.0*dt*dt/sigma)*exp(-dt*dt/sigma); + } + else + { + phi_p[index] = amplitude/2.0*(R-dt)/R* + exp( - sqr( (R - radius - dt)/ sigma ) ) + + amplitude/2.0*(R+dt)/R* + exp( - sqr( (R - radius + dt)/ sigma ) ); + } } } } |