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