src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

fp t23,
   t24,
   t29,
   t21,
   t22,
   t28,
   t25,
   t20,
   t19;
      t23 = x*x;
      t24 = t23*t23;
      t29 = RATIONAL(-1.0,24.0)*t23+RATIONAL(1.0,24.0)*t24;
      t21 = RATIONAL(2.0,3.0);
      t22 = RATIONAL(-1.0,6.0);
      t28 = t21*t23+t22*t24;
      t25 = x*t23;
      t20 = RATIONAL(1.0,12.0);
      t19 = RATIONAL(-1.0,12.0);
      coeff_I_m2 = t20*x+t19*t25+t29;
      coeff_I_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(1.0,6.0)*t25+t28;
      coeff_I_0 = RATIONAL(-5.0,4.0)*t23+RATIONAL(1.0,1.0)+RATIONAL(1.0,4.0)*
t24;
      coeff_I_p1 = t22*t25+t21*x+t28;
      coeff_I_p2 = t20*t25+t19*x+t29;