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