From 5f565bb6944f8896ca8ea80931be2e2b78344589 Mon Sep 17 00:00:00 2001 From: allen Date: Wed, 21 Aug 2002 15:44:37 +0000 Subject: Fixed error setting solution at special case of r=0 Closes Cactus/1085 Also made the gaussian in IDScalarWaveCXX consistent with gaussians in the other ID thorns, hence new testsuite results. git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWave/trunk@59 f5a6acaf-da7d-456b-b0a8-35edbc60b392 --- src/InitialData.F77 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/InitialData.F77 b/src/InitialData.F77 index 8a8da75..08de276 100644 --- a/src/InitialData.F77 +++ b/src/InitialData.F77 @@ -65,7 +65,7 @@ do k=1, cctk_lsh(3) do j=1, cctk_lsh(2) do i=1, cctk_lsh(1) - + xp=x(i,j,k) yp=y(i,j,k) zp=z(i,j,k) @@ -75,7 +75,7 @@ phi(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) if (rp .eq. 0.0) then - phi_p(i,j,k) = (1.0 - 2.0*dt**2/sigma)*exp(-dt**2/sigma) + phi_p(i,j,k) = amplitude*(1.0 - 2.0*dt**2/sigma)*exp(-dt**2/sigma) else phi_p(i,j,k) = amplitude/2.0*(rp-dt)/rp* & exp( - ( (rp-radius-dt)/sigma)**2 ) -- cgit v1.2.3