aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxy.compute.c
blob: d5b8bd96edcd69a362823cc517a2db1e0292e3ab (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 t521;
fp t527;
fp t541;
fp t496;
fp t522;
fp t513;
fp t560;
fp t518;
fp t515;
fp t517;
fp t559;
fp t542;
fp t498;
fp t511;
fp t558;
fp t524;
fp t516;
fp t523;
fp t510;
fp t557;
fp t545;
fp t500;
fp t504;
fp t556;
fp t512;
fp t499;
fp t555;
fp t507;
fp t508;
fp t554;
fp t492;
fp t506;
fp t553;
fp t546;
fp t503;
fp t552;
fp t493;
fp t505;
fp t551;
fp t495;
fp t489;
fp t550;
fp t509;
fp t491;
fp t549;
fp t548;
fp t547;
fp t526;
fp t514;
fp t544;
fp t525;
fp t519;
fp t540;
fp t486;
fp t494;
fp t543;
fp t539;
fp t538;
fp t537;
fp t536;
fp t535;
fp t534;
fp t533;
fp t532;
fp t531;
fp t530;
fp t529;
fp t490;
fp t528;
fp t501;
fp t497;
      t521 = RATIONAL(1.0,2.0);
      t527 = z*z;
      t541 = t521*t527;
      t496 = y*t541;
      t522 = RATIONAL(-1.0,2.0);
      t513 = t522*y;
      t560 = t496+t513;
      t518 = RATIONAL(-1.0,4.0);
      t515 = t518*t527;
      t517 = RATIONAL(1.0,4.0);
      t559 = t515+t517;
      t542 = t522*t527;
      t498 = y*t542;
      t511 = t521*y;
      t558 = t498+t511;
      t524 = RATIONAL(-1.0,8.0);
      t516 = t524*z;
      t523 = RATIONAL(1.0,8.0);
      t510 = t523*t527;
      t557 = t516+t510;
      t545 = z*t518;
      t500 = y*t545;
      t504 = y*t515;
      t556 = t500+t504;
      t512 = t517*t527;
      t499 = y*t512;
      t555 = t499+t500;
      t507 = t523*z;
      t508 = t524*t527;
      t554 = t507+t508;
      t492 = x*t545;
      t506 = x*t515;
      t553 = t492+t506;
      t546 = z*t517;
      t503 = y*t546;
      t552 = t499+t503;
      t493 = x*t512;
      t505 = x*t546;
      t551 = t493+t505;
      t495 = z*t513;
      t489 = x*t496;
      t550 = x*t495+t489;
      t509 = t521*x;
      t491 = z*t509;
      t549 = y*t491+t489;
      t548 = t504+t503;
      t547 = t518+t512;
      t526 = y*x;
      t514 = t522*x;
      t544 = y*t527;
      t525 = z*t526;
      t519 = RATIONAL(2.0,1.0);
      t540 = x*t544;
      t486 = t519*t540;
      t494 = RATIONAL(-2.0,1.0)*t526;
      t543 = t494+t486;
      t539 = (1.0-t527)*x;
      t538 = y-t544;
      t537 = -t525-t540;
      t536 = t525-t540;
      t535 = x*t541-t540;
      t534 = t495-t540;
      t533 = t516+t508+t549;
      t532 = t505+t506+t550;
      t531 = t507+t510+t549;
      t530 = t492+t493+t550;
      t529 = t514+t526+t535;
      t490 = x*t542;
      t528 = t490+t526-t540+t509;
      t501 = z*t514;
      t497 = z*t511;
      coeffs_dxy->coeff_m1_m1_m1 = t532+t548+t557;
      coeffs_dxy->coeff_0_m1_m1 = t501+t525+t535;
      coeffs_dxy->coeff_p1_m1_m1 = t532+t554+t555;
      coeffs_dxy->coeff_m1_0_m1 = t525+t496+t534;
      coeffs_dxy->coeff_0_0_m1 = t486+z*t494;
      coeffs_dxy->coeff_p1_0_m1 = t497+t498+t536;
      coeffs_dxy->coeff_m1_p1_m1 = t530+t548+t554;
      coeffs_dxy->coeff_0_p1_m1 = t491+t490+t536;
      coeffs_dxy->coeff_p1_p1_m1 = t530+t555+t557;
      coeffs_dxy->coeff_m1_m1_0 = t529+t559+t560;
      coeffs_dxy->coeff_0_m1_0 = t539+t543;
      coeffs_dxy->coeff_p1_m1_0 = t529+t547+t558;
      coeffs_dxy->coeff_m1_0_0 = t538+t543;
      coeffs_dxy->coeff_0_0_0 = (RATIONAL(-4.0,1.0)*t527+RATIONAL(4.0,1.0))*
t526;
      coeffs_dxy->coeff_p1_0_0 = -t538+t543;
      coeffs_dxy->coeff_m1_p1_0 = t528+t547+t560;
      coeffs_dxy->coeff_0_p1_0 = -t539+t543;
      coeffs_dxy->coeff_p1_p1_0 = t528+t558+t559;
      coeffs_dxy->coeff_m1_m1_p1 = t531+t553+t556;
      coeffs_dxy->coeff_0_m1_p1 = t491-t525+t535;
      coeffs_dxy->coeff_p1_m1_p1 = t533+t552+t553;
      coeffs_dxy->coeff_m1_0_p1 = t497+t496+t537;
      coeffs_dxy->coeff_0_0_p1 = t486+t519*t525;
      coeffs_dxy->coeff_p1_0_p1 = -t525+t498+t534;
      coeffs_dxy->coeff_m1_p1_p1 = t533+t551+t556;
      coeffs_dxy->coeff_0_p1_p1 = t490+t501+t537;
      coeffs_dxy->coeff_p1_p1_p1 = t531+t551+t552;