aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
blob: fef7036fd016001a0d67c7c5b049c26b13ba59d8 (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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
fp t374;
fp t375;
fp t376;
fp t388;
fp t362;
fp t358;
fp t369;
fp t361;
fp t373;
fp t356;
fp t370;
fp t387;
fp t386;
fp t357;
fp t348;
fp t363;
fp t385;
fp t372;
fp t355;
fp t384;
fp t354;
fp t345;
fp t383;
fp t367;
fp t382;
fp t360;
fp t359;
fp t381;
fp t380;
fp t379;
fp t378;
fp t371;
fp t368;
fp t366;
fp t365;
fp t353;
fp t352;
fp t351;
fp t350;
fp t346;
fp t344;
fp t343;
      t374 = x*y;
      t375 = y*y;
      t376 = t375*y;
      t388 = x*t376;
      t362 = RATIONAL(-1.0,2.0);
      t358 = t362*y;
      t369 = RATIONAL(-1.0,3.0);
      t361 = t369*y;
      t373 = RATIONAL(-1.0,6.0);
      t356 = t373*y;
      t370 = RATIONAL(1.0,3.0);
      t387 = t370*y;
      t386 = x*t375;
      t357 = t362*t376;
      t348 = x*t357;
      t363 = RATIONAL(1.0,2.0);
      t385 = t363*t374+t348;
      t372 = RATIONAL(1.0,6.0);
      t355 = t372*t376;
      t384 = t355+t356;
      t354 = t363*t376;
      t345 = x*t354;
      t383 = t345+x*t358;
      t367 = RATIONAL(3.0,2.0);
      t382 = t367*t386-t375;
      t360 = t373*t376;
      t359 = t363*t375;
      t381 = y+t357+t359;
      t380 = t386-x;
      t379 = t359+t361+t360;
      t378 = t354+t358+RATIONAL(1.0,1.0)-t375;
      t371 = RATIONAL(3.0,1.0);
      t368 = RATIONAL(-3.0,1.0);
      t366 = RATIONAL(-2.0,1.0);
      t365 = RATIONAL(-3.0,2.0);
      t353 = t367*t388;
      t352 = t365*t388;
      t351 = x*t355;
      t350 = t362*t386;
      t346 = t365*t386;
      t344 = x*t360;
      t343 = x*t359;
      coeffs_dxx->coeff_m1_m1 = t350+t351+x*t387+t379;
      coeffs_dxx->coeff_0_m1 = -t374+t370*t376+t348+RATIONAL(2.0,3.0)*y+t382;
      coeffs_dxx->coeff_p1_m1 = t345+t374+t346+t379;
      coeffs_dxx->coeff_p2_m1 = x*t361+t344+t343;
      coeffs_dxx->coeff_m1_0 = t378+t380+t385;
      coeffs_dxx->coeff_0_0 = -t376+y+RATIONAL(2.0,1.0)*t375+t366+t353+(t365*y+
t371+t368*t375)*x;
      coeffs_dxx->coeff_p1_0 = t352+(t371*t375+t368+t367*y)*x+t378;
      coeffs_dxx->coeff_p2_0 = -t380+t383;
      coeffs_dxx->coeff_m1_p1 = -t374+t350+t345+t381;
      coeffs_dxx->coeff_0_p1 = t352+t376+(t371*x+t366)*y+t382;
      coeffs_dxx->coeff_p1_p1 = t353+t346+t368*t374+t381;
      coeffs_dxx->coeff_p2_p1 = t374+t343+t348;
      coeffs_dxx->coeff_m1_p2 = t372*t374+t344+t384;
      coeffs_dxx->coeff_0_p2 = t369*t376+t387+t383;
      coeffs_dxx->coeff_p1_p2 = t384+t385;
      coeffs_dxx->coeff_p2_p2 = t351+x*t356;