aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c
blob: bc709b8053a3c081c3803a62ac21128e28fcbb8f (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
fp t35;
fp t41;
fp t38;
fp t29;
fp t36;
fp t40;
fp t37;
fp t39;
fp t34;
fp t33;
fp t30;
fp t28;
fp t27;
      t35 = RATIONAL(1.0,2.0);
      t41 = x*t35;
      t38 = y*y;
      t29 = t38*t41;
      t36 = RATIONAL(1.0,4.0);
      t40 = t29+t36*t38;
      t37 = RATIONAL(-1.0,4.0);
      t39 = t29+t37*t38;
      t34 = RATIONAL(-1.0,2.0);
      t33 = t36*y;
      t30 = t37*y;
      t28 = t34*x*y;
      t27 = y*t41;
      coeffs_dx->coeff_m1_m1 = t33+t28+t39;
      coeffs_dx->coeff_0_m1 = (y-t38)*x;
      coeffs_dx->coeff_p1_m1 = t28+t30+t40;
      coeffs_dx->coeff_m1_0 = x+t34+(t35-x)*t38;
      coeffs_dx->coeff_0_0 = (RATIONAL(-2.0,1.0)+RATIONAL(2.0,1.0)*t38)*x;
      coeffs_dx->coeff_p1_0 = x+t35+(t34-x)*t38;
      coeffs_dx->coeff_m1_p1 = t30+t27+t39;
      coeffs_dx->coeff_0_p1 = (-y-t38)*x;
      coeffs_dx->coeff_p1_p1 = t27+t33+t40;