aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyy.compute.c
blob: 1033b59a3796f3db8b4d7ca9c6eb7575d4d1e4b7 (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 t680;
fp t658;
fp t679;
fp t660;
fp t678;
fp t664;
fp t677;
fp t663;
fp t676;
fp t654;
fp t655;
fp t662;
fp t657;
fp t640;
fp t675;
fp t659;
fp t674;
fp t650;
fp t673;
fp t672;
fp t671;
fp t670;
fp t669;
fp t668;
fp t667;
fp t666;
fp t665;
fp t652;
fp t651;
fp t648;
fp t647;
fp t645;
fp t644;
fp t642;
fp t641;
fp t639;
fp t638;
fp t637;
fp t636;
fp t635;
fp t634;
fp t633;
fp t632;
fp t631;
      t680 = x*z;
      t658 = RATIONAL(1.0,2.0);
      t679 = t658*z;
      t660 = RATIONAL(-1.0,2.0);
      t678 = t660*z;
      t664 = x*x;
      t677 = t664*z;
      t663 = z*z;
      t676 = x*t663;
      t654 = t658*t663;
      t655 = t658*t664;
      t662 = t664*t663;
      t657 = RATIONAL(1.0,4.0);
      t640 = t657*t662;
      t675 = t640+t657*t680;
      t659 = RATIONAL(-1.0,4.0);
      t674 = t640+t659*t680;
      t650 = t660*t662;
      t673 = x*t679+t650;
      t672 = t650+x*t678;
      t671 = t650+t654;
      t670 = t650+t655;
      t669 = t662-t663;
      t668 = -t664+t662;
      t667 = t664+t663;
      t666 = t677-z;
      t665 = t676-x;
      t652 = t660*t676;
      t651 = t660*t677;
      t648 = t659*t676;
      t647 = t659*t677;
      t645 = x*t654;
      t644 = z*t655;
      t642 = t657*t676;
      t641 = t657*t677;
      t639 = t662+RATIONAL(1.0,1.0)-t667;
      t638 = t678+t644+t671;
      t637 = t658*x+t652+t670;
      t636 = t651+t679+t671;
      t635 = t645+t660*x+t670;
      t634 = t647+t642+t674;
      t633 = t648+t641+t674;
      t632 = t647+t648+t675;
      t631 = t641+t642+t675;
      coeffs_dyy->coeff_m1_m1_m1 = t632;
      coeffs_dyy->coeff_0_m1_m1 = t638;
      coeffs_dyy->coeff_p1_m1_m1 = t634;
      coeffs_dyy->coeff_m1_0_m1 = t645+t644+t672;
      coeffs_dyy->coeff_0_0_m1 = -t666+t669;
      coeffs_dyy->coeff_p1_0_m1 = t644+t652+t673;
      coeffs_dyy->coeff_m1_p1_m1 = t632;
      coeffs_dyy->coeff_0_p1_m1 = t638;
      coeffs_dyy->coeff_p1_p1_m1 = t634;
      coeffs_dyy->coeff_m1_m1_0 = t635;
      coeffs_dyy->coeff_0_m1_0 = t639;
      coeffs_dyy->coeff_p1_m1_0 = t637;
      coeffs_dyy->coeff_m1_0_0 = -t665+t668;
      coeffs_dyy->coeff_0_0_0 = (1.0+t662)*RATIONAL(-2.0,1.0)+t667*RATIONAL(2.0
,1.0);
      coeffs_dyy->coeff_p1_0_0 = t665+t668;
      coeffs_dyy->coeff_m1_p1_0 = t635;
      coeffs_dyy->coeff_0_p1_0 = t639;
      coeffs_dyy->coeff_p1_p1_0 = t637;
      coeffs_dyy->coeff_m1_m1_p1 = t633;
      coeffs_dyy->coeff_0_m1_p1 = t636;
      coeffs_dyy->coeff_p1_m1_p1 = t631;
      coeffs_dyy->coeff_m1_0_p1 = t651+t645+t673;
      coeffs_dyy->coeff_0_0_p1 = t666+t669;
      coeffs_dyy->coeff_p1_0_p1 = t651+t652+t672;
      coeffs_dyy->coeff_m1_p1_p1 = t633;
      coeffs_dyy->coeff_0_p1_p1 = t636;
      coeffs_dyy->coeff_p1_p1_p1 = t631;