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


1
2
3
4
5
6
7
8

fp t13,
   t12;
      t13 = x*x;
      t12 = RATIONAL(-1.0,2.0);
      coeff_dx_m1 = x+t12*t13+RATIONAL(-1.0,3.0);
      coeff_dx_0 = t12+RATIONAL(-2.0,1.0)*x+RATIONAL(3.0,2.0)*t13;
      coeff_dx_p1 = RATIONAL(-3.0,2.0)*t13+RATIONAL(1.0,1.0)+x;
      coeff_dx_p2 = RATIONAL(1.0,2.0)*t13+RATIONAL(-1.0,6.0);