aboutsummaryrefslogtreecommitdiff
path: root/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
blob: a6f79fdcfccbb8fda33ddb7e7a01be4c13798631 (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
fp t493;
fp t494;
fp t508;
fp t481;
fp t507;
fp t488;
fp t477;
fp t495;
fp t506;
fp t489;
fp t505;
fp t492;
fp t480;
fp t482;
fp t476;
fp t473;
fp t491;
fp t478;
fp t475;
fp t479;
fp t486;
fp t504;
fp t503;
fp t472;
fp t502;
fp t462;
fp t501;
fp t500;
fp t499;
fp t498;
fp t497;
fp t490;
fp t487;
fp t485;
fp t484;
fp t471;
fp t469;
fp t468;
fp t466;
fp t465;
fp t464;
fp t463;
      t493 = x*y;
      t494 = x*x;
      t508 = t494*y;
      t481 = RATIONAL(-1.0,2.0);
      t507 = y*t481;
      t488 = RATIONAL(-1.0,3.0);
      t477 = t488*x;
      t495 = x*t494;
      t506 = t495*y;
      t489 = RATIONAL(1.0,3.0);
      t505 = t489*x;
      t492 = RATIONAL(-1.0,6.0);
      t480 = t492*x;
      t482 = RATIONAL(1.0,2.0);
      t476 = t482*t494;
      t473 = t482*t495;
      t491 = RATIONAL(1.0,6.0);
      t478 = t491*t495;
      t475 = t492*t495;
      t479 = t481*t495;
      t486 = RATIONAL(3.0,2.0);
      t504 = -t494+t486*t508;
      t503 = t480+t478;
      t472 = y*t473;
      t502 = x*t507+t472;
      t462 = y*t479;
      t501 = t462+t482*t493;
      t500 = x+t479+t476;
      t499 = t508-y;
      t498 = t475+t477+t476;
      t497 = -t494+RATIONAL(1.0,1.0)+t473+t481*x;
      t490 = RATIONAL(3.0,1.0);
      t487 = RATIONAL(-3.0,1.0);
      t485 = RATIONAL(-2.0,1.0);
      t484 = RATIONAL(-3.0,2.0);
      t471 = y*t475;
      t469 = t486*t506;
      t468 = t484*t506;
      t466 = t484*t508;
      t465 = t494*t507;
      t464 = y*t478;
      t463 = y*t476;
      coeffs_dyy->coeff_m1_m1 = t464+t465+y*t505+t498;
      coeffs_dyy->coeff_0_m1 = t497+t499+t501;
      coeffs_dyy->coeff_p1_m1 = t465-t493+t472+t500;
      coeffs_dyy->coeff_p2_m1 = t491*t493+t471+t503;
      coeffs_dyy->coeff_m1_0 = RATIONAL(2.0,3.0)*x-t493+t462+t489*t495+t504;
      coeffs_dyy->coeff_0_0 = x+RATIONAL(2.0,1.0)*t494+t485+t469-t495+(t484*x+
t490+t487*t494)*y;
      coeffs_dyy->coeff_p1_0 = t495+t468+(t485+t490*y)*x+t504;
      coeffs_dyy->coeff_p2_0 = t505+t488*t495+t502;
      coeffs_dyy->coeff_m1_p1 = t472+t493+t466+t498;
      coeffs_dyy->coeff_0_p1 = t468+(t490*t494+t487+t486*x)*y+t497;
      coeffs_dyy->coeff_p1_p1 = t466+t469+t487*t493+t500;
      coeffs_dyy->coeff_p2_p1 = t501+t503;
      coeffs_dyy->coeff_m1_p2 = y*t477+t471+t463;
      coeffs_dyy->coeff_0_p2 = -t499+t502;
      coeffs_dyy->coeff_p1_p2 = t463+t493+t462;
      coeffs_dyy->coeff_p2_p2 = y*t480+t464;