aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c
blob: 56cfd488b95f71f2f85b18f55f3e36cca9953c6d (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
fp t435;
fp t424;
fp t439;
fp t409;
fp t460;
fp t438;
fp t416;
fp t459;
fp t451;
fp t429;
fp t445;
fp t436;
fp t403;
fp t458;
fp t423;
fp t408;
fp t457;
fp t456;
fp t437;
fp t455;
fp t454;
fp t427;
fp t453;
fp t405;
fp t452;
fp t421;
fp t422;
fp t390;
fp t450;
fp t449;
fp t448;
fp t434;
fp t447;
fp t446;
fp t425;
fp t444;
fp t443;
fp t426;
fp t391;
fp t442;
fp t441;
fp t420;
fp t440;
fp t433;
fp t432;
fp t431;
fp t430;
fp t428;
fp t394;
fp t389;
      t435 = RATIONAL(-1.0,6.0);
      t424 = RATIONAL(1.0,2.0);
      t439 = y*y;
      t409 = t424*t439;
      t460 = t435+t409;
      t438 = x*x;
      t416 = t424*t438;
      t459 = t416+t435;
      t451 = t439+t438;
      t429 = RATIONAL(-1.0,3.0);
      t445 = t439*t438;
      t436 = RATIONAL(3.0,4.0);
      t403 = t436*t445;
      t458 = t429+t403;
      t423 = RATIONAL(-1.0,2.0);
      t408 = t423*t438;
      t457 = (t429+t408)*y;
      t456 = t423+RATIONAL(-9.0,4.0)*t445;
      t437 = x*y;
      t455 = t438*y;
      t454 = y+RATIONAL(9.0,4.0)*t445;
      t427 = RATIONAL(3.0,2.0);
      t453 = t427*t439;
      t405 = t423*t439;
      t452 = t427*t438;
      t421 = RATIONAL(1.0,4.0);
      t422 = RATIONAL(-3.0,4.0);
      t390 = t422*t445;
      t450 = t421*t439+t390;
      t449 = t459*y;
      t448 = t460*x;
      t434 = RATIONAL(-1.0,4.0);
      t447 = RATIONAL(1.0,18.0)+t434*t445;
      t446 = t421*t438+t390;
      t425 = RATIONAL(-3.0,2.0);
      t444 = x+t437+t425*x*t439;
      t443 = t437+(t429+t405)*x;
      t426 = RATIONAL(-2.0,1.0);
      t391 = t426*t437;
      t442 = t391+(t423+t452)*y;
      t441 = t391+(t453+t423)*x;
      t420 = RATIONAL(1.0,12.0);
      t440 = t420+t403+t451*t434;
      t433 = RATIONAL(2.0,3.0);
      t432 = RATIONAL(1.0,6.0);
      t431 = RATIONAL(3.0,1.0);
      t430 = RATIONAL(1.0,3.0);
      t428 = RATIONAL(-3.0,1.0);
      t394 = t421*t445;
      t389 = t425*t455;
      coeffs_dxy->coeff_m1_m1 = t394+RATIONAL(1.0,9.0)+t451*t432+t443+t457;
      coeffs_dxy->coeff_0_m1 = t408+t432+(t433+t439)*x+t442+t450;
      coeffs_dxy->coeff_p1_m1 = t405+t389+y+t416+t443+t458;
      coeffs_dxy->coeff_p2_m1 = t420*t439+t435*t438+t447+t449;
      coeffs_dxy->coeff_m1_0 = t405+t432+(t433+t438)*y+t441+t446;
      coeffs_dxy->coeff_0_0 = t428*t455+t421+t451*t422+(t428*t439+1.0+RATIONAL(
4.0,1.0)*y)*x+t454;
      coeffs_dxy->coeff_p1_0 = t436*t438+t453+(t431*t438+t426)*y+t441+t456;
      coeffs_dxy->coeff_p2_0 = (-t438+t430)*y+t440;
      coeffs_dxy->coeff_m1_p1 = t408+t409+t444+t457+t458;
      coeffs_dxy->coeff_0_p1 = t452+t436*t439+(t426+t431*t439)*x+t442+t456;
      coeffs_dxy->coeff_p1_p1 = t389+RATIONAL(1.0,1.0)+t451*t425+t444+t454;
      coeffs_dxy->coeff_p2_p1 = t449+t450+t459;
      coeffs_dxy->coeff_m1_p2 = t420*t438+t435*t439+t447+t448;
      coeffs_dxy->coeff_0_p2 = (-t439+t430)*x+t440;
      coeffs_dxy->coeff_p1_p2 = t446+t448+t460;
      coeffs_dxy->coeff_p2_p2 = RATIONAL(1.0,36.0)+t394+t451*RATIONAL(-1.0,12.0
);