From 266e41ac2faa66e868afb82d510f58fe551c4d81 Mon Sep 17 00:00:00 2001
From: ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b>
Date: Sun, 9 Sep 2001 11:47:21 +0000
Subject: I forgot add following files ...

But it looks like there is bug ... I'm working on it.

ryoji


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiOddBrillBH/trunk@34 b6f3ac56-194f-0410-8878-cdf6079d7f1b
---
 src/ckerr.x            |  4 ++++
 src/kerr.x             |  4 ++++
 src/kij_sergio.x       | 48 ++++++++++++++++++++++++++++++++++++++++++++++++
 src/ksq_sergio.x       | 49 +++++++++++++++++++++++++++++++++++++++++++++++++
 test/test_axioddbh.par |  2 +-
 5 files changed, 106 insertions(+), 1 deletion(-)
 create mode 100644 src/ckerr.x
 create mode 100644 src/kerr.x
 create mode 100644 src/kij_sergio.x
 create mode 100644 src/ksq_sergio.x

diff --git a/src/ckerr.x b/src/ckerr.x
new file mode 100644
index 0000000..b6dbed5
--- /dev/null
+++ b/src/ckerr.x
@@ -0,0 +1,4 @@
+        mass = sqrt(2+sqrt(4+byJ**2))
+        a = byJ/mass
+        rBL = sqrt(mass**2-a**2)/2*exp(eta(i,j,k))*(1+(mass+a)/(sqrt(mass**2-a**2
+     &  )*exp(eta(i,j,k))))*(1+(mass-a)/(sqrt(mass**2-a**2)*exp(eta(i,j,k))))
diff --git a/src/kerr.x b/src/kerr.x
new file mode 100644
index 0000000..57f01b8
--- /dev/null
+++ b/src/kerr.x
@@ -0,0 +1,4 @@
+        mass = sqrt(2+sqrt(4+byJ**2))
+        a = byJ/mass
+        rBL = sqrt(mass**2-a**2)/2*exp(etagrd(i))*(1+(mass+a)/(sqrt(mass**2-a**2
+     &  )*exp(etagrd(i))))*(1+(mass-a)/(sqrt(mass**2-a**2)*exp(etagrd(i))))
diff --git a/src/kij_sergio.x b/src/kij_sergio.x
new file mode 100644
index 0000000..bb73eb5
--- /dev/null
+++ b/src/kij_sergio.x
@@ -0,0 +1,48 @@
+        o1 = a**2
+        o2 = a*o1
+        o3 = -2.00000000000000d0*mass*rBL
+        o4 = rBL**2
+        o5 = o1 + o3 + o4
+        o6 = sqrt(o5)
+        o7 = cos(phi(i,j,k))
+        o8 = cos(q(i,j,k))
+        o9 = o8**2
+        o10 = o1*o9
+        o11 = o10 + o4
+        o12 = o11**2
+        o13 = 1/o12
+        o14 = -3.00000000000000d0*eta(i,j,k)
+        o15 = exp(o14)
+        o16 = sin(phi(i,j,k))
+        o17 = sin(q(i,j,k))
+        o18 = o17**2
+        o19 = -o1
+        o20 = mass**2
+        o21 = o19 + o20
+        o22 = sqrt(o21)
+        o23 = 1/o4
+        o24 = 2.50000000000000d-1*o21*o23
+        o25 = 1.00000000000000d0 + o24
+        o26 = o1 + o4
+        o27 = 2.00000000000000d0*o26*o4
+        o28 = o19 + o4
+        o29 = o11*o28
+        o30 = o27 + o29
+        o31 = -2.00000000000000d0*eta(i,j,k)
+        o32 = exp(o31)
+        o33 = o7**2
+        o34 = o16**2
+        o35 = o17*o18
+        kxx(i,j,k) = -(a*mass*o13*o16*o18*o22*o25*o30*o32*o7) + 4.000000
+     &  0000000d0*mass*o13*o15*o16*o18*o2*o6*o7*o9*rBL
+        kxy(i,j,k) = 5.0000000000000d-1*a*mass*o13*o18*o22*o25*o30*o32*o
+     &  33 - 5.0000000000000d-1*a*mass*o13*o18*o22*o25*o30*o32*o34 - 2.0
+     &  0000000000000d0*mass*o13*o15*o18*o2*o33*o6*o9*rBL + 2.0000000000
+     &  0000d0*mass*o13*o15*o18*o2*o34*o6*o9*rBL
+        kxz(i,j,k) = -5.0000000000000d-1*a*mass*o13*o16*o17*o22*o25*o30*
+     &  o32*o8 - 2.00000000000000d0*mass*o13*o15*o16*o2*o35*o6*o8*rBL
+        kyy(i,j,k) = a*mass*o13*o16*o18*o22*o25*o30*o32*o7 - 4.000000000
+     &  0000d0*mass*o13*o15*o16*o18*o2*o6*o7*o9*rBL
+        kyz(i,j,k) = 5.0000000000000d-1*a*mass*o13*o17*o22*o25*o30*o32*o
+     &  7*o8 + 2.00000000000000d0*mass*o13*o15*o2*o35*o6*o7*o8*rBL
+        kzz(i,j,k) = 0
diff --git a/src/ksq_sergio.x b/src/ksq_sergio.x
new file mode 100644
index 0000000..8a6180b
--- /dev/null
+++ b/src/ksq_sergio.x
@@ -0,0 +1,49 @@
+        o1 = a**2
+        o2 = -2.00000000000000d0*mass*rBL
+        o3 = rBL**2
+        o4 = o1 + o2 + o3
+        o5 = sqrt(o4)
+        o6 = cos(qgrd(j))
+        o7 = o6**2
+        o8 = o1*o7
+        o9 = o3 + o8
+        o10 = o9**2
+        o11 = 1/o10
+        o12 = 1.0!psi2dv(eta,q)
+        o13 = o12**2
+        o14 = 1/o13
+        o15 = sin(qgrd(j))
+        o16 = o15**2
+        o17 = -o1
+        o18 = mass**2
+        o19 = o17 + o18
+        o20 = sqrt(o19)
+        o21 = 1/o3
+        o22 = 2.50000000000000d-1*o19*o21
+        o23 = 1.00000000000000d0 + o22
+        o24 = o1 + o3
+        o25 = 2.00000000000000d0*o24*o3
+        o26 = o17 + o3
+        o27 = o26*o9
+        o28 = o25 + o27
+        o29 = o23**2
+        o30 = o10**2
+        o31 = 1/o30
+        o32 = o28**2
+        o33 = o13**2
+        o34 = o33**2
+        o35 = o33*o34
+        o36 = 1/o35
+        o37 = o1**2
+        o38 = o16**2
+c        exc33(i,j) = 0
+c        exc32(i,j) = -2.00000000000000d0*a*mass*o1*o11*o14*o15*o16*o5*o6
+c     &  *rBL
+c        exc31(i,j) = 5.0000000000000d-1*a*mass*o11*o14*o16*o20*o23*o28*e
+c     &  xp(eta)
+c        exc22(i,j) = 0
+c        exc21(i,j) = 0
+c        exc11(i,j) = 0
+        ksq(i,j) = 8.0000000000000d0*o1*o18*o3*o31*o36*o37*o38*o4*o7 + 5
+     &  .0000000000000d-1*o1*o16*o18*o19*o29*o31*o32*o36*exp(2.000000000
+     &  00000d0*etagrd(i))
diff --git a/test/test_axioddbh.par b/test/test_axioddbh.par
index 7090fb5..b58374a 100755
--- a/test/test_axioddbh.par
+++ b/test/test_axioddbh.par
@@ -2,7 +2,7 @@
 #!DESC "Initial data for axisymmetric distored black hole"
 ######################################################################
 
-ActiveThorns = "time ADMconstraints pugh pughreduce pughslab pughinterp cartgrid3d einstein ADM IDAxiOddBrillBH ioascii iobasic ioutil"
+ActiveThorns = "time boundary ADMconstraints pugh pughreduce pughslab pughinterp cartgrid3d einstein ADM IDAxiOddBrillBH ioascii iobasic ioutil"
 
 # GENERAL 
 
-- 
cgit v1.2.3