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


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

fp t71,
   t72,
   t76,
   t75,
   t70,
   t67;
      t71 = x*x;
      t72 = x*t71;
      t76 = RATIONAL(-1.0,12.0)*x+RATIONAL(1.0,6.0)*t72;
      t75 = RATIONAL(-2.0,3.0)*t72+RATIONAL(4.0,3.0)*x;
      t70 = t71*t71;
      t67 = RATIONAL(-1.0,8.0)*t71;
      coeff_dx_m2 = t67+RATIONAL(-1.0,24.0)*t70+RATIONAL(1.0,20.0)+t76;
      coeff_dx_m1 = RATIONAL(-1.0,2.0)+t67+RATIONAL(5.0,24.0)*t70+t75;
      coeff_dx_0 = RATIONAL(5.0,4.0)*t71+t72+RATIONAL(-5.0,12.0)*t70+RATIONAL(
-5.0,2.0)*x+RATIONAL(-1.0,3.0);
      coeff_dx_p1 = RATIONAL(5.0,12.0)*t70+RATIONAL(-7.0,4.0)*t71+RATIONAL(1.0,
1.0)+t75;
      coeff_dx_p2 = RATIONAL(-1.0,4.0)+RATIONAL(7.0,8.0)*t71+RATIONAL(-5.0,24.0
)*t70+t76;
      coeff_dx_p3 = RATIONAL(1.0,24.0)*t70+RATIONAL(1.0,30.0)+t67;