src/Lagrange-tensor-product/3d.coeffs/3d.cube.order1.smooth0/coeffs-I.compute.c


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

fp t2;
fp t8;
fp t1;
fp t7;
fp t6;
fp t4;
fp t5;
      t2 = y*z;
      t8 = z-t2;
      t1 = x*t2;
      t7 = t2-t1;
      t6 = x*z-t1;
      t4 = x*y;
      t5 = x-t4-t6;
      coeffs_I->coeff_0_0_0 = RATIONAL(1.0,1.0)-y-t5-t8;
      coeffs_I->coeff_p1_0_0 = t5;
      coeffs_I->coeff_0_p1_0 = y-t4-t7;
      coeffs_I->coeff_p1_p1_0 = -t1+t4;
      coeffs_I->coeff_0_0_p1 = -t6+t8;
      coeffs_I->coeff_p1_0_p1 = t6;
      coeffs_I->coeff_0_p1_p1 = t7;
      coeffs_I->coeff_p1_p1_p1 = t1;