aboutsummaryrefslogtreecommitdiff
path: root/src/psi_1st_deriv.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/psi_1st_deriv.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/psi_1st_deriv.x')
-rw-r--r--src/psi_1st_deriv.x20
1 files changed, 20 insertions, 0 deletions
diff --git a/src/psi_1st_deriv.x b/src/psi_1st_deriv.x
new file mode 100644
index 0000000..c7725c5
--- /dev/null
+++ b/src/psi_1st_deriv.x
@@ -0,0 +1,20 @@
+ o1 = 5.0000000000000d-1*eta(i,j,k)
+ o2 = exp(o1)
+ o3 = psi3d(i,j,k)
+ o4 = 1/o3
+ o5 = cos(phi(i,j,k))
+ o6 = cos(q(i,j,k))
+ o7 = dqpsi3d(i,j,k)
+ o8 = -1.50000000000000d0*eta(i,j,k)
+ o9 = exp(o8)
+ o10 = dphipsi3d(i,j,k)
+ o11 = sin(phi(i,j,k))
+ o12 = sin(q(i,j,K))
+ o13 = 1/o12
+ o14 = detapsi3d(i,j,k)
+ psix(i,j,k) = o2*o4*(-(o10*o11*o13*o9) + o12*o14*o5*o9 - 5.00000
+ & 00000000d-1*o12*o3*o5*o9 + o5*o6*o7*o9)
+ psiy(i,j,k) = o2*o4*(o11*o12*o14*o9 - 5.0000000000000d-1*o11*o12
+ & *o3*o9 + 1.00000000000000d0*o10*o13*o5*o9 + o11*o6*o7*o9)
+ psiz(i,j,k) = o2*o4*(o14*o6*o9 - 5.0000000000000d-1*o3*o6*o9 - o
+ & 12*o7*o9)