You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
1680 lines
33 KiB
1680 lines
33 KiB
Shat7(1,z6(1,1,1,1,1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - 5.D0*Csum1(ep)
|
|
& - 5.D0*Csum2(ep)
|
|
& - 10.D0*Csum11(ep)
|
|
& - 20.D0*Csum12(ep)
|
|
& - 10.D0*Csum22(ep)
|
|
& - 10.D0*Csum111(ep)
|
|
& - 10.D0*Csum222(ep)
|
|
& - 30.D0*Csum112(ep)
|
|
& - 30.D0*Csum122(ep)
|
|
& - 5.D0*Csum1111(ep)
|
|
& - 5.D0*Csum2222(ep)
|
|
& - 20.D0*Csum1112(ep)
|
|
& - 30.D0*Csum1122(ep)
|
|
& - 20.D0*Csum1222(ep)
|
|
& - Csum11111(ep)
|
|
& - Csum22222(ep)
|
|
& - 5.D0*Csum11112(ep)
|
|
& - 5.D0*Csum12222(ep)
|
|
& - 10.D0*Csum11122(ep)
|
|
& - 10.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,1,1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - 5.D0*Csum1(ep)
|
|
& - 5.D0*Csum2(ep)
|
|
& - 10.D0*Csum11(ep)
|
|
& - 20.D0*Csum12(ep)
|
|
& - 10.D0*Csum22(ep)
|
|
& - 10.D0*Csum111(ep)
|
|
& - 10.D0*Csum222(ep)
|
|
& - 30.D0*Csum112(ep)
|
|
& - 30.D0*Csum122(ep)
|
|
& - 5.D0*Csum1111(ep)
|
|
& - 5.D0*Csum2222(ep)
|
|
& - 20.D0*Csum1112(ep)
|
|
& - 30.D0*Csum1122(ep)
|
|
& - 20.D0*Csum1222(ep)
|
|
& - Csum11111(ep)
|
|
& - Csum22222(ep)
|
|
& - 5.D0*Csum11112(ep)
|
|
& - 5.D0*Csum12222(ep)
|
|
& - 10.D0*Csum11122(ep)
|
|
& - 10.D0*Csum11222(ep)
|
|
& + Cv(cc111111 + C124,ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,1,1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - 5.D0*Csum1(ep)
|
|
& - 5.D0*Csum2(ep)
|
|
& - 10.D0*Csum11(ep)
|
|
& - 20.D0*Csum12(ep)
|
|
& - 10.D0*Csum22(ep)
|
|
& - 10.D0*Csum111(ep)
|
|
& - 10.D0*Csum222(ep)
|
|
& - 30.D0*Csum112(ep)
|
|
& - 30.D0*Csum122(ep)
|
|
& - 5.D0*Csum1111(ep)
|
|
& - 5.D0*Csum2222(ep)
|
|
& - 20.D0*Csum1112(ep)
|
|
& - 30.D0*Csum1122(ep)
|
|
& - 20.D0*Csum1222(ep)
|
|
& - Csum11111(ep)
|
|
& - Csum22222(ep)
|
|
& - 5.D0*Csum11112(ep)
|
|
& - 5.D0*Csum12222(ep)
|
|
& - 10.D0*Csum11122(ep)
|
|
& - 10.D0*Csum11222(ep)
|
|
& + Cv(cc111111 + C123,ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,1,1,2),ep)=
|
|
& + Csum1(ep)
|
|
& + 4.D0*Csum11(ep)
|
|
& + 4.D0*Csum12(ep)
|
|
& + 6.D0*Csum111(ep)
|
|
& + 12.D0*Csum112(ep)
|
|
& + 6.D0*Csum122(ep)
|
|
& + 4.D0*Csum1111(ep)
|
|
& + 12.D0*Csum1112(ep)
|
|
& + 12.D0*Csum1122(ep)
|
|
& + 4.D0*Csum1222(ep)
|
|
& + Csum11111(ep)
|
|
& + 4.D0*Csum11112(ep)
|
|
& + Csum12222(ep)
|
|
& + 6.D0*Csum11122(ep)
|
|
& + 4.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,1,1,2),ep)=
|
|
& + Csum1(ep)
|
|
& + 4.D0*Csum11(ep)
|
|
& + 4.D0*Csum12(ep)
|
|
& + 6.D0*Csum111(ep)
|
|
& + 12.D0*Csum112(ep)
|
|
& + 6.D0*Csum122(ep)
|
|
& + 4.D0*Csum1111(ep)
|
|
& + 12.D0*Csum1112(ep)
|
|
& + 12.D0*Csum1122(ep)
|
|
& + 4.D0*Csum1222(ep)
|
|
& + Csum11111(ep)
|
|
& + 4.D0*Csum11112(ep)
|
|
& + Csum12222(ep)
|
|
& + 6.D0*Csum11122(ep)
|
|
& + 4.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,1,1,2),ep)=
|
|
& + Csum1(ep)
|
|
& + 4.D0*Csum11(ep)
|
|
& + 4.D0*Csum12(ep)
|
|
& + 6.D0*Csum111(ep)
|
|
& + 12.D0*Csum112(ep)
|
|
& + 6.D0*Csum122(ep)
|
|
& + 4.D0*Csum1111(ep)
|
|
& + 12.D0*Csum1112(ep)
|
|
& + 12.D0*Csum1122(ep)
|
|
& + 4.D0*Csum1222(ep)
|
|
& + Csum11111(ep)
|
|
& + 4.D0*Csum11112(ep)
|
|
& + Csum12222(ep)
|
|
& + 6.D0*Csum11122(ep)
|
|
& + 4.D0*Csum11222(ep)
|
|
& + Cv(cc111112 + C123,ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,1,1,3),ep)=
|
|
& + Csum2(ep)
|
|
& + 4.D0*Csum12(ep)
|
|
& + 4.D0*Csum22(ep)
|
|
& + 6.D0*Csum222(ep)
|
|
& + 6.D0*Csum112(ep)
|
|
& + 12.D0*Csum122(ep)
|
|
& + 4.D0*Csum2222(ep)
|
|
& + 4.D0*Csum1112(ep)
|
|
& + 12.D0*Csum1122(ep)
|
|
& + 12.D0*Csum1222(ep)
|
|
& + Csum22222(ep)
|
|
& + Csum11112(ep)
|
|
& + 4.D0*Csum12222(ep)
|
|
& + 4.D0*Csum11122(ep)
|
|
& + 6.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,1,1,3),ep)=
|
|
& + Csum2(ep)
|
|
& + 4.D0*Csum12(ep)
|
|
& + 4.D0*Csum22(ep)
|
|
& + 6.D0*Csum222(ep)
|
|
& + 6.D0*Csum112(ep)
|
|
& + 12.D0*Csum122(ep)
|
|
& + 4.D0*Csum2222(ep)
|
|
& + 4.D0*Csum1112(ep)
|
|
& + 12.D0*Csum1122(ep)
|
|
& + 12.D0*Csum1222(ep)
|
|
& + Csum22222(ep)
|
|
& + Csum11112(ep)
|
|
& + 4.D0*Csum12222(ep)
|
|
& + 4.D0*Csum11122(ep)
|
|
& + 6.D0*Csum11222(ep)
|
|
& + Cv(cc111112 + C124,ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,1,1,3),ep)=
|
|
& + Csum2(ep)
|
|
& + 4.D0*Csum12(ep)
|
|
& + 4.D0*Csum22(ep)
|
|
& + 6.D0*Csum222(ep)
|
|
& + 6.D0*Csum112(ep)
|
|
& + 12.D0*Csum122(ep)
|
|
& + 4.D0*Csum2222(ep)
|
|
& + 4.D0*Csum1112(ep)
|
|
& + 12.D0*Csum1122(ep)
|
|
& + 12.D0*Csum1222(ep)
|
|
& + Csum22222(ep)
|
|
& + Csum11112(ep)
|
|
& + 4.D0*Csum12222(ep)
|
|
& + 4.D0*Csum11122(ep)
|
|
& + 6.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,1,2,2),ep)=
|
|
& - Csum11(ep)
|
|
& - 3.D0*Csum111(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum1111(ep)
|
|
& - 6.D0*Csum1112(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - Csum11111(ep)
|
|
& - 3.D0*Csum11112(ep)
|
|
& - 3.D0*Csum11122(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,1,2,2),ep)=
|
|
& - Csum11(ep)
|
|
& - 3.D0*Csum111(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum1111(ep)
|
|
& - 6.D0*Csum1112(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - Csum11111(ep)
|
|
& - 3.D0*Csum11112(ep)
|
|
& - 3.D0*Csum11122(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,1,2,2),ep)=
|
|
& - Csum11(ep)
|
|
& - 3.D0*Csum111(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum1111(ep)
|
|
& - 6.D0*Csum1112(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - Csum11111(ep)
|
|
& - 3.D0*Csum11112(ep)
|
|
& - 3.D0*Csum11122(ep)
|
|
& - Csum11222(ep)
|
|
& + Cv(cc111122 + C123,ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,1,2,3),ep)=
|
|
& - Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum1112(ep)
|
|
& - 6.D0*Csum1122(ep)
|
|
& - 3.D0*Csum1222(ep)
|
|
& - Csum11112(ep)
|
|
& - Csum12222(ep)
|
|
& - 3.D0*Csum11122(ep)
|
|
& - 3.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,1,2,3),ep)=
|
|
& - Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum1112(ep)
|
|
& - 6.D0*Csum1122(ep)
|
|
& - 3.D0*Csum1222(ep)
|
|
& - Csum11112(ep)
|
|
& - Csum12222(ep)
|
|
& - 3.D0*Csum11122(ep)
|
|
& - 3.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,1,2,3),ep)=
|
|
& - Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum1112(ep)
|
|
& - 6.D0*Csum1122(ep)
|
|
& - 3.D0*Csum1222(ep)
|
|
& - Csum11112(ep)
|
|
& - Csum12222(ep)
|
|
& - 3.D0*Csum11122(ep)
|
|
& - 3.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,1,3,3),ep)=
|
|
& - Csum22(ep)
|
|
& - 3.D0*Csum222(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum2222(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - 6.D0*Csum1222(ep)
|
|
& - Csum22222(ep)
|
|
& - 3.D0*Csum12222(ep)
|
|
& - Csum11122(ep)
|
|
& - 3.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,1,3,3),ep)=
|
|
& - Csum22(ep)
|
|
& - 3.D0*Csum222(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum2222(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - 6.D0*Csum1222(ep)
|
|
& - Csum22222(ep)
|
|
& - 3.D0*Csum12222(ep)
|
|
& - Csum11122(ep)
|
|
& - 3.D0*Csum11222(ep)
|
|
& + Cv(cc111122 + C124,ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,1,3,3),ep)=
|
|
& - Csum22(ep)
|
|
& - 3.D0*Csum222(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum2222(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - 6.D0*Csum1222(ep)
|
|
& - Csum22222(ep)
|
|
& - 3.D0*Csum12222(ep)
|
|
& - Csum11122(ep)
|
|
& - 3.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,2,2,2),ep)=
|
|
& + Csum111(ep)
|
|
& + 2.D0*Csum1111(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + Csum11111(ep)
|
|
& + 2.D0*Csum11112(ep)
|
|
& + Csum11122(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,2,2,2),ep)=
|
|
& + Csum111(ep)
|
|
& + 2.D0*Csum1111(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + Csum11111(ep)
|
|
& + 2.D0*Csum11112(ep)
|
|
& + Csum11122(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,2,2,2),ep)=
|
|
& + Csum111(ep)
|
|
& + 2.D0*Csum1111(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + Csum11111(ep)
|
|
& + 2.D0*Csum11112(ep)
|
|
& + Csum11122(ep)
|
|
& + Cv(cc111222 + C123,ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,2,2,3),ep)=
|
|
& + Csum112(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + Csum11112(ep)
|
|
& + 2.D0*Csum11122(ep)
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,2,2,3),ep)=
|
|
& + Csum112(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + Csum11112(ep)
|
|
& + 2.D0*Csum11122(ep)
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,2,2,3),ep)=
|
|
& + Csum112(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + Csum11112(ep)
|
|
& + 2.D0*Csum11122(ep)
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,2,3,3),ep)=
|
|
& + Csum122(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum12222(ep)
|
|
& + Csum11122(ep)
|
|
& + 2.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,2,3,3),ep)=
|
|
& + Csum122(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum12222(ep)
|
|
& + Csum11122(ep)
|
|
& + 2.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,2,3,3),ep)=
|
|
& + Csum122(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum12222(ep)
|
|
& + Csum11122(ep)
|
|
& + 2.D0*Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,1,3,3,3),ep)=
|
|
& + Csum222(ep)
|
|
& + 2.D0*Csum2222(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum22222(ep)
|
|
& + 2.D0*Csum12222(ep)
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,1,3,3,3),ep)=
|
|
& + Csum222(ep)
|
|
& + 2.D0*Csum2222(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum22222(ep)
|
|
& + 2.D0*Csum12222(ep)
|
|
& + Csum11222(ep)
|
|
& + Cv(cc111222 + C124,ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,1,3,3,3),ep)=
|
|
& + Csum222(ep)
|
|
& + 2.D0*Csum2222(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum22222(ep)
|
|
& + 2.D0*Csum12222(ep)
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,2,2,2,2),ep)=
|
|
& - Csum1111(ep)
|
|
& - Csum11111(ep)
|
|
& - Csum11112(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,2,2,2,2),ep)=
|
|
& - Csum1111(ep)
|
|
& - Csum11111(ep)
|
|
& - Csum11112(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,2,2,2,2),ep)=
|
|
& - Csum1111(ep)
|
|
& - Csum11111(ep)
|
|
& - Csum11112(ep)
|
|
& + Cv(cc112222 + C123,ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,2,2,2,3),ep)=
|
|
& - Csum1112(ep)
|
|
& - Csum11112(ep)
|
|
& - Csum11122(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,2,2,2,3),ep)=
|
|
& - Csum1112(ep)
|
|
& - Csum11112(ep)
|
|
& - Csum11122(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,2,2,2,3),ep)=
|
|
& - Csum1112(ep)
|
|
& - Csum11112(ep)
|
|
& - Csum11122(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,2,2,3,3),ep)=
|
|
& - Csum1122(ep)
|
|
& - Csum11122(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,2,2,3,3),ep)=
|
|
& - Csum1122(ep)
|
|
& - Csum11122(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,2,2,3,3),ep)=
|
|
& - Csum1122(ep)
|
|
& - Csum11122(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,2,3,3,3),ep)=
|
|
& - Csum1222(ep)
|
|
& - Csum12222(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,2,3,3,3),ep)=
|
|
& - Csum1222(ep)
|
|
& - Csum12222(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,2,3,3,3),ep)=
|
|
& - Csum1222(ep)
|
|
& - Csum12222(ep)
|
|
& - Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,1,3,3,3,3),ep)=
|
|
& - Csum2222(ep)
|
|
& - Csum22222(ep)
|
|
& - Csum12222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,1,3,3,3,3),ep)=
|
|
& - Csum2222(ep)
|
|
& - Csum22222(ep)
|
|
& - Csum12222(ep)
|
|
& + Cv(cc112222 + C124,ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,1,3,3,3,3),ep)=
|
|
& - Csum2222(ep)
|
|
& - Csum22222(ep)
|
|
& - Csum12222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,2,2,2,2,2),ep)=
|
|
& + Csum11111(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,2,2,2,2,2),ep)=
|
|
& + Csum11111(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,2,2,2,2,2),ep)=
|
|
& + Csum11111(ep)
|
|
& + Cv(cc122222 + C123,ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,2,2,2,2,3),ep)=
|
|
& + Csum11112(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,2,2,2,2,3),ep)=
|
|
& + Csum11112(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,2,2,2,2,3),ep)=
|
|
& + Csum11112(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,2,2,2,3,3),ep)=
|
|
& + Csum11122(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,2,2,2,3,3),ep)=
|
|
& + Csum11122(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,2,2,2,3,3),ep)=
|
|
& + Csum11122(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,2,2,3,3,3),ep)=
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,2,2,3,3,3),ep)=
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,2,2,3,3,3),ep)=
|
|
& + Csum11222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,2,3,3,3,3),ep)=
|
|
& + Csum12222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,2,3,3,3,3),ep)=
|
|
& + Csum12222(ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,2,3,3,3,3),ep)=
|
|
& + Csum12222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(1,3,3,3,3,3),ep)=
|
|
& + Csum22222(ep)
|
|
&
|
|
|
|
Shat7(2,z6(1,3,3,3,3,3),ep)=
|
|
& + Csum22222(ep)
|
|
& + Cv(cc122222 + C124,ep)
|
|
&
|
|
|
|
Shat7(3,z6(1,3,3,3,3,3),ep)=
|
|
& + Csum22222(ep)
|
|
&
|
|
|
|
Shat7(1,z6(2,2,2,2,2,2),ep)=
|
|
& + Cv(cc111111 + C134,ep)
|
|
& - Cv(cc111111 + C234,ep)
|
|
&
|
|
|
|
Shat7(2,z6(2,2,2,2,2,2),ep)=
|
|
& - Cv(cc111111 + C234,ep)
|
|
&
|
|
|
|
Shat7(3,z6(2,2,2,2,2,2),ep)=
|
|
& - Cv(cc111111 + C234,ep)
|
|
& + Cv(cc222222 + C123,ep)
|
|
&
|
|
|
|
Shat7(1,z6(2,2,2,2,2,3),ep)=
|
|
& + Cv(cc111112 + C134,ep)
|
|
& - Cv(cc111112 + C234,ep)
|
|
&
|
|
|
|
Shat7(2,z6(2,2,2,2,2,3),ep)=
|
|
& - Cv(cc111112 + C234,ep)
|
|
&
|
|
|
|
Shat7(3,z6(2,2,2,2,2,3),ep)=
|
|
& - Cv(cc111112 + C234,ep)
|
|
&
|
|
|
|
Shat7(1,z6(2,2,2,2,3,3),ep)=
|
|
& + Cv(cc111122 + C134,ep)
|
|
& - Cv(cc111122 + C234,ep)
|
|
&
|
|
|
|
Shat7(2,z6(2,2,2,2,3,3),ep)=
|
|
& - Cv(cc111122 + C234,ep)
|
|
&
|
|
|
|
Shat7(3,z6(2,2,2,2,3,3),ep)=
|
|
& - Cv(cc111122 + C234,ep)
|
|
&
|
|
|
|
Shat7(1,z6(2,2,2,3,3,3),ep)=
|
|
& + Cv(cc111222 + C134,ep)
|
|
& - Cv(cc111222 + C234,ep)
|
|
&
|
|
|
|
Shat7(2,z6(2,2,2,3,3,3),ep)=
|
|
& - Cv(cc111222 + C234,ep)
|
|
&
|
|
|
|
Shat7(3,z6(2,2,2,3,3,3),ep)=
|
|
& - Cv(cc111222 + C234,ep)
|
|
&
|
|
|
|
Shat7(1,z6(2,2,3,3,3,3),ep)=
|
|
& + Cv(cc112222 + C134,ep)
|
|
& - Cv(cc112222 + C234,ep)
|
|
&
|
|
|
|
Shat7(2,z6(2,2,3,3,3,3),ep)=
|
|
& - Cv(cc112222 + C234,ep)
|
|
&
|
|
|
|
Shat7(3,z6(2,2,3,3,3,3),ep)=
|
|
& - Cv(cc112222 + C234,ep)
|
|
&
|
|
|
|
Shat7(1,z6(2,3,3,3,3,3),ep)=
|
|
& + Cv(cc122222 + C134,ep)
|
|
& - Cv(cc122222 + C234,ep)
|
|
&
|
|
|
|
Shat7(2,z6(2,3,3,3,3,3),ep)=
|
|
& - Cv(cc122222 + C234,ep)
|
|
&
|
|
|
|
Shat7(3,z6(2,3,3,3,3,3),ep)=
|
|
& - Cv(cc122222 + C234,ep)
|
|
&
|
|
|
|
Shat7(1,z6(3,3,3,3,3,3),ep)=
|
|
& + Cv(cc222222 + C134,ep)
|
|
& - Cv(cc222222 + C234,ep)
|
|
&
|
|
|
|
Shat7(2,z6(3,3,3,3,3,3),ep)=
|
|
& + Cv(cc222222 + C124,ep)
|
|
& - Cv(cc222222 + C234,ep)
|
|
&
|
|
|
|
Shat7(3,z6(3,3,3,3,3,3),ep)=
|
|
& - Cv(cc222222 + C234,ep)
|
|
|
|
|
|
Shat6(1,z5(1,1,1,1,1),ep)=
|
|
& + Csum0(ep)
|
|
& + 4.D0*Csum1(ep)
|
|
& + 4.D0*Csum2(ep)
|
|
& + 4.D0*Csum222(ep)
|
|
& + 4.D0*Csum111(ep)
|
|
& + 6.D0*Csum11(ep)
|
|
& + 6.D0*Csum22(ep)
|
|
& + 12.D0*Csum12(ep)
|
|
& + 12.D0*Csum112(ep)
|
|
& + 12.D0*Csum122(ep)
|
|
& + 6.D0*Csum1122(ep)
|
|
& + 4.D0*Csum1112(ep)
|
|
& + 4.D0*Csum1222(ep)
|
|
& + Csum1111(ep)
|
|
& + Csum2222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,1,1,1),ep)=
|
|
& + Csum0(ep)
|
|
& + 4.D0*Csum1(ep)
|
|
& + 4.D0*Csum2(ep)
|
|
& + 4.D0*Csum222(ep)
|
|
& + 4.D0*Csum111(ep)
|
|
& + 6.D0*Csum11(ep)
|
|
& + 6.D0*Csum22(ep)
|
|
& + 12.D0*Csum12(ep)
|
|
& + 12.D0*Csum112(ep)
|
|
& + 12.D0*Csum122(ep)
|
|
& + 6.D0*Csum1122(ep)
|
|
& + 4.D0*Csum1112(ep)
|
|
& + 4.D0*Csum1222(ep)
|
|
& + Csum1111(ep)
|
|
& + Csum2222(ep)
|
|
& + Cv(cc11111 + C124,ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,1,1,1),ep)=
|
|
& + Csum0(ep)
|
|
& + 4.D0*Csum1(ep)
|
|
& + 4.D0*Csum2(ep)
|
|
& + 4.D0*Csum222(ep)
|
|
& + 4.D0*Csum111(ep)
|
|
& + 6.D0*Csum11(ep)
|
|
& + 6.D0*Csum22(ep)
|
|
& + 12.D0*Csum12(ep)
|
|
& + 12.D0*Csum112(ep)
|
|
& + 12.D0*Csum122(ep)
|
|
& + 6.D0*Csum1122(ep)
|
|
& + 4.D0*Csum1112(ep)
|
|
& + 4.D0*Csum1222(ep)
|
|
& + Csum1111(ep)
|
|
& + Csum2222(ep)
|
|
& + Cv(cc11111 + C123,ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,1,1,2),ep)=
|
|
& - Csum1(ep)
|
|
& - 3.D0*Csum111(ep)
|
|
& - 3.D0*Csum11(ep)
|
|
& - 3.D0*Csum12(ep)
|
|
& - 6.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - 3.D0*Csum1112(ep)
|
|
& - Csum1222(ep)
|
|
& - Csum1111(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,1,1,2),ep)=
|
|
& - Csum1(ep)
|
|
& - 3.D0*Csum111(ep)
|
|
& - 3.D0*Csum11(ep)
|
|
& - 3.D0*Csum12(ep)
|
|
& - 6.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - 3.D0*Csum1112(ep)
|
|
& - Csum1222(ep)
|
|
& - Csum1111(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,1,1,2),ep)=
|
|
& - Csum1(ep)
|
|
& - 3.D0*Csum111(ep)
|
|
& - 3.D0*Csum11(ep)
|
|
& - 3.D0*Csum12(ep)
|
|
& - 6.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - 3.D0*Csum1112(ep)
|
|
& - Csum1222(ep)
|
|
& - Csum1111(ep)
|
|
& + Cv(cc11112 + C123,ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,1,1,3),ep)=
|
|
& - Csum2(ep)
|
|
& - 3.D0*Csum222(ep)
|
|
& - 3.D0*Csum22(ep)
|
|
& - 3.D0*Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 6.D0*Csum122(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - Csum1112(ep)
|
|
& - 3.D0*Csum1222(ep)
|
|
& - Csum2222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,1,1,3),ep)=
|
|
& - Csum2(ep)
|
|
& - 3.D0*Csum222(ep)
|
|
& - 3.D0*Csum22(ep)
|
|
& - 3.D0*Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 6.D0*Csum122(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - Csum1112(ep)
|
|
& - 3.D0*Csum1222(ep)
|
|
& - Csum2222(ep)
|
|
& + Cv(cc11112 + C124,ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,1,1,3),ep)=
|
|
& - Csum2(ep)
|
|
& - 3.D0*Csum222(ep)
|
|
& - 3.D0*Csum22(ep)
|
|
& - 3.D0*Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 6.D0*Csum122(ep)
|
|
& - 3.D0*Csum1122(ep)
|
|
& - Csum1112(ep)
|
|
& - 3.D0*Csum1222(ep)
|
|
& - Csum2222(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,1,2,2),ep)=
|
|
& + 2.D0*Csum111(ep)
|
|
& + Csum11(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + Csum1122(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + Csum1111(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,1,2,2),ep)=
|
|
& + 2.D0*Csum111(ep)
|
|
& + Csum11(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + Csum1122(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + Csum1111(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,1,2,2),ep)=
|
|
& + 2.D0*Csum111(ep)
|
|
& + Csum11(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + Csum1122(ep)
|
|
& + 2.D0*Csum1112(ep)
|
|
& + Csum1111(ep)
|
|
& + Cv(cc11122 + C123,ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,1,2,3),ep)=
|
|
& + Csum12(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + Csum1112(ep)
|
|
& + Csum1222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,1,2,3),ep)=
|
|
& + Csum12(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + Csum1112(ep)
|
|
& + Csum1222(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,1,2,3),ep)=
|
|
& + Csum12(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
& + 2.D0*Csum1122(ep)
|
|
& + Csum1112(ep)
|
|
& + Csum1222(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,1,3,3),ep)=
|
|
& + 2.D0*Csum222(ep)
|
|
& + Csum22(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
& + Csum1122(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum2222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,1,3,3),ep)=
|
|
& + 2.D0*Csum222(ep)
|
|
& + Csum22(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
& + Csum1122(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum2222(ep)
|
|
& + Cv(cc11122 + C124,ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,1,3,3),ep)=
|
|
& + 2.D0*Csum222(ep)
|
|
& + Csum22(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
& + Csum1122(ep)
|
|
& + 2.D0*Csum1222(ep)
|
|
& + Csum2222(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,2,2,2),ep)=
|
|
& - Csum111(ep)
|
|
& - Csum1112(ep)
|
|
& - Csum1111(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,2,2,2),ep)=
|
|
& - Csum111(ep)
|
|
& - Csum1112(ep)
|
|
& - Csum1111(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,2,2,2),ep)=
|
|
& - Csum111(ep)
|
|
& - Csum1112(ep)
|
|
& - Csum1111(ep)
|
|
& + Cv(cc11222 + C123,ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,2,2,3),ep)=
|
|
& - Csum112(ep)
|
|
& - Csum1122(ep)
|
|
& - Csum1112(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,2,2,3),ep)=
|
|
& - Csum112(ep)
|
|
& - Csum1122(ep)
|
|
& - Csum1112(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,2,2,3),ep)=
|
|
& - Csum112(ep)
|
|
& - Csum1122(ep)
|
|
& - Csum1112(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,2,3,3),ep)=
|
|
& - Csum122(ep)
|
|
& - Csum1122(ep)
|
|
& - Csum1222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,2,3,3),ep)=
|
|
& - Csum122(ep)
|
|
& - Csum1122(ep)
|
|
& - Csum1222(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,2,3,3),ep)=
|
|
& - Csum122(ep)
|
|
& - Csum1122(ep)
|
|
& - Csum1222(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,1,3,3,3),ep)=
|
|
& - Csum222(ep)
|
|
& - Csum1222(ep)
|
|
& - Csum2222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,1,3,3,3),ep)=
|
|
& - Csum222(ep)
|
|
& - Csum1222(ep)
|
|
& - Csum2222(ep)
|
|
& + Cv(cc11222 + C124,ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,1,3,3,3),ep)=
|
|
& - Csum222(ep)
|
|
& - Csum1222(ep)
|
|
& - Csum2222(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,2,2,2,2),ep)=
|
|
& + Csum1111(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,2,2,2,2),ep)=
|
|
& + Csum1111(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,2,2,2,2),ep)=
|
|
& + Csum1111(ep)
|
|
& + Cv(cc12222 + C123,ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,2,2,2,3),ep)=
|
|
& + Csum1112(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,2,2,2,3),ep)=
|
|
& + Csum1112(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,2,2,2,3),ep)=
|
|
& + Csum1112(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,2,2,3,3),ep)=
|
|
& + Csum1122(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,2,2,3,3),ep)=
|
|
& + Csum1122(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,2,2,3,3),ep)=
|
|
& + Csum1122(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,2,3,3,3),ep)=
|
|
& + Csum1222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,2,3,3,3),ep)=
|
|
& + Csum1222(ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,2,3,3,3),ep)=
|
|
& + Csum1222(ep)
|
|
&
|
|
|
|
Shat6(1,z5(1,3,3,3,3),ep)=
|
|
& + Csum2222(ep)
|
|
&
|
|
|
|
Shat6(2,z5(1,3,3,3,3),ep)=
|
|
& + Csum2222(ep)
|
|
& + Cv(cc12222 + C124,ep)
|
|
&
|
|
|
|
Shat6(3,z5(1,3,3,3,3),ep)=
|
|
& + Csum2222(ep)
|
|
&
|
|
|
|
Shat6(1,z5(2,2,2,2,2),ep)=
|
|
& + Cv(cc11111 + C134,ep)
|
|
& - Cv(cc11111 + C234,ep)
|
|
&
|
|
|
|
Shat6(2,z5(2,2,2,2,2),ep)=
|
|
& - Cv(cc11111 + C234,ep)
|
|
&
|
|
|
|
Shat6(3,z5(2,2,2,2,2),ep)=
|
|
& - Cv(cc11111 + C234,ep)
|
|
& + Cv(cc22222 + C123,ep)
|
|
&
|
|
|
|
Shat6(1,z5(2,2,2,2,3),ep)=
|
|
& + Cv(cc11112 + C134,ep)
|
|
& - Cv(cc11112 + C234,ep)
|
|
&
|
|
|
|
Shat6(2,z5(2,2,2,2,3),ep)=
|
|
& - Cv(cc11112 + C234,ep)
|
|
&
|
|
|
|
Shat6(3,z5(2,2,2,2,3),ep)=
|
|
& - Cv(cc11112 + C234,ep)
|
|
&
|
|
|
|
Shat6(1,z5(2,2,2,3,3),ep)=
|
|
& + Cv(cc11122 + C134,ep)
|
|
& - Cv(cc11122 + C234,ep)
|
|
&
|
|
|
|
Shat6(2,z5(2,2,2,3,3),ep)=
|
|
& - Cv(cc11122 + C234,ep)
|
|
&
|
|
|
|
Shat6(3,z5(2,2,2,3,3),ep)=
|
|
& - Cv(cc11122 + C234,ep)
|
|
&
|
|
|
|
Shat6(1,z5(2,2,3,3,3),ep)=
|
|
& + Cv(cc11222 + C134,ep)
|
|
& - Cv(cc11222 + C234,ep)
|
|
&
|
|
|
|
Shat6(2,z5(2,2,3,3,3),ep)=
|
|
& - Cv(cc11222 + C234,ep)
|
|
&
|
|
|
|
Shat6(3,z5(2,2,3,3,3),ep)=
|
|
& - Cv(cc11222 + C234,ep)
|
|
&
|
|
|
|
Shat6(1,z5(2,3,3,3,3),ep)=
|
|
& + Cv(cc12222 + C134,ep)
|
|
& - Cv(cc12222 + C234,ep)
|
|
&
|
|
|
|
Shat6(2,z5(2,3,3,3,3),ep)=
|
|
& - Cv(cc12222 + C234,ep)
|
|
&
|
|
|
|
Shat6(3,z5(2,3,3,3,3),ep)=
|
|
& - Cv(cc12222 + C234,ep)
|
|
&
|
|
|
|
Shat6(1,z5(3,3,3,3,3),ep)=
|
|
& + Cv(cc22222 + C134,ep)
|
|
& - Cv(cc22222 + C234,ep)
|
|
&
|
|
|
|
Shat6(2,z5(3,3,3,3,3),ep)=
|
|
& + Cv(cc22222 + C124,ep)
|
|
& - Cv(cc22222 + C234,ep)
|
|
&
|
|
|
|
Shat6(3,z5(3,3,3,3,3),ep)=
|
|
& - Cv(cc22222 + C234,ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,1,1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - 3.D0*Csum1(ep)
|
|
& - 3.D0*Csum2(ep)
|
|
& - Csum222(ep)
|
|
& - Csum111(ep)
|
|
& - 3.D0*Csum11(ep)
|
|
& - 3.D0*Csum22(ep)
|
|
& - 6.D0*Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,1,1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - 3.D0*Csum1(ep)
|
|
& - 3.D0*Csum2(ep)
|
|
& - Csum222(ep)
|
|
& - Csum111(ep)
|
|
& - 3.D0*Csum11(ep)
|
|
& - 3.D0*Csum22(ep)
|
|
& - 6.D0*Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& + Cv(cc1111 + C124,ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,1,1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - 3.D0*Csum1(ep)
|
|
& - 3.D0*Csum2(ep)
|
|
& - Csum222(ep)
|
|
& - Csum111(ep)
|
|
& - 3.D0*Csum11(ep)
|
|
& - 3.D0*Csum22(ep)
|
|
& - 6.D0*Csum12(ep)
|
|
& - 3.D0*Csum112(ep)
|
|
& - 3.D0*Csum122(ep)
|
|
& + Cv(cc1111 + C123,ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,1,1,2),ep)=
|
|
& + Csum1(ep)
|
|
& + Csum111(ep)
|
|
& + 2.D0*Csum11(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + Csum122(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,1,1,2),ep)=
|
|
& + Csum1(ep)
|
|
& + Csum111(ep)
|
|
& + 2.D0*Csum11(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + Csum122(ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,1,1,2),ep)=
|
|
& + Csum1(ep)
|
|
& + Csum111(ep)
|
|
& + 2.D0*Csum11(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + 2.D0*Csum112(ep)
|
|
& + Csum122(ep)
|
|
& + Cv(cc1112 + C123,ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,1,1,3),ep)=
|
|
& + Csum2(ep)
|
|
& + Csum222(ep)
|
|
& + 2.D0*Csum22(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + Csum112(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,1,1,3),ep)=
|
|
& + Csum2(ep)
|
|
& + Csum222(ep)
|
|
& + 2.D0*Csum22(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + Csum112(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
& + Cv(cc1112 + C124,ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,1,1,3),ep)=
|
|
& + Csum2(ep)
|
|
& + Csum222(ep)
|
|
& + 2.D0*Csum22(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + Csum112(ep)
|
|
& + 2.D0*Csum122(ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,1,2,2),ep)=
|
|
& - Csum111(ep)
|
|
& - Csum11(ep)
|
|
& - Csum112(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,1,2,2),ep)=
|
|
& - Csum111(ep)
|
|
& - Csum11(ep)
|
|
& - Csum112(ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,1,2,2),ep)=
|
|
& - Csum111(ep)
|
|
& - Csum11(ep)
|
|
& - Csum112(ep)
|
|
& + Cv(cc1122 + C123,ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,1,2,3),ep)=
|
|
& - Csum12(ep)
|
|
& - Csum112(ep)
|
|
& - Csum122(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,1,2,3),ep)=
|
|
& - Csum12(ep)
|
|
& - Csum112(ep)
|
|
& - Csum122(ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,1,2,3),ep)=
|
|
& - Csum12(ep)
|
|
& - Csum112(ep)
|
|
& - Csum122(ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,1,3,3),ep)=
|
|
& - Csum222(ep)
|
|
& - Csum22(ep)
|
|
& - Csum122(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,1,3,3),ep)=
|
|
& - Csum222(ep)
|
|
& - Csum22(ep)
|
|
& - Csum122(ep)
|
|
& + Cv(cc1122 + C124,ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,1,3,3),ep)=
|
|
& - Csum222(ep)
|
|
& - Csum22(ep)
|
|
& - Csum122(ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,2,2,2),ep)=
|
|
& + Csum111(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,2,2,2),ep)=
|
|
& + Csum111(ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,2,2,2),ep)=
|
|
& + Csum111(ep)
|
|
& + Cv(cc1222 + C123,ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,2,2,3),ep)=
|
|
& + Csum112(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,2,2,3),ep)=
|
|
& + Csum112(ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,2,2,3),ep)=
|
|
& + Csum112(ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,2,3,3),ep)=
|
|
& + Csum122(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,2,3,3),ep)=
|
|
& + Csum122(ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,2,3,3),ep)=
|
|
& + Csum122(ep)
|
|
&
|
|
|
|
Shat5(1,z4(1,3,3,3),ep)=
|
|
& + Csum222(ep)
|
|
&
|
|
|
|
Shat5(2,z4(1,3,3,3),ep)=
|
|
& + Csum222(ep)
|
|
& + Cv(cc1222 + C124,ep)
|
|
&
|
|
|
|
Shat5(3,z4(1,3,3,3),ep)=
|
|
& + Csum222(ep)
|
|
&
|
|
|
|
Shat5(1,z4(2,2,2,2),ep)=
|
|
& + Cv(cc1111 + C134,ep)
|
|
& - Cv(cc1111 + C234,ep)
|
|
&
|
|
|
|
Shat5(2,z4(2,2,2,2),ep)=
|
|
& - Cv(cc1111 + C234,ep)
|
|
&
|
|
|
|
Shat5(3,z4(2,2,2,2),ep)=
|
|
& - Cv(cc1111 + C234,ep)
|
|
& + Cv(cc2222 + C123,ep)
|
|
&
|
|
|
|
Shat5(1,z4(2,2,2,3),ep)=
|
|
& + Cv(cc1112 + C134,ep)
|
|
& - Cv(cc1112 + C234,ep)
|
|
&
|
|
|
|
Shat5(2,z4(2,2,2,3),ep)=
|
|
& - Cv(cc1112 + C234,ep)
|
|
&
|
|
|
|
Shat5(3,z4(2,2,2,3),ep)=
|
|
& - Cv(cc1112 + C234,ep)
|
|
&
|
|
|
|
Shat5(1,z4(2,2,3,3),ep)=
|
|
& + Cv(cc1122 + C134,ep)
|
|
& - Cv(cc1122 + C234,ep)
|
|
&
|
|
|
|
Shat5(2,z4(2,2,3,3),ep)=
|
|
& - Cv(cc1122 + C234,ep)
|
|
&
|
|
|
|
Shat5(3,z4(2,2,3,3),ep)=
|
|
& - Cv(cc1122 + C234,ep)
|
|
&
|
|
|
|
Shat5(1,z4(2,3,3,3),ep)=
|
|
& + Cv(cc1222 + C134,ep)
|
|
& - Cv(cc1222 + C234,ep)
|
|
&
|
|
|
|
Shat5(2,z4(2,3,3,3),ep)=
|
|
& - Cv(cc1222 + C234,ep)
|
|
&
|
|
|
|
Shat5(3,z4(2,3,3,3),ep)=
|
|
& - Cv(cc1222 + C234,ep)
|
|
&
|
|
|
|
Shat5(1,z4(3,3,3,3),ep)=
|
|
& + Cv(cc2222 + C134,ep)
|
|
& - Cv(cc2222 + C234,ep)
|
|
&
|
|
|
|
Shat5(2,z4(3,3,3,3),ep)=
|
|
& + Cv(cc2222 + C124,ep)
|
|
& - Cv(cc2222 + C234,ep)
|
|
&
|
|
|
|
Shat5(3,z4(3,3,3,3),ep)=
|
|
& - Cv(cc2222 + C234,ep)
|
|
&
|
|
|
|
Shat4(1,z3(1,1,1),ep)=
|
|
& + Csum0(ep)
|
|
& + 2.D0*Csum1(ep)
|
|
& + 2.D0*Csum2(ep)
|
|
& + Csum11(ep)
|
|
& + Csum22(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
&
|
|
|
|
Shat4(2,z3(1,1,1),ep)=
|
|
& + Csum0(ep)
|
|
& + 2.D0*Csum1(ep)
|
|
& + 2.D0*Csum2(ep)
|
|
& + Csum11(ep)
|
|
& + Csum22(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + Cv(cc111 + C124,ep)
|
|
&
|
|
|
|
Shat4(3,z3(1,1,1),ep)=
|
|
& + Csum0(ep)
|
|
& + 2.D0*Csum1(ep)
|
|
& + 2.D0*Csum2(ep)
|
|
& + Csum11(ep)
|
|
& + Csum22(ep)
|
|
& + 2.D0*Csum12(ep)
|
|
& + Cv(cc111 + C123,ep)
|
|
&
|
|
|
|
Shat4(1,z3(2,2,2),ep)=
|
|
& + Cv(cc111 + C134,ep)
|
|
& - Cv(cc111 + C234,ep)
|
|
&
|
|
|
|
Shat4(2,z3(2,2,2),ep)=
|
|
& - Cv(cc111 + C234,ep)
|
|
&
|
|
|
|
Shat4(3,z3(2,2,2),ep)=
|
|
& - Cv(cc111 + C234,ep)
|
|
& + Cv(cc222 + C123,ep)
|
|
&
|
|
|
|
Shat4(1,z3(3,3,3),ep)=
|
|
& + Cv(cc222 + C134,ep)
|
|
& - Cv(cc222 + C234,ep)
|
|
&
|
|
|
|
Shat4(2,z3(3,3,3),ep)=
|
|
& + Cv(cc222 + C124,ep)
|
|
& - Cv(cc222 + C234,ep)
|
|
&
|
|
|
|
Shat4(3,z3(3,3,3),ep)=
|
|
& - Cv(cc222 + C234,ep)
|
|
&
|
|
|
|
Shat4(1,z3(1,1,2),ep)=
|
|
& - Csum1(ep)
|
|
& - Csum11(ep)
|
|
& - Csum12(ep)
|
|
&
|
|
|
|
Shat4(2,z3(1,1,2),ep)=
|
|
& - Csum1(ep)
|
|
& - Csum11(ep)
|
|
& - Csum12(ep)
|
|
&
|
|
|
|
Shat4(3,z3(1,1,2),ep)=
|
|
& - Csum1(ep)
|
|
& - Csum11(ep)
|
|
& - Csum12(ep)
|
|
& + Cv(cc112 + C123,ep)
|
|
&
|
|
|
|
Shat4(1,z3(1,1,3),ep)=
|
|
& - Csum2(ep)
|
|
& - Csum22(ep)
|
|
& - Csum12(ep)
|
|
&
|
|
|
|
Shat4(2,z3(1,1,3),ep)=
|
|
& - Csum2(ep)
|
|
& - Csum22(ep)
|
|
& - Csum12(ep)
|
|
& + Cv(cc112 + C124,ep)
|
|
&
|
|
|
|
Shat4(3,z3(1,1,3),ep)=
|
|
& - Csum2(ep)
|
|
& - Csum22(ep)
|
|
& - Csum12(ep)
|
|
&
|
|
|
|
Shat4(1,z3(1,2,2),ep)=
|
|
& + Csum11(ep)
|
|
&
|
|
|
|
Shat4(2,z3(1,2,2),ep)=
|
|
& + Csum11(ep)
|
|
&
|
|
|
|
Shat4(3,z3(1,2,2),ep)=
|
|
& + Csum11(ep)
|
|
& + Cv(cc122 + C123,ep)
|
|
&
|
|
|
|
Shat4(1,z3(1,3,3),ep)=
|
|
& + Csum22(ep)
|
|
&
|
|
|
|
Shat4(2,z3(1,3,3),ep)=
|
|
& + Csum22(ep)
|
|
& + Cv(cc122 + C124,ep)
|
|
&
|
|
|
|
Shat4(3,z3(1,3,3),ep)=
|
|
& + Csum22(ep)
|
|
&
|
|
|
|
Shat4(1,z3(2,2,3),ep)=
|
|
& + Cv(cc112 + C134,ep)
|
|
& - Cv(cc112 + C234,ep)
|
|
&
|
|
|
|
Shat4(2,z3(2,2,3),ep)=
|
|
& - Cv(cc112 + C234,ep)
|
|
&
|
|
|
|
Shat4(3,z3(2,2,3),ep)=
|
|
& - Cv(cc112 + C234,ep)
|
|
&
|
|
|
|
Shat4(1,z3(2,3,3),ep)=
|
|
& + Cv(cc122 + C134,ep)
|
|
& - Cv(cc122 + C234,ep)
|
|
&
|
|
|
|
Shat4(2,z3(2,3,3),ep)=
|
|
& - Cv(cc122 + C234,ep)
|
|
&
|
|
|
|
Shat4(3,z3(2,3,3),ep)=
|
|
& - Cv(cc122 + C234,ep)
|
|
&
|
|
|
|
Shat4(1,z3(1,2,3),ep)=
|
|
& + Csum12(ep)
|
|
&
|
|
|
|
Shat4(2,z3(1,2,3),ep)=
|
|
& + Csum12(ep)
|
|
&
|
|
|
|
Shat4(3,z3(1,2,3),ep)=
|
|
& + Csum12(ep)
|
|
&
|
|
|
|
Shat3(1,z2(1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - Csum1(ep)
|
|
& - Csum2(ep)
|
|
&
|
|
|
|
Shat3(2,z2(1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - Csum1(ep)
|
|
& - Csum2(ep)
|
|
& + Cv(cc11 + C124,ep)
|
|
&
|
|
|
|
Shat3(3,z2(1,1),ep)=
|
|
& - Csum0(ep)
|
|
& - Csum1(ep)
|
|
& - Csum2(ep)
|
|
& + Cv(cc11 + C123,ep)
|
|
&
|
|
|
|
Shat3(1,z2(2,2),ep)=
|
|
& + Cv(cc11 + C134,ep)
|
|
& - Cv(cc11 + C234,ep)
|
|
&
|
|
|
|
Shat3(2,z2(2,2),ep)=
|
|
& - Cv(cc11 + C234,ep)
|
|
&
|
|
|
|
Shat3(3,z2(2,2),ep)=
|
|
& - Cv(cc11 + C234,ep)
|
|
& + Cv(cc22 + C123,ep)
|
|
&
|
|
|
|
Shat3(1,z2(3,3),ep)=
|
|
& + Cv(cc22 + C134,ep)
|
|
& - Cv(cc22 + C234,ep)
|
|
&
|
|
|
|
Shat3(2,z2(3,3),ep)=
|
|
& + Cv(cc22 + C124,ep)
|
|
& - Cv(cc22 + C234,ep)
|
|
&
|
|
|
|
Shat3(3,z2(3,3),ep)=
|
|
& - Cv(cc22 + C234,ep)
|
|
&
|
|
|
|
Shat3(1,z2(1,2),ep)=
|
|
& + Csum1(ep)
|
|
&
|
|
|
|
Shat3(2,z2(1,2),ep)=
|
|
& + Csum1(ep)
|
|
&
|
|
|
|
Shat3(3,z2(1,2),ep)=
|
|
& + Csum1(ep)
|
|
& + Cv(cc12 + C123,ep)
|
|
&
|
|
|
|
Shat3(1,z2(1,3),ep)=
|
|
& + Csum2(ep)
|
|
&
|
|
|
|
Shat3(2,z2(1,3),ep)=
|
|
& + Csum2(ep)
|
|
& + Cv(cc12 + C124,ep)
|
|
&
|
|
|
|
Shat3(3,z2(1,3),ep)=
|
|
& + Csum2(ep)
|
|
&
|
|
|
|
Shat3(1,z2(2,3),ep)=
|
|
& + Cv(cc12 + C134,ep)
|
|
& - Cv(cc12 + C234,ep)
|
|
&
|
|
|
|
Shat3(2,z2(2,3),ep)=
|
|
& - Cv(cc12 + C234,ep)
|
|
&
|
|
|
|
Shat3(3,z2(2,3),ep)=
|
|
& - Cv(cc12 + C234,ep)
|
|
&
|
|
|
|
Shat2(1,1,ep)=
|
|
& + Csum0(ep)
|
|
&
|
|
|
|
Shat2(2,1,ep)=
|
|
& + Csum0(ep)
|
|
& + Cv(cc1 + C124,ep)
|
|
&
|
|
|
|
Shat2(3,1,ep)=
|
|
& + Csum0(ep)
|
|
& + Cv(cc1 + C123,ep)
|
|
&
|
|
|
|
Shat2(1,2,ep)=
|
|
& + Cv(cc1 + C134,ep)
|
|
& - Cv(cc1 + C234,ep)
|
|
&
|
|
|
|
Shat2(2,2,ep)=
|
|
& - Cv(cc1 + C234,ep)
|
|
&
|
|
|
|
Shat2(3,2,ep)=
|
|
& - Cv(cc1 + C234,ep)
|
|
& + Cv(cc2 + C123,ep)
|
|
&
|
|
|
|
Shat2(1,3,ep)=
|
|
& + Cv(cc2 + C134,ep)
|
|
& - Cv(cc2 + C234,ep)
|
|
&
|
|
|
|
Shat2(2,3,ep)=
|
|
& + Cv(cc2 + C124,ep)
|
|
& - Cv(cc2 + C234,ep)
|
|
&
|
|
|
|
Shat2(3,3,ep)=
|
|
& - Cv(cc2 + C234,ep)
|
|
&
|
|
|
|
Shat1(1,ep) =
|
|
& + Cv(cc0 + C134,ep)
|
|
& - Cv(cc0 + C234,ep)
|
|
&
|
|
|
|
Shat1(2,ep) =
|
|
& + Cv(cc0 + C124,ep)
|
|
& - Cv(cc0 + C234,ep)
|
|
&
|
|
|
|
Shat1(3,ep) =
|
|
& + Cv(cc0 + C123,ep)
|
|
& - Cv(cc0 + C234,ep)
|
|
&
|
|
|