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