aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c
blob: 8ed883817cb4fb923dd302d82d78d062e91280d9 (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
165
166
167
168
169
170
fp t715;
fp t722;
fp t736;
fp t689;
fp t682;
fp t713;
fp t754;
fp t716;
fp t753;
fp t705;
fp t691;
fp t683;
fp t714;
fp t704;
fp t685;
fp t752;
fp t717;
fp t709;
fp t693;
fp t751;
fp t718;
fp t710;
fp t699;
fp t684;
fp t750;
fp t749;
fp t686;
fp t703;
fp t748;
fp t708;
fp t747;
fp t746;
fp t690;
fp t711;
fp t745;
fp t694;
fp t744;
fp t743;
fp t735;
fp t697;
fp t742;
fp t698;
fp t706;
fp t741;
fp t721;
fp t737;
fp t733;
fp t740;
fp t739;
fp t738;
fp t720;
fp t712;
fp t681;
fp t701;
fp t734;
fp t732;
fp t731;
fp t730;
fp t729;
fp t687;
fp t688;
fp t728;
fp t696;
fp t727;
fp t695;
fp t726;
fp t725;
fp t724;
fp t723;
fp t700;
fp t692;
      t715 = RATIONAL(1.0,2.0);
      t722 = x*x;
      t736 = t715*t722;
      t689 = z*t736;
      t682 = y*t689;
      t713 = RATIONAL(1.0,8.0);
      t754 = t682+t713*t722;
      t716 = RATIONAL(-1.0,8.0);
      t753 = t682+t716*t722;
      t705 = t715*z;
      t691 = x*t705;
      t683 = y*t691;
      t714 = RATIONAL(1.0,4.0);
      t704 = t714*t722;
      t685 = z*t704;
      t752 = t683+t685;
      t717 = RATIONAL(-1.0,4.0);
      t709 = t717*t722;
      t693 = z*t709;
      t751 = t683+t693;
      t718 = RATIONAL(-1.0,2.0);
      t710 = t718*z;
      t699 = x*t710;
      t684 = y*t699;
      t750 = t684+t685;
      t749 = t684+t693;
      t686 = y*t704;
      t703 = t713*x;
      t748 = t686+t703;
      t708 = t716*x;
      t747 = t686+t708;
      t746 = t689+t710;
      t690 = y*t736;
      t711 = t718*y;
      t745 = t690+t711;
      t694 = y*t709;
      t744 = t694+t703;
      t743 = t694+t708;
      t735 = t718*t722;
      t697 = z*t735;
      t742 = t697+t705;
      t698 = y*t735;
      t706 = t715*y;
      t741 = t698+t706;
      t721 = y*z;
      t737 = t722*z;
      t733 = y*t737;
      t740 = t721-t733;
      t739 = t714*x;
      t738 = t717*x;
      t720 = x*t721;
      t712 = RATIONAL(2.0,1.0);
      t681 = t712*t733;
      t701 = RATIONAL(-2.0,1.0)*t721;
      t734 = t681+t701;
      t732 = (-1.0+t722)*y;
      t731 = t737-z;
      t730 = -t720-t733;
      t729 = t720-t733;
      t687 = z*t739;
      t688 = y*t739;
      t728 = t687+t688+t754;
      t696 = y*t738;
      t727 = t687+t696+t753;
      t695 = z*t738;
      t726 = t688+t695+t753;
      t725 = t695+t696+t754;
      t724 = t704+t717+t740;
      t723 = t709+t714+t740;
      t700 = x*t711;
      t692 = x*t706;
      coeffs_dyz->coeff_m1_m1_m1 = t728+t743+t749;
      coeffs_dyz->coeff_0_m1_m1 = t723+t745+t746;
      coeffs_dyz->coeff_p1_m1_m1 = t725+t744+t751;
      coeffs_dyz->coeff_m1_0_m1 = t700+t690+t729;
      coeffs_dyz->coeff_0_0_m1 = -t732+t734;
      coeffs_dyz->coeff_p1_0_m1 = t692+t690+t730;
      coeffs_dyz->coeff_m1_p1_m1 = t726+t744+t750;
      coeffs_dyz->coeff_0_p1_m1 = t724+t742+t745;
      coeffs_dyz->coeff_p1_p1_m1 = t727+t743+t752;
      coeffs_dyz->coeff_m1_m1_0 = t689+t699+t729;
      coeffs_dyz->coeff_0_m1_0 = -t731+t734;
      coeffs_dyz->coeff_p1_m1_0 = t691+t689+t730;
      coeffs_dyz->coeff_m1_0_0 = t681+x*t701;
      coeffs_dyz->coeff_0_0_0 = (RATIONAL(4.0,1.0)+RATIONAL(-4.0,1.0)*t722)*
t721;
      coeffs_dyz->coeff_p1_0_0 = t712*t720+t681;
      coeffs_dyz->coeff_m1_p1_0 = t691+t697+t729;
      coeffs_dyz->coeff_0_p1_0 = t731+t734;
      coeffs_dyz->coeff_p1_p1_0 = t699+t697+t730;
      coeffs_dyz->coeff_m1_m1_p1 = t727+t748+t749;
      coeffs_dyz->coeff_0_m1_p1 = t724+t741+t746;
      coeffs_dyz->coeff_p1_m1_p1 = t726+t747+t751;
      coeffs_dyz->coeff_m1_0_p1 = t692+t698+t729;
      coeffs_dyz->coeff_0_0_p1 = t732+t734;
      coeffs_dyz->coeff_p1_0_p1 = t700+t698+t730;
      coeffs_dyz->coeff_m1_p1_p1 = t725+t747+t750;
      coeffs_dyz->coeff_0_p1_p1 = t723+t741+t742;
      coeffs_dyz->coeff_p1_p1_p1 = t728+t748+t752;