source: CIVL/examples/fortran/cv2018/equivalence/harmonic/c_function.c

main
Last change on this file was ea777aa, checked in by Alex Wilton <awilton@…>, 3 years ago

Moved examples, include, build_default.properties, common.xml, and README out from dev.civl.com into the root of the repo.

git-svn-id: ​svn://vsl.cis.udel.edu/civl/trunk@5704 fb995dde-84ed-4084-dfe6-e5aef3e2452c
  • Property mode set to 100644
File size: 3.9 KB
RevLine 
[d5096aa]1#include <math.h>
2
3void 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}
Note: See TracBrowser for help on using the repository browser.