aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxx.compute.c
blob: e8ab5c4d2c4dc2341b5a789764243cff9c4c44a7 (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
117
118
119
120
fp t465;
fp t485;
fp t464;
fp t484;
fp t467;
fp t483;
fp t468;
fp t482;
fp t460;
fp t481;
fp t461;
fp t480;
fp t466;
fp t479;
fp t478;
fp t477;
fp t473;
fp t448;
fp t476;
fp t459;
fp t458;
fp t452;
fp t453;
fp t475;
fp t474;
fp t444;
fp t472;
fp t471;
fp t470;
fp t469;
fp t455;
fp t454;
fp t451;
fp t450;
fp t447;
fp t446;
fp t445;
fp t443;
fp t442;
fp t441;
fp t440;
fp t439;
fp t438;
fp t437;
fp t436;
fp t435;
      t465 = RATIONAL(-1.0,2.0);
      t485 = t465*y;
      t464 = RATIONAL(1.0,2.0);
      t484 = t464*y;
      t467 = z*z;
      t483 = y*t467;
      t468 = y*y;
      t482 = t468*z;
      t460 = RATIONAL(1.0,4.0);
      t481 = y*t460;
      t461 = RATIONAL(-1.0,4.0);
      t480 = z*t461;
      t466 = t467*t468;
      t479 = -t468+t466;
      t478 = t466-t467;
      t477 = t468+t467;
      t473 = t460*t468;
      t448 = t467*t473;
      t476 = z*t473+t448;
      t459 = t464*t467;
      t458 = t464*t468;
      t452 = z*t458;
      t453 = t465*t466;
      t475 = t452+t453;
      t474 = t448+t468*t480;
      t444 = t465*t483;
      t472 = t444+t453;
      t471 = t453+t465*t482;
      t470 = t482-z;
      t469 = t483-y;
      t455 = t461*t483;
      t454 = y*t459;
      t451 = z*t485;
      t450 = y*t480;
      t447 = t467*t481;
      t446 = z*t484;
      t445 = z*t481;
      t443 = t466+RATIONAL(1.0,1.0)-t477;
      t442 = t464*z+t459+t471;
      t441 = t458+t484+t472;
      t440 = t465*z+t459+t475;
      t439 = t485+t454+t453+t458;
      t438 = t447+t445+t476;
      t437 = t455+t450+t476;
      t436 = t455+t445+t474;
      t435 = t450+t447+t474;
      coeffs_dxx->coeff_m1_m1_m1 = t436;
      coeffs_dxx->coeff_0_m1_m1 = t451+t454+t475;
      coeffs_dxx->coeff_p1_m1_m1 = t436;
      coeffs_dxx->coeff_m1_0_m1 = t440;
      coeffs_dxx->coeff_0_0_m1 = -t470+t478;
      coeffs_dxx->coeff_p1_0_m1 = t440;
      coeffs_dxx->coeff_m1_p1_m1 = t435;
      coeffs_dxx->coeff_0_p1_m1 = t452+t446+t472;
      coeffs_dxx->coeff_p1_p1_m1 = t435;
      coeffs_dxx->coeff_m1_m1_0 = t439;
      coeffs_dxx->coeff_0_m1_0 = -t469+t479;
      coeffs_dxx->coeff_p1_m1_0 = t439;
      coeffs_dxx->coeff_m1_0_0 = t443;
      coeffs_dxx->coeff_0_0_0 = (t466+1.0)*RATIONAL(-2.0,1.0)+t477*RATIONAL(2.0
,1.0);
      coeffs_dxx->coeff_p1_0_0 = t443;
      coeffs_dxx->coeff_m1_p1_0 = t441;
      coeffs_dxx->coeff_0_p1_0 = t469+t479;
      coeffs_dxx->coeff_p1_p1_0 = t441;
      coeffs_dxx->coeff_m1_m1_p1 = t437;
      coeffs_dxx->coeff_0_m1_p1 = t446+t454+t471;
      coeffs_dxx->coeff_p1_m1_p1 = t437;
      coeffs_dxx->coeff_m1_0_p1 = t442;
      coeffs_dxx->coeff_0_0_p1 = t470+t478;
      coeffs_dxx->coeff_p1_0_p1 = t442;
      coeffs_dxx->coeff_m1_p1_p1 = t438;
      coeffs_dxx->coeff_0_p1_p1 = t444+t451+t471;
      coeffs_dxx->coeff_p1_p1_p1 = t438;