PLaSK library
Loading...
Searching...
No Matches
diffusion3d-eval-shb.ipp
Go to the documentation of this file.
1K(e.i00, e.i00) += (1.0/1225.0)*e.X*e.Y*(100*P[e.n00].c00*dG.c00 + 100*P[e.n00].c11*dG.c11 + 30*P[e.n01].c00*dG.c00 + 30*P[e.n01].c11*dG.c11 + 30*P[e.n10].c00*dG.c00 + 30*P[e.n10].c11*dG.c11 + 9*P[e.n11].c00*dG.c00 + 9*P[e.n11].c11*dG.c11);
2K(e.i00, e.i01) += (1.0/14700.0)*e.X*(e.Y*e.Y)*(150*P[e.n00].c00*dG.c00 + 150*P[e.n00].c11*dG.c11 + 70*P[e.n01].c00*dG.c00 + 70*P[e.n01].c11*dG.c11 + 45*P[e.n10].c00*dG.c00 + 45*P[e.n10].c11*dG.c11 + 21*P[e.n11].c00*dG.c00 + 21*P[e.n11].c11*dG.c11);
3K(e.i01, e.i01) += (1.0/29400.0)*e.X*(e.Y*e.Y*e.Y)*(50*P[e.n00].c00*dG.c00 + 50*P[e.n00].c11*dG.c11 + 30*P[e.n01].c00*dG.c00 + 30*P[e.n01].c11*dG.c11 + 15*P[e.n10].c00*dG.c00 + 15*P[e.n10].c11*dG.c11 + 9*P[e.n11].c00*dG.c00 + 9*P[e.n11].c11*dG.c11);
4K(e.i00, e.i10) += (1.0/14700.0)*(e.X*e.X)*e.Y*(150*P[e.n00].c00*dG.c00 + 150*P[e.n00].c11*dG.c11 + 45*P[e.n01].c00*dG.c00 + 45*P[e.n01].c11*dG.c11 + 70*P[e.n10].c00*dG.c00 + 70*P[e.n10].c11*dG.c11 + 21*P[e.n11].c00*dG.c00 + 21*P[e.n11].c11*dG.c11);
5K(e.i01, e.i10) += (e.X*e.X)*(e.Y*e.Y)*((1.0/784.0)*P[e.n00].c00*dG.c00 + (1.0/784.0)*P[e.n00].c11*dG.c11 + (1.0/1680.0)*P[e.n01].c00*dG.c00 + (1.0/1680.0)*P[e.n01].c11*dG.c11 + (1.0/1680.0)*P[e.n10].c00*dG.c00 + (1.0/1680.0)*P[e.n10].c11*dG.c11 + (1.0/3600.0)*P[e.n11].c00*dG.c00 + (1.0/3600.0)*P[e.n11].c11*dG.c11);
6K(e.i10, e.i10) += (1.0/29400.0)*(e.X*e.X*e.X)*e.Y*(50*P[e.n00].c00*dG.c00 + 50*P[e.n00].c11*dG.c11 + 15*P[e.n01].c00*dG.c00 + 15*P[e.n01].c11*dG.c11 + 30*P[e.n10].c00*dG.c00 + 30*P[e.n10].c11*dG.c11 + 9*P[e.n11].c00*dG.c00 + 9*P[e.n11].c11*dG.c11);
7K(e.i00, e.i02) += (9.0/4900.0)*e.X*e.Y*(10*P[e.n00].c00*dG.c00 + 10*P[e.n00].c11*dG.c11 + 10*P[e.n01].c00*dG.c00 + 10*P[e.n01].c11*dG.c11 + 3*P[e.n10].c00*dG.c00 + 3*P[e.n10].c11*dG.c11 + 3*P[e.n11].c00*dG.c00 + 3*P[e.n11].c11*dG.c11);
8K(e.i01, e.i02) += (1.0/14700.0)*e.X*(e.Y*e.Y)*(60*P[e.n00].c00*dG.c00 + 60*P[e.n00].c11*dG.c11 + 70*P[e.n01].c00*dG.c00 + 70*P[e.n01].c11*dG.c11 + 18*P[e.n10].c00*dG.c00 + 18*P[e.n10].c11*dG.c11 + 21*P[e.n11].c00*dG.c00 + 21*P[e.n11].c11*dG.c11);
9K(e.i10, e.i02) += (3.0/19600.0)*(e.X*e.X)*e.Y*(15*P[e.n00].c00*dG.c00 + 15*P[e.n00].c11*dG.c11 + 15*P[e.n01].c00*dG.c00 + 15*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
10K(e.i02, e.i02) += (1.0/1225.0)*e.X*e.Y*(30*P[e.n00].c00*dG.c00 + 30*P[e.n00].c11*dG.c11 + 100*P[e.n01].c00*dG.c00 + 100*P[e.n01].c11*dG.c11 + 9*P[e.n10].c00*dG.c00 + 9*P[e.n10].c11*dG.c11 + 30*P[e.n11].c00*dG.c00 + 30*P[e.n11].c11*dG.c11);
11K(e.i00, e.i03) += -1.0/14700.0*e.X*(e.Y*e.Y)*(70*P[e.n00].c00*dG.c00 + 70*P[e.n00].c11*dG.c11 + 60*P[e.n01].c00*dG.c00 + 60*P[e.n01].c11*dG.c11 + 21*P[e.n10].c00*dG.c00 + 21*P[e.n10].c11*dG.c11 + 18*P[e.n11].c00*dG.c00 + 18*P[e.n11].c11*dG.c11);
12K(e.i01, e.i03) += -1.0/9800.0*e.X*(e.Y*e.Y*e.Y)*(10*P[e.n00].c00*dG.c00 + 10*P[e.n00].c11*dG.c11 + 10*P[e.n01].c00*dG.c00 + 10*P[e.n01].c11*dG.c11 + 3*P[e.n10].c00*dG.c00 + 3*P[e.n10].c11*dG.c11 + 3*P[e.n11].c00*dG.c00 + 3*P[e.n11].c11*dG.c11);
13K(e.i10, e.i03) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/1680.0*P[e.n00].c00*dG.c00 - 1.0/1680.0*P[e.n00].c11*dG.c11 - 1.0/1960.0*P[e.n01].c00*dG.c00 - 1.0/1960.0*P[e.n01].c11*dG.c11 - 1.0/3600.0*P[e.n10].c00*dG.c00 - 1.0/3600.0*P[e.n10].c11*dG.c11 - 1.0/4200.0*P[e.n11].c00*dG.c00 - 1.0/4200.0*P[e.n11].c11*dG.c11);
14K(e.i02, e.i03) += -1.0/14700.0*e.X*(e.Y*e.Y)*(70*P[e.n00].c00*dG.c00 + 70*P[e.n00].c11*dG.c11 + 150*P[e.n01].c00*dG.c00 + 150*P[e.n01].c11*dG.c11 + 21*P[e.n10].c00*dG.c00 + 21*P[e.n10].c11*dG.c11 + 45*P[e.n11].c00*dG.c00 + 45*P[e.n11].c11*dG.c11);
15K(e.i03, e.i03) += (1.0/29400.0)*e.X*(e.Y*e.Y*e.Y)*(30*P[e.n00].c00*dG.c00 + 30*P[e.n00].c11*dG.c11 + 50*P[e.n01].c00*dG.c00 + 50*P[e.n01].c11*dG.c11 + 9*P[e.n10].c00*dG.c00 + 9*P[e.n10].c11*dG.c11 + 15*P[e.n11].c00*dG.c00 + 15*P[e.n11].c11*dG.c11);
16K(e.i00, e.i12) += (3.0/19600.0)*(e.X*e.X)*e.Y*(15*P[e.n00].c00*dG.c00 + 15*P[e.n00].c11*dG.c11 + 15*P[e.n01].c00*dG.c00 + 15*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
17K(e.i01, e.i12) += (e.X*e.X)*(e.Y*e.Y)*((1.0/1960.0)*P[e.n00].c00*dG.c00 + (1.0/1960.0)*P[e.n00].c11*dG.c11 + (1.0/1680.0)*P[e.n01].c00*dG.c00 + (1.0/1680.0)*P[e.n01].c11*dG.c11 + (1.0/4200.0)*P[e.n10].c00*dG.c00 + (1.0/4200.0)*P[e.n10].c11*dG.c11 + (1.0/3600.0)*P[e.n11].c00*dG.c00 + (1.0/3600.0)*P[e.n11].c11*dG.c11);
18K(e.i10, e.i12) += (3.0/39200.0)*(e.X*e.X*e.X)*e.Y*(5*P[e.n00].c00*dG.c00 + 5*P[e.n00].c11*dG.c11 + 5*P[e.n01].c00*dG.c00 + 5*P[e.n01].c11*dG.c11 + 3*P[e.n10].c00*dG.c00 + 3*P[e.n10].c11*dG.c11 + 3*P[e.n11].c00*dG.c00 + 3*P[e.n11].c11*dG.c11);
19K(e.i02, e.i12) += (1.0/14700.0)*(e.X*e.X)*e.Y*(45*P[e.n00].c00*dG.c00 + 45*P[e.n00].c11*dG.c11 + 150*P[e.n01].c00*dG.c00 + 150*P[e.n01].c11*dG.c11 + 21*P[e.n10].c00*dG.c00 + 21*P[e.n10].c11*dG.c11 + 70*P[e.n11].c00*dG.c00 + 70*P[e.n11].c11*dG.c11);
20K(e.i03, e.i12) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/1680.0*P[e.n00].c00*dG.c00 - 1.0/1680.0*P[e.n00].c11*dG.c11 - 1.0/784.0*P[e.n01].c00*dG.c00 - 1.0/784.0*P[e.n01].c11*dG.c11 - 1.0/3600.0*P[e.n10].c00*dG.c00 - 1.0/3600.0*P[e.n10].c11*dG.c11 - 1.0/1680.0*P[e.n11].c00*dG.c00 - 1.0/1680.0*P[e.n11].c11*dG.c11);
21K(e.i12, e.i12) += (1.0/29400.0)*(e.X*e.X*e.X)*e.Y*(15*P[e.n00].c00*dG.c00 + 15*P[e.n00].c11*dG.c11 + 50*P[e.n01].c00*dG.c00 + 50*P[e.n01].c11*dG.c11 + 9*P[e.n10].c00*dG.c00 + 9*P[e.n10].c11*dG.c11 + 30*P[e.n11].c00*dG.c00 + 30*P[e.n11].c11*dG.c11);
22K(e.i00, e.i20) += (9.0/4900.0)*e.X*e.Y*(10*P[e.n00].c00*dG.c00 + 10*P[e.n00].c11*dG.c11 + 3*P[e.n01].c00*dG.c00 + 3*P[e.n01].c11*dG.c11 + 10*P[e.n10].c00*dG.c00 + 10*P[e.n10].c11*dG.c11 + 3*P[e.n11].c00*dG.c00 + 3*P[e.n11].c11*dG.c11);
23K(e.i01, e.i20) += (3.0/19600.0)*e.X*(e.Y*e.Y)*(15*P[e.n00].c00*dG.c00 + 15*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 15*P[e.n10].c00*dG.c00 + 15*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
24K(e.i10, e.i20) += (1.0/14700.0)*(e.X*e.X)*e.Y*(60*P[e.n00].c00*dG.c00 + 60*P[e.n00].c11*dG.c11 + 18*P[e.n01].c00*dG.c00 + 18*P[e.n01].c11*dG.c11 + 70*P[e.n10].c00*dG.c00 + 70*P[e.n10].c11*dG.c11 + 21*P[e.n11].c00*dG.c00 + 21*P[e.n11].c11*dG.c11);
25K(e.i02, e.i20) += (81.0/19600.0)*e.X*e.Y*(P[e.n00].c00*dG.c00 + P[e.n00].c11*dG.c11 + P[e.n01].c00*dG.c00 + P[e.n01].c11*dG.c11 + P[e.n10].c00*dG.c00 + P[e.n10].c11*dG.c11 + P[e.n11].c00*dG.c00 + P[e.n11].c11*dG.c11);
26K(e.i03, e.i20) += -3.0/19600.0*e.X*(e.Y*e.Y)*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 6*P[e.n01].c00*dG.c00 + 6*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 6*P[e.n11].c00*dG.c00 + 6*P[e.n11].c11*dG.c11);
27K(e.i12, e.i20) += (3.0/19600.0)*(e.X*e.X)*e.Y*(6*P[e.n00].c00*dG.c00 + 6*P[e.n00].c11*dG.c11 + 6*P[e.n01].c00*dG.c00 + 6*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
28K(e.i20, e.i20) += (1.0/1225.0)*e.X*e.Y*(30*P[e.n00].c00*dG.c00 + 30*P[e.n00].c11*dG.c11 + 9*P[e.n01].c00*dG.c00 + 9*P[e.n01].c11*dG.c11 + 100*P[e.n10].c00*dG.c00 + 100*P[e.n10].c11*dG.c11 + 30*P[e.n11].c00*dG.c00 + 30*P[e.n11].c11*dG.c11);
29K(e.i00, e.i21) += (3.0/19600.0)*e.X*(e.Y*e.Y)*(15*P[e.n00].c00*dG.c00 + 15*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 15*P[e.n10].c00*dG.c00 + 15*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
30K(e.i01, e.i21) += (3.0/39200.0)*e.X*(e.Y*e.Y*e.Y)*(5*P[e.n00].c00*dG.c00 + 5*P[e.n00].c11*dG.c11 + 3*P[e.n01].c00*dG.c00 + 3*P[e.n01].c11*dG.c11 + 5*P[e.n10].c00*dG.c00 + 5*P[e.n10].c11*dG.c11 + 3*P[e.n11].c00*dG.c00 + 3*P[e.n11].c11*dG.c11);
31K(e.i10, e.i21) += (e.X*e.X)*(e.Y*e.Y)*((1.0/1960.0)*P[e.n00].c00*dG.c00 + (1.0/1960.0)*P[e.n00].c11*dG.c11 + (1.0/4200.0)*P[e.n01].c00*dG.c00 + (1.0/4200.0)*P[e.n01].c11*dG.c11 + (1.0/1680.0)*P[e.n10].c00*dG.c00 + (1.0/1680.0)*P[e.n10].c11*dG.c11 + (1.0/3600.0)*P[e.n11].c00*dG.c00 + (1.0/3600.0)*P[e.n11].c11*dG.c11);
32K(e.i02, e.i21) += (3.0/19600.0)*e.X*(e.Y*e.Y)*(6*P[e.n00].c00*dG.c00 + 6*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 6*P[e.n10].c00*dG.c00 + 6*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
33K(e.i03, e.i21) += -9.0/39200.0*e.X*(e.Y*e.Y*e.Y)*(P[e.n00].c00*dG.c00 + P[e.n00].c11*dG.c11 + P[e.n01].c00*dG.c00 + P[e.n01].c11*dG.c11 + P[e.n10].c00*dG.c00 + P[e.n10].c11*dG.c11 + P[e.n11].c00*dG.c00 + P[e.n11].c11*dG.c11);
34K(e.i12, e.i21) += (e.X*e.X)*(e.Y*e.Y)*((1.0/4900.0)*P[e.n00].c00*dG.c00 + (1.0/4900.0)*P[e.n00].c11*dG.c11 + (1.0/4200.0)*P[e.n01].c00*dG.c00 + (1.0/4200.0)*P[e.n01].c11*dG.c11 + (1.0/4200.0)*P[e.n10].c00*dG.c00 + (1.0/4200.0)*P[e.n10].c11*dG.c11 + (1.0/3600.0)*P[e.n11].c00*dG.c00 + (1.0/3600.0)*P[e.n11].c11*dG.c11);
35K(e.i20, e.i21) += (1.0/14700.0)*e.X*(e.Y*e.Y)*(45*P[e.n00].c00*dG.c00 + 45*P[e.n00].c11*dG.c11 + 21*P[e.n01].c00*dG.c00 + 21*P[e.n01].c11*dG.c11 + 150*P[e.n10].c00*dG.c00 + 150*P[e.n10].c11*dG.c11 + 70*P[e.n11].c00*dG.c00 + 70*P[e.n11].c11*dG.c11);
36K(e.i21, e.i21) += (1.0/29400.0)*e.X*(e.Y*e.Y*e.Y)*(15*P[e.n00].c00*dG.c00 + 15*P[e.n00].c11*dG.c11 + 9*P[e.n01].c00*dG.c00 + 9*P[e.n01].c11*dG.c11 + 50*P[e.n10].c00*dG.c00 + 50*P[e.n10].c11*dG.c11 + 30*P[e.n11].c00*dG.c00 + 30*P[e.n11].c11*dG.c11);
37K(e.i00, e.i30) += -1.0/14700.0*(e.X*e.X)*e.Y*(70*P[e.n00].c00*dG.c00 + 70*P[e.n00].c11*dG.c11 + 21*P[e.n01].c00*dG.c00 + 21*P[e.n01].c11*dG.c11 + 60*P[e.n10].c00*dG.c00 + 60*P[e.n10].c11*dG.c11 + 18*P[e.n11].c00*dG.c00 + 18*P[e.n11].c11*dG.c11);
38K(e.i01, e.i30) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/1680.0*P[e.n00].c00*dG.c00 - 1.0/1680.0*P[e.n00].c11*dG.c11 - 1.0/3600.0*P[e.n01].c00*dG.c00 - 1.0/3600.0*P[e.n01].c11*dG.c11 - 1.0/1960.0*P[e.n10].c00*dG.c00 - 1.0/1960.0*P[e.n10].c11*dG.c11 - 1.0/4200.0*P[e.n11].c00*dG.c00 - 1.0/4200.0*P[e.n11].c11*dG.c11);
39K(e.i10, e.i30) += -1.0/9800.0*(e.X*e.X*e.X)*e.Y*(10*P[e.n00].c00*dG.c00 + 10*P[e.n00].c11*dG.c11 + 3*P[e.n01].c00*dG.c00 + 3*P[e.n01].c11*dG.c11 + 10*P[e.n10].c00*dG.c00 + 10*P[e.n10].c11*dG.c11 + 3*P[e.n11].c00*dG.c00 + 3*P[e.n11].c11*dG.c11);
40K(e.i02, e.i30) += -3.0/19600.0*(e.X*e.X)*e.Y*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 6*P[e.n10].c00*dG.c00 + 6*P[e.n10].c11*dG.c11 + 6*P[e.n11].c00*dG.c00 + 6*P[e.n11].c11*dG.c11);
41K(e.i03, e.i30) += (e.X*e.X)*(e.Y*e.Y)*((1.0/3600.0)*P[e.n00].c00*dG.c00 + (1.0/3600.0)*P[e.n00].c11*dG.c11 + (1.0/4200.0)*P[e.n01].c00*dG.c00 + (1.0/4200.0)*P[e.n01].c11*dG.c11 + (1.0/4200.0)*P[e.n10].c00*dG.c00 + (1.0/4200.0)*P[e.n10].c11*dG.c11 + (1.0/4900.0)*P[e.n11].c00*dG.c00 + (1.0/4900.0)*P[e.n11].c11*dG.c11);
42K(e.i12, e.i30) += -9.0/39200.0*(e.X*e.X*e.X)*e.Y*(P[e.n00].c00*dG.c00 + P[e.n00].c11*dG.c11 + P[e.n01].c00*dG.c00 + P[e.n01].c11*dG.c11 + P[e.n10].c00*dG.c00 + P[e.n10].c11*dG.c11 + P[e.n11].c00*dG.c00 + P[e.n11].c11*dG.c11);
43K(e.i20, e.i30) += -1.0/14700.0*(e.X*e.X)*e.Y*(70*P[e.n00].c00*dG.c00 + 70*P[e.n00].c11*dG.c11 + 21*P[e.n01].c00*dG.c00 + 21*P[e.n01].c11*dG.c11 + 150*P[e.n10].c00*dG.c00 + 150*P[e.n10].c11*dG.c11 + 45*P[e.n11].c00*dG.c00 + 45*P[e.n11].c11*dG.c11);
44K(e.i21, e.i30) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/1680.0*P[e.n00].c00*dG.c00 - 1.0/1680.0*P[e.n00].c11*dG.c11 - 1.0/3600.0*P[e.n01].c00*dG.c00 - 1.0/3600.0*P[e.n01].c11*dG.c11 - 1.0/784.0*P[e.n10].c00*dG.c00 - 1.0/784.0*P[e.n10].c11*dG.c11 - 1.0/1680.0*P[e.n11].c00*dG.c00 - 1.0/1680.0*P[e.n11].c11*dG.c11);
45K(e.i30, e.i30) += (1.0/29400.0)*(e.X*e.X*e.X)*e.Y*(30*P[e.n00].c00*dG.c00 + 30*P[e.n00].c11*dG.c11 + 9*P[e.n01].c00*dG.c00 + 9*P[e.n01].c11*dG.c11 + 50*P[e.n10].c00*dG.c00 + 50*P[e.n10].c11*dG.c11 + 15*P[e.n11].c00*dG.c00 + 15*P[e.n11].c11*dG.c11);
46K(e.i00, e.i22) += (81.0/19600.0)*e.X*e.Y*(P[e.n00].c00*dG.c00 + P[e.n00].c11*dG.c11 + P[e.n01].c00*dG.c00 + P[e.n01].c11*dG.c11 + P[e.n10].c00*dG.c00 + P[e.n10].c11*dG.c11 + P[e.n11].c00*dG.c00 + P[e.n11].c11*dG.c11);
47K(e.i01, e.i22) += (3.0/19600.0)*e.X*(e.Y*e.Y)*(6*P[e.n00].c00*dG.c00 + 6*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 6*P[e.n10].c00*dG.c00 + 6*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
48K(e.i10, e.i22) += (3.0/19600.0)*(e.X*e.X)*e.Y*(6*P[e.n00].c00*dG.c00 + 6*P[e.n00].c11*dG.c11 + 6*P[e.n01].c00*dG.c00 + 6*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 7*P[e.n11].c00*dG.c00 + 7*P[e.n11].c11*dG.c11);
49K(e.i02, e.i22) += (9.0/4900.0)*e.X*e.Y*(3*P[e.n00].c00*dG.c00 + 3*P[e.n00].c11*dG.c11 + 10*P[e.n01].c00*dG.c00 + 10*P[e.n01].c11*dG.c11 + 3*P[e.n10].c00*dG.c00 + 3*P[e.n10].c11*dG.c11 + 10*P[e.n11].c00*dG.c00 + 10*P[e.n11].c11*dG.c11);
50K(e.i03, e.i22) += -3.0/19600.0*e.X*(e.Y*e.Y)*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 15*P[e.n01].c00*dG.c00 + 15*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 15*P[e.n11].c00*dG.c00 + 15*P[e.n11].c11*dG.c11);
51K(e.i12, e.i22) += (1.0/14700.0)*(e.X*e.X)*e.Y*(18*P[e.n00].c00*dG.c00 + 18*P[e.n00].c11*dG.c11 + 60*P[e.n01].c00*dG.c00 + 60*P[e.n01].c11*dG.c11 + 21*P[e.n10].c00*dG.c00 + 21*P[e.n10].c11*dG.c11 + 70*P[e.n11].c00*dG.c00 + 70*P[e.n11].c11*dG.c11);
52K(e.i20, e.i22) += (9.0/4900.0)*e.X*e.Y*(3*P[e.n00].c00*dG.c00 + 3*P[e.n00].c11*dG.c11 + 3*P[e.n01].c00*dG.c00 + 3*P[e.n01].c11*dG.c11 + 10*P[e.n10].c00*dG.c00 + 10*P[e.n10].c11*dG.c11 + 10*P[e.n11].c00*dG.c00 + 10*P[e.n11].c11*dG.c11);
53K(e.i21, e.i22) += (1.0/14700.0)*e.X*(e.Y*e.Y)*(18*P[e.n00].c00*dG.c00 + 18*P[e.n00].c11*dG.c11 + 21*P[e.n01].c00*dG.c00 + 21*P[e.n01].c11*dG.c11 + 60*P[e.n10].c00*dG.c00 + 60*P[e.n10].c11*dG.c11 + 70*P[e.n11].c00*dG.c00 + 70*P[e.n11].c11*dG.c11);
54K(e.i30, e.i22) += -3.0/19600.0*(e.X*e.X)*e.Y*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 15*P[e.n10].c00*dG.c00 + 15*P[e.n10].c11*dG.c11 + 15*P[e.n11].c00*dG.c00 + 15*P[e.n11].c11*dG.c11);
55K(e.i22, e.i22) += (1.0/1225.0)*e.X*e.Y*(9*P[e.n00].c00*dG.c00 + 9*P[e.n00].c11*dG.c11 + 30*P[e.n01].c00*dG.c00 + 30*P[e.n01].c11*dG.c11 + 30*P[e.n10].c00*dG.c00 + 30*P[e.n10].c11*dG.c11 + 100*P[e.n11].c00*dG.c00 + 100*P[e.n11].c11*dG.c11);
56K(e.i00, e.i23) += -3.0/19600.0*e.X*(e.Y*e.Y)*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 6*P[e.n01].c00*dG.c00 + 6*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 6*P[e.n11].c00*dG.c00 + 6*P[e.n11].c11*dG.c11);
57K(e.i01, e.i23) += -9.0/39200.0*e.X*(e.Y*e.Y*e.Y)*(P[e.n00].c00*dG.c00 + P[e.n00].c11*dG.c11 + P[e.n01].c00*dG.c00 + P[e.n01].c11*dG.c11 + P[e.n10].c00*dG.c00 + P[e.n10].c11*dG.c11 + P[e.n11].c00*dG.c00 + P[e.n11].c11*dG.c11);
58K(e.i10, e.i23) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/4200.0*P[e.n00].c00*dG.c00 - 1.0/4200.0*P[e.n00].c11*dG.c11 - 1.0/4900.0*P[e.n01].c00*dG.c00 - 1.0/4900.0*P[e.n01].c11*dG.c11 - 1.0/3600.0*P[e.n10].c00*dG.c00 - 1.0/3600.0*P[e.n10].c11*dG.c11 - 1.0/4200.0*P[e.n11].c00*dG.c00 - 1.0/4200.0*P[e.n11].c11*dG.c11);
59K(e.i02, e.i23) += -3.0/19600.0*e.X*(e.Y*e.Y)*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 15*P[e.n01].c00*dG.c00 + 15*P[e.n01].c11*dG.c11 + 7*P[e.n10].c00*dG.c00 + 7*P[e.n10].c11*dG.c11 + 15*P[e.n11].c00*dG.c00 + 15*P[e.n11].c11*dG.c11);
60K(e.i03, e.i23) += (3.0/39200.0)*e.X*(e.Y*e.Y*e.Y)*(3*P[e.n00].c00*dG.c00 + 3*P[e.n00].c11*dG.c11 + 5*P[e.n01].c00*dG.c00 + 5*P[e.n01].c11*dG.c11 + 3*P[e.n10].c00*dG.c00 + 3*P[e.n10].c11*dG.c11 + 5*P[e.n11].c00*dG.c00 + 5*P[e.n11].c11*dG.c11);
61K(e.i12, e.i23) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/4200.0*P[e.n00].c00*dG.c00 - 1.0/4200.0*P[e.n00].c11*dG.c11 - 1.0/1960.0*P[e.n01].c00*dG.c00 - 1.0/1960.0*P[e.n01].c11*dG.c11 - 1.0/3600.0*P[e.n10].c00*dG.c00 - 1.0/3600.0*P[e.n10].c11*dG.c11 - 1.0/1680.0*P[e.n11].c00*dG.c00 - 1.0/1680.0*P[e.n11].c11*dG.c11);
62K(e.i20, e.i23) += -1.0/14700.0*e.X*(e.Y*e.Y)*(21*P[e.n00].c00*dG.c00 + 21*P[e.n00].c11*dG.c11 + 18*P[e.n01].c00*dG.c00 + 18*P[e.n01].c11*dG.c11 + 70*P[e.n10].c00*dG.c00 + 70*P[e.n10].c11*dG.c11 + 60*P[e.n11].c00*dG.c00 + 60*P[e.n11].c11*dG.c11);
63K(e.i21, e.i23) += -1.0/9800.0*e.X*(e.Y*e.Y*e.Y)*(3*P[e.n00].c00*dG.c00 + 3*P[e.n00].c11*dG.c11 + 3*P[e.n01].c00*dG.c00 + 3*P[e.n01].c11*dG.c11 + 10*P[e.n10].c00*dG.c00 + 10*P[e.n10].c11*dG.c11 + 10*P[e.n11].c00*dG.c00 + 10*P[e.n11].c11*dG.c11);
64K(e.i30, e.i23) += (e.X*e.X)*(e.Y*e.Y)*((1.0/3600.0)*P[e.n00].c00*dG.c00 + (1.0/3600.0)*P[e.n00].c11*dG.c11 + (1.0/4200.0)*P[e.n01].c00*dG.c00 + (1.0/4200.0)*P[e.n01].c11*dG.c11 + (1.0/1680.0)*P[e.n10].c00*dG.c00 + (1.0/1680.0)*P[e.n10].c11*dG.c11 + (1.0/1960.0)*P[e.n11].c00*dG.c00 + (1.0/1960.0)*P[e.n11].c11*dG.c11);
65K(e.i22, e.i23) += -1.0/14700.0*e.X*(e.Y*e.Y)*(21*P[e.n00].c00*dG.c00 + 21*P[e.n00].c11*dG.c11 + 45*P[e.n01].c00*dG.c00 + 45*P[e.n01].c11*dG.c11 + 70*P[e.n10].c00*dG.c00 + 70*P[e.n10].c11*dG.c11 + 150*P[e.n11].c00*dG.c00 + 150*P[e.n11].c11*dG.c11);
66K(e.i23, e.i23) += (1.0/29400.0)*e.X*(e.Y*e.Y*e.Y)*(9*P[e.n00].c00*dG.c00 + 9*P[e.n00].c11*dG.c11 + 15*P[e.n01].c00*dG.c00 + 15*P[e.n01].c11*dG.c11 + 30*P[e.n10].c00*dG.c00 + 30*P[e.n10].c11*dG.c11 + 50*P[e.n11].c00*dG.c00 + 50*P[e.n11].c11*dG.c11);
67K(e.i00, e.i32) += -3.0/19600.0*(e.X*e.X)*e.Y*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 6*P[e.n10].c00*dG.c00 + 6*P[e.n10].c11*dG.c11 + 6*P[e.n11].c00*dG.c00 + 6*P[e.n11].c11*dG.c11);
68K(e.i01, e.i32) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/4200.0*P[e.n00].c00*dG.c00 - 1.0/4200.0*P[e.n00].c11*dG.c11 - 1.0/3600.0*P[e.n01].c00*dG.c00 - 1.0/3600.0*P[e.n01].c11*dG.c11 - 1.0/4900.0*P[e.n10].c00*dG.c00 - 1.0/4900.0*P[e.n10].c11*dG.c11 - 1.0/4200.0*P[e.n11].c00*dG.c00 - 1.0/4200.0*P[e.n11].c11*dG.c11);
69K(e.i10, e.i32) += -9.0/39200.0*(e.X*e.X*e.X)*e.Y*(P[e.n00].c00*dG.c00 + P[e.n00].c11*dG.c11 + P[e.n01].c00*dG.c00 + P[e.n01].c11*dG.c11 + P[e.n10].c00*dG.c00 + P[e.n10].c11*dG.c11 + P[e.n11].c00*dG.c00 + P[e.n11].c11*dG.c11);
70K(e.i02, e.i32) += -1.0/14700.0*(e.X*e.X)*e.Y*(21*P[e.n00].c00*dG.c00 + 21*P[e.n00].c11*dG.c11 + 70*P[e.n01].c00*dG.c00 + 70*P[e.n01].c11*dG.c11 + 18*P[e.n10].c00*dG.c00 + 18*P[e.n10].c11*dG.c11 + 60*P[e.n11].c00*dG.c00 + 60*P[e.n11].c11*dG.c11);
71K(e.i03, e.i32) += (e.X*e.X)*(e.Y*e.Y)*((1.0/3600.0)*P[e.n00].c00*dG.c00 + (1.0/3600.0)*P[e.n00].c11*dG.c11 + (1.0/1680.0)*P[e.n01].c00*dG.c00 + (1.0/1680.0)*P[e.n01].c11*dG.c11 + (1.0/4200.0)*P[e.n10].c00*dG.c00 + (1.0/4200.0)*P[e.n10].c11*dG.c11 + (1.0/1960.0)*P[e.n11].c00*dG.c00 + (1.0/1960.0)*P[e.n11].c11*dG.c11);
72K(e.i12, e.i32) += -1.0/9800.0*(e.X*e.X*e.X)*e.Y*(3*P[e.n00].c00*dG.c00 + 3*P[e.n00].c11*dG.c11 + 10*P[e.n01].c00*dG.c00 + 10*P[e.n01].c11*dG.c11 + 3*P[e.n10].c00*dG.c00 + 3*P[e.n10].c11*dG.c11 + 10*P[e.n11].c00*dG.c00 + 10*P[e.n11].c11*dG.c11);
73K(e.i20, e.i32) += -3.0/19600.0*(e.X*e.X)*e.Y*(7*P[e.n00].c00*dG.c00 + 7*P[e.n00].c11*dG.c11 + 7*P[e.n01].c00*dG.c00 + 7*P[e.n01].c11*dG.c11 + 15*P[e.n10].c00*dG.c00 + 15*P[e.n10].c11*dG.c11 + 15*P[e.n11].c00*dG.c00 + 15*P[e.n11].c11*dG.c11);
74K(e.i21, e.i32) += (e.X*e.X)*(e.Y*e.Y)*(-1.0/4200.0*P[e.n00].c00*dG.c00 - 1.0/4200.0*P[e.n00].c11*dG.c11 - 1.0/3600.0*P[e.n01].c00*dG.c00 - 1.0/3600.0*P[e.n01].c11*dG.c11 - 1.0/1960.0*P[e.n10].c00*dG.c00 - 1.0/1960.0*P[e.n10].c11*dG.c11 - 1.0/1680.0*P[e.n11].c00*dG.c00 - 1.0/1680.0*P[e.n11].c11*dG.c11);
75K(e.i30, e.i32) += (3.0/39200.0)*(e.X*e.X*e.X)*e.Y*(3*P[e.n00].c00*dG.c00 + 3*P[e.n00].c11*dG.c11 + 3*P[e.n01].c00*dG.c00 + 3*P[e.n01].c11*dG.c11 + 5*P[e.n10].c00*dG.c00 + 5*P[e.n10].c11*dG.c11 + 5*P[e.n11].c00*dG.c00 + 5*P[e.n11].c11*dG.c11);
76K(e.i22, e.i32) += -1.0/14700.0*(e.X*e.X)*e.Y*(21*P[e.n00].c00*dG.c00 + 21*P[e.n00].c11*dG.c11 + 70*P[e.n01].c00*dG.c00 + 70*P[e.n01].c11*dG.c11 + 45*P[e.n10].c00*dG.c00 + 45*P[e.n10].c11*dG.c11 + 150*P[e.n11].c00*dG.c00 + 150*P[e.n11].c11*dG.c11);
77K(e.i23, e.i32) += (e.X*e.X)*(e.Y*e.Y)*((1.0/3600.0)*P[e.n00].c00*dG.c00 + (1.0/3600.0)*P[e.n00].c11*dG.c11 + (1.0/1680.0)*P[e.n01].c00*dG.c00 + (1.0/1680.0)*P[e.n01].c11*dG.c11 + (1.0/1680.0)*P[e.n10].c00*dG.c00 + (1.0/1680.0)*P[e.n10].c11*dG.c11 + (1.0/784.0)*P[e.n11].c00*dG.c00 + (1.0/784.0)*P[e.n11].c11*dG.c11);
78K(e.i32, e.i32) += (1.0/29400.0)*(e.X*e.X*e.X)*e.Y*(9*P[e.n00].c00*dG.c00 + 9*P[e.n00].c11*dG.c11 + 30*P[e.n01].c00*dG.c00 + 30*P[e.n01].c11*dG.c11 + 15*P[e.n10].c00*dG.c00 + 15*P[e.n10].c11*dG.c11 + 50*P[e.n11].c00*dG.c00 + 50*P[e.n11].c11*dG.c11);
79F[e.i00] += (1.0/400.0)*e.X*e.Y*(49*Ug*P[e.n00].c00*dG.c00 + 49*Ug*P[e.n00].c11*dG.c11 + 21*Ug*P[e.n01].c00*dG.c00 + 21*Ug*P[e.n01].c11*dG.c11 + 21*Ug*P[e.n10].c00*dG.c00 + 21*Ug*P[e.n10].c11*dG.c11 + 9*Ug*P[e.n11].c00*dG.c00 + 9*Ug*P[e.n11].c11*dG.c11 - 49*G.c00*P[e.n00].c00 - 21*G.c00*P[e.n01].c00 - 21*G.c00*P[e.n10].c00 - 9*G.c00*P[e.n11].c00 - 49*G.c11*P[e.n00].c11 - 21*G.c11*P[e.n01].c11 - 21*G.c11*P[e.n10].c11 - 9*G.c11*P[e.n11].c11);
80F[e.i01] += (1.0/1200.0)*e.X*(e.Y*e.Y)*(21*Ug*P[e.n00].c00*dG.c00 + 21*Ug*P[e.n00].c11*dG.c11 + 14*Ug*P[e.n01].c00*dG.c00 + 14*Ug*P[e.n01].c11*dG.c11 + 9*Ug*P[e.n10].c00*dG.c00 + 9*Ug*P[e.n10].c11*dG.c11 + 6*Ug*P[e.n11].c00*dG.c00 + 6*Ug*P[e.n11].c11*dG.c11 - 21*G.c00*P[e.n00].c00 - 14*G.c00*P[e.n01].c00 - 9*G.c00*P[e.n10].c00 - 6*G.c00*P[e.n11].c00 - 21*G.c11*P[e.n00].c11 - 14*G.c11*P[e.n01].c11 - 9*G.c11*P[e.n10].c11 - 6*G.c11*P[e.n11].c11);
81F[e.i10] += (1.0/1200.0)*(e.X*e.X)*e.Y*(21*Ug*P[e.n00].c00*dG.c00 + 21*Ug*P[e.n00].c11*dG.c11 + 9*Ug*P[e.n01].c00*dG.c00 + 9*Ug*P[e.n01].c11*dG.c11 + 14*Ug*P[e.n10].c00*dG.c00 + 14*Ug*P[e.n10].c11*dG.c11 + 6*Ug*P[e.n11].c00*dG.c00 + 6*Ug*P[e.n11].c11*dG.c11 - 21*G.c00*P[e.n00].c00 - 9*G.c00*P[e.n01].c00 - 14*G.c00*P[e.n10].c00 - 6*G.c00*P[e.n11].c00 - 21*G.c11*P[e.n00].c11 - 9*G.c11*P[e.n01].c11 - 14*G.c11*P[e.n10].c11 - 6*G.c11*P[e.n11].c11);
82F[e.i02] += (1.0/400.0)*e.X*e.Y*(21*Ug*P[e.n00].c00*dG.c00 + 21*Ug*P[e.n00].c11*dG.c11 + 49*Ug*P[e.n01].c00*dG.c00 + 49*Ug*P[e.n01].c11*dG.c11 + 9*Ug*P[e.n10].c00*dG.c00 + 9*Ug*P[e.n10].c11*dG.c11 + 21*Ug*P[e.n11].c00*dG.c00 + 21*Ug*P[e.n11].c11*dG.c11 - 21*G.c00*P[e.n00].c00 - 49*G.c00*P[e.n01].c00 - 9*G.c00*P[e.n10].c00 - 21*G.c00*P[e.n11].c00 - 21*G.c11*P[e.n00].c11 - 49*G.c11*P[e.n01].c11 - 9*G.c11*P[e.n10].c11 - 21*G.c11*P[e.n11].c11);
83F[e.i03] += (1.0/1200.0)*e.X*(e.Y*e.Y)*(-14*Ug*P[e.n00].c00*dG.c00 - 14*Ug*P[e.n00].c11*dG.c11 - 21*Ug*P[e.n01].c00*dG.c00 - 21*Ug*P[e.n01].c11*dG.c11 - 6*Ug*P[e.n10].c00*dG.c00 - 6*Ug*P[e.n10].c11*dG.c11 - 9*Ug*P[e.n11].c00*dG.c00 - 9*Ug*P[e.n11].c11*dG.c11 + 14*G.c00*P[e.n00].c00 + 21*G.c00*P[e.n01].c00 + 6*G.c00*P[e.n10].c00 + 9*G.c00*P[e.n11].c00 + 14*G.c11*P[e.n00].c11 + 21*G.c11*P[e.n01].c11 + 6*G.c11*P[e.n10].c11 + 9*G.c11*P[e.n11].c11);
84F[e.i12] += (1.0/1200.0)*(e.X*e.X)*e.Y*(9*Ug*P[e.n00].c00*dG.c00 + 9*Ug*P[e.n00].c11*dG.c11 + 21*Ug*P[e.n01].c00*dG.c00 + 21*Ug*P[e.n01].c11*dG.c11 + 6*Ug*P[e.n10].c00*dG.c00 + 6*Ug*P[e.n10].c11*dG.c11 + 14*Ug*P[e.n11].c00*dG.c00 + 14*Ug*P[e.n11].c11*dG.c11 - 9*G.c00*P[e.n00].c00 - 21*G.c00*P[e.n01].c00 - 6*G.c00*P[e.n10].c00 - 14*G.c00*P[e.n11].c00 - 9*G.c11*P[e.n00].c11 - 21*G.c11*P[e.n01].c11 - 6*G.c11*P[e.n10].c11 - 14*G.c11*P[e.n11].c11);
85F[e.i20] += (1.0/400.0)*e.X*e.Y*(21*Ug*P[e.n00].c00*dG.c00 + 21*Ug*P[e.n00].c11*dG.c11 + 9*Ug*P[e.n01].c00*dG.c00 + 9*Ug*P[e.n01].c11*dG.c11 + 49*Ug*P[e.n10].c00*dG.c00 + 49*Ug*P[e.n10].c11*dG.c11 + 21*Ug*P[e.n11].c00*dG.c00 + 21*Ug*P[e.n11].c11*dG.c11 - 21*G.c00*P[e.n00].c00 - 9*G.c00*P[e.n01].c00 - 49*G.c00*P[e.n10].c00 - 21*G.c00*P[e.n11].c00 - 21*G.c11*P[e.n00].c11 - 9*G.c11*P[e.n01].c11 - 49*G.c11*P[e.n10].c11 - 21*G.c11*P[e.n11].c11);
86F[e.i21] += (1.0/1200.0)*e.X*(e.Y*e.Y)*(9*Ug*P[e.n00].c00*dG.c00 + 9*Ug*P[e.n00].c11*dG.c11 + 6*Ug*P[e.n01].c00*dG.c00 + 6*Ug*P[e.n01].c11*dG.c11 + 21*Ug*P[e.n10].c00*dG.c00 + 21*Ug*P[e.n10].c11*dG.c11 + 14*Ug*P[e.n11].c00*dG.c00 + 14*Ug*P[e.n11].c11*dG.c11 - 9*G.c00*P[e.n00].c00 - 6*G.c00*P[e.n01].c00 - 21*G.c00*P[e.n10].c00 - 14*G.c00*P[e.n11].c00 - 9*G.c11*P[e.n00].c11 - 6*G.c11*P[e.n01].c11 - 21*G.c11*P[e.n10].c11 - 14*G.c11*P[e.n11].c11);
87F[e.i30] += (1.0/1200.0)*(e.X*e.X)*e.Y*(-14*Ug*P[e.n00].c00*dG.c00 - 14*Ug*P[e.n00].c11*dG.c11 - 6*Ug*P[e.n01].c00*dG.c00 - 6*Ug*P[e.n01].c11*dG.c11 - 21*Ug*P[e.n10].c00*dG.c00 - 21*Ug*P[e.n10].c11*dG.c11 - 9*Ug*P[e.n11].c00*dG.c00 - 9*Ug*P[e.n11].c11*dG.c11 + 14*G.c00*P[e.n00].c00 + 6*G.c00*P[e.n01].c00 + 21*G.c00*P[e.n10].c00 + 9*G.c00*P[e.n11].c00 + 14*G.c11*P[e.n00].c11 + 6*G.c11*P[e.n01].c11 + 21*G.c11*P[e.n10].c11 + 9*G.c11*P[e.n11].c11);
88F[e.i22] += (1.0/400.0)*e.X*e.Y*(9*Ug*P[e.n00].c00*dG.c00 + 9*Ug*P[e.n00].c11*dG.c11 + 21*Ug*P[e.n01].c00*dG.c00 + 21*Ug*P[e.n01].c11*dG.c11 + 21*Ug*P[e.n10].c00*dG.c00 + 21*Ug*P[e.n10].c11*dG.c11 + 49*Ug*P[e.n11].c00*dG.c00 + 49*Ug*P[e.n11].c11*dG.c11 - 9*G.c00*P[e.n00].c00 - 21*G.c00*P[e.n01].c00 - 21*G.c00*P[e.n10].c00 - 49*G.c00*P[e.n11].c00 - 9*G.c11*P[e.n00].c11 - 21*G.c11*P[e.n01].c11 - 21*G.c11*P[e.n10].c11 - 49*G.c11*P[e.n11].c11);
89F[e.i23] += (1.0/1200.0)*e.X*(e.Y*e.Y)*(-6*Ug*P[e.n00].c00*dG.c00 - 6*Ug*P[e.n00].c11*dG.c11 - 9*Ug*P[e.n01].c00*dG.c00 - 9*Ug*P[e.n01].c11*dG.c11 - 14*Ug*P[e.n10].c00*dG.c00 - 14*Ug*P[e.n10].c11*dG.c11 - 21*Ug*P[e.n11].c00*dG.c00 - 21*Ug*P[e.n11].c11*dG.c11 + 6*G.c00*P[e.n00].c00 + 9*G.c00*P[e.n01].c00 + 14*G.c00*P[e.n10].c00 + 21*G.c00*P[e.n11].c00 + 6*G.c11*P[e.n00].c11 + 9*G.c11*P[e.n01].c11 + 14*G.c11*P[e.n10].c11 + 21*G.c11*P[e.n11].c11);
90F[e.i32] += (1.0/1200.0)*(e.X*e.X)*e.Y*(-6*Ug*P[e.n00].c00*dG.c00 - 6*Ug*P[e.n00].c11*dG.c11 - 14*Ug*P[e.n01].c00*dG.c00 - 14*Ug*P[e.n01].c11*dG.c11 - 9*Ug*P[e.n10].c00*dG.c00 - 9*Ug*P[e.n10].c11*dG.c11 - 21*Ug*P[e.n11].c00*dG.c00 - 21*Ug*P[e.n11].c11*dG.c11 + 6*G.c00*P[e.n00].c00 + 14*G.c00*P[e.n01].c00 + 9*G.c00*P[e.n10].c00 + 21*G.c00*P[e.n11].c00 + 6*G.c11*P[e.n00].c11 + 14*G.c11*P[e.n01].c11 + 9*G.c11*P[e.n10].c11 + 21*G.c11*P[e.n11].c11);