src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c


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 t50;
fp t56;
fp t53;
fp t43;
fp t51;
fp t55;
fp t52;
fp t54;
fp t49;
fp t47;
fp t45;
fp t44;
fp t42;
      t50 = RATIONAL(1.0,2.0);
      t56 = y*t50;
      t53 = x*x;
      t43 = t53*t56;
      t51 = RATIONAL(1.0,4.0);
      t55 = t43+t51*x;
      t52 = RATIONAL(-1.0,4.0);
      t54 = t43+t52*x;
      t49 = RATIONAL(-1.0,2.0);
      t47 = t52*t53;
      t45 = t51*t53;
      t44 = x*t56;
      t42 = t49*x*y;
      coeffs_dy->coeff_m1_m1 = t42+t47+t55;
      coeffs_dy->coeff_0_m1 = t49+y+(-y+t50)*t53;
      coeffs_dy->coeff_p1_m1 = t44+t47+t54;
      coeffs_dy->coeff_m1_0 = (x-t53)*y;
      coeffs_dy->coeff_0_0 = (RATIONAL(-2.0,1.0)+RATIONAL(2.0,1.0)*t53)*y;
      coeffs_dy->coeff_p1_0 = (-x-t53)*y;
      coeffs_dy->coeff_m1_p1 = t45+t42+t54;
      coeffs_dy->coeff_0_p1 = t50+y+(t49-y)*t53;
      coeffs_dy->coeff_p1_p1 = t44+t45+t55;