| 1 | #include <math.h>
|
|---|
| 2 |
|
|---|
| 3 | void F_HARMONIC(double _alpha, double _beta, double _gamma, double _beta2, double _beta3, double *_f) {
|
|---|
| 4 | double t1,t2,t3,t4,t5,t6,t7,t8,t9;
|
|---|
| 5 | double t10,t11,t12,t13,t14,t15,t16,t17,t18,t19;
|
|---|
| 6 | double t20,t21,t22,t23,t24,t25,t26,t27,t28,t29;
|
|---|
| 7 | double t30,t31,t32,t33,t34,t35,t36,t37,t38,t39;
|
|---|
| 8 | double t40,t41,t42,t43,t44,t45,t46,t47,t48,t49;
|
|---|
| 9 | double t50,t51,t52,t53,t54,t55;
|
|---|
| 10 | double t11_2,t31_2,t41_2,t30_2,t37_2,t44_2,t10_2,t49_2,t28_2,t44_3,t46_2,t53_2,t55_2;
|
|---|
| 11 | double _f_0,_f_1,_f_2,_f_3,_f_4,_f_5,_f_6,_f_7,_f_8;
|
|---|
| 12 |
|
|---|
| 13 | t1 = cos(_beta);
|
|---|
| 14 | t2 = cos(_alpha);
|
|---|
| 15 | t3 = sin(_beta2);
|
|---|
| 16 | t4 = t2 - t3;
|
|---|
| 17 | t5 = t2 + t3;
|
|---|
| 18 | t6 = t1*t1;
|
|---|
| 19 | t7 = t6 + 0.1e1*_beta3;
|
|---|
| 20 | t8 = cos(_gamma);
|
|---|
| 21 | t9 = sin(_gamma);
|
|---|
| 22 | t10 = t2*t2;
|
|---|
| 23 | t11 = 0.2e1 * t2;
|
|---|
| 24 | t12 = t10 + (-t11 - t3) * t3;
|
|---|
| 25 | t11_2 = t10 + (t11 - t3) * t3;
|
|---|
| 26 | t13 = 0.1e1 + (0.6e1 + t6) * t6;
|
|---|
| 27 | t14 = t9*t9;
|
|---|
| 28 | t15 = t8*t8;
|
|---|
| 29 | t16 = t15*t15 + t14*t14;
|
|---|
| 30 | t17 = t15 - t14;
|
|---|
| 31 | t18 = sqrt(0.5e1);
|
|---|
| 32 | t19 = sin(_beta);
|
|---|
| 33 | t20 = sqrt(0.3e1);
|
|---|
| 34 | t21 = sqrt(0.7e1);
|
|---|
| 35 | t22 = sqrt(0.14e2);
|
|---|
| 36 | t23 = sqrt(0.30e2);
|
|---|
| 37 | t24 = t8 - t9;
|
|---|
| 38 | t25 = t8 + t9;
|
|---|
| 39 | t26 = t19*t19;
|
|---|
| 40 | t27 = t26*t26;
|
|---|
| 41 | t28 = sqrt(0.2e1);
|
|---|
| 42 | t29 = sqrt(0.6e1);
|
|---|
| 43 | t30 = sqrt(0.42e2);
|
|---|
| 44 | t31 = 0.2e1 * t23;
|
|---|
| 45 | t32 = t18 * (0.4e1 * t20 + t28 * t29) * t30 + t31 * t20 * t21;
|
|---|
| 46 | t33 = t15 - t14 / 0.3e1;
|
|---|
| 47 | t34 = t15 - 0.3e1 * t14;
|
|---|
| 48 | t35 = t6 + 0.3e1;
|
|---|
| 49 | t36 = t12 * t11_2;
|
|---|
| 50 | t37 = t29 * t20;
|
|---|
| 51 | t38 = t28 + t37;
|
|---|
| 52 | t39 = -0.3e1 / 0.16e2;
|
|---|
| 53 | t40 = t15 + t14;
|
|---|
| 54 | t31_2 = t20 * (t6 * t18 * t29 + t31 * (t6 - t26 * t40 / 0.2e1)) * t28 + (0.12e2 * t18 + t29 * t23) * t6;
|
|---|
| 55 | t41 = t40 * t6;
|
|---|
| 56 | t42 = t41 - t15 - t14 - 0.2e1 * t26;
|
|---|
| 57 | t43 = sqrt(0.10e2);
|
|---|
| 58 | t44 = sqrt(0.15e2);
|
|---|
| 59 | t45 = t20 * t28;
|
|---|
| 60 | t46 = t45 * t42;
|
|---|
| 61 | t41_2 = t46 * t43 + (t41 - t15 - t14 - 0.8e1 * t26) * t44;
|
|---|
| 62 | t47 = t1 - 0.1e1;
|
|---|
| 63 | t48 = t1 + 0.1e1;
|
|---|
| 64 | t30_2 = _alpha*t41_2 * sqrt(0.105e3) - 0.30e2 * t30 * t26 * t29 + 0.15e2 * t40 * t47 * t48 * t21;
|
|---|
| 65 | t49 = _beta*t3*t3;
|
|---|
| 66 | t50 = t10*t10 + (-0.6e1 * t10 + t49) * t49*_beta2;
|
|---|
| 67 | t51 = 0.4e1 * t2 * t3;
|
|---|
| 68 | t37_2 = (t45 * (t6 - 0.3e1 / 0.8e1 * t26 * t40) + t29 * t6 / 0.2e1) * t44 + (-t37 * t26 * t40 / 0.2e1 + 0.3e1 / 0.2e1 * t28 * (t6 - 0.5e1 / 0.4e1 * t26 * t40)) * t18;
|
|---|
| 69 | t44_2 = sqrt(0.70e2);
|
|---|
| 70 | t52 = sqrt(0.21e2);
|
|---|
| 71 | t10_2 = t51 * (t10 - t49)*_beta3;
|
|---|
| 72 | t49_2 = t50 * t1;
|
|---|
| 73 | t53 = t18 / 0.1680e4*_beta2;
|
|---|
| 74 | t28_2 = t53 * ((0.6e1 * t28 * t44_2 * t43 + 0.15e2 * t52 * t29) * t48 * t47 * t40 + t46 * t44_2 * t23 + t41_2 * sqrt(0.210e3));
|
|---|
| 75 | t44_3 = t1 * t21 / 0.6e1*_beta*_gamma;
|
|---|
| 76 | t46_2 = -t26 * t40 / 0.4e1;
|
|---|
| 77 | t53_2 = t26 * t25 * t24 * t31_2 * t21 / 0.24e2 - t53 * t30_2 * (t36 * t7 * t25 * t24 - 0.16e2 * t2 * t8 * t3 * t9 * t4 * t5 * t1);
|
|---|
| 78 | t54 = t19 * (t1 * t32 * t8 * t26 * t34 * t21 / 0.96e2 - t38 * ((-0.36e2 * t6 - 0.12e2) * t33 * t9 * t5 * t4 * t3 * t2 + t36 * t8 * t34 * t1 * t35) * t18 / 0.16e2);
|
|---|
| 79 | t55 = 0.2e1 * t8;
|
|---|
| 80 | t55_2 = (-0.384e3 * t8 * t9 * t17 * t7 * t5 * t4 * t3 * t2 * t1 + t36 * t13 * (0.6e1 * t16 - 0.36e2 * t14 * t15)) * t18 / 0.48e2 + t27 * t20 * t21 * t22 * t23 * (t15 + (-t55 - t9) * t9) * (t15 + (t55 - t9) * t9) / 0.48e2;
|
|---|
| 81 |
|
|---|
| 82 | _f_0 = ((0.24e2 * t16 - 0.144e3 * t14 * t15) * t7 * t5 * t4 * t3 * t2 * t1 + 0.6e1 * t9 * t12 * t11_2 * t13 * t8 * t17) * t18 / 0.12e2 + t27 * t8 * t9 * t20 * t21 * t22 * t23 * t24 * t25 / 0.12e2;
|
|---|
| 83 | _f_1 = t19 * (t9 * t1 * t26 * t32 * t33 * t21 / 0.32e2 + t39 * ((0.4e1 * t6 + 0.4e1 / 0.3e1) * t34 * t8 * t5 * t4 * t3 * t2 + t36 * t9 * t33 * t1 * t35) * t38 * t18);
|
|---|
| 84 | _f_2 = t21 * t9 * t31_2 * t26 * t8 / 0.12e2 - t30_2 * (t7 * t50 * t9 * t8 + t51 * t24 * t25 * t4 * t5 * t1) * t18 / 0.840e3;
|
|---|
| 85 | _f_3 = t19 * t44_3 * t9 * t37_2 + t28_2 * (t49_2 * t9 + t10_2 * t8);
|
|---|
| 86 | _f_4 = -t36 * (-t40 * t47 * t41_2 * t48 * t52 / 0.2e1 + t26 * (t23 * t42 * t20 + 0.6e1 * t43 * t40 * t47 * t48) * t22) * t18 / 0.420e3 + t21 * (t46_2 * t45 * t29 * (t6 + t46_2) + (t6 - 0.3e1 / 0.2e1 * t26 * t40) * t6);
|
|---|
| 87 | _f_5 = t19 * (t44_3 * t37_2 * t8 + t28_2 * (t49_2 * t8 - t10_2 * t9));
|
|---|
| 88 | _f_6 = t53_2;
|
|---|
| 89 | _f_7 = t54*1;
|
|---|
| 90 | _f_8 = t55_2;
|
|---|
| 91 |
|
|---|
| 92 | *_f = _f_0 + _f_1 + _f_2 + _f_3 + _f_4 + _f_5 + _f_6 + _f_7 + _f_8;
|
|---|
| 93 | } |
|---|
