diff options
| author | rhaas <rhaas@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2014-01-13 17:16:03 +0000 |
|---|---|---|
| committer | rhaas <rhaas@5c0f84ea-6048-4d6e-bfa4-55cd5f2e0dd7> | 2014-01-13 17:16:03 +0000 |
| commit | 1e7726df7ee9e189088a3869de2364e839003c45 (patch) | |
| tree | 362bea531943fd62c52078ac79f3da4b8b01e8cc | |
| parent | 202a2eb9ed7ac09b291fce2df16ba9d0377e38a1 (diff) | |
add m=2 ftp coefficient for teukolsky wavessvn
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
| -rw-r--r-- | src/teukwaves.F77 | 3 |
1 files changed, 2 insertions, 1 deletions
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 |