aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeffs-I.compute.c
blob: 18caffa96340d4ac655b9523e04b090c3a1b6bba (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
fp t38;
fp t41;
fp t65;
fp t40;
fp t39;
fp t63;
fp t48;
fp t49;
fp t36;
fp t68;
fp t29;
fp t67;
fp t66;
fp t32;
fp t27;
fp t53;
fp t7;
fp t51;
fp t64;
fp t58;
fp t62;
fp t35;
fp t33;
fp t55;
fp t8;
fp t61;
fp t37;
fp t57;
fp t60;
fp t59;
fp t56;
fp t54;
fp t52;
fp t31;
fp t50;
fp t26;
fp t28;
fp t47;
fp t45;
fp t34;
fp t30;
fp t25;
fp t24;
fp t23;
fp t22;
fp t21;
fp t20;
fp t19;
fp t17;
fp t15;
fp t14;
fp t13;
fp t12;
fp t11;
fp t10;
fp t4;
      t38 = y*y;
      t41 = y*t38;
      t65 = -t41-y;
      t40 = x*x;
      t39 = x*t40;
      t63 = t39+x;
      t48 = t63*t38;
      t49 = t65*t40;
      t36 = RATIONAL(1.0,4.0);
      t68 = y*t36;
      t29 = RATIONAL(-1.0,4.0);
      t67 = t29*x;
      t66 = t29*y;
      t32 = RATIONAL(5.0,4.0);
      t27 = RATIONAL(3.0,4.0);
      t53 = t27*t41;
      t7 = t39*t53;
      t51 = t40*t38;
      t64 = t32*t51+t7;
      t58 = t39*t41;
      t62 = RATIONAL(-9.0,4.0)*t58+RATIONAL(-5.0,1.0)*t51;
      t35 = RATIONAL(-5.0,2.0);
      t33 = RATIONAL(-3.0,4.0);
      t55 = t33*t41;
      t8 = t39*t55;
      t61 = t8+t35*t51;
      t37 = RATIONAL(-1.0,2.0);
      t57 = t29*t39;
      t60 = t41*t57+t37*t51;
      t59 = t40*t41;
      t56 = t33*t39;
      t54 = t27*t39;
      t52 = t39*t38;
      t31 = RATIONAL(2.0,1.0);
      t50 = t31*t51+x*t66+t7;
      t26 = RATIONAL(-5.0,4.0);
      t28 = RATIONAL(1.0,2.0);
      t47 = t26*t38+t28;
      t45 = t28+t26*t40;
      t34 = RATIONAL(3.0,2.0);
      t30 = RATIONAL(-3.0,2.0);
      t25 = RATIONAL(3.0,1.0);
      t24 = RATIONAL(-3.0,1.0);
      t23 = RATIONAL(15.0,4.0);
      t22 = RATIONAL(-15.0,4.0);
      t21 = x*t53;
      t20 = x*t68;
      t19 = y*t54;
      t17 = x*t55;
      t15 = RATIONAL(9.0,4.0)*t58;
      t14 = y*t56;
      t13 = t41*t67;
      t12 = t36*x*t41;
      t11 = t39*t68;
      t10 = t36*t58;
      t4 = y*t57;
      coeffs_I->coeff_m1_m1 = t20+t12+t11+t10+t51+(-t49+t48)*t37;
      coeffs_I->coeff_0_m1 = t38+t14+t34*t52+t61-t65*(t32*t40+t37);
      coeffs_I->coeff_p1_m1 = t13+t19+t49+(t30*t39+t28*x)*t38+t50;
      coeffs_I->coeff_p2_m1 = t28*t52+t4-t36*t49+t60;
      coeffs_I->coeff_m1_0 = t17+t40+t34*t59+t61+t63*(t32*t38+t37);
      coeffs_I->coeff_0_0 = t15+RATIONAL(1.0,1.0)+(t35+t22*t41)*t40+(t41+t39)*
t34+(t35+t22*t39+RATIONAL(25.0,4.0)*t40)*t38;
      coeffs_I->coeff_p1_0 = t21+(t31+t25*t41)*t40+(t30+t23*t38)*t39+t47*x+t62;
      coeffs_I->coeff_p2_0 = (t37+t55)*t40+t47*t39+t64;
      coeffs_I->coeff_m1_p1 = t4+t21+(t30*t41+t28*y)*t40-t48+t50;
      coeffs_I->coeff_0_p1 = t19+(t23*t40+t30)*t41+(t25*t39+t31)*t38+t45*y+t62;
      coeffs_I->coeff_p1_p1 = t20+t17+t14+t15+(t24*t41+y)*t40+(RATIONAL(4.0,1.0
)*t40+x+t24*t39)*t38;
      coeffs_I->coeff_p2_p1 = t8+t11+t52+(t53+t66-t38)*t40;
      coeffs_I->coeff_m1_p2 = t28*t59+t13+t36*t48+t60;
      coeffs_I->coeff_0_p2 = t45*t41+(t56+t37)*t38+t64;
      coeffs_I->coeff_p1_p2 = t8+t59+t12+(t54+t67-t40)*t38;
      coeffs_I->coeff_p2_p2 = t10+t36*t51+(t59+t52)*t29;