aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c
blob: 25d4963eb621922e7074fbd38a9fc36461ddb27e (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
fp t127;
fp t135;
fp t137;
fp t172;
fp t129;
fp t136;
fp t138;
fp t150;
fp t177;
fp t134;
fp t151;
fp t126;
fp t152;
fp t176;
fp t130;
fp t123;
fp t157;
fp t93;
fp t175;
fp t133;
fp t154;
fp t128;
fp t119;
fp t99;
fp t174;
fp t156;
fp t91;
fp t132;
fp t165;
fp t173;
fp t171;
fp t170;
fp t169;
fp t168;
fp t167;
fp t125;
fp t166;
fp t164;
fp t160;
fp t163;
fp t121;
fp t155;
fp t162;
fp t161;
fp t124;
fp t159;
fp t120;
fp t158;
fp t153;
fp t149;
fp t148;
fp t101;
fp t147;
fp t115;
fp t146;
fp t145;
fp t144;
fp t143;
fp t106;
fp t142;
fp t141;
fp t131;
fp t122;
fp t118;
fp t116;
fp t114;
fp t113;
fp t109;
fp t87;
fp t84;
fp t83;
      t127 = RATIONAL(-1.0,2.0);
      t135 = x*x;
      t137 = t135*x;
      t172 = x*y;
      t129 = RATIONAL(-1.0,4.0);
      t136 = y*y;
      t138 = t136*y;
      t150 = t129*t138;
      t177 = t127*t172+t137*t150;
      t134 = RATIONAL(-1.0,6.0);
      t151 = t134*t135;
      t126 = RATIONAL(-1.0,12.0);
      t152 = t126*t135;
      t176 = y*t151+t138*t152;
      t130 = RATIONAL(-1.0,3.0);
      t123 = RATIONAL(1.0,12.0);
      t157 = t123*t137;
      t93 = t138*t157;
      t175 = t93+t130*t172;
      t133 = RATIONAL(1.0,4.0);
      t154 = t133*t138;
      t128 = RATIONAL(1.0,2.0);
      t119 = t128*t135;
      t99 = y*t119;
      t174 = t99+t137*t154;
      t156 = t126*t137;
      t91 = t138*t156;
      t132 = RATIONAL(1.0,6.0);
      t165 = y*t132;
      t173 = t91+x*t165;
      t171 = x*t138;
      t170 = t134*x;
      t169 = x*t136;
      t168 = y*t137;
      t167 = t133*y;
      t125 = RATIONAL(1.0,18.0);
      t166 = y*t125;
      t164 = t126*t169+t136*t157;
      t160 = t123*t138;
      t163 = y*t152+t135*t160;
      t121 = RATIONAL(-1.0,36.0);
      t155 = t121*t137;
      t162 = t138*t155+x*t166;
      t161 = t91+y*t170;
      t124 = RATIONAL(1.0,36.0);
      t159 = t124*t137;
      t120 = t128*t136;
      t158 = t135*t136;
      t153 = t129*t137;
      t149 = t133*t136;
      t148 = t134*t136;
      t101 = t135*t149;
      t147 = t136*t156+x*t148+t101;
      t115 = x*t120;
      t146 = t136*t153+t101+t115;
      t145 = t93+t126*t171+y*t156+t123*t172;
      t144 = t132*t135+t134;
      t143 = t134+t132*t136;
      t106 = t127*t158;
      t142 = t119+t106+(t154+t129*y)*t135;
      t141 = t137*t149+t129*t169+t106+t120;
      t131 = RATIONAL(1.0,3.0);
      t122 = RATIONAL(-1.0,18.0);
      t118 = t132*t171;
      t116 = t137*t165;
      t114 = t138*t170;
      t113 = t135*t150;
      t109 = t134*t168;
      t87 = x*t160;
      t84 = y*t157;
      t83 = t138*t159;
      coeffs_I->coeff_m1_m1 = t137*t166+t83+(t125*t138+RATIONAL(1.0,9.0)*y)*x+
t147+t176;
      coeffs_I->coeff_0_m1 = t87+t109+t144*t138+(t131*t135+t130)*y+t141+t173;
      coeffs_I->coeff_p1_m1 = t114+t116+t146+t175+t176;
      coeffs_I->coeff_p2_m1 = t122*t168+t124*t171+t162+t164;
      coeffs_I->coeff_m1_0 = t84+t114+t143*t137+(t130+t131*t136)*x+t142+t173;
      coeffs_I->coeff_0_0 = -t136+RATIONAL(1.0,1.0)+t115-t135+y*t153+t158+(t138
+t137)*t128+(t167+t150)*x+(t135*t138+t137*t136+x+y)*t127+t174;
      coeffs_I->coeff_p1_0 = (t128*t138-t136+1.0)*x+(t127+t120+t167)*t137+t142+
t177;
      coeffs_I->coeff_p2_0 = (t132+t148)*t137+t143*x+t145;
      coeffs_I->coeff_m1_p1 = t99+t118+t113+t109+t147+t175;
      coeffs_I->coeff_0_p1 = (t128*t137-t135+1.0)*y+(t127+t133*x+t119)*t138+
t141+t177;
      coeffs_I->coeff_p1_p1 = t127*t168+t113+(y+t127*t138)*x+t146+t174;
      coeffs_I->coeff_p2_p1 = t116+t87+t161+t164;
      coeffs_I->coeff_m1_p2 = y*t159+t122*t171+t162+t163;
      coeffs_I->coeff_0_p2 = (t132+t151)*t138+t144*y+t145;
      coeffs_I->coeff_p1_p2 = t84+t118+t161+t163;
      coeffs_I->coeff_p2_p2 = y*t155+t83+(t121*t138+t124*y)*x;