aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeff-I.compute.c
blob: 09905397690a8d2b55af5d1cda625d84af3cfe92 (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
fp t39,
   t46,
   t44,
   t36,
   t80,
   t43,
   t33,
   t42,
   t20,
   t79,
   t32,
   t71,
   t21,
   t78,
   t38,
   t24,
   t40,
   t26,
   t77,
   t45,
   t30,
   t41,
   t76,
   t27,
   t75,
   t34,
   t74,
   t70,
   t23,
   t73,
   t29,
   t47,
   t31,
   t72,
   t25,
   t69,
   t68,
   t37,
   t67,
   t35,
   t66,
   t28,
   t65,
   t64,
   t63,
   t62,
   t61,
   t60,
   t59,
   t58,
   t57,
   t56,
   t22,
   t55,
   t54,
   t53,
   t52,
   t51,
   t50,
   t49,
   t19,
   t48;
      t39 = RATIONAL(1.0,27.0);
      t46 = x*x;
      t44 = RATIONAL(-1.0,9.0);
      t36 = t44*t46;
      t80 = t39+t36;
      t43 = RATIONAL(-1.0,18.0);
      t33 = t43*x;
      t42 = RATIONAL(1.0,12.0);
      t20 = t42*x*z;
      t79 = t33+t20;
      t32 = t43*y;
      t71 = y*t42;
      t21 = x*t71;
      t78 = t32+t21;
      t38 = RATIONAL(-1.0,12.0);
      t24 = t38*x*y;
      t40 = RATIONAL(1.0,18.0);
      t26 = t40*y;
      t77 = t24+t26;
      t45 = z*z;
      t30 = t40*t45;
      t41 = RATIONAL(-2.0,27.0);
      t76 = t30+t41;
      t27 = t40*x;
      t75 = t26+t27;
      t34 = t44*t45;
      t74 = t39+t34;
      t70 = z*t38;
      t23 = x*t70;
      t73 = t23+t27;
      t29 = t40*t46;
      t47 = y*y;
      t31 = t40*t47;
      t72 = t29+t31;
      t25 = t40*z;
      t69 = t30+t25;
      t68 = t31+t34;
      t37 = RATIONAL(4.0,27.0);
      t67 = t37+t36;
      t35 = t44*t47;
      t66 = t29+t35;
      t28 = t43*z;
      t65 = t30+t28;
      t64 = t29+t39;
      t63 = t21+t72;
      t62 = t39+t69;
      t61 = t31+t65;
      t60 = t30+t35+t67;
      t59 = t37+t34+t66;
      t58 = t67+t68;
      t57 = t31+t36+t62;
      t56 = t29+t33+t23+t69;
      t22 = y*t70;
      t55 = t22+t26+t61;
      t54 = t32+t24+t27+t72;
      t53 = t33+t64+t68;
      t52 = t35+t64+t65;
      t51 = t29+t41+t55;
      t50 = t41+t31+t56;
      t49 = t54+t76;
      t19 = z*t71;
      t48 = t19+t20+t63+t76;
      coeff_I_m1_m1_m1 = t28+t32+t33+t48;
      coeff_I_0_m1_m1 = t32+t19+t61+t80;
      coeff_I_p1_m1_m1 = t28+t23+t19+t49;
      coeff_I_m1_0_m1 = t52+t79;
      coeff_I_0_0_m1 = t28+t60;
      coeff_I_p1_0_m1 = t52+t73;
      coeff_I_m1_p1_m1 = t24+t51+t79;
      coeff_I_0_p1_m1 = t55+t80;
      coeff_I_p1_p1_m1 = t21+t51+t73;
      coeff_I_m1_m1_0 = t53+t78;
      coeff_I_0_m1_0 = t32+t58;
      coeff_I_p1_m1_0 = t54+t74;
      coeff_I_m1_0_0 = t33+t59;
      coeff_I_0_0_0 = t35+t34+RATIONAL(7.0,27.0)+t36;
      coeff_I_p1_0_0 = t27+t59;
      coeff_I_m1_p1_0 = t53+t77;
      coeff_I_0_p1_0 = t26+t58;
      coeff_I_p1_p1_0 = t63+t74+t75;
      coeff_I_m1_m1_p1 = t22+t50+t78;
      coeff_I_0_m1_p1 = t32+t22+t57;
      coeff_I_p1_m1_p1 = t22+t25+t20+t49;
      coeff_I_m1_0_p1 = t35+t39+t56;
      coeff_I_0_0_p1 = t25+t60;
      coeff_I_p1_0_p1 = t27+t20+t62+t66;
      coeff_I_m1_p1_p1 = t19+t50+t77;
      coeff_I_0_p1_p1 = t26+t19+t57;
      coeff_I_p1_p1_p1 = t25+t48+t75;