aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c
blob: 9ed02d4b25956d0802c1857e09b1bdee5f4a0845 (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
fp t233;
fp t234;
fp t228;
fp t205;
fp t230;
fp t213;
fp t258;
fp t229;
fp t232;
fp t240;
fp t192;
fp t224;
fp t208;
fp t257;
fp t207;
fp t256;
fp t227;
fp t231;
fp t242;
fp t255;
fp t254;
fp t214;
fp t201;
fp t218;
fp t206;
fp t219;
fp t253;
fp t252;
fp t251;
fp t209;
fp t250;
fp t217;
fp t216;
fp t249;
fp t225;
fp t248;
fp t247;
fp t246;
fp t245;
fp t223;
fp t182;
fp t244;
fp t243;
fp t241;
fp t239;
fp t198;
fp t238;
fp t237;
fp t236;
fp t226;
fp t222;
fp t221;
fp t220;
fp t215;
fp t211;
fp t193;
fp t184;
fp t180;
fp t179;
      t233 = y*y;
      t234 = t233*y;
      t228 = RATIONAL(1.0,6.0);
      t205 = t228*t234;
      t230 = RATIONAL(-1.0,6.0);
      t213 = t230*y;
      t258 = t213+t205;
      t229 = RATIONAL(1.0,4.0);
      t232 = x*x;
      t240 = t229*t232;
      t192 = t234*t240;
      t224 = RATIONAL(-1.0,3.0);
      t208 = t224*y;
      t257 = t192+t208;
      t207 = t230*t234;
      t256 = (t208+t207)*x;
      t227 = RATIONAL(3.0,4.0);
      t231 = x*y;
      t242 = t232*t234;
      t255 = t231+t227*t242;
      t254 = t232*y;
      t214 = RATIONAL(1.0,12.0);
      t201 = t214*y;
      t218 = RATIONAL(-1.0,2.0);
      t206 = t218*y;
      t219 = RATIONAL(1.0,2.0);
      t253 = t219*t234;
      t252 = t218*t232;
      t251 = t219*t232;
      t209 = t219*t233;
      t250 = t218*t234;
      t217 = RATIONAL(-1.0,12.0);
      t216 = RATIONAL(1.0,18.0);
      t249 = t217*t242+t216*y;
      t225 = RATIONAL(-3.0,4.0);
      t248 = t225*t242+t206;
      t247 = t258*x;
      t246 = t232*t233;
      t245 = (t217+t240)*t233;
      t223 = RATIONAL(-1.0,4.0);
      t182 = t223*t242;
      t244 = t182+y*t240;
      t243 = t227*t246+t223*t233;
      t241 = t182+t214*t234;
      t239 = (1.0+t253+t206)*x;
      t198 = x*t209;
      t238 = t198+t209+t225*t246;
      t237 = t198+t223*t246+t230*t233;
      t236 = t217*t234+t201+t192+t223*t254;
      t226 = RATIONAL(1.0,3.0);
      t222 = RATIONAL(3.0,2.0);
      t221 = RATIONAL(-2.0,1.0);
      t220 = RATIONAL(-3.0,2.0);
      t215 = RATIONAL(1.0,36.0);
      t211 = t228*y;
      t193 = y*t251;
      t184 = y*t252;
      t180 = x*t250;
      t179 = t214*t242;
      coeffs_dx->coeff_m1_m1 = t179+t216*t234+(RATIONAL(1.0,9.0)+t228*t232)*y+
t237+t256;
      coeffs_dx->coeff_0_m1 = t184+t211+(RATIONAL(2.0,3.0)*y-t233+t226*t234)*x+
t241+t243;
      coeffs_dx->coeff_p1_m1 = t193+t207+t238+t256+t257;
      coeffs_dx->coeff_p2_m1 = t230*t254+t215*t234+t245+t249;
      coeffs_dx->coeff_m1_0 = t252+t211+t207+t224+(-x+t251+t226)*t233+t239+t244
;
      coeffs_dx->coeff_0_0 = t223*t234+t209+t218+t229*y+(t225*y+t222+t220*t233)
*t232+(-t234+t221+RATIONAL(2.0,1.0)*t233)*x+t255;
      coeffs_dx->coeff_p1_0 = RATIONAL(1.0,1.0)+t253-x*t233-t233+(t222*t233+
t220+t227*y)*t232+t239+t248;
      coeffs_dx->coeff_p2_0 = t228*t233+t230+(t219+t218*t233)*t232+t236;
      coeffs_dx->coeff_m1_p1 = t205+t180+t184+t231+t237+t257;
      coeffs_dx->coeff_0_p1 = t222*t254+t229*t234+(-t233+t234+t221*y)*x+t243+
t248;
      coeffs_dx->coeff_p1_p1 = t220*t254+t180+t250+y+t238+t255;
      coeffs_dx->coeff_p2_p1 = t193+t213+t241+t245;
      coeffs_dx->coeff_m1_p2 = t232*t201+RATIONAL(-1.0,18.0)*t234+t247+t249;
      coeffs_dx->coeff_0_p2 = (t224*t234+t226*y)*x+t236;
      coeffs_dx->coeff_p1_p2 = t244+t247+t258;
      coeffs_dx->coeff_p2_p2 = t179+RATIONAL(-1.0,36.0)*t234+(t215+t217*t232)*y
;