aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c
blob: b189aa69a556a3e9c2204e2b0ef8ec2746c3d552 (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
fp t103,
   t107,
   t106,
   t102,
   t113,
   t112,
   t96,
   t101,
   t100,
   t111,
   t105,
   t99,
   t98,
   t97;
      t103 = x*x;
      t107 = x*t103;
      t106 = t107*t107;
      t102 = t103*t103;
      t113 = RATIONAL(-1.0,144.0)*t102+RATIONAL(1.0,720.0)*t106+RATIONAL(1.0,
180.0)*t103;
      t112 = RATIONAL(-1.0,120.0)*t106+RATIONAL(-3.0,40.0)*t103+RATIONAL(1.0,
12.0)*t102;
      t96 = RATIONAL(-13.0,48.0);
      t101 = RATIONAL(3.0,4.0);
      t100 = RATIONAL(1.0,48.0);
      t111 = t96*t102+t100*t106+t101*t103;
      t105 = t107*t103;
      t99 = RATIONAL(-1.0,60.0);
      t98 = RATIONAL(1.0,60.0);
      t97 = RATIONAL(-1.0,48.0);
      coeff_I_m3 = t99*x+t100*t107+RATIONAL(-1.0,240.0)*t105+t113;
      coeff_I_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t98*t105+t112;
      coeff_I_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t97*t105+t111;
      coeff_I_0 = RATIONAL(7.0,18.0)*t102+RATIONAL(-49.0,36.0)*t103+RATIONAL(
1.0,1.0)+RATIONAL(-1.0,36.0)*t106;
      coeff_I_p1 = t96*t107+t101*x+t100*t105+t111;
      coeff_I_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t99*t105+t112;
      coeff_I_p3 = t97*t107+t98*x+RATIONAL(1.0,240.0)*t105+t113;