From dfb59643fb29687a3f28321a8061afc7e105d123 Mon Sep 17 00:00:00 2001
From: ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b>
Date: Wed, 9 Mar 2005 14:55:29 +0000
Subject: sorry, i missed to commit this file.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiOddBrillBH/trunk@52 b6f3ac56-194f-0410-8878-cdf6079d7f1b
---
 src/qfunc_even.x | 26 ++++++++++++++++++++++++++
 1 file changed, 26 insertions(+)
 create mode 100644 src/qfunc_even.x

diff --git a/src/qfunc_even.x b/src/qfunc_even.x
new file mode 100644
index 0000000..d351f25
--- /dev/null
+++ b/src/qfunc_even.x
@@ -0,0 +1,26 @@
+        o1 = -xi
+        o2 = cos(qgrd(j))
+        o3 = xi*o2
+        o4 = 1.d0 + o1 + o3
+        o5 = -eta0_me
+        o6 = etagrd(i) + o5
+        o7 = o6**2
+        o8 = sigma_me**2
+        o9 = o8**2
+        o10 = 1/o9
+        o11 = 1/o8
+        o12 = -(o11*o7)
+        o13 = exp(o12)
+        o14 = eta0_me + etagrd(i)
+        o15 = o14**2
+        o16 = -(o11*o15)
+        o17 = exp(o16)
+        o18 = sin(qgrd(j))
+        o19 = o18**n_me
+        o20 = o2**2
+        o21 = o13 + o17
+        qfetaeta(i,j) = amp_me*o19*o4*(-2.d0*o11*o13 - 2.d0*o11*o17 + 4.d0*o10*o15*o
+     &  17 + 4.d0*o10*o13*o7)
+        qfqq(i,j) = -(amp_me*xi*n_me*o19*o2*o21) - amp_me*xi*(1.d0 + n_me)*o19*o2*o21 -
+     &   amp_me*n_me*o19*o21*o4 + amp_me*(-1.d0 + n_me)*n_me*o18**(-2 + n)*o20*o21*o
+     &  4
-- 
cgit v1.2.3