This commit was generated by cvs2svn to compensate for changes in r2,

which included commits to RCS files with non-trunk default branches.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/RotatingDBHIVP/trunk@3 535fb057-194f-0410-b5a5-c63992f15602

1 files changed, 111 insertions, 0 deletions

diff --git a/src/ksq_axi.x b/src/ksq_axi.x
new file mode 100644
index 0000000..57ad7ec
--- /dev/null
+++ b/src/ksq_axi.x
@@ -0,0 +1,111 @@
+        o1 = cos(q)
+        o2 = 2.00000000000000d0*q
+        o3 = cos(o2)
+        o4 = dnnnnngtil
+        o5 = sin(q)
+        o6 = o5**2
+        o7 = o5*o6
+        o8 = dnnngtil
+        o9 = dngtil
+        o10 = dnnnnnnngtil
+        o11 = o6**2
+        o12 = 4.0000000000000d0*q
+        o13 = cos(o12)
+        o14 = dnnnngtil
+        o15 = dnngtil
+        o16 = o1**2
+        o17 = dnnnnnngtil
+        o18 = o11*o6
+        o19 = gtil
+        o20 = byJ**2
+        o21 = o14**2
+        o22 = o3**2
+        o23 = o13**2
+        o24 = o15**2
+        o25 = o4**2
+        o26 = o8**2
+        o27 = o9**2
+        o28 = o16**2
+        o29 = o17**2
+        o30 = o10**2
+        o31 = o11**2
+        o32 = o31*o6
+        o33 = o19**2
+c        exc33(i,j,k) = 0
+c        exc32(i,j,k) = 1.33333333333333d0*o1*o10*o11*o5 - 2.083333333333
+c     &  33d-2*o1*(3.5200000000000d2 - 9.2800000000000d2*o3)*o4*o7 + 1.66
+c     &  666666666667d-1*o1*(1.15000000000000d2 - 4.1500000000000d2*o3)*o
+c     &  7*o8 - 1.56250000000000d0*o1*(8.0000000000000d0 - 4.000000000000
+c     &  0d1*o3)*o7*o9
+c        exc31(i,j,k) = -8.0000000000000d0*o11*o16*o17 + 1.33333333333333
+c     &  d0*o17*o18 + 3.00000000000000d0*byJ*o6 - 2.08333333333333d-1*o15*(
+c     &  -1.11000000000000d2 - 5.8100000000000d2*o13 - 2.68000000000000d2
+c     &  *o3)*o6 - 5.2083333333333d-3*o14*(6.7200000000000d2 + 6.49600000
+c     &  00000d3*o13 + 2.04800000000000d3*o3)*o6 + o19*(5.0000000000000d1
+c     &  *o11 + 7.5000000000000d2*o11*o16 - 1.25000000000000d2*o18 - 2.00
+c     &  000000000000d2*o16*o6)
+c        exc22(i,j,k) = 0
+c        exc21(i,j,k) = 0
+c        exc11(i,j,k) = 0
+        t1 = -9.6000000000000d1*byJ*o11*o16*o17 + 1.12000000000000d2*o11*o
+     &  14*o16*o17 + 1.08266666666667d3*o11*o13*o14*o16*o17 - 7.40000000
+     &  00000d2*o11*o15*o16*o17 - 3.8733333333333d3*o11*o13*o15*o16*o17 
+     &  + 1.60000000000000d1*byJ*o17*o18 - 1.86666666666667d1*o14*o17*o18 
+     &  - 1.80444444444444d2*o13*o14*o17*o18 + 1.23333333333333d2*o15*o1
+     &  7*o18 + 6.4555555555556d2*o13*o15*o17*o18 + 6.0000000000000d2*byJ*
+     &  o11*o19 - 7.0000000000000d2*o11*o14*o19 - 6.7666666666667d3*o11*
+     &  o13*o14*o19 + 4.6250000000000d3*o11*o15*o19 + 2.42083333333333d4
+     &  *o11*o13*o15*o19 + 9.0000000000000d3*byJ*o11*o16*o19 - 1.050000000
+     &  00000d4*o11*o14*o16*o19 - 1.01500000000000d5*o11*o13*o14*o16*o19
+     &   + 6.9375000000000d4*o11*o15*o16*o19 + 3.6312500000000d5*o11*o13
+     &  *o15*o16*o19 - 1.50000000000000d3*byJ*o18*o19 + 1.75000000000000d3
+     &  *o14*o18*o19 + 1.69166666666667d4*o13*o14*o18*o19 - 1.1562500000
+     &  0000d4*o15*o18*o19 - 6.0520833333333d4*o13*o15*o18*o19 - 2.66666
+     &  666666667d3*o16*o17*o18*o19 + 1.07555555555556d2*o11*o16*o25
+        t2 = 7.4755555555556d2*o11*o16*o22*o25 + 7.3472222222222d2*o11*o
+     &  16*o26 + 9.5680555555556d3*o11*o16*o22*o26 + 3.12500000000000d2*
+     &  o11*o16*o27 + 7.8125000000000d3*o11*o16*o22*o27 + 6.400000000000
+     &  0d3*o11*o17*o19*o28 - 2.40000000000000d4*o17*o18*o19*o28 + 1.280
+     &  00000000000d2*o18*o28*o29 + 3.4133333333333d2*o11*o14*o16*o17*o3
+     &   - 1.78666666666667d3*o11*o15*o16*o17*o3 - 5.6888888888889d1*o14
+     &  *o17*o18*o3 + 2.97777777777778d2*o15*o17*o18*o3 - 2.133333333333
+     &  33d3*o11*o14*o19*o3 + 1.11666666666667d4*o11*o15*o19*o3 - 3.2000
+     &  000000000d4*o11*o14*o16*o19*o3 + 1.67500000000000d5*o11*o15*o16*
+     &  o19*o3 + 5.3333333333333d3*o14*o18*o19*o3 - 2.79166666666667d4*o
+     &  15*o18*o19*o3 - 5.6711111111111d2*o11*o16*o25*o3 - 5.30277777777
+     &  78d3*o11*o16*o26*o3 - 3.12500000000000d3*o11*o16*o27*o3 + 2.6666
+     &  6666666667d2*o17*o19*o31 + 8.0000000000000d3*o16*o17*o19*o31 - 4
+     &  .2666666666667d1*o16*o29*o31 + 3.5555555555556d0*o16*o30*o31 - 6
+     &  .6666666666667d2*o17*o19*o32 + 3.5555555555556d0*o29*o32 - 4.000
+     &  0000000000d4*o11*o16*o33
+        t3 = 5.0000000000000d3*o18*o33 + 2.50000000000000d5*o16*o18*o33 
+     &  - 6.0000000000000d5*o11*o28*o33 + 1.12500000000000d6*o18*o28*o33
+     &   - 2.50000000000000d4*o31*o33 - 3.7500000000000d5*o16*o31*o33 + 
+     &  3.12500000000000d4*o32*o33 - 3.9111111111111d1*o10*o16*o18*o4 + 
+     &  1.03111111111111d2*o10*o16*o18*o3*o4 - 4.2000000000000d1*byJ*o14*o
+     &  6 - 4.0600000000000d2*byJ*o13*o14*o6 + 2.77500000000000d2*byJ*o15*o6
+     &   + 1.45250000000000d3*byJ*o13*o15*o6 - 3.2375000000000d2*o14*o15*o
+     &  6 - 4.8241666666667d3*o13*o14*o15*o6 - 2.40000000000000d3*byJ*o16*
+     &  o19*o6 + 2.80000000000000d3*o14*o16*o19*o6 + 2.70666666666667d4*
+     &  o13*o14*o16*o19*o6 - 1.85000000000000d4*o15*o16*o19*o6 - 9.68333
+     &  33333333d4*o13*o15*o16*o19*o6 + 1.80000000000000d1*o20*o6 + 2.45
+     &  000000000000d1*o21*o6 + 4.7366666666667d2*o13*o21*o6 - 2.3822222
+     &  2222222d3*o14*o15*o22*o6 + 2.27555555555556d2*o21*o22*o6 - 1.638
+     &  09722222222d4*o14*o15*o23*o6 + 2.28938888888889d3*o21*o23*o6
+        t4 = 1.06953125000000d3*o24*o6 + 1.11963541666667d4*o13*o24*o6 +
+     &   6.2347222222222d3*o22*o24*o6 + 2.93021701388889d4*o23*o24*o6 - 
+     &  1.28000000000000d2*byJ*o14*o3*o6 + 6.7000000000000d2*byJ*o15*o3*o6 -
+     &   1.76833333333333d3*o14*o15*o3*o6 - 1.27205555555556d4*o13*o14*o
+     &  15*o3*o6 + 8.5333333333333d3*o14*o16*o19*o3*o6 - 4.4666666666667
+     &  d4*o15*o16*o19*o3*o6 + 1.49333333333333d2*o21*o3*o6 + 1.44355555
+     &  555556d3*o13*o21*o3*o6 + 5.1645833333333d3*o24*o3*o6 + 2.7032638
+     &  8888889d4*o13*o24*o3*o6 + 8.0000000000000d4*o28*o33*o6 + 1.02222
+     &  222222222d2*o10*o16*o18*o8 - 3.6888888888889d2*o10*o16*o18*o3*o8
+     &   - 5.6222222222222d2*o11*o16*o4*o8 - 5.3488888888889d3*o11*o16*o
+     &  22*o4*o8 + 3.5111111111111d3*o11*o16*o3*o4*o8 - 6.6666666666667d
+     &  1*o10*o16*o18*o9 + 3.3333333333333d2*o10*o16*o18*o3*o9 + 3.66666
+     &  66666667d2*o11*o16*o4*o9 + 4.8333333333333d3*o11*o16*o22*o4*o9 -
+     &   2.80000000000000d3*o11*o16*o3*o4*o9 - 9.5833333333333d2*o11*o16
+     &  *o8*o9 - 1.72916666666667d4*o11*o16*o22*o8*o9 + 8.2500000000000d
+     &  3*o11*o16*o3*o8*o9
+        ksq(i,j,k) = t1 + t2 + t3 + t4