src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.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

fp t72;
fp t68;
fp t74;
fp t69;
fp t73;
fp t71;
fp t70;
fp t65;
fp t64;
fp t63;
      t72 = x*y;
      t68 = RATIONAL(-1.0,2.0);
      t74 = t68*x+t72;
      t69 = RATIONAL(1.0,2.0);
      t73 = t69*x+t72;
      t71 = RATIONAL(-1.0,4.0);
      t70 = RATIONAL(1.0,4.0);
      t65 = t68*y;
      t64 = t69*y;
      t63 = RATIONAL(-2.0,1.0)*t72;
      coeffs_dxy->coeff_m1_m1 = t65+t70+t74;
      coeffs_dxy->coeff_0_m1 = t63+x;
      coeffs_dxy->coeff_p1_m1 = t71+t64+t74;
      coeffs_dxy->coeff_m1_0 = t63+y;
      coeffs_dxy->coeff_0_0 = RATIONAL(4.0,1.0)*t72;
      coeffs_dxy->coeff_p1_0 = -y+t63;
      coeffs_dxy->coeff_m1_p1 = t65+t71+t73;
      coeffs_dxy->coeff_0_p1 = t63-x;
      coeffs_dxy->coeff_p1_p1 = t70+t64+t73;