From b68b294b343f555efeaa085fe08ba03aa6d9cd11 Mon Sep 17 00:00:00 2001 From: schnetter Date: Sat, 19 Mar 2005 15:29:53 +0000 Subject: Use sqrt instead of **0.5 git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@88 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7 --- src/planewaves.F77 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/planewaves.F77 b/src/planewaves.F77 index 4a40688..e33ece2 100644 --- a/src/planewaves.F77 +++ b/src/planewaves.F77 @@ -108,12 +108,12 @@ c precalc kx = 2*pi*sin(the)*cos(phi)/wave ky = 2*pi*sin(the)*sin(phi)/wave kz = 2*pi*cos(the)/wave - w = (kx*kx+ky*ky+kz*kz)**0.5 + w = sqrt(kx*kx+ky*ky+kz*kz) kxp = 2*pi*sin(the)*cos(phi)/wavep kyp = 2*pi*sin(the)*sin(phi)/wavep kzp = 2*pi*cos(the)/wavep - wp = (kxp*kxp+kyp*kyp+kzp*kzp)**0.5 + wp = sqrt(kxp*kxp+kyp*kyp+kzp*kzp) c *************** plane waves ******************** do k=1,cctk_lsh(3) -- cgit v1.2.3