aboutsummaryrefslogtreecommitdiff
path: root/src/kij_axi.x
diff options
context:
space:
mode:
authorryoji <ryoji@535fb057-194f-0410-b5a5-c63992f15602>1999-12-01 13:34:57 +0000
committerryoji <ryoji@535fb057-194f-0410-b5a5-c63992f15602>1999-12-01 13:34:57 +0000
commit5d27c1507b8402c8be8ec6433bbc5483db0203ec (patch)
treeccef534f16a8818747f1faafb1d946e9fad79d45 /src/kij_axi.x
parent4f2663df48d28c15194cde1de8e03a279534fb0b (diff)
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
Diffstat (limited to 'src/kij_axi.x')
-rw-r--r--src/kij_axi.x62
1 files changed, 62 insertions, 0 deletions
diff --git a/src/kij_axi.x b/src/kij_axi.x
new file mode 100644
index 0000000..e91a8e0
--- /dev/null
+++ b/src/kij_axi.x
@@ -0,0 +1,62 @@
+ o1 = cos(phi)
+ o2 = cos(q)
+ o3 = o2**2
+ o4 = -3.00000000000000d0*eta
+ o5 = exp(o4)
+ o6 = sin(phi)
+ o7 = sin(q)
+ o8 = o7**2
+ o9 = dnnnnngtil
+ o10 = -3.6000000000000d1*o9
+ o11 = dnnnnnnngtil
+ o12 = -4.0000000000000d0*o11*o8
+ o13 = 1.16000000000000d2*o8*o9
+ o14 = dnnngtil
+ o15 = -4.1500000000000d2*o8
+ o16 = 1.50000000000000d2 + o15
+ o17 = o14*o16
+ o18 = dngtil
+ o19 = 5.0000000000000d0*o8
+ o20 = -2.00000000000000d0 + o19
+ o21 = 7.5000000000000d1*o18*o20
+ o22 = o10 + o12 + o13 + o17 + o21
+ o23 = dnnnnnngtil
+ o24 = 3.00000000000000d0*o23
+ o25 = 2.00000000000000d0*q
+ o26 = cos(o25)
+ o27 = 4.0000000000000d0*o23*o26
+ o28 = 4.0000000000000d0*q
+ o29 = cos(o28)
+ o30 = -7.0000000000000d0*o23*o29
+ o31 = dnnnngtil
+ o32 = 2.10000000000000d1*o31
+ o33 = 6.4000000000000d1*o26*o31
+ o34 = 2.03000000000000d2*o29*o31
+ o35 = o24 + o27 + o30 + o32 + o33 + o34
+ o36 = -4.0000000000000d0*o35
+ o37 = 2.68000000000000d2*o26
+ o38 = 5.8100000000000d2*o29
+ o39 = 1.11000000000000d2 + o37 + o38
+ o40 = dnngtil
+ o41 = 5.0000000000000d0*o39*o40
+ o42 = 2.00000000000000d1*o26
+ o43 = 3.5000000000000d1*o29
+ o44 = 9.0000000000000d0 + o42 + o43
+ o45 = gtil
+ o46 = -7.5000000000000d1*o44*o45
+ o47 = o36 + o41 + o46
+ o48 = o1**2
+ o49 = o6**2
+ o50 = o7*o8
+ hxx(i,j,k) = 6.6666666666667d-1*o1*o22*o3*o5*o6*o8 - 8.333333333
+ & 3333d-2*o1*o47*o5*o6*o8
+ hxy(i,j,k) = -3.3333333333333d-1*o22*o3*o48*o5*o8 + 4.1666666666
+ & 667d-2*o47*o48*o5*o8 + 3.3333333333333d-1*o22*o3*o49*o5*o8 - 4.1
+ & 666666666667d-2*o47*o49*o5*o8
+ hxz(i,j,k) = -3.3333333333333d-1*o2*o22*o5*o50*o6 - 4.1666666666
+ & 667d-2*o2*o47*o5*o6*o7
+ hyy(i,j,k) = -6.6666666666667d-1*o1*o22*o3*o5*o6*o8 + 8.33333333
+ & 33333d-2*o1*o47*o5*o6*o8
+ hyz(i,j,k) = 3.3333333333333d-1*o1*o2*o22*o5*o50 + 4.16666666666
+ & 67d-2*o1*o2*o47*o5*o7
+ hzz(i,j,k) = 0