aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxx.compute.c
blob: b04274827653ae6e6af14eb12cb007ec75b1a246 (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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
fp t480;
fp t458;
fp t479;
fp t460;
fp t478;
fp t464;
fp t477;
fp t463;
fp t476;
fp t454;
fp t455;
fp t462;
fp t457;
fp t440;
fp t475;
fp t459;
fp t474;
fp t450;
fp t473;
fp t472;
fp t471;
fp t470;
fp t469;
fp t468;
fp t467;
fp t466;
fp t465;
fp t452;
fp t451;
fp t448;
fp t447;
fp t445;
fp t444;
fp t442;
fp t441;
fp t439;
fp t438;
fp t437;
fp t436;
fp t435;
fp t434;
fp t433;
fp t432;
fp t431;
      t480 = y*z;
      t458 = RATIONAL(1.0,2.0);
      t479 = t458*z;
      t460 = RATIONAL(-1.0,2.0);
      t478 = t460*z;
      t464 = y*y;
      t477 = t464*z;
      t463 = z*z;
      t476 = y*t463;
      t454 = t458*t463;
      t455 = t458*t464;
      t462 = t464*t463;
      t457 = RATIONAL(1.0,4.0);
      t440 = t457*t462;
      t475 = t440+t457*t480;
      t459 = RATIONAL(-1.0,4.0);
      t474 = t440+t459*t480;
      t450 = t460*t462;
      t473 = y*t479+t450;
      t472 = t450+y*t478;
      t471 = t450+t454;
      t470 = t450+t455;
      t469 = t462-t463;
      t468 = t462-t464;
      t467 = t464+t463;
      t466 = y-t476;
      t465 = z-t477;
      t452 = t460*t476;
      t451 = t460*t477;
      t448 = t459*t476;
      t447 = t459*t477;
      t445 = y*t454;
      t444 = z*t455;
      t442 = t457*t476;
      t441 = t457*t477;
      t439 = RATIONAL(1.0,1.0)+t462-t467;
      t438 = t458*y+t452+t470;
      t437 = t478+t444+t471;
      t436 = t479+t451+t471;
      t435 = t445+t460*y+t470;
      t434 = t447+t442+t474;
      t433 = t447+t448+t475;
      t432 = t442+t441+t475;
      t431 = t441+t448+t474;
      coeffs_dxx->coeff_m1_m1_m1 = t433;
      coeffs_dxx->coeff_0_m1_m1 = t444+t445+t472;
      coeffs_dxx->coeff_p1_m1_m1 = t433;
      coeffs_dxx->coeff_m1_0_m1 = t437;
      coeffs_dxx->coeff_0_0_m1 = t465+t469;
      coeffs_dxx->coeff_p1_0_m1 = t437;
      coeffs_dxx->coeff_m1_p1_m1 = t434;
      coeffs_dxx->coeff_0_p1_m1 = t444+t452+t473;
      coeffs_dxx->coeff_p1_p1_m1 = t434;
      coeffs_dxx->coeff_m1_m1_0 = t435;
      coeffs_dxx->coeff_0_m1_0 = t466+t468;
      coeffs_dxx->coeff_p1_m1_0 = t435;
      coeffs_dxx->coeff_m1_0_0 = t439;
      coeffs_dxx->coeff_0_0_0 = (t462+1.0)*RATIONAL(-2.0,1.0)+t467*RATIONAL(2.0
,1.0);
      coeffs_dxx->coeff_p1_0_0 = t439;
      coeffs_dxx->coeff_m1_p1_0 = t438;
      coeffs_dxx->coeff_0_p1_0 = -t466+t468;
      coeffs_dxx->coeff_p1_p1_0 = t438;
      coeffs_dxx->coeff_m1_m1_p1 = t431;
      coeffs_dxx->coeff_0_m1_p1 = t451+t445+t473;
      coeffs_dxx->coeff_p1_m1_p1 = t431;
      coeffs_dxx->coeff_m1_0_p1 = t436;
      coeffs_dxx->coeff_0_0_p1 = -t465+t469;
      coeffs_dxx->coeff_p1_0_p1 = t436;
      coeffs_dxx->coeff_m1_p1_p1 = t432;
      coeffs_dxx->coeff_0_p1_p1 = t452+t451+t472;
      coeffs_dxx->coeff_p1_p1_p1 = t432;