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


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

fp t36,
   t37,
   t34,
   t40,
   t35,
   t39;
      t36 = x*x;
      t37 = x*t36;
      t34 = RATIONAL(-1.0,12.0);
      t40 = t34*x+RATIONAL(1.0,6.0)*t37;
      t35 = RATIONAL(-2.0,3.0);
      t39 = t35*t37+RATIONAL(4.0,3.0)*x;
      coeff_dx_m2 = RATIONAL(-1.0,4.0)*t36+RATIONAL(1.0,12.0)+t40;
      coeff_dx_m1 = t35+RATIONAL(1.0,2.0)*t36+t39;
      coeff_dx_0 = t37+RATIONAL(-5.0,2.0)*x;
      coeff_dx_p1 = RATIONAL(2.0,3.0)+RATIONAL(-1.0,2.0)*t36+t39;
      coeff_dx_p2 = t34+RATIONAL(1.0,4.0)*t36+t40;