aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dzz.compute.c
blob: 49a4403c56b2f17d2dcd53928fcb9c3979015e7e (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
fp t804;
fp t782;
fp t803;
fp t784;
fp t802;
fp t788;
fp t801;
fp t787;
fp t800;
fp t778;
fp t779;
fp t786;
fp t781;
fp t764;
fp t799;
fp t783;
fp t798;
fp t774;
fp t797;
fp t796;
fp t795;
fp t794;
fp t793;
fp t792;
fp t791;
fp t790;
fp t789;
fp t776;
fp t775;
fp t772;
fp t771;
fp t769;
fp t768;
fp t766;
fp t765;
fp t763;
fp t762;
fp t761;
fp t760;
fp t759;
fp t758;
fp t757;
fp t756;
fp t755;
      t804 = x*y;
      t782 = RATIONAL(1.0,2.0);
      t803 = t782*y;
      t784 = RATIONAL(-1.0,2.0);
      t802 = t784*y;
      t788 = x*x;
      t801 = t788*y;
      t787 = y*y;
      t800 = x*t787;
      t778 = t782*t787;
      t779 = t782*t788;
      t786 = t788*t787;
      t781 = RATIONAL(1.0,4.0);
      t764 = t781*t786;
      t799 = t764+t781*t804;
      t783 = RATIONAL(-1.0,4.0);
      t798 = t764+t783*t804;
      t774 = t784*t786;
      t797 = x*t803+t774;
      t796 = t774+x*t802;
      t795 = t774+t778;
      t794 = t774+t779;
      t793 = -t787+t786;
      t792 = -t788+t786;
      t791 = t788+t787;
      t790 = t801-y;
      t789 = t800-x;
      t776 = t784*t800;
      t775 = t784*t801;
      t772 = t783*t800;
      t771 = t783*t801;
      t769 = x*t778;
      t768 = y*t779;
      t766 = t781*t800;
      t765 = t781*t801;
      t763 = t786+RATIONAL(1.0,1.0)-t791;
      t762 = t782*x+t776+t794;
      t761 = t769+t784*x+t794;
      t760 = t803+t775+t795;
      t759 = t768+t802+t795;
      t758 = t765+t772+t798;
      t757 = t771+t772+t799;
      t756 = t766+t765+t799;
      t755 = t766+t771+t798;
      coeffs_dzz->coeff_m1_m1_m1 = t757;
      coeffs_dzz->coeff_0_m1_m1 = t759;
      coeffs_dzz->coeff_p1_m1_m1 = t755;
      coeffs_dzz->coeff_m1_0_m1 = t761;
      coeffs_dzz->coeff_0_0_m1 = t763;
      coeffs_dzz->coeff_p1_0_m1 = t762;
      coeffs_dzz->coeff_m1_p1_m1 = t758;
      coeffs_dzz->coeff_0_p1_m1 = t760;
      coeffs_dzz->coeff_p1_p1_m1 = t756;
      coeffs_dzz->coeff_m1_m1_0 = t768+t769+t796;
      coeffs_dzz->coeff_0_m1_0 = -t790+t793;
      coeffs_dzz->coeff_p1_m1_0 = t776+t768+t797;
      coeffs_dzz->coeff_m1_0_0 = -t789+t792;
      coeffs_dzz->coeff_0_0_0 = (t786+1.0)*RATIONAL(-2.0,1.0)+t791*RATIONAL(2.0
,1.0);
      coeffs_dzz->coeff_p1_0_0 = t789+t792;
      coeffs_dzz->coeff_m1_p1_0 = t769+t775+t797;
      coeffs_dzz->coeff_0_p1_0 = t790+t793;
      coeffs_dzz->coeff_p1_p1_0 = t775+t776+t796;
      coeffs_dzz->coeff_m1_m1_p1 = t757;
      coeffs_dzz->coeff_0_m1_p1 = t759;
      coeffs_dzz->coeff_p1_m1_p1 = t755;
      coeffs_dzz->coeff_m1_0_p1 = t761;
      coeffs_dzz->coeff_0_0_p1 = t763;
      coeffs_dzz->coeff_p1_0_p1 = t762;
      coeffs_dzz->coeff_m1_p1_p1 = t758;
      coeffs_dzz->coeff_0_p1_p1 = t760;
      coeffs_dzz->coeff_p1_p1_p1 = t756;