aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Hermite/2d.coeffs/2d.cube.order2/coeff-I.compute.c
blob: b18719b67bb0a3b5b09f7a822d705cc2f2f81ca0 (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
fp t41,
   t40,
   t68,
   t38,
   t39,
   t67,
   t48,
   t49,
   t36,
   t69,
   t28,
   t66,
   t35,
   t34,
   t56,
   t7,
   t52,
   t65,
   t60,
   t37,
   t64,
   t61,
   t63,
   t27,
   t55,
   t21,
   t32,
   t51,
   t62,
   t59,
   t58,
   t57,
   t54,
   t53,
   t31,
   t50,
   t25,
   t29,
   t47,
   t46,
   t33,
   t30,
   t26,
   t24,
   t23,
   t22,
   t20,
   t19,
   t18,
   t16,
   t12,
   t11,
   t10,
   t9,
   t8,
   t6,
   t4;
      t41 = x*x;
      t40 = t41*x;
      t68 = t40+x;
      t38 = y*y;
      t39 = y*t38;
      t67 = y+t39;
      t48 = t68*t38;
      t49 = t67*t41;
      t36 = RATIONAL(1.0,4.0);
      t69 = x*t36;
      t28 = RATIONAL(-1.0,4.0);
      t66 = t28*y;
      t35 = RATIONAL(-5.0,2.0);
      t34 = RATIONAL(-3.0,4.0);
      t56 = t34*t40;
      t7 = t39*t56;
      t52 = t41*t38;
      t65 = t7+t35*t52;
      t60 = t28*t39;
      t37 = RATIONAL(-1.0,2.0);
      t64 = t40*t60+t37*t52;
      t61 = t39*t40;
      t63 = RATIONAL(-9.0,4.0)*t61+RATIONAL(-5.0,1.0)*t52;
      t27 = RATIONAL(3.0,4.0);
      t55 = t27*t39;
      t21 = t40*t55;
      t32 = RATIONAL(5.0,4.0);
      t51 = t32*t38;
      t62 = t41*t51+t21;
      t59 = t36*t40;
      t58 = t34*t39;
      t57 = t27*t40;
      t54 = t41*t39;
      t53 = t40*t38;
      t31 = RATIONAL(2.0,1.0);
      t50 = t21+t31*t52+x*t66;
      t25 = RATIONAL(-5.0,4.0);
      t29 = RATIONAL(1.0,2.0);
      t47 = t25*t38+t29;
      t46 = t29+t25*t41;
      t33 = RATIONAL(3.0,2.0);
      t30 = RATIONAL(-3.0,2.0);
      t26 = RATIONAL(3.0,1.0);
      t24 = RATIONAL(-3.0,1.0);
      t23 = RATIONAL(15.0,4.0);
      t22 = RATIONAL(-15.0,4.0);
      t20 = y*t69;
      t19 = x*t60;
      t18 = x*t55;
      t16 = t40*t66;
      t12 = RATIONAL(9.0,4.0)*t61;
      t11 = x*t58;
      t10 = t39*t69;
      t9 = y*t59;
      t8 = t39*t59;
      t6 = y*t56;
      t4 = y*t57;
      coeff_I_m1_m1 = t20+t10+t9+t8+t52+(t49+t48)*t37;
      coeff_I_0_m1 = t38+t6+t33*t53+t65+t67*(t32*t41+t37);
      coeff_I_p1_m1 = t19+t4-t49+(t30*t40+t29*x)*t38+t50;
      coeff_I_p2_m1 = t29*t53+t16+t36*t49+t64;
      coeff_I_m1_0 = t11+t41+t33*t54+t65+t68*(t51+t37);
      coeff_I_0_0 = t12+RATIONAL(1.0,1.0)+(t35+t22*t39)*t41+(t39+t40)*t33+(t35+
t22*t40+RATIONAL(25.0,4.0)*t41)*t38;
      coeff_I_p1_0 = t18+(t31+t26*t39)*t41+(t30+t23*t38)*t40+t47*x+t63;
      coeff_I_p2_0 = (t37+t58)*t41+t47*t40+t62;
      coeff_I_m1_p1 = t16+t18+(t30*t39+t29*y)*t41-t48+t50;
      coeff_I_0_p1 = t4+(t23*t41+t30)*t39+(t26*t40+t31)*t38+t46*y+t63;
      coeff_I_p1_p1 = t20+t11+t6+t12+(t24*t39+y)*t41+(RATIONAL(4.0,1.0)*t41+x+
t24*t40)*t38;
      coeff_I_p2_p1 = t7+t9+t53+(t55+t66-t38)*t41;
      coeff_I_m1_p2 = t29*t54+t19+t36*t48+t64;
      coeff_I_0_p2 = t46*t39+(t56+t37)*t38+t62;
      coeff_I_p1_p2 = t7+t54+t10+(t57+t28*x-t41)*t38;
      coeff_I_p2_p2 = t8+t36*t52+(t54+t53)*t28;