diff options
Diffstat (limited to 'src/bhbrill3d.x')
| -rw-r--r-- | src/bhbrill3d.x | 126 |
1 files changed, 126 insertions, 0 deletions
diff --git a/src/bhbrill3d.x b/src/bhbrill3d.x new file mode 100644 index 0000000..01a731b --- /dev/null +++ b/src/bhbrill3d.x @@ -0,0 +1,126 @@ + o1 = deta**2 + o2 = 1/o1 + o3 = dq**2 + o4 = 1/o3 + o5 = 1/dq + o6 = tan(q) + o7 = 1/o6 + o8 = dphi**2 + o9 = 1/o8 + o10 = eta**2 + o11 = sigma**2 + o12 = 1/o11 + o13 = -(o10*o12) + o14 = -2.00000000000000d0*eta*eta0*o12 + o15 = eta0**2 + o16 = -(o12*o15) + o17 = o13 + o14 + o16 + o18 = exp(o17) + o19 = sin(q) + o20 = o19**n + o21 = 2.00000000000000d0*amp*o18*o20 + o22 = cos(phi) + o23 = o22**2 + o24 = 2.00000000000000d0*amp*c*o18*o20*o23 + o25 = 2.00000000000000d0*eta*eta0*o12 + o26 = o13 + o16 + o25 + o27 = exp(o26) + o28 = 2.00000000000000d0*amp*o20*o27 + o29 = 2.00000000000000d0*amp*c*o20*o23*o27 + o30 = o21 + o24 + o28 + o29 + o31 = exp(o30) + o32 = o19**2 + o33 = 1/o32 + o34 = 1.00000000000000d0*o31*o33*o9 + o35 = 1/dphi + o36 = o13 + o14 + o16 + o21 + o24 + o28 + o29 + o37 = exp(o36) + o38 = sin(phi) + o39 = -2.00000000000000d0 + n + o40 = o19**o39 + o41 = o13 + o16 + o21 + o24 + o25 + o28 + o29 + o42 = exp(o41) + o43 = cos(q) + o44 = o43**2 + o45 = n**2 + o46 = o38**2 + o47 = o11**2 + o48 = 1/o47 + o49 = amp**2 + o50 = c**2 + o51 = -2.00000000000000d0*o10*o12 + o52 = -2.00000000000000d0*o12*o15 + o53 = o21 + o24 + o28 + o29 + o51 + o52 + o54 = exp(o53) + o55 = 2.00000000000000d0*n + o56 = -2.00000000000000d0 + o55 + o57 = o19**o56 + o58 = -4.0000000000000d0*eta*eta0*o12 + o59 = o21 + o24 + o28 + o29 + o51 + o52 + o58 + o60 = exp(o59) + o61 = 4.0000000000000d0*eta*eta0*o12 + o62 = o21 + o24 + o28 + o29 + o51 + o52 + o61 + o63 = exp(o62) + o64 = o6**2 + o65 = 1/o64 + o66 = 5.0000000000000d-1*eta + o67 = cosh(o66) + Ae(i,j,k) = o2 + Aw(i,j,k) = o2 + An(i,j,k) = o4 + 5.0000000000000d-1*o5*o7 + As(i,j,k) = o4 - 5.0000000000000d-1*o5*o7 + Ab(i,j,k) = o34 - 2.00000000000000d0*amp*c*o22*o35*o37*o38*o40 - + & 2.00000000000000d0*amp*c*o22*o35*o38*o40*o42 + Aq(i,j,k) = o34 + 2.00000000000000d0*amp*c*o22*o35*o37*o38*o40 + + & 2.00000000000000d0*amp*c*o22*o35*o38*o40*o42 + Ac(i,j,k) = -1.25000000000000d-1 - 2.00000000000000d0*o2 - 2.500 + & 00000000000d-1*amp*n*o18*o20 - 5.0000000000000d-1*amp*o12*o18*o2 + & 0 - 2.50000000000000d-1*amp*c*n*o18*o20*o23 - 5.0000000000000d-1 + & *amp*c*o12*o18*o20*o23 - 2.50000000000000d-1*amp*n*o20*o27 - 5.0 + & 000000000000d-1*amp*o12*o20*o27 - 2.50000000000000d-1*amp*c*n*o2 + & 0*o23*o27 - 5.0000000000000d-1*amp*c*o12*o20*o23*o27 - 1.2500000 + & 0000000d-1*o33 - 2.00000000000000d0*o4 - amp*c*o23*o37*o40 - amp + & *c*o23*o40*o42 - 2.50000000000000d-1*amp*n*o18*o40*o44 - 2.50000 + & 000000000d-1*amp*c*n*o18*o23*o40*o44 - 2.50000000000000d-1*amp*n + & *o27*o40*o44 - 2.50000000000000d-1*amp*c*n*o23*o27*o40*o44 + 2.5 + & 0000000000000d-1*amp*o18*o40*o44*o45 + 2.50000000000000d-1*amp*c + & *o18*o23*o40*o44*o45 + 2.50000000000000d-1*amp*o27*o40*o44*o45 + + & 2.50000000000000d-1*amp*c*o23*o27*o40*o44*o45 + amp*c*o37*o40*o + & 46 + amp*c*o40*o42*o46 + 2.00000000000000d0*amp*eta*eta0*o18*o20 + & *o48 + amp*o10*o18*o20*o48 + amp*o15*o18*o20*o48 + 2.00000000000 + & 000d0*amp*c*eta*eta0*o18*o20*o23*o48 + amp*c*o10*o18*o20*o23*o48 + & + amp*c*o15*o18*o20*o23*o48 - 2.00000000000000d0*amp*eta*eta0*o + & 20*o27*o48 + amp*o10*o20*o27*o48 + amp*o15*o20*o27*o48 - 2.00000 + & 000000000d0*amp*c*eta*eta0*o20*o23*o27*o48 + amp*c*o10*o20*o23*o + & 27*o48 + amp*c*o15*o20*o23*o27*o48 + 6.0000000000000d0*o23*o46*o + & 49*o50*o54*o57 + 3.00000000000000d0*o23*o46*o49*o50*o57*o60 + 3. + & 00000000000000d0*o23*o46*o49*o50*o57*o63 + 1.25000000000000d-1*o + & 65 - 2.00000000000000d0*o31*o33*o9 + t1 = -2.50000000000000d-1*o67 + 5.0000000000000d-1*amp*n*o18*o20 + & *o67 + amp*o12*o18*o20*o67 + 5.0000000000000d-1*amp*c*n*o18*o20* + & o23*o67 + amp*c*o12*o18*o20*o23*o67 + 5.0000000000000d-1*amp*n*o + & 20*o27*o67 + amp*o12*o20*o27*o67 + 5.0000000000000d-1*amp*c*n*o2 + & 0*o23*o27*o67 + amp*c*o12*o20*o23*o27*o67 + 2.50000000000000d-1* + & o33*o67 + 2.00000000000000d0*amp*c*o23*o37*o40*o67 + 2.000000000 + & 00000d0*amp*c*o23*o40*o42*o67 + 5.0000000000000d-1*amp*n*o18*o40 + & *o44*o67 + 5.0000000000000d-1*amp*c*n*o18*o23*o40*o44*o67 + 5.00 + & 00000000000d-1*amp*n*o27*o40*o44*o67 + 5.0000000000000d-1*amp*c* + & n*o23*o27*o40*o44*o67 - 5.0000000000000d-1*amp*o18*o40*o44*o45*o + & 67 - 5.0000000000000d-1*amp*c*o18*o23*o40*o44*o45*o67 - 5.000000 + & 0000000d-1*amp*o27*o40*o44*o45*o67 + t2 = -5.0000000000000d-1*amp*c*o23*o27*o40*o44*o45*o67 - 2.00000 + & 000000000d0*amp*c*o37*o40*o46*o67 - 2.00000000000000d0*amp*c*o40 + & *o42*o46*o67 - 4.0000000000000d0*amp*eta*eta0*o18*o20*o48*o67 - + & 2.00000000000000d0*amp*o10*o18*o20*o48*o67 - 2.00000000000000d0* + & amp*o15*o18*o20*o48*o67 - 4.0000000000000d0*amp*c*eta*eta0*o18*o + & 20*o23*o48*o67 - 2.00000000000000d0*amp*c*o10*o18*o20*o23*o48*o6 + & 7 - 2.00000000000000d0*amp*c*o15*o18*o20*o23*o48*o67 + 4.0000000 + & 000000d0*amp*eta*eta0*o20*o27*o48*o67 - 2.00000000000000d0*amp*o + & 10*o20*o27*o48*o67 - 2.00000000000000d0*amp*o15*o20*o27*o48*o67 + & + 4.0000000000000d0*amp*c*eta*eta0*o20*o23*o27*o48*o67 - 2.00000 + & 000000000d0*amp*c*o10*o20*o23*o27*o48*o67 - 2.00000000000000d0*a + & mp*c*o15*o20*o23*o27*o48*o67 - 1.20000000000000d1*o23*o46*o49*o5 + & 0*o54*o57*o67 - 6.0000000000000d0*o23*o46*o49*o50*o57*o60*o67 - + & 6.0000000000000d0*o23*o46*o49*o50*o57*o63*o67 - 2.50000000000000 + & d-1*o65*o67 + Rhs(i,j,k) = t1 + t2 |