aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c
blob: a9bdd0380c18d9ec315adb52be50b8a9a76a854a (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
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
fp t17;
fp t15;
fp t18;
fp t19;
fp t4;
fp t16;
fp t24;
fp t26;
fp t25;
fp t14;
fp t12;
fp t13;
fp t20;
fp t9;
fp t23;
fp t11;
fp t22;
fp t21;
fp t8;
fp t6;
fp t5;
fp t3;
      t17 = x*x;
      t15 = RATIONAL(1.0,4.0);
      t18 = y*y;
      t19 = t15*t18;
      t4 = t17*t19;
      t16 = RATIONAL(-1.0,4.0);
      t24 = y*t16;
      t26 = t4+t17*t24;
      t25 = y*t15;
      t14 = RATIONAL(1.0,2.0);
      t12 = t14*t18;
      t13 = RATIONAL(-1.0,2.0);
      t20 = t13*t18;
      t9 = t17*t20;
      t23 = t12+t9;
      t11 = t14*t17;
      t22 = t11+t9;
      t21 = t4+t17*t25;
      t8 = x*t19;
      t6 = x*t24;
      t5 = x*t25;
      t3 = t16*x*t18;
      coeffs_I->coeff_m1_m1 = t3+t5+t26;
      coeffs_I->coeff_0_m1 = (t13+t11)*y+t23;
      coeffs_I->coeff_p1_m1 = t8+t6+t26;
      coeffs_I->coeff_m1_0 = (t12+t13)*x+t22;
      coeffs_I->coeff_0_0 = -t18+t17*t18+RATIONAL(1.0,1.0)-t17;
      coeffs_I->coeff_p1_0 = (t20+t14)*x+t22;
      coeffs_I->coeff_m1_p1 = t6+t3+t21;
      coeffs_I->coeff_0_p1 = (t13*t17+t14)*y+t23;
      coeffs_I->coeff_p1_p1 = t5+t8+t21;