diff options
| author | ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b> | 2005-03-09 14:55:29 +0000 |
|---|---|---|
| committer | ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b> | 2005-03-09 14:55:29 +0000 |
| commit | dfb59643fb29687a3f28321a8061afc7e105d123 (patch) | |
| tree | 394d441f48bd38d0f7ff0ad3e849b272a627aabe | |
| parent | 59c381790dc39733e20107641dca2b772ca49080 (diff) | |
sorry, i missed to commit this file.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiOddBrillBH/trunk@52 b6f3ac56-194f-0410-8878-cdf6079d7f1b
| -rw-r--r-- | src/qfunc_even.x | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/src/qfunc_even.x b/src/qfunc_even.x new file mode 100644 index 0000000..d351f25 --- /dev/null +++ b/src/qfunc_even.x @@ -0,0 +1,26 @@ + o1 = -xi + o2 = cos(qgrd(j)) + o3 = xi*o2 + o4 = 1.d0 + o1 + o3 + o5 = -eta0_me + o6 = etagrd(i) + o5 + o7 = o6**2 + o8 = sigma_me**2 + o9 = o8**2 + o10 = 1/o9 + o11 = 1/o8 + o12 = -(o11*o7) + o13 = exp(o12) + o14 = eta0_me + etagrd(i) + o15 = o14**2 + o16 = -(o11*o15) + o17 = exp(o16) + o18 = sin(qgrd(j)) + o19 = o18**n_me + o20 = o2**2 + o21 = o13 + o17 + qfetaeta(i,j) = amp_me*o19*o4*(-2.d0*o11*o13 - 2.d0*o11*o17 + 4.d0*o10*o15*o + & 17 + 4.d0*o10*o13*o7) + qfqq(i,j) = -(amp_me*xi*n_me*o19*o2*o21) - amp_me*xi*(1.d0 + n_me)*o19*o2*o21 - + & amp_me*n_me*o19*o21*o4 + amp_me*(-1.d0 + n_me)*n_me*o18**(-2 + n)*o20*o21*o + & 4 |