aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c
blob: a17d58128cbbf5bdedb40b1ad2526dc2d5579808 (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
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
fp t716;
fp t720;
fp t714;
fp t737;
fp t694;
fp t682;
fp t756;
fp t711;
fp t740;
fp t693;
fp t710;
fp t742;
fp t696;
fp t755;
fp t715;
fp t754;
fp t705;
fp t753;
fp t691;
fp t752;
fp t717;
fp t738;
fp t692;
fp t681;
fp t708;
fp t751;
fp t699;
fp t750;
fp t704;
fp t739;
fp t684;
fp t680;
fp t749;
fp t697;
fp t701;
fp t688;
fp t748;
fp t695;
fp t747;
fp t746;
fp t698;
fp t686;
fp t745;
fp t735;
fp t687;
fp t744;
fp t743;
fp t719;
fp t741;
fp t733;
fp t679;
fp t690;
fp t736;
fp t734;
fp t732;
fp t731;
fp t730;
fp t729;
fp t728;
fp t727;
fp t726;
fp t725;
fp t724;
fp t723;
fp t722;
fp t721;
fp t718;
fp t689;
      t716 = RATIONAL(1.0,8.0);
      t720 = x*x;
      t714 = RATIONAL(1.0,2.0);
      t737 = t714*t720;
      t694 = y*t737;
      t682 = z*t694;
      t756 = t682+t716*t720;
      t711 = RATIONAL(-1.0,4.0);
      t740 = x*t711;
      t693 = y*t740;
      t710 = RATIONAL(1.0,4.0);
      t742 = x*t710;
      t696 = z*t742;
      t755 = t693+t696;
      t715 = RATIONAL(-1.0,8.0);
      t754 = t715*t720+t682;
      t705 = t711*t720;
      t753 = t705+t710;
      t691 = z*t740;
      t752 = t693+t691;
      t717 = RATIONAL(-1.0,2.0);
      t738 = x*t717;
      t692 = y*t738;
      t681 = z*t692;
      t708 = t715*x;
      t751 = t681+t708;
      t699 = y*t742;
      t750 = t696+t699;
      t704 = t716*x;
      t739 = x*t714;
      t684 = z*t739;
      t680 = y*t684;
      t749 = t704+t680;
      t697 = y*t705;
      t701 = t710*t720;
      t688 = z*t701;
      t748 = t697+t688;
      t695 = z*t737;
      t747 = t717*z+t695;
      t746 = t711+t701;
      t698 = z*t705;
      t686 = y*t701;
      t745 = t698+t686;
      t735 = t717*t720;
      t687 = z*t735;
      t744 = t714*z+t687;
      t743 = t699+t691;
      t719 = y*z;
      t741 = y*t720;
      t733 = RATIONAL(2.0,1.0)*t719;
      t679 = t720*t733;
      t690 = RATIONAL(-2.0,1.0)*t719;
      t736 = t679+t690;
      t734 = z*t741;
      t732 = (t720-1.0)*z;
      t731 = -y+t741;
      t730 = -t734+z*t738;
      t729 = -t734+t684;
      t728 = t694-t734;
      t727 = y*t735-t734;
      t726 = t704+t681+t754;
      t725 = t686+t688+t756;
      t724 = t697+t698+t756;
      t723 = t680+t708+t754;
      t722 = t719+t714*y+t727;
      t721 = t719+t717*y+t728;
      t718 = x*t719;
      t689 = y*t739;
      coeffs_dyz->coeff_m1_m1_m1 = t724+t750+t751;
      coeffs_dyz->coeff_0_m1_m1 = t721+t747+t753;
      coeffs_dyz->coeff_p1_m1_m1 = t724+t749+t752;
      coeffs_dyz->coeff_m1_0_m1 = t692+t718+t728;
      coeffs_dyz->coeff_0_0_m1 = -t731+t736;
      coeffs_dyz->coeff_p1_0_m1 = t689-t718+t728;
      coeffs_dyz->coeff_m1_p1_m1 = t726+t743+t748;
      coeffs_dyz->coeff_0_p1_m1 = t721+t744+t746;
      coeffs_dyz->coeff_p1_p1_m1 = t723+t748+t755;
      coeffs_dyz->coeff_m1_m1_0 = t695+t718+t730;
      coeffs_dyz->coeff_0_m1_0 = -t732+t736;
      coeffs_dyz->coeff_p1_m1_0 = t695-t718+t729;
      coeffs_dyz->coeff_m1_0_0 = x*t690+t679;
      coeffs_dyz->coeff_0_0_0 = (RATIONAL(4.0,1.0)+RATIONAL(-4.0,1.0)*t720)*
t719;
      coeffs_dyz->coeff_p1_0_0 = x*t733+t679;
      coeffs_dyz->coeff_m1_p1_0 = t687+t718+t729;
      coeffs_dyz->coeff_0_p1_0 = t732+t736;
      coeffs_dyz->coeff_p1_p1_0 = t687-t718+t730;
      coeffs_dyz->coeff_m1_m1_p1 = t726+t745+t755;
      coeffs_dyz->coeff_0_m1_p1 = t722+t746+t747;
      coeffs_dyz->coeff_p1_m1_p1 = t723+t743+t745;
      coeffs_dyz->coeff_m1_0_p1 = t689+t718+t727;
      coeffs_dyz->coeff_0_0_p1 = t731+t736;
      coeffs_dyz->coeff_p1_0_p1 = t692-t718+t727;
      coeffs_dyz->coeff_m1_p1_p1 = t725+t751+t752;
      coeffs_dyz->coeff_0_p1_p1 = t722+t744+t753;
      coeffs_dyz->coeff_p1_p1_p1 = t725+t749+t750;