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/IDAxiOddBrillBH/trunk@3 b6f3ac56-194f-0410-8878-cdf6079d7f1b

1 files changed, 95 insertions, 0 deletions

diff --git a/src/ksq_axi_rdbh.x b/src/ksq_axi_rdbh.x
new file mode 100644
index 0000000..a590721
--- /dev/null
+++ b/src/ksq_axi_rdbh.x
@@ -0,0 +1,95 @@
+        o1 = cos(qgrd(j))
+        o2 = 2.00000000000000d0*qgrd(j)
+        o3 = cos(o2)
+        o4 = dnnnnngtil
+        o5 = sin(qgrd(j))
+        o6 = o5**2
+        o7 = o5*o6
+        o8 = dnnngtil
+        o9 = dngtil
+        o10 = dnnnnnnngtil
+        o11 = o6**2
+        o12 = 4.0000000000000d0*qgrd(j)
+        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
+        t1 = 2.88000000000000d1*byJ*o11*o16*o17 + 1.00800000000000d1*o11*o
+     &  14*o16*o17 + 9.7440000000000d1*o11*o13*o14*o16*o17 - 6.660000000
+     &  0000d1*o11*o15*o16*o17 - 3.4860000000000d2*o11*o13*o15*o16*o17 -
+     &   4.8000000000000d0*byJ*o17*o18 - 1.68000000000000d0*o14*o17*o18 - 
+     &  1.62400000000000d1*o13*o14*o17*o18 + 1.11000000000000d1*o15*o17*
+     &  o18 + 5.8100000000000d1*o13*o15*o17*o18 - 1.80000000000000d2*byJ*o
+     &  11*o19 - 6.3000000000000d1*o11*o14*o19 - 6.0900000000000d2*o11*o
+     &  13*o14*o19 + 4.1625000000000d2*o11*o15*o19 + 2.17875000000000d3*
+     &  o11*o13*o15*o19 - 2.70000000000000d3*byJ*o11*o16*o19 - 9.450000000
+     &  0000d2*o11*o14*o16*o19 - 9.1350000000000d3*o11*o13*o14*o16*o19 +
+     &   6.2437500000000d3*o11*o15*o16*o19 + 3.2681250000000d4*o11*o13*o
+     &  15*o16*o19 + 4.5000000000000d2*byJ*o18*o19 + 1.57500000000000d2*o1
+     &  4*o18*o19 + 1.52250000000000d3*o13*o14*o18*o19 - 1.0406250000000
+     &  0d3*o15*o18*o19 - 5.4468750000000d3*o13*o15*o18*o19 - 2.40000000
+     &  000000d2*o16*o17*o18*o19 + 9.6800000000000d0*o11*o16*o25
+        t2 = 6.7280000000000d1*o11*o16*o22*o25 + 6.6125000000000d1*o11*o
+     &  16*o26 + 8.6112500000000d2*o11*o16*o22*o26 + 2.81250000000000d1*
+     &  o11*o16*o27 + 7.0312500000000d2*o11*o16*o22*o27 + 5.760000000000
+     &  0d2*o11*o17*o19*o28 - 2.16000000000000d3*o17*o18*o19*o28 + 1.152
+     &  00000000000d1*o18*o28*o29 + 3.07200000000000d1*o11*o14*o16*o17*o
+     &  3 - 1.60800000000000d2*o11*o15*o16*o17*o3 - 5.1200000000000d0*o1
+     &  4*o17*o18*o3 + 2.68000000000000d1*o15*o17*o18*o3 - 1.92000000000
+     &  000d2*o11*o14*o19*o3 + 1.00500000000000d3*o11*o15*o19*o3 - 2.880
+     &  00000000000d3*o11*o14*o16*o19*o3 + 1.50750000000000d4*o11*o15*o1
+     &  6*o19*o3 + 4.8000000000000d2*o14*o18*o19*o3 - 2.51250000000000d3
+     &  *o15*o18*o19*o3 - 5.1040000000000d1*o11*o16*o25*o3 - 4.772500000
+     &  0000d2*o11*o16*o26*o3 - 2.81250000000000d2*o11*o16*o27*o3 + 2.40
+     &  000000000000d1*o17*o19*o31 + 7.2000000000000d2*o16*o17*o19*o31 -
+     &   3.8400000000000d0*o16*o29*o31 + 3.2000000000000d-1*o16*o30*o31 
+     &  - 6.0000000000000d1*o17*o19*o32 + 3.2000000000000d-1*o29*o32 - 3
+     &  .6000000000000d3*o11*o16*o33
+        t3 = 4.5000000000000d2*o18*o33 + 2.25000000000000d4*o16*o18*o33 
+     &  - 5.4000000000000d4*o11*o28*o33 + 1.01250000000000d5*o18*o28*o33
+     &   - 2.25000000000000d3*o31*o33 - 3.3750000000000d4*o16*o31*o33 + 
+     &  2.81250000000000d3*o32*o33 - 3.5200000000000d0*o10*o16*o18*o4 + 
+     &  9.2800000000000d0*o10*o16*o18*o3*o4 + 1.26000000000000d1*byJ*o14*o
+     &  6 + 1.21800000000000d2*byJ*o13*o14*o6 - 8.3250000000000d1*byJ*o15*o6
+     &   - 4.3575000000000d2*byJ*o13*o15*o6 - 2.91375000000000d1*o14*o15*o
+     &  6 - 4.3417500000000d2*o13*o14*o15*o6 + 7.2000000000000d2*byJ*o16*o
+     &  19*o6 + 2.52000000000000d2*o14*o16*o19*o6 + 2.43600000000000d3*o
+     &  13*o14*o16*o19*o6 - 1.66500000000000d3*o15*o16*o19*o6 - 8.715000
+     &  0000000d3*o13*o15*o16*o19*o6 + 1.80000000000000d1*o20*o6 + 2.205
+     &  00000000000d0*o21*o6 + 4.2630000000000d1*o13*o21*o6 - 2.14400000
+     &  000000d2*o14*o15*o22*o6 + 2.04800000000000d1*o21*o22*o6 - 1.4742
+     &  8750000000d3*o14*o15*o23*o6 + 2.06045000000000d2*o21*o23*o6
+        t4 = 9.6257812500000d1*o24*o6 + 1.00767187500000d3*o13*o24*o6 + 
+     &  5.6112500000000d2*o22*o24*o6 + 2.63719531250000d3*o23*o24*o6 + 3
+     &  .8400000000000d1*byJ*o14*o3*o6 - 2.01000000000000d2*byJ*o15*o3*o6 - 
+     &  1.59150000000000d2*o14*o15*o3*o6 - 1.14485000000000d3*o13*o14*o1
+     &  5*o3*o6 + 7.6800000000000d2*o14*o16*o19*o3*o6 - 4.0200000000000d
+     &  3*o15*o16*o19*o3*o6 + 1.34400000000000d1*o21*o3*o6 + 1.299200000
+     &  00000d2*o13*o21*o3*o6 + 4.6481250000000d2*o24*o3*o6 + 2.43293750
+     &  000000d3*o13*o24*o3*o6 + 7.2000000000000d3*o28*o33*o6 + 9.200000
+     &  0000000d0*o10*o16*o18*o8 - 3.3200000000000d1*o10*o16*o18*o3*o8 -
+     &   5.0600000000000d1*o11*o16*o4*o8 - 4.8140000000000d2*o11*o16*o22
+     &  *o4*o8 + 3.16000000000000d2*o11*o16*o3*o4*o8 - 6.0000000000000d0
+     &  *o10*o16*o18*o9 + 3.00000000000000d1*o10*o16*o18*o3*o9 + 3.30000
+     &  00000000d1*o11*o16*o4*o9 + 4.3500000000000d2*o11*o16*o22*o4*o9 -
+     &   2.52000000000000d2*o11*o16*o3*o4*o9 - 8.6250000000000d1*o11*o16
+     &  *o8*o9 - 1.55625000000000d3*o11*o16*o22*o8*o9 + 7.4250000000000d
+     &  2*o11*o16*o3*o8*o9
+        ksq(i,j) = t1 + t2 + t3 + t4