aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c
blob: 48bfb8ac6049f60aab7931baa113cfbfcf756f7b (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 t591;
fp t598;
fp t612;
fp t567;
fp t558;
fp t589;
fp t630;
fp t592;
fp t629;
fp t581;
fp t566;
fp t559;
fp t590;
fp t580;
fp t561;
fp t628;
fp t593;
fp t585;
fp t569;
fp t627;
fp t594;
fp t586;
fp t574;
fp t560;
fp t626;
fp t625;
fp t563;
fp t579;
fp t624;
fp t584;
fp t623;
fp t565;
fp t622;
fp t587;
fp t621;
fp t571;
fp t620;
fp t619;
fp t611;
fp t573;
fp t618;
fp t575;
fp t582;
fp t617;
fp t597;
fp t613;
fp t609;
fp t616;
fp t615;
fp t614;
fp t596;
fp t588;
fp t557;
fp t577;
fp t610;
fp t608;
fp t607;
fp t606;
fp t605;
fp t562;
fp t564;
fp t604;
fp t572;
fp t603;
fp t570;
fp t602;
fp t601;
fp t600;
fp t599;
fp t576;
fp t568;
      t591 = RATIONAL(1.0,2.0);
      t598 = y*y;
      t612 = t591*t598;
      t567 = x*t612;
      t558 = z*t567;
      t589 = RATIONAL(1.0,8.0);
      t630 = t558+t589*t598;
      t592 = RATIONAL(-1.0,8.0);
      t629 = t558+t592*t598;
      t581 = t591*z;
      t566 = y*t581;
      t559 = x*t566;
      t590 = RATIONAL(1.0,4.0);
      t580 = t590*t598;
      t561 = z*t580;
      t628 = t559+t561;
      t593 = RATIONAL(-1.0,4.0);
      t585 = t593*t598;
      t569 = z*t585;
      t627 = t559+t569;
      t594 = RATIONAL(-1.0,2.0);
      t586 = t594*z;
      t574 = y*t586;
      t560 = x*t574;
      t626 = t560+t561;
      t625 = t560+t569;
      t563 = x*t580;
      t579 = t589*y;
      t624 = t563+t579;
      t584 = t592*y;
      t623 = t563+t584;
      t565 = z*t612;
      t622 = t565+t586;
      t587 = t594*x;
      t621 = t567+t587;
      t571 = x*t585;
      t620 = t571+t579;
      t619 = t571+t584;
      t611 = t594*t598;
      t573 = z*t611;
      t618 = t573+t581;
      t575 = x*t611;
      t582 = t591*x;
      t617 = t575+t582;
      t597 = x*z;
      t613 = x*t598;
      t609 = z*t613;
      t616 = t597-t609;
      t615 = t590*y;
      t614 = t593*y;
      t596 = y*t597;
      t588 = RATIONAL(2.0,1.0);
      t557 = t588*t609;
      t577 = RATIONAL(-2.0,1.0)*t597;
      t610 = t557+t577;
      t608 = x-t613;
      t607 = (1.0-t598)*z;
      t606 = -t596-t609;
      t605 = t596-t609;
      t562 = z*t615;
      t564 = x*t615;
      t604 = t562+t564+t630;
      t572 = x*t614;
      t603 = t562+t572+t629;
      t570 = z*t614;
      t602 = t564+t570+t629;
      t601 = t570+t572+t630;
      t600 = t580+t593+t616;
      t599 = t585+t590+t616;
      t576 = y*t587;
      t568 = y*t582;
      coeffs_dxz->coeff_m1_m1_m1 = t604+t619+t625;
      coeffs_dxz->coeff_0_m1_m1 = t567+t576+t605;
      coeffs_dxz->coeff_p1_m1_m1 = t602+t620+t626;
      coeffs_dxz->coeff_m1_0_m1 = t599+t621+t622;
      coeffs_dxz->coeff_0_0_m1 = t608+t610;
      coeffs_dxz->coeff_p1_0_m1 = t600+t618+t621;
      coeffs_dxz->coeff_m1_p1_m1 = t601+t620+t627;
      coeffs_dxz->coeff_0_p1_m1 = t567+t568+t606;
      coeffs_dxz->coeff_p1_p1_m1 = t603+t619+t628;
      coeffs_dxz->coeff_m1_m1_0 = t565+t574+t605;
      coeffs_dxz->coeff_0_m1_0 = y*t577+t557;
      coeffs_dxz->coeff_p1_m1_0 = t573+t566+t605;
      coeffs_dxz->coeff_m1_0_0 = t607+t610;
      coeffs_dxz->coeff_0_0_0 = (RATIONAL(4.0,1.0)+RATIONAL(-4.0,1.0)*t598)*
t597;
      coeffs_dxz->coeff_p1_0_0 = -t607+t610;
      coeffs_dxz->coeff_m1_p1_0 = t566+t565+t606;
      coeffs_dxz->coeff_0_p1_0 = t557+t588*t596;
      coeffs_dxz->coeff_p1_p1_0 = t573+t574+t606;
      coeffs_dxz->coeff_m1_m1_p1 = t603+t624+t625;
      coeffs_dxz->coeff_0_m1_p1 = t568+t575+t605;
      coeffs_dxz->coeff_p1_m1_p1 = t601+t623+t626;
      coeffs_dxz->coeff_m1_0_p1 = t600+t617+t622;
      coeffs_dxz->coeff_0_0_p1 = -t608+t610;
      coeffs_dxz->coeff_p1_0_p1 = t599+t617+t618;
      coeffs_dxz->coeff_m1_p1_p1 = t602+t623+t627;
      coeffs_dxz->coeff_0_p1_p1 = t575+t576+t606;
      coeffs_dxz->coeff_p1_p1_p1 = t604+t624+t628;