From 1e7726df7ee9e189088a3869de2364e839003c45 Mon Sep 17 00:00:00 2001 From: rhaas Date: Mon, 13 Jan 2014 17:16:03 +0000 Subject: add m=2 ftp coefficient for teukolsky waves The theta-phi coefficient was not set at all and contained random data. This commit implements the coefficient as found in Eq. 7 of PRD 26:745 (1982) which is the paper referenced in the docs. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDLinearWaves/trunk@122 5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7 --- src/teukwaves.F77 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/teukwaves.F77 b/src/teukwaves.F77 index a7a365f..8e19e5a 100644 --- a/src/teukwaves.F77 +++ b/src/teukwaves.F77 @@ -358,7 +358,7 @@ c and outgoing to have the same sign. sinp = yp/(xp**2+yp**2)**0.5 cosp = xp/(xp**2+yp**2)**0.5 -c mvalue +c mvalue (Eq. 7 of PRD 26:745 (1982)) if (m.eq.0) then frr = 2-3*sint**2 frt = -3*sint*cost @@ -393,6 +393,7 @@ c mvalue frp = sint*(1-2*sinp**2) ftt1 = (1+cost**2)*2*sinp*cosp ftt2 = -2*sinp*cosp + ftp = cost*(1-2*cosp**2) fpp1 = -ftt1 fpp2 = cost**2*2*sinp*cosp drt = 8*sint*sinp*cosp -- cgit v1.2.3