aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c
blob: 017bfb0105b60092d0fa572316a84f7862321ec9 (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
fp t305;
fp t290;
fp t313;
fp t315;
fp t310;
fp t314;
fp t321;
fp t275;
fp t341;
fp t311;
fp t331;
fp t304;
fp t340;
fp t335;
fp t300;
fp t287;
fp t308;
fp t339;
fp t291;
fp t338;
fp t299;
fp t286;
fp t337;
fp t293;
fp t309;
fp t289;
fp t336;
fp t334;
fp t333;
fp t332;
fp t330;
fp t298;
fp t328;
fp t297;
fp t329;
fp t327;
fp t306;
fp t326;
fp t268;
fp t325;
fp t324;
fp t295;
fp t323;
fp t322;
fp t320;
fp t261;
fp t319;
fp t312;
fp t318;
fp t317;
fp t307;
fp t303;
fp t302;
fp t301;
fp t296;
fp t283;
fp t277;
fp t266;
fp t265;
fp t264;
fp t260;
      t305 = RATIONAL(-1.0,3.0);
      t290 = t305*x;
      t313 = x*x;
      t315 = t313*x;
      t310 = RATIONAL(1.0,4.0);
      t314 = y*y;
      t321 = t310*t314;
      t275 = t315*t321;
      t341 = t290+t275;
      t311 = RATIONAL(-1.0,6.0);
      t331 = t313*t314;
      t304 = RATIONAL(-1.0,4.0);
      t340 = t304*t331+t311*t313;
      t335 = t314*t315;
      t300 = RATIONAL(1.0,2.0);
      t287 = t300*t313;
      t308 = RATIONAL(3.0,4.0);
      t339 = t287+t308*t335;
      t291 = t311*x;
      t338 = x*t314;
      t299 = RATIONAL(-1.0,2.0);
      t286 = t299*x;
      t337 = t299*t314;
      t293 = t311*t315;
      t309 = RATIONAL(1.0,6.0);
      t289 = t309*t315;
      t336 = t300*t314;
      t334 = t299*t315;
      t333 = t300*t315;
      t332 = t314*t304;
      t330 = (t291+t289)*y;
      t298 = RATIONAL(-1.0,12.0);
      t328 = t298*t314;
      t297 = RATIONAL(1.0,18.0);
      t329 = t297*x+t315*t328;
      t327 = (t321+t298)*t313;
      t306 = RATIONAL(-3.0,4.0);
      t326 = t286+t306*t335;
      t268 = t315*t332;
      t325 = t291+t268;
      t324 = t304*t313+t308*t331;
      t295 = RATIONAL(1.0,12.0);
      t323 = t314*t295;
      t322 = t268+t309*x;
      t320 = (1.0+t286+t333)*y;
      t261 = y*t287;
      t319 = t261+(t290+t293)*y;
      t312 = x*y;
      t318 = t312+t261+y*t334;
      t317 = t275+t298*t315+(t295+t332)*x;
      t307 = RATIONAL(1.0,3.0);
      t303 = RATIONAL(3.0,2.0);
      t302 = RATIONAL(-2.0,1.0);
      t301 = RATIONAL(-3.0,2.0);
      t296 = RATIONAL(1.0,36.0);
      t283 = t295*t315;
      t277 = x*t337;
      t266 = x*t321;
      t265 = t306*t331;
      t264 = x*t336;
      t260 = t315*t323;
      coeffs_dy->coeff_m1_m1 = t260+t297*t315+(RATIONAL(1.0,9.0)+t309*t314)*x+
t319+t340;
      coeffs_dy->coeff_0_m1 = t305+t293+t266+t337+(t336-y+t307)*t313+t320+t322;
      coeffs_dy->coeff_p1_m1 = t277+t289+t318+t340+t341;
      coeffs_dy->coeff_p2_m1 = x*t323+RATIONAL(-1.0,18.0)*t315+t329+t330;
      coeffs_dy->coeff_m1_0 = t283+t277+(t307*t315+RATIONAL(2.0,3.0)*x-t313)*y+
t322+t324;
      coeffs_dy->coeff_0_0 = t299+t304*t315+t312+t310*x+(t303+t306*x+t301*t313)
*t314+(-t315+t302+RATIONAL(2.0,1.0)*t313)*y+t339;
      coeffs_dy->coeff_p1_0 = t303*t338+t310*t315+(t315+t302*x-t313)*y+t324+
t326;
      coeffs_dy->coeff_p2_0 = (t305*t315+t307*x)*y+t317;
      coeffs_dy->coeff_m1_p1 = t287+t265+t264+t293+t319+t341;
      coeffs_dy->coeff_0_p1 = t333+RATIONAL(1.0,1.0)-t313*y-t313+(t301+t308*x+
t303*t313)*t314+t320+t326;
      coeffs_dy->coeff_p1_p1 = x+t265+t334+t301*t338+t318+t339;
      coeffs_dy->coeff_p2_p1 = t289+t266+t325+t330;
      coeffs_dy->coeff_m1_p2 = t311*t338+t296*t315+t327+t329;
      coeffs_dy->coeff_0_p2 = t336+t311+(t337+t309)*t313+t317;
      coeffs_dy->coeff_p1_p2 = t283+t264+t325+t327;
      coeffs_dy->coeff_p2_p2 = t260+RATIONAL(-1.0,36.0)*t315+(t328+t296)*x;