aboutsummaryrefslogtreecommitdiff
path: root/src/psi_1st_deriv.x
diff options
context:
space:
mode:
authorallen <allen@0a4070d5-58f5-498f-b6c0-2693e757fa0f>1999-11-01 11:27:16 +0000
committerallen <allen@0a4070d5-58f5-498f-b6c0-2693e757fa0f>1999-11-01 11:27:16 +0000
commitcab79658b8e19947c3b4058f75961f01a055ef2e (patch)
tree0532e9dac09c72f6e1aa2942186ac65dba2b5e66 /src/psi_1st_deriv.x
parentf6ebcf616a79e8a6ed1683187dc8d58e33980c21 (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/IDAxiBrillBH/trunk@3 0a4070d5-58f5-498f-b6c0-2693e757fa0f
Diffstat (limited to 'src/psi_1st_deriv.x')
-rw-r--r--src/psi_1st_deriv.x18
1 files changed, 18 insertions, 0 deletions
diff --git a/src/psi_1st_deriv.x b/src/psi_1st_deriv.x
new file mode 100644
index 0000000..3c3335e
--- /dev/null
+++ b/src/psi_1st_deriv.x
@@ -0,0 +1,18 @@
+ o1 = 5.0000000000000d-1*eta(i,j,k)
+ o2 = exp(o1)
+ o3 = psi2dv(i,j,k)
+ o4 = 1/o3
+ o5 = cos(phi(i,j,k))
+ o6 = cos(q(i,j,k))
+ o7 = dqpsi2dv(i,j,k)
+ o8 = -1.50000000000000d0*eta(i,j,k)
+ o9 = exp(o8)
+ o10 = detapsi2dv(i,j,k)
+ o11 = sin(q(i,j,k))
+ o12 = sin(phi(i,j,k))
+ psix(i,j,k) = o2*o4*(o10*o11*o5*o9 - 5.0000000000000d-1*o11*o3*o
+ & 5*o9 + o5*o6*o7*o9)
+ psiy(i,j,k) = o2*o4*(o10*o11*o12*o9 - 5.0000000000000d-1*o11*o12
+ & *o3*o9 + o12*o6*o7*o9)
+ psiz(i,j,k) = o2*o4*(o10*o6*o9 - 5.0000000000000d-1*o3*o6*o9 - o
+ & 11*o7*o9)