$ set Noon %SET-W-NOTSET, error modifying TO49AX$DKB100: -SET-E-INVDEV, device is invalid for requested operation FORM version 2.3 Dec 10 1997 i i,j; i ab1,ab2; i mu,nu,al,be; i b1,b2,b3,b4; i a1,a2,a3,a4,a5,a6; i aa1,aa2,aa3,aa4,aa5,aa6,aa7; i bb1,bb2,bb3,bb4,bb5,bb6,bb7; * v pp,qp,Qp,Qm,k,kb,BQ; * s n; s zi; s I0; s PZ; s XX; s PZs; s Com; s gexp; s GT,T; s selv; s Xpi; s selb; s n1,n2; s [2-3]; s [5-6]; s [2+3]; s [5+6]; s eq,egt; s Qu,Qd,srt; s denk,denkb; s one,aux,int; s ups,x1,x2,z; s rl1,rl2,rl3; s sw,swC,Gw,GwC; s gramd,restore; s l1,l2,l,lp,lv; s vv,vvs,vvc,vvq; s sel0,seld,seldd; s O,Os,Oc,Oq,Op,Om; s fws,fwt,fwsC,fwtC; s um,dm,ums,dms,QmS,QpS; s omx1,x12,x2z,ux; s Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2; s Pws,Pwt,PwsC,PwtC,PwsM2,PwtM2; s [fws-fwt],[fwsC-fwtC],[Sigmas-Sigmat]; s Sigmas,SigmasC,Sigmat,SigmatC,SigmasM2; s [1+X],[x2-z],[ups+X+x1],[x1-x2],[1-x1]: s Xps,Xpis,yc,ys,omx1s,omx1c,omx1q,x12s,x12c,x2zq,x2zc,x2zs,uxs; s y,s,rs,sp,t,tp,u,up,kp,km,zp,zm,tM,uM,X,zX,Xp,Xph,tuM,tdM,Y,zXi; s cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts,sst; * f sl,sq,slC,sqC; f sum1_gg,sum1_zz,sum2; f wsub_gg,wsub_gz,wsub_zz; * cf es; cf fun; cf Rep; cf redf; cf kbf,qpf; cf ae,au,ad; cf tv1g,tv1gC; cf tv1z,tv1zC; cf Mods,Re,Im; cf tv0,tv1,tv0C,tv1C; cf ve,vu,vd,ss,cs,csi; cf vvvg,vvvgC,wg,wgC,wgM2; cf vvvz,vvvzC,wz,wzC,wzM2; cf I_uu_21,I_uu_22,I_uu_23,I_uu_24; cf I_dd_21,I_dd_22,I_dd_23,I_dd_24; cf I_nn_21,I_nn_22,I_nn_23,I_nn_24; cf rs1wg,rs2wg,rt1wg,rt2wg,rwwg; cf rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz; cf I_un_21,I_un_22,I_un_23,I_un_24,I_un_25; cf I_dn_21,I_dn_22,I_dn_23,I_dn_24,I_dn_25; cf I_u_11,I_u_12,I_d_11,I_d_12,I_n_11,I_n_12; cf C00,B00t,B00s,B00g,C00C,B00tC,B00sC,B00gC; * .global * *---- * #procedure wtensor(compo) * * First W_{mu,nu} is constructed, gg or gz or zz * g w_'compo'= wsub_'compo'(mu,nu)*int; * id wsub_gg(mu?,nu?)= sum_(l1,1,4,sum1_gg(l1,mu,nu)); * id wsub_zz(mu?,nu?)= sum_(l1,5,9,sum1_zz(l1,mu,nu)); * id wsub_gz(mu?,nu?)= sum_(l1,1,4,sum1_zz(l1,mu,nu))+ sum_(l1,5,9,sum1_gg(l1,mu,nu)); * id sum1_gg(l?,mu?,nu?)= sum_(l2,1,4,sum2(l,l2,mu,nu)); id sum1_zz(l?,mu?,nu?)= sum_(l2,5,9,sum2(l,l2,mu,nu)); * id sum2(l?,lp?,mu?,nu?)= -sqC(i,l,nu)*(-i_*g_(i,k)+um*gi_(i))* sq(i,lp,mu)*(-i_*g_(i,kb)-dm*gi_(i))* slC(j,l,nu)*g_(j,pp)*sl(j,lp,mu)* (g_(j,pp)-g_(j,Qp)); * #call diagrams{dummy} * trace4,i; trace4,j; contract; * .store * #call kinematics{dummy} #call lin{'compo'} * save wrl$grp04:[giampiero]outputampgg_B.sav wwww_'compo'; *save wrl$grp04:[giampiero]outputampgz_B.sav wwww_'compo'; *save wrl$grp04:[giampiero]outputampzz_B.sav wwww_'compo'; * #endprocedure * *--------------------------------------------------------------- * #procedure diagrams(dummy) * * The diagrams are constructed * id sl(i?,1,mu?)= i_/8*Gt* (tv0(mu,a1,a5)+Gs*(tv1g(mu,a1,a5)-Sigmas*tv0(mu,a1,a5)))* es(T)*g_(i,a1)*g6_(i)*vvs/Pws/Pwt/s^2; id slC(i?,1,nu?)= i_/8*GtC* (tv0C(nu,b1,a6)+GsC*(tv1gC(nu,b1,a6)-SigmasC*tv0C(nu,b1,a6) ))* es(T)*g_(i,b1)*g6_(i)*vvs/PwsC/PwtC/s^2; * id sl(i?,5,mu?)= i_/32*Gt* (tv0(mu,aa1,aa5)+Gs*(tv1z(mu,aa1,aa5)-Sigmas*tv0(mu,aa1,aa5)))* GT/PZ*g_(i,aa1)*g6_(i)*vvs/Pws/Pwt/s^3; id slC(i?,5,nu?)= i_/32*GtC* (tv0C(nu,bb1,aa6)+GsC*(tv1zC(nu,bb1,aa6)-SigmasC*tv0C(nu,bb1,aa6) ))* GT/PZ*g_(i,bb1)*g6_(i)*vvs/PwsC/PwtC/s^3; * al sl(i?,2,mu?)= i_/(8*y*s)*Gs*es(T)* g_(i,mu)*(g_(i,pp)+g_(i,Qm))*g_(i,a2)*g6_(i)*vv/Pws/s; al slC(i?,2,nu?)= i_/(8*y*s)*GsC*es(T)* g_(i,b2)*g6_(i)*(g_(i,pp)+g_(i,Qm))*g_(i,nu)*vv/PwsC/s; * al sl(i?,6,mu?)= -i_/(8*y*s)*Gs*GT*csi(T)/16/PZ* g_(i,mu)*(ve(T)*g6_(i)+ae(T)*g7_(i))* (g_(i,pp)+g_(i,Qm))*g_(i,aa2)*g6_(i)*vv/Pws/s^2; al slC(i?,6,nu?)= -i_/(8*y*s)*GsC*GT*csi(T)/16/PZ* g_(i,bb2)*g6_(i)*(g_(i,pp)+g_(i,Qm))* g_(i,nu)*(ve(T)*g6_(i)+ae(T)*g7_(i))*vv/PwsC/s^2; * al sl(i?,3,mu?)= i_/8*Qu*Gt*es(T)* g_(i,a3)*g6_(i)*vv/Pwt/s/denk; al slC(i?,3,nu?)= i_/8*Qu*GtC*es(T)* g_(i,b3)*g6_(i)*vv/PwtC/s/denk; * al sl(i?,7,mu?)= -i_/8*Gt*GT*csi(T)/16/PZ* g_(i,aa3)*g6_(i)*vv/Pwt/s^2/denk; al slC(i?,7,nu?)= -i_/8*GtC*GT*csi(T)/16/PZ* g_(i,bb3)*g6_(i)*vv/PwtC/s^2/denk; * al sl(i?,4,mu?)= -i_/8*Qd*Gt*es(T)* g_(i,a4)*g6_(i)*vv/Pwt/s/denkb; al slC(i?,4,nu?)= -i_/8*Qd*GtC*es(T)* g_(i,b4)*g6_(i)*vv/PwtC/s/denkb; * al sl(i?,8,mu?)= -i_/8*Gt*GT*csi(T)/16/PZ* g_(i,aa4)*g6_(i)*vv/Pwt/s^2/denkb; al slC(i?,8,nu?)= -i_/8*GtC*GT*csi(T)/16/PZ* g_(i,bb4)*g6_(i)*vv/PwtC/s^2/denkb; * al sl(i?,9,mu?)= i_/32/km*Gs*GT*csi(T)/PZ* g_(i,aa7)*g6_(i)*(g_(i,Qm)-g_(i,qp))* g_(i,mu)*g6_(i)*vv/Pws/s^2/4; al slC(i?,9,nu?)= i_/32/km*GsC*GT*csi(T)/PZ* g_(i,nu)*g6_(i)*(g_(i,Qm)-g_(i,qp))* g_(i,bb7)*g6_(i)*vv/PwsC/s^2/4; * id sq(i?,1,mu?)= g_(i,a5)*g6_(i); id sqC(i?,1,nu?)= -g_(i,a6)*g6_(i); * al sq(i?,5,mu?)= g_(i,aa5)*g6_(i); al sqC(i?,5,nu?)= -g_(i,aa6)*g6_(i); * al sq(i?,2,mu?)= g_(i,a2)*g6_(i); al sqC(i?,2,nu?)= -g_(i,b2)*g6_(i); * al sq(i?,6,mu?)= g_(i,aa2)*g6_(i); al sqC(i?,6,nu?)= -g_(i,bb2)*g6_(i); * al sq(i?,3,mu?)= g_(i,mu)*(g_(i,Qm)-g_(i,k)-i_*um*gi_(i))*g_(i,a3)*g6_(i); al sqC(i?,3,nu?)= g_(i,b3)*g6_(i)*(-(g_(i,Qm)-g_(i,k))+i_*um*gi_(i))*g_(i,n u); * al sq(i?,7,mu?)= g_(i,mu)*(vu(T)*g6_(i)+au(T)*g7_(i))* (g_(i,Qm)-g_(i,k)-i_*um*gi_(i))*g_(i,aa3)*g6_(i); al sqC(i?,7,nu?)= g_(i,bb3)*g6_(i)*(-(g_(i,Qm)-g_(i,k))+i_*um*gi_(i))* g_(i,nu)*(vu(T)*g6_(i)+au(T)*g7_(i)); * al sq(i?,4,mu?)= g_(i,a4)*g6_(i)*(g_(i,Qm)-g_(i,kb)+i_*dm*gi_(i))*g_(i,mu); al sqC(i?,4,nu?)= g_(i,nu)*(-(g_(i,Qm)-g_(i,kb))-i_*dm*gi_(i))*g_(i,b4)*g6_ (i); * al sq(i?,8,mu?)= g_(i,aa4)*g6_(i)*(g_(i,Qm)-g_(i,kb)+i_*dm*gi_(i))* g_(i,mu)*(vd(T)*g6_(i)+ad(T)*g7_(i)) ; al sqC(i?,8,nu?)= g_(i,nu)*(vd(T)*g6_(i)+ad(T)*g7_(i))* (-(g_(i,Qm)-g_(i,kb))-i_*dm*gi_(i))*g_(i,bb4)*g6_(i); * al sq(i?,9,mu?)= g_(i,aa7)*g6_(i); al sqC(i?,9,nu?)= -g_(i,bb7)*g6_(i); * #endprocedure * .store * #procedure kinematics(dummy) * * Mostly kinematics * g ww_'compo'= w_'compo'; *nwrite statistics; * id Qu= 2/3; al Qd= -1/3; id srt^-2= 1/2; * id pp.pp= 0; al qp.qp= 0; al k.k= -ums; al kb.kb= -dms; al Qm.Qm= QmS; id um^2= ums; al dm^2= dms; * *--------------------------------------------------------------------- id Gs= gexp*Gs; al GsC= gexp*GsC; id Gt= gexp*Gt; al GtC= gexp*GtC; *--------------------------------------------------------------------- * id Gs*GsC= GsM2; al Gt*GtC= GtM2; * id Pws^-1*PwsC^-1= 1/PwsM2; id Pwt^-1*PwtC^-1= 1/PwtM2; id Pws^-1= PwsC/PwsM2; id PwsC^-1= Pws/PwsM2; id Pwt^-1= PwtC/PwtM2; id PwtC^-1= Pwt/PwtM2; * * Triple-vertex is substituted * id tv0(mu?,al?,be?)= d_(mu,al)*(Qm(be)-Qp(be))+ d_(al,be)*(2*Qp(mu)+Qm(mu))- d_(mu,be)*(2*Qm(al)+Qp(al)); al tv0C(mu?,al?,be?)= d_(mu,al)*(Qm(be)-Qp(be))+ d_(al,be)*(2*Qp(mu)+Qm(mu))- d_(mu,be)*(2*Qm(al)+Qp(al)); id tv1g(mu?,al?,be?)= sum_(lv,1,10,(vvvg(lv,mu,al,be))); id tv1gC(mu?,al?,be?)= sum_(lv,1,10,(vvvgC(lv,mu,al,be))); id tv1z(mu?,al?,be?)= sum_(lv,1,10,(vvvz(lv,mu,al,be))); id tv1zC(mu?,al?,be?)= sum_(lv,1,10,(vvvzC(lv,mu,al,be))); * * Form factors are introduced * id vvvg(1,mu?,al?,be?)= 4*wg(1)/rs*(d_(al,be)*Qp(mu)+ d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al))); id vvvg(2,mu?,al?,be?)= 2*wg(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(a l))); id vvvg(3,mu?,al?,be?)= 2*wg(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)*(Qm(al)+Qp(a l))); id vvvg(4,mu?,al?,be?)= 2*wg(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu); id vvvg(5,mu?,al?,be?)= wg(5)/rs*e_(Qm,mu,al,be); id vvvg(6,mu?,al?,be?)= wg(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be)); id vvvg(7,mu?,al?,be?)= wg(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu)); id vvvg(8,mu?,al?,be?)= 2*wg(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al)); id vvvg(9,mu?,al?,be?)= 2*wg(9)/rs*d_(al,be)*Qm(mu); id vvvg(10,mu?,al?,be?)= -wg(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu); * id vvvgC(1,mu?,al?,be?)= 4*wgC(1)/rs*(d_(al,be)*Qp(mu)+ d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al))); id vvvgC(2,mu?,al?,be?)= 2*wgC(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)* (Qm(al)+Qp(al))); id vvvgC(3,mu?,al?,be?)= 2*wgC(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)* (Qm(al)+Qp(al))); id vvvgC(4,mu?,al?,be?)= 2*wgC(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu); id vvvgC(5,mu?,al?,be?)= -wgC(5)/rs*e_(Qm,mu,al,be); id vvvgC(6,mu?,al?,be?)= -wgC(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be)); id vvvgC(7,mu?,al?,be?)= -wgC(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu)) ; id vvvgC(8,mu?,al?,be?)= 2*wgC(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al)); id vvvgC(9,mu?,al?,be?)= 2*wgC(9)/rs*d_(al,be)*Qm(mu); id vvvgC(10,mu?,al?,be?)= wgC(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu); * id vvvz(1,mu?,al?,be?)= 4*wz(1)/rs*(d_(al,be)*Qp(mu)+ d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al))); id vvvz(2,mu?,al?,be?)= 2*wz(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(a l))); id vvvz(3,mu?,al?,be?)= 2*wz(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)*(Qm(al)+Qp(a l))); id vvvz(4,mu?,al?,be?)= 2*wz(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu); id vvvz(5,mu?,al?,be?)= wz(5)/rs*e_(Qm,mu,al,be); id vvvz(6,mu?,al?,be?)= wz(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be)); id vvvz(7,mu?,al?,be?)= wz(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu)); id vvvz(8,mu?,al?,be?)= 2*wz(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al)); id vvvz(9,mu?,al?,be?)= 2*wz(9)/rs*d_(al,be)*Qm(mu); id vvvz(10,mu?,al?,be?)= -wz(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu); * id vvvzC(1,mu?,al?,be?)= 4*wzC(1)/rs*(d_(al,be)*Qp(mu)+ d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al))); id vvvzC(2,mu?,al?,be?)= 2*wzC(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)* (Qm(al)+Qp(al))); id vvvzC(3,mu?,al?,be?)= 2*wzC(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)* (Qm(al)+Qp(al))); id vvvzC(4,mu?,al?,be?)= 2*wzC(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu); id vvvzC(5,mu?,al?,be?)= -wzC(5)/rs*e_(Qm,mu,al,be); id vvvzC(6,mu?,al?,be?)= -wzC(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be)); id vvvzC(7,mu?,al?,be?)= -wzC(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu)) ; id vvvzC(8,mu?,al?,be?)= 2*wzC(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al)); id vvvzC(9,mu?,al?,be?)= 2*wzC(9)/rs*d_(al,be)*Qm(mu); id vvvzC(10,mu?,al?,be?)= wzC(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu); * contract; * id Sigmas= 0; al SigmasC= 0; *id Sigmas*SigmasC= SigmasM2; * id BQ(mu?)= pp(mu)-qp(mu)+Qm(mu); id Qp(mu?)= pp(mu)-qp(mu); * id e_(a1?,Qp,a2?,a3?)= e_(a1,pp,a2,a3)-e_(a1,qp,a2,a3); id e_(a1?,a2?,Qp,a3?)= e_(a1,a2,pp,a3)-e_(a1,a2,qp,a3); id e_(a1?,a2?,a3?,Qp)= e_(a1,a2,a3,pp)-e_(a1,a2,a3,qp); * id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; * .store * #endprocedure * #procedure lin(compo) * * Linear independent invariants are introduced * g www_'compo'= ww_'compo'*selb; *nwrite statistics; * contract; * id k.k= -ums; al kb.kb= -dms; * id denk^-1*denkb^-1= -1/2/(Qp.Qm)*(1/denk+1/denkb); * id pp.pp= 0; al qp.qp= 0; al k.k= -ums; al kb.kb= -dms; al Qm.Qm= QmS; * id selb*denkb^-2= 1/denkb^2; id selb*denkb^-1= 1/denkb; id selb= 0; * id denkb^-1= -1/uM; * id e_(k,a1?,a2?,a3?)= e_(pp,a1,a2,a3)-e_(qp,a1,a2,a3)+e_(Qm,a1,a2,a3)- e_(kb,a1,a2,a3); id k(mu?)= pp(mu)-qp(mu)+Qm(mu)-kb(mu); id QpS= -2*pp.qp; id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; * id pp.Qm= -1/2*(QmS+y*s); al k.kb= 1/2*(-sp+ums+dms); al Qm.k= 1/2*(QmS+t-ums); al pp.kb= 1/2*(tp-dms); al Qm.kb= 1/2*(QmS+u-dms); al pp.k= 1/2*(up-ums); al pp.qp= 1/2*kp; al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(dms-zp); al k.qp= 1/2*(ums-zm); * id sp= x2*y*s; al km= z*y*s; al km^-1= 1/(z*y*s); al zm= -(x2-x1)*y*s; id QmS= X*y*s; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * .store * g wwww_'compo'= www_'compo'; *nwrite statistics; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id Qp.Qm^-1= -2/(zX*y*s); ratio,z,zX,aux; id aux= 1+2*X; id aux^-1= 1/2/Xph; ratio,X,Xph,aux; id aux= 1/2; al aux^-1= 2; * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * .store * #endprocedure * #procedure integration(compo) * * Integration starts * id kbf(mu?)*kbf(nu?)*int= I_dd_21*pp(mu)*pp(nu)+I_dd_22*Qm(mu)*Qm(nu)+ I_dd_23*(pp(mu)*Qm(nu)+pp(nu)*Qm(mu))+ I_dd_24*d_(mu,nu); id qpf(mu?)*qpf(nu?)*int= I_nn_21*pp(mu)*pp(nu)+I_nn_22*Qm(mu)*Qm(nu)+ I_nn_23*(pp(mu)*Qm(nu)+pp(nu)*Qm(mu))+ I_nn_24*d_(mu,nu); id kbf(mu?)*qpf(nu?)*int= I_dn_21*pp(mu)*pp(nu)+I_dn_22*Qm(mu)*Qm(nu)+ I_dn_23*pp(mu)*Qm(nu)+I_dn_24*pp(nu)*Qm(mu)+ I_dn_25*d_(mu,nu); id kbf(mu?)*int= I_d_11*pp(mu)+I_d_12*Qm(mu); id qpf(mu?)*int= I_n_11*pp(mu)+I_n_12*Qm(mu); id int= I0; * .sort * id I_d_11 = + I0*dms*Xp^-1*y^-1*s^-1 * ( 1 ) + I0*Xp^-2 * ( - 2*ups + 2*x1 - 2*x2 ) + I0*Xp^-1 * ( - 1 + ups - 2*x1 + 2*x2 ) + I0 * ( 1 ); al I_n_11 = + I0*Xp^-2 * ( 2*x2 - 2*z ) + I0*Xp^-1 * ( - 1 - 2*x2 + z ) + I0 * ( 1 ); al I_d_12 = + I0*Xp^-1 * ( ups - x1 + x2 ) + I0 * ( 1 ); al I_n_12 = + I0*Xp^-1 * ( - x2 + z ) + I0 * ( 1 ); * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * .sort * id I_dd_21 = + I0*dms*Xp^-3*y^-1*s^-1 * ( - 6*ups + 6*x1 - 6*x2 ) + I0*dms*Xp^-2*y^-1*s^-1 * ( - 2 + 4*ups - 6*x1 + 6*x2 ) + I0*dms*Xp^-1*y^-1*s^-1 * ( 2 ) + I0*Xp^-4 * ( - 12*ups*x1 + 12*ups*x2 + 6*ups^2 - 12*x1*x2 + 6*x1^ 2 + 6*x2^2 ) + I0*Xp^-3 * ( 18*ups*x1 - 18*ups*x2 + 6*ups - 6*ups^2 + 24*x1*x2 - 6* x1 - 12*x1^2 + 6*x2 - 12*x2^2 ) + I0*Xp^-2 * ( 1 - 6*ups*x1 + 6*ups*x2 - 8*ups + ups^2 - 12*x1*x2 + 12 *x1 + 6*x1^2 - 12*x2 + 6*x2^2 ) + I0*Xp^-1 * ( - 2 + 2*ups - 6*x1 + 6*x2 ) + I0 * ( 1 ); al I_nn_21 = + I0*Xp^-4 * ( - 12*x2*z + 6*x2^2 + 6*z^2 ) + I0*Xp^-3 * ( 18*x2*z - 6*x2 - 12*x2^2 + 6*z - 6*z^2 ) + I0*Xp^-2 * ( 1 - 6*x2*z + 12*x2 + 6*x2^2 - 8*z + z^2 ) + I0*Xp^-1 * ( - 2 - 6*x2 + 2*z ) + I0 * ( 1 ); al I_dn_21 = + I0*ums*Xp^-2*y^-1*s^-1 * ( 1 ) + I0*ums*Xp^-1*y^-1*s^-1 * ( - 1 ) + I0*dms*Xp^-3*y^-1*s^-1 * ( 3*x2 - 3*z ) + I0*dms*Xp^-2*y^-1*s^-1 * ( - 2 - 3*x2 + 2*z ) + I0*dms*Xp^-1*y^-1*s^-1 * ( 2 ) + I0*Xp^-4 * ( - 6*ups*x2 + 6*ups*z + 6*x1*x2 - 6*x1*z + 6*x2*z - 6 * x2^2 ) + I0*Xp^-3 * ( 9*ups*x2 - 6*ups*z + 3*ups - 12*x1*x2 + 9*x1*z - 3*x1 - 9*x2*z + 12*x2^2 + 3*z ) + I0*Xp^-2 * ( 2 - 3*ups*x2 + ups*z - 4*ups + 6*x1*x2 - 3*x1*z + 5*x 1 + 3*x2*z - 6*x2^2 - 4*z ) + I0*Xp^-1 * ( - 3 + ups - 2*x1 + z ) + I0 * ( 1 ); * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * .sort * id I_dd_22 = + I0*Xp^-2 * ( - 2*ups*x1 + 2*ups*x2 + ups^2 - 2*x1*x2 + x1^2 + x2^ 2 ) + I0*Xp^-1 * ( 2*ups - 2*x1 + 2*x2 ) + I0 * ( 1 ); al I_nn_22 = + I0*Xp^-2 * ( - 2*x2*z + x2^2 + z^2 ) + I0*Xp^-1 * ( - 2*x2 + 2*z ) + I0 * ( 1 ); al I_dn_22 = + I0*Xp^-2 * ( - ups*x2 + ups*z + x1*x2 - x1*z + x2*z - x2^2 ) + I0*Xp^-1 * ( ups - x1 + z ) + I0 * ( 1 ); * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * .sort * id I_dd_23 = + I0*dms*Xp^-2*y^-1*s^-1 * ( 2*ups - 2*x1 + 2*x2 ) + I0*dms*Xp^-1*y^-1*s^-1 * ( 1 ) + I0*Xp^-3 * ( 6*ups*x1 - 6*ups*x2 - 3*ups^2 + 6*x1*x2 - 3*x1^2 - 3* x2^2 ) + I0*Xp^-2 * ( - 4*ups*x1 + 4*ups*x2 - 4*ups + ups^2 - 6*x1*x2 + 4* x1 + 3*x1^2 - 4*x2 + 3*x2^2 ) + I0*Xp^-1 * ( - 1 + 2*ups - 4*x1 + 4*x2 ) + I0 * ( 1 ); al I_nn_23 = + I0*Xp^-3 * ( 6*x2*z - 3*x2^2 - 3*z^2 ) + I0*Xp^-2 * ( - 4*x2*z + 4*x2 + 3*x2^2 - 4*z + z^2 ) + I0*Xp^-1 * ( - 1 - 4*x2 + 2*z ) + I0 * ( 1 ); al I_dn_23 = + I0*ums*Xp^-1*y^-1*s^-1 * ( - 1/2 ) + I0*dms*Xp^-2*y^-1*s^-1 * ( - 3/2*x2 + 3/2*z ) + I0*dms*Xp^-1*y^-1*s^-1 * ( 3/2 ) + I0*Xp^-3 * ( 3*ups*x2 - 3*ups*z - 3*x1*x2 + 3*x1*z - 3*x2*z + 3*x2 ^2 ) + I0*Xp^-2 * ( - 3/2*ups*x2 + ups*z - 5/2*ups + 3*x1*x2 - 5/2*x1*z + 5/2*x1 + 5/2*x2*z - x2 - 3*x2^2 - 3/2*z ) + I0*Xp^-1 * ( - 3/2 + ups - 2*x1 + x2 + z ) + I0 * ( 1 ); * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * .sort * id I_dd_24 = + I0*dms*Xp^-1 * ( 1/2*ups - 1/2*x1 + 1/2*x2 ) + I0*Xp^-2*y*s * ( ups*x1 - ups*x2 - 1/2*ups^2 + x1*x2 - 1/2*x1^2 - 1/ 2*x2^2 ) + I0*Xp^-1*y*s * ( - 1/2*ups*x1 + 1/2*ups*x2 - 1/2*ups - x1*x2 + 1/ 2* x1 + 1/2*x1^2 - 1/2*x2 + 1/2*x2^2 ) + I0*y*s * ( - 1/2*x1 + 1/2*x2 ); id I_nn_24 = + I0*Xp^-2*y*s * ( x2*z - 1/2*x2^2 - 1/2*z^2 ) + I0*Xp^-1*y*s * ( - 1/2*x2*z + 1/2*x2 + 1/2*x2^2 - 1/2*z ) + I0*y*s * ( - 1/2*x2 ); id I_dn_24 = + I0*ums*Xp^-1*y^-1*s^-1 * ( - 1/2 ) + I0*dms*Xp^-2*y^-1*s^-1 * ( - 1/2*x2 + 1/2*z ) + I0*dms*Xp^-1*y^-1*s^-1 * ( 1/2 ) + I0*Xp^-3 * ( 3*ups*x2 - 3*ups*z - 3*x1*x2 + 3*x1*z - 3*x2*z + 3*x2 ^2 ) + I0*Xp^-2 * ( - 5/2*ups*x2 + ups*z - 3/2*ups + 3*x1*x2 - 3/2*x1*z + 3/2*x1 + 3/2*x2*z + x2 - 3*x2^2 - 5/2*z ) + I0*Xp^-1 * ( - 3/2 + ups - x1 - x2 + z ) + I0 * ( 1 ); * * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * .sort * id I_dn_25 = + I0*ums * ( - 1/4 ) + I0*dms*Xp^-1 * ( - 1/4*x2 + 1/4*z ) + I0*dms * ( 1/4 ) + I0*Xp^-2*y*s * ( 1/2*ups*x2 - 1/2*ups*z - 1/2*x1*x2 + 1/2*x1*z - 1 /2 *x2*z + 1/2*x2^2 ) + I0*Xp^-1*y*s * ( - 1/4*ups*x2 - 1/4*ups + 1/2*x1*x2 - 1/4*x1*z + 1/ 4*x1 + 1/4*x2*z - 1/2*x2^2 - 1/4*z ) + I0*y*s * ( - 1/4 ); * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id ums= 0; id dms*uM^-2= aux; id dms= 0; id aux= dms*uM^-2; * id uM^-1= -1/(y*s)/tdM; * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * ratio,z,zX,aux; id aux= 1+2*X; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * #endprocedure * #procedure totalgg(dummy) * * Now the full stuff * load wrl$grp04:[giampiero]outputampgg_B.sav; * g T_gg= PwsM2*PwtM2*wwww_gg/(eq*I0); *nwrite statistics; * id es(T)*es(T)= eq; * id wg(n?)= wg(n)*rs; al wgC(n?)= wgC(n)*rs; * id wg(4)= X*wg(4); al wg(7)= X*wg(7); al wg(8)= X*wg(8); al wg(10)= X*wg(10); id wgC(4)= X*wgC(4); al wgC(7)= X*wgC(7); al wgC(8)= X*wgC(8); al wgC(10)= X*wgC(10); * id wg(n?)*wgC(n?)= wgM2(n); * id GtM2*GsM2*wg(2)*wgC(1)= GtM2*GsM2*(-wgC(2)*wg(1)+ 2*Re(wg(2)*wgC(1))); al GtM2*GsM2*wg(3)*wgC(1)= GtM2*GsM2*(-wgC(3)*wg(1)+ 2*Re(wg(3)*wgC(1))); al GtM2*GsM2*wg(4)*wgC(1)= GtM2*GsM2*(-wgC(4)*wg(1)+ 2*Re(wg(4)*wgC(1))); al GtM2*GsM2*wg(5)*wgC(1)= GtM2*GsM2*(-wgC(5)*wg(1)+ 2*Re(wg(5)*wgC(1))); al GtM2*GsM2*wg(6)*wgC(1)= GtM2*GsM2*(-wgC(6)*wg(1)+ 2*Re(wg(6)*wgC(1))); al GtM2*GsM2*wg(7)*wgC(1)= GtM2*GsM2*(-wgC(7)*wg(1)+ 2*Re(wg(7)*wgC(1))); al GtM2*GsM2*wg(8)*wgC(1)= GtM2*GsM2*(-wgC(8)*wg(1)+ 2*Re(wg(8)*wgC(1))); al GtM2*GsM2*wg(9)*wgC(1)= GtM2*GsM2*(-wgC(9)*wg(1)+ 2*Re(wg(9)*wgC(1))); al GtM2*GsM2*wg(10)*wgC(1)= GtM2*GsM2*(-wgC(10)*wg(1)+ 2*Re(wg(10)*wgC(1))); * al GtM2*GsM2*wg(3)*wgC(2)= GtM2*GsM2*(-wgC(3)*wg(2)+ 2*Re(wg(3)*wgC(2))); al GtM2*GsM2*wg(4)*wgC(2)= GtM2*GsM2*(-wgC(4)*wg(2)+ 2*Re(wg(4)*wgC(2))); al GtM2*GsM2*wg(5)*wgC(2)= GtM2*GsM2*(-wgC(5)*wg(2)+ 2*Re(wg(5)*wgC(2))); al GtM2*GsM2*wg(6)*wgC(2)= GtM2*GsM2*(-wgC(6)*wg(2)+ 2*Re(wg(6)*wgC(2))); al GtM2*GsM2*wg(7)*wgC(2)= GtM2*GsM2*(-wgC(7)*wg(2)+ 2*Re(wg(7)*wgC(2))); al GtM2*GsM2*wg(8)*wgC(2)= GtM2*GsM2*(-wgC(8)*wg(2)+ 2*Re(wg(8)*wgC(2))); al GtM2*GsM2*wg(9)*wgC(2)= GtM2*GsM2*(-wgC(9)*wg(2)+ 2*Re(wg(9)*wgC(2))); al GtM2*GsM2*wg(10)*wgC(2)= GtM2*GsM2*(-wgC(10)*wg(2)+ 2*Re(wg(10)*wgC(2))); * al GtM2*GsM2*wg(4)*wgC(3)= GtM2*GsM2*(-wgC(4)*wg(3)+ 2*Re(wg(4)*wgC(3))); al GtM2*GsM2*wg(5)*wgC(3)= GtM2*GsM2*(-wgC(5)*wg(3)+ 2*Re(wg(5)*wgC(3))); al GtM2*GsM2*wg(6)*wgC(3)= GtM2*GsM2*(-wgC(6)*wg(3)+ 2*Re(wg(6)*wgC(3))); al GtM2*GsM2*wg(7)*wgC(3)= GtM2*GsM2*(-wgC(7)*wg(3)+ 2*Re(wg(7)*wgC(3))); al GtM2*GsM2*wg(8)*wgC(3)= GtM2*GsM2*(-wgC(8)*wg(3)+ 2*Re(wg(8)*wgC(3))); al GtM2*GsM2*wg(9)*wgC(3)= GtM2*GsM2*(-wgC(9)*wg(3)+ 2*Re(wg(9)*wgC(3))); al GtM2*GsM2*wg(10)*wgC(3)= GtM2*GsM2*(-wgC(10)*wg(3)+ 2*Re(wg(10)*wgC(3))); * al GtM2*GsM2*wg(5)*wgC(4)= GtM2*GsM2*(-wgC(5)*wg(4)+ 2*Re(wg(5)*wgC(4))); al GtM2*GsM2*wg(6)*wgC(4)= GtM2*GsM2*(-wgC(6)*wg(4)+ 2*Re(wg(6)*wgC(4))); al GtM2*GsM2*wg(7)*wgC(4)= GtM2*GsM2*(-wgC(7)*wg(4)+ 2*Re(wg(7)*wgC(4))); al GtM2*GsM2*wg(8)*wgC(4)= GtM2*GsM2*(-wgC(8)*wg(4)+ 2*Re(wg(8)*wgC(4))); al GtM2*GsM2*wg(9)*wgC(4)= GtM2*GsM2*(-wgC(9)*wg(4)+ 2*Re(wg(9)*wgC(4))); al GtM2*GsM2*wg(10)*wgC(4)= GtM2*GsM2*(-wgC(10)*wg(4)+ 2*Re(wg(10)*wgC(4))); * al GtM2*GsM2*wg(6)*wgC(5)= GtM2*GsM2*(-wgC(6)*wg(5)+ 2*Re(wg(6)*wgC(5))); al GtM2*GsM2*wg(7)*wgC(5)= GtM2*GsM2*(-wgC(7)*wg(5)+ 2*Re(wg(7)*wgC(5))); al GtM2*GsM2*wg(8)*wgC(5)= GtM2*GsM2*(-wgC(8)*wg(5)+ 2*Re(wg(8)*wgC(5))); al GtM2*GsM2*wg(9)*wgC(5)= GtM2*GsM2*(-wgC(9)*wg(5)+ 2*Re(wg(9)*wgC(5))); al GtM2*GsM2*wg(10)*wgC(5)= GtM2*GsM2*(-wgC(10)*wg(5)+ 2*Re(wg(10)*wgC(5))); * al GtM2*GsM2*wg(7)*wgC(6)= GtM2*GsM2*(-wgC(7)*wg(6)+ 2*Re(wg(7)*wgC(6))); al GtM2*GsM2*wg(8)*wgC(6)= GtM2*GsM2*(-wgC(8)*wg(6)+ 2*Re(wg(8)*wgC(6))); al GtM2*GsM2*wg(9)*wgC(6)= GtM2*GsM2*(-wgC(9)*wg(6)+ 2*Re(wg(9)*wgC(6))); al GtM2*GsM2*wg(10)*wgC(6)= GtM2*GsM2*(-wgC(10)*wg(6)+ 2*Re(wg(10)*wgC(6))); * al GtM2*GsM2*wg(8)*wgC(7)= GtM2*GsM2*(-wgC(8)*wg(7)+ 2*Re(wg(8)*wgC(7))); al GtM2*GsM2*wg(9)*wgC(7)= GtM2*GsM2*(-wgC(9)*wg(7)+ 2*Re(wg(9)*wgC(7))); al GtM2*GsM2*wg(10)*wgC(7)= GtM2*GsM2*(-wgC(10)*wg(7)+ 2*Re(wg(10)*wgC(7))); * al GtM2*GsM2*wg(9)*wgC(8)= GtM2*GsM2*(-wgC(9)*wg(8)+ 2*Re(wg(9)*wgC(8))); al GtM2*GsM2*wg(10)*wgC(8)= GtM2*GsM2*(-wgC(10)*wg(8)+ 2*Re(wg(10)*wgC(8))); * al GtM2*GsM2*wg(10)*wgC(9)= GtM2*GsM2*(-wgC(10)*wg(9)+ 2*Re(wg(10)*wgC(9))); * id wg(n?)*Gs*GtM2*PwsC= GtM2*(-wgC(n)*GsC*Pws+2*Re(wg(n)*Gs*PwsC)); id wg(n?)*Gt*GsM2*PwtC= GsM2*(-wgC(n)*GtC*Pwt+2*Re(wg(n)*Gt*PwtC)); id wg(n?)*Gs*GtM2= GtM2*(-GsC*wgC(n)+2*Re(Gs*wg(n))); id wg(n?)*Gt*GsM2= GsM2*(-GtC*wgC(n)+2*Re(Gt*wg(n))); * id GsM2*GtC*PwtC*wgC(n?)= GsM2*(-Gt*Pwt*wg(n)+2*Re(Gt*Pwt*wg(n))); id GtM2*GsC*PwsC*wgC(n?)= GtM2*(-Gs*Pws*wg(n)+2*Re(Gs*Pws*wg(n))); * id GsM2*GtC*PwtC= GsM2*(-Gt*Pwt+2*Re(Gt*Pwt)); id GtM2*GsC*PwsC= GtM2*(-Gs*Pws+2*Re(Gs*Pws)); * id GsC*PwsC*Gt*Pwt= -Gs*Pws*GtC*PwtC+2*Re(GsC*Gt*PwsC*Pwt); *id GsC*Gt*Pws*PwtC= -Gs*GtC*PwsC*Pwt+2*Re(Gs*GtC*PwsC*Pwt); id GsC*Gt*Pws*PwtC= Re(Gs*GtC*PwsC*Pwt)+i_*Im(Gs*GtC*PwsC*Pwt); id Gs*GtC*PwsC*Pwt= Re(Gs*GtC*PwsC*Pwt)-i_*Im(Gs*GtC*PwsC*Pwt); id GsC*Gt*PwsC*PwtC= -Gs*GtC*Pws*Pwt+2*Re(Gs*GtC*Pws*Pwt); * id Gs*GtC*Pwt= -GsC*Gt*PwtC+2*Re(Gs*GtC*Pwt); * id GtM2*PwsC= GtM2*(-Pws+2*Re(Pws)); id GsM2*PwtC= GsM2*(-Pwt+2*Re(Pwt)); * id Re(wgC(2)*wg(1))= Re(wg(2)*wgC(1)); al Re(wgC(3)*wg(1))= Re(wg(3)*wgC(1)); al Re(wgC(4)*wg(1))= Re(wg(4)*wgC(1)); al Re(wgC(5)*wg(1))= Re(wg(5)*wgC(1)); al Re(wgC(6)*wg(1))= Re(wg(6)*wgC(1)); al Re(wgC(7)*wg(1))= Re(wg(7)*wgC(1)); al Re(wgC(8)*wg(1))= Re(wg(8)*wgC(1)); al Re(wgC(9)*wg(1))= Re(wg(9)*wgC(1)); al Re(wgC(10)*wg(1))= Re(wg(10)*wgC(1)); * id Re(wgC(3)*wg(2))= Re(wg(3)*wgC(2)); al Re(wgC(4)*wg(2))= Re(wg(4)*wgC(2)); al Re(wgC(5)*wg(2))= Re(wg(5)*wgC(2)); al Re(wgC(6)*wg(2))= Re(wg(6)*wgC(2)); al Re(wgC(7)*wg(2))= Re(wg(7)*wgC(2)); al Re(wgC(8)*wg(2))= Re(wg(8)*wgC(2)); al Re(wgC(9)*wg(2))= Re(wg(9)*wgC(2)); al Re(wgC(10)*wg(2))= Re(wg(10)*wgC(2)); * id Re(wgC(4)*wg(3))= Re(wg(4)*wgC(3)); al Re(wgC(5)*wg(3))= Re(wg(5)*wgC(3)); al Re(wgC(6)*wg(3))= Re(wg(6)*wgC(3)); al Re(wgC(7)*wg(3))= Re(wg(7)*wgC(3)); al Re(wgC(8)*wg(3))= Re(wg(8)*wgC(3)); al Re(wgC(9)*wg(3))= Re(wg(9)*wgC(3)); al Re(wgC(10)*wg(3))= Re(wg(10)*wgC(3)); * id Re(wgC(5)*wg(4))= Re(wg(5)*wgC(4)); al Re(wgC(6)*wg(4))= Re(wg(6)*wgC(4)); al Re(wgC(7)*wg(4))= Re(wg(7)*wgC(4)); al Re(wgC(8)*wg(4))= Re(wg(8)*wgC(4)); al Re(wgC(9)*wg(4))= Re(wg(9)*wgC(4)); al Re(wgC(10)*wg(4))= Re(wg(10)*wgC(4)); * id Re(wgC(6)*wg(5))= Re(wg(6)*wgC(5)); al Re(wgC(7)*wg(5))= Re(wg(7)*wgC(5)); al Re(wgC(8)*wg(5))= Re(wg(8)*wgC(5)); al Re(wgC(9)*wg(5))= Re(wg(9)*wgC(5)); al Re(wgC(10)*wg(5))= Re(wg(10)*wgC(5)); * id Re(wgC(7)*wg(6))= Re(wg(7)*wgC(6)); al Re(wgC(8)*wg(6))= Re(wg(8)*wgC(6)); al Re(wgC(9)*wg(6))= Re(wg(9)*wgC(6)); al Re(wgC(10)*wg(6))= Re(wg(10)*wgC(6)); * id Re(wgC(8)*wg(7))= Re(wg(8)*wgC(7)); al Re(wgC(9)*wg(7))= Re(wg(9)*wgC(7)); al Re(wgC(10)*wg(7))= Re(wg(10)*wgC(7)); * id Re(wgC(9)*wg(8))= Re(wg(9)*wgC(8)); al Re(wgC(10)*wg(8))= Re(wg(10)*wgC(8)); * id Re(wgC(10)*wg(9))= Re(wg(10)*wgC(9)); * .store * #endprocedure * #procedure contractgg(dummy) * g Wd_gg= T_gg*d_(mu,nu); g Wpp_gg= T_gg*pp(mu)*pp(nu); * id pp.pp= 0; al qp.qp= 0; al Qm.Qm= QmS; al kb.kb= -dms; id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; * id pp.Qm= -1/2*(QmS+y*s); al k.kb= 1/2*(-sp+ums+dms); al Qm.k= 1/2*(QmS+t-ums); al pp.kb= 1/2*(tp-dms); al Qm.kb= 1/2*(QmS+u-dms); al pp.k= 1/2*(up-ums); al pp.qp= 1/2*kp; al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(dms-zp); al k.qp= 1/2*(ums-zm); * id sp= x2*y*s; al km= z*y*s; al km^-1= 1/(z*y*s); al zm= -(x2-x1)*y*s; id QmS= X*y*s; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id kb(mu?)= kbf(mu); id qp(mu?)= qpf(mu); id e_(kb,qp,nu?,al?)= kbf(ab1)*qpf(ab2)*e_(ab1,ab2,nu,al); al e_(kb,mu?,nu?,al?)= kbf(ab1)*e_(ab1,mu,nu,al); al e_(qp,mu?,nu?,al?)= qpf(ab2)*e_(ab2,mu,nu,al); * #call integration{dummy} * id rs^-2= -1/(X*y*s); id pp.pp= 0; al qp.qp= 0; al Qm.Qm= QmS; id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; * id pp.Qm= -1/2*(QmS+y*s); al k.kb= 1/2*(-sp+ums+dms); al Qm.k= 1/2*(QmS+t-ums); al pp.kb= 1/2*(tp-dms); al Qm.kb= 1/2*(QmS+u-dms); al pp.k= 1/2*(up-ums); al pp.qp= 1/2*kp; al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(dms-zp); al k.qp= 1/2*(ums-zm); * id sp= x2*y*s; al km= z*y*s; al km^-1= 1/(z*y*s); al zm= -(x2-x1)*y*s; id QmS= X*y*s; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * .store * g GG10= (-1/2*Wd_gg-2*X/Xp^2*Wpp_gg/(y*s))*sel0*s; g GG20= (-1/2*Wd_gg-6*X/Xp^2*Wpp_gg/(y*s))*sel0*s; g GG11= (-1/2*Wd_gg-2*X/Xp^2*Wpp_gg/(y*s))*seld*s; g GG21= (-1/2*Wd_gg-6*X/Xp^2*Wpp_gg/(y*s))*seld*s; g GG12= (-1/2*Wd_gg-2*X/Xp^2*Wpp_gg/(y*s))*seldd*s; g GG22= (-1/2*Wd_gg-6*X/Xp^2*Wpp_gg/(y*s))*seldd*s; * ratio,tdM,ups,aux; id aux= dms/(y*s); * id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * id tdM^-2*seldd= 1; id seldd= 0; * id tdM^-2*seld= 0; id tdM^-1*seld= 1; id seld= 0; * id tdM^-2*sel0= 0; id tdM^-1*sel0= 0; * id sel0= 1; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * *-- id Re(wg(n?)*wgC(2))= 1/2*(Re(wg(n)*wgC([2+3]))+Re(wg(n)*wgC([2-3]))); id Re(wg(n?)*wgC(3))= 1/2*(Re(wg(n)*wgC([2+3]))-Re(wg(n)*wgC([2-3]))); id Re(wg(2)*wgC(2))= wgM2(2); id Re(wg(3)*wgC(3))= wgM2(3); * id Re(wg(2)*wgC(n?))= 1/2*(Re(wg([2+3])*wgC(n))+Re(wg([2-3])*wgC(n))); id Re(wg(3)*wgC(n?))= 1/2*(Re(wg([2+3])*wgC(n))-Re(wg([2-3])*wgC(n))); id Re(wg(2)*wgC(2))= wgM2(2); id Re(wg(3)*wgC(3))= wgM2(3); * id Re(wg(n?)*wgC(5))= 1/2*(Re(wg(n)*wgC([5+6]))+Re(wg(n)*wgC([5-6]))); id Re(wg(n?)*wgC(6))= 1/2*(Re(wg(n)*wgC([5+6]))-Re(wg(n)*wgC([5-6]))); id Re(wg(5)*wgC(5))= wgM2(5); id Re(wg(6)*wgC(6))= wgM2(6); * id Re(wg(5)*wgC(n?))= 1/2*(Re(wg([5+6])*wgC(n))+Re(wg([5-6])*wgC(n))); id Re(wg(6)*wgC(n?))= 1/2*(Re(wg([5+6])*wgC(n))-Re(wg([5-6])*wgC(n))); id Re(wg(5)*wgC(5))= wgM2(5); id Re(wg(6)*wgC(6))= wgM2(6); * id Re(wg(2)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)+Re(wg([2-3])*X*Y)); id Re(wg(3)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)-Re(wg([2-3])*X*Y)); id Re(wg(5)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)+Re(wg([5-6])*X*Y)); id Re(wg(6)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)-Re(wg([5-6])*X*Y)); * id Re(wg(2)*X?)= 1/2*(Re(wg([2+3])*X)+Re(wg([2-3])*X)); id Re(wg(3)*X?)= 1/2*(Re(wg([2+3])*X)-Re(wg([2-3])*X)); id Re(wg(5)*X?)= 1/2*(Re(wg([5+6])*X)+Re(wg([5-6])*X)); id Re(wg(6)*X?)= 1/2*(Re(wg([5+6])*X)-Re(wg([5-6])*X)); * id Re(wg([2+3])*wgC([2-3]))= Re(wg([2-3])*wgC([2+3])); id Re(wg([2+3])*wgC([5-6]))= Re(wg([5-6])*wgC([2+3])); id Re(wg([2-3])*wgC([5+6]))= Re(wg([5+6])*wgC([2-3])); id Re(wg([5+6])*wgC([5-6]))= Re(wg([5-6])*wgC([5+6])); * id Re(wg([2+3])*wgC([2+3]))= wgM2(2)+wgM2(3)+2*Re(wg(2)*wgC(3)); id Re(wg([2-3])*wgC([2-3]))= wgM2(2)+wgM2(3)-2*Re(wg(2)*wgC(3)); id Re(wg([5+6])*wgC([5+6]))= wgM2(5)+wgM2(6)+2*Re(wg(5)*wgC(6)); id Re(wg([5-6])*wgC([5-6]))= wgM2(5)+wgM2(6)-2*Re(wg(5)*wgC(6)); *-- id X= [1+X]-1; id x1= 1-[1-x1]; id x2= 1-[1-x1]-[x1-x2]; id z= 1-[1-x1]-[x1-x2]-[x2-z]; id ups= [ups+X+x1]-[1+X]+[1-x1]; id Xp^-1= 1/[1+X]; ratio,X,[1+X],one; id one= 1; al one^-1= 1; * id Re(wg(X?)*wgC(Y?))= Rep(X,Y); id GsM2*GtM2= GstM2; id [1+X] = Xp; al [1+X]^-1= Xpi; id [1-x1]= omx1; al [x1-x2]= x12; al [x2-z]= x2z; al [ups+X+x1]= ux; al zX^-1= zXi; id Xpi^2= Xpi-aux; id aux= XX; * id Xp^2= Xps; id Xpi^2= Xpis; id y^3= yc; id y^2= ys; * id omx1= O; id O^4= Oq; id O^3= Oc; id O^2= Os; id x12^3= x12c; id x12^2= x12s; id x2z^4= x2zq; id x2z^3= x2zc; id x2z^2= x2zs; id ux^2= uxs; id vv^2= vvs; id vv^3= vvc; id vvs*vv= vvc; id vvs^2= vvq; id vvc*vv= vvq; * id Re(Gs*GtC*Pwt*PwsC)= rl1; id Re(Gs*GtC*Pwt)= rl2; id Re(Pws)= rl3; id Re(wg(n?)*Gs*PwsC)= rs1wg(n); id Re(wg(n?)*Gs)= rs2wg(n); id Re(wg(n?)*Gt*PwtC)= rt1wg(n); id Re(wg(n?)*Gt)= rt2wg(n); id Rep(n1?,n2?)= rwwg(n1,n2); id Re(wg(n1?)*wgC(n2?))= rwwg(n1,n2); id wgM2(n?)= rwwg(n,n); * id fun?([2+3],n?)= fun(11,n); al fun?([2-3],n?)= fun(12,n); al fun?([5+6],n?)= fun(13,n); al fun?([5-6],n?)= fun(14,n); * id fun?(n?,[2+3])= fun(n,11); al fun?(n?,[2-3])= fun(n,12); al fun?(n?,[5+6])= fun(n,13); al fun?(n?,[5-6])= fun(n,14); * id fun?([2+3])= fun(11); al fun?([2-3])= fun(12); al fun?([5+6])= fun(13); al fun?([5-6])= fun(14); * * .store * *----g^2-------------------- * g GG10s= GG10; g GG20s= GG20; g GG11s= GG11; g GG21s= GG21; g GG12s= GG12; g GG22s= GG22; * id gexp^4= 0; id gexp^3= 0; id gexp^2= 1; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wg,wgC,wgM2,rs1wg,rs2wg, rt1wg,rt2wg,rwwg,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, * ux,uxs,x2z,x2zs,x2zc,x2zq,x12,x12s,x12c,zXi,XX; zXi,XX,vv,vvs,vvc,vvq; format fortran; print; * .store * *----g^3-------------------- * g GG10c= GG10; g GG20c= GG20; g GG11c= GG11; g GG21c= GG21; g GG12c= GG12; g GG22c= GG22; * id gexp^4= 0; id gexp^3= 1; id gexp^2= 0; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wg,wgC,wgM2,rs1wg,rs2wg, rt1wg,rt2wg,rwwg,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, * ux,uxs,x2z,x2zs,x2zc,x2zq,x12,x12s,x12c,zXi,XX; zXi,XX,vv,vvs,vvc,vvq; format fortran; print; * .store * *----g^4-------------------- * g GG10q= GG10/GstM2; g GG20q= GG20/GstM2; g GG11q= GG11/GstM2; g GG21q= GG21/GstM2; g GG12q= GG12/GstM2; g GG22q= GG22/GstM2; * id gexp^4= 1; id gexp^3= 0; id gexp^2= 0; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wg,wgC,wgM2,rs1wg,rs2wg, rt1wg,rt2wg,rwwg,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, * ux,uxs,x2z,x2zs,x2zc,x2zq,x12,x12s,x12c,zXi,XX; zXi,XX,vv,vvs,vvc,vvq; format fortran; print; * #endprocedure * #procedure totalgz(dummy) * * Now the full stuff * load wrl$grp04:[giampiero]outputampgz_B.sav; * g T_gz= PwsM2*PwtM2*wwww_gz/(egt*I0); *nwrite statistics; * id es(T)*GT= egt; * id wg(n?)= wg(n)*rs; al wgC(n?)= wgC(n)*rs; id wz(n?)= wz(n)*rs; al wzC(n?)= wzC(n)*rs; * id wg(4)= X*wg(4); al wg(7)= X*wg(7); al wg(8)= X*wg(8); al wg(10)= X*wg(10); id wgC(4)= X*wgC(4); al wgC(7)= X*wgC(7); al wgC(8)= X*wgC(8); al wgC(10)= X*wgC(10); * id wz(4)= X*wz(4); al wz(7)= X*wz(7); al wz(8)= X*wz(8); al wz(10)= X*wz(10); id wzC(4)= X*wzC(4); al wzC(7)= X*wzC(7); al wzC(8)= X*wzC(8); al wzC(10)= X*wzC(10); * id wg(n?)*wgC(n?)= wgM2(n); id wz(n?)*wzC(n?)= wzM2(n); * id GtM2*GsM2*wg(2)*wzC(1)= GtM2*GsM2*(-wgC(2)*wz(1)+ 2*Re(wg(2)*wzC(1))); al GtM2*GsM2*wg(3)*wzC(1)= GtM2*GsM2*(-wgC(3)*wz(1)+ 2*Re(wg(3)*wzC(1))); al GtM2*GsM2*wg(4)*wzC(1)= GtM2*GsM2*(-wgC(4)*wz(1)+ 2*Re(wg(4)*wzC(1))); al GtM2*GsM2*wg(5)*wzC(1)= GtM2*GsM2*(-wgC(5)*wz(1)+ 2*Re(wg(5)*wzC(1))); al GtM2*GsM2*wg(6)*wzC(1)= GtM2*GsM2*(-wgC(6)*wz(1)+ 2*Re(wg(6)*wzC(1))); al GtM2*GsM2*wg(7)*wzC(1)= GtM2*GsM2*(-wgC(7)*wz(1)+ 2*Re(wg(7)*wzC(1))); al GtM2*GsM2*wg(8)*wzC(1)= GtM2*GsM2*(-wgC(8)*wz(1)+ 2*Re(wg(8)*wzC(1))); al GtM2*GsM2*wg(9)*wzC(1)= GtM2*GsM2*(-wgC(9)*wz(1)+ 2*Re(wg(9)*wzC(1))); al GtM2*GsM2*wg(10)*wzC(1)= GtM2*GsM2*(-wgC(10)*wz(1)+ 2*Re(wg(10)*wzC(1))); * al GtM2*GsM2*wg(3)*wzC(2)= GtM2*GsM2*(-wgC(3)*wz(2)+ 2*Re(wg(3)*wzC(2))); al GtM2*GsM2*wg(4)*wzC(2)= GtM2*GsM2*(-wgC(4)*wz(2)+ 2*Re(wg(4)*wzC(2))); al GtM2*GsM2*wg(5)*wzC(2)= GtM2*GsM2*(-wgC(5)*wz(2)+ 2*Re(wg(5)*wzC(2))); al GtM2*GsM2*wg(6)*wzC(2)= GtM2*GsM2*(-wgC(6)*wz(2)+ 2*Re(wg(6)*wzC(2))); al GtM2*GsM2*wg(7)*wzC(2)= GtM2*GsM2*(-wgC(7)*wz(2)+ 2*Re(wg(7)*wzC(2))); al GtM2*GsM2*wg(8)*wzC(2)= GtM2*GsM2*(-wgC(8)*wz(2)+ 2*Re(wg(8)*wzC(2))); al GtM2*GsM2*wg(9)*wzC(2)= GtM2*GsM2*(-wgC(9)*wz(2)+ 2*Re(wg(9)*wzC(2))); al GtM2*GsM2*wg(10)*wzC(2)= GtM2*GsM2*(-wgC(10)*wz(2)+ 2*Re(wg(10)*wzC(2))); * al GtM2*GsM2*wg(4)*wzC(3)= GtM2*GsM2*(-wgC(4)*wz(3)+ 2*Re(wg(4)*wzC(3))); al GtM2*GsM2*wg(5)*wzC(3)= GtM2*GsM2*(-wgC(5)*wz(3)+ 2*Re(wg(5)*wzC(3))); al GtM2*GsM2*wg(6)*wzC(3)= GtM2*GsM2*(-wgC(6)*wz(3)+ 2*Re(wg(6)*wzC(3))); al GtM2*GsM2*wg(7)*wzC(3)= GtM2*GsM2*(-wgC(7)*wz(3)+ 2*Re(wg(7)*wzC(3))); al GtM2*GsM2*wg(8)*wzC(3)= GtM2*GsM2*(-wgC(8)*wz(3)+ 2*Re(wg(8)*wzC(3))); al GtM2*GsM2*wg(9)*wzC(3)= GtM2*GsM2*(-wgC(9)*wz(3)+ 2*Re(wg(9)*wzC(3))); al GtM2*GsM2*wg(10)*wzC(3)= GtM2*GsM2*(-wgC(10)*wz(3)+ 2*Re(wg(10)*wzC(3))); * al GtM2*GsM2*wg(5)*wzC(4)= GtM2*GsM2*(-wgC(5)*wz(4)+ 2*Re(wg(5)*wzC(4))); al GtM2*GsM2*wg(6)*wzC(4)= GtM2*GsM2*(-wgC(6)*wz(4)+ 2*Re(wg(6)*wzC(4))); al GtM2*GsM2*wg(7)*wzC(4)= GtM2*GsM2*(-wgC(7)*wz(4)+ 2*Re(wg(7)*wzC(4))); al GtM2*GsM2*wg(8)*wzC(4)= GtM2*GsM2*(-wgC(8)*wz(4)+ 2*Re(wg(8)*wzC(4))); al GtM2*GsM2*wg(9)*wzC(4)= GtM2*GsM2*(-wgC(9)*wz(4)+ 2*Re(wg(9)*wzC(4))); al GtM2*GsM2*wg(10)*wzC(4)= GtM2*GsM2*(-wgC(10)*wz(4)+ 2*Re(wg(10)*wzC(4))); * al GtM2*GsM2*wg(6)*wzC(5)= GtM2*GsM2*(-wgC(6)*wz(5)+ 2*Re(wg(6)*wzC(5))); al GtM2*GsM2*wg(7)*wzC(5)= GtM2*GsM2*(-wgC(7)*wz(5)+ 2*Re(wg(7)*wzC(5))); al GtM2*GsM2*wg(8)*wzC(5)= GtM2*GsM2*(-wgC(8)*wz(5)+ 2*Re(wg(8)*wzC(5))); al GtM2*GsM2*wg(9)*wzC(5)= GtM2*GsM2*(-wgC(9)*wz(5)+ 2*Re(wg(9)*wzC(5))); al GtM2*GsM2*wg(10)*wzC(5)= GtM2*GsM2*(-wgC(10)*wz(5)+ 2*Re(wg(10)*wzC(5))); * al GtM2*GsM2*wg(7)*wzC(6)= GtM2*GsM2*(-wgC(7)*wz(6)+ 2*Re(wg(7)*wzC(6))); al GtM2*GsM2*wg(8)*wzC(6)= GtM2*GsM2*(-wgC(8)*wz(6)+ 2*Re(wg(8)*wzC(6))); al GtM2*GsM2*wg(9)*wzC(6)= GtM2*GsM2*(-wgC(9)*wz(6)+ 2*Re(wg(9)*wzC(6))); al GtM2*GsM2*wg(10)*wzC(6)= GtM2*GsM2*(-wgC(10)*wz(6)+ 2*Re(wg(10)*wzC(6))); * al GtM2*GsM2*wg(8)*wzC(7)= GtM2*GsM2*(-wgC(8)*wz(7)+ 2*Re(wg(8)*wzC(7))); al GtM2*GsM2*wg(9)*wzC(7)= GtM2*GsM2*(-wgC(9)*wz(7)+ 2*Re(wg(9)*wzC(7))); al GtM2*GsM2*wg(10)*wzC(7)= GtM2*GsM2*(-wgC(10)*wz(7)+ 2*Re(wg(10)*wzC(7))); * al GtM2*GsM2*wg(9)*wzC(8)= GtM2*GsM2*(-wgC(9)*wz(8)+ 2*Re(wg(9)*wzC(8))); al GtM2*GsM2*wg(10)*wzC(8)= GtM2*GsM2*(-wgC(10)*wz(8)+ 2*Re(wg(10)*wzC(8))); * al GtM2*GsM2*wg(10)*wzC(9)= GtM2*GsM2*(-wgC(10)*wz(9)+ 2*Re(wg(10)*wzC(9))); * id GtM2*GsM2*wz(2)*wgC(1)= GtM2*GsM2*(-wzC(2)*wg(1)+ 2*Re(wz(2)*wgC(1))); al GtM2*GsM2*wz(3)*wgC(1)= GtM2*GsM2*(-wzC(3)*wg(1)+ 2*Re(wz(3)*wgC(1))); al GtM2*GsM2*wz(4)*wgC(1)= GtM2*GsM2*(-wzC(4)*wg(1)+ 2*Re(wz(4)*wgC(1))); al GtM2*GsM2*wz(5)*wgC(1)= GtM2*GsM2*(-wzC(5)*wg(1)+ 2*Re(wz(5)*wgC(1))); al GtM2*GsM2*wz(6)*wgC(1)= GtM2*GsM2*(-wzC(6)*wg(1)+ 2*Re(wz(6)*wgC(1))); al GtM2*GsM2*wz(7)*wgC(1)= GtM2*GsM2*(-wzC(7)*wg(1)+ 2*Re(wz(7)*wgC(1))); al GtM2*GsM2*wz(8)*wgC(1)= GtM2*GsM2*(-wzC(8)*wg(1)+ 2*Re(wz(8)*wgC(1))); al GtM2*GsM2*wz(9)*wgC(1)= GtM2*GsM2*(-wzC(9)*wg(1)+ 2*Re(wz(9)*wgC(1))); al GtM2*GsM2*wz(10)*wgC(1)= GtM2*GsM2*(-wzC(10)*wg(1)+ 2*Re(wz(10)*wgC(1))); * al GtM2*GsM2*wz(3)*wgC(2)= GtM2*GsM2*(-wzC(3)*wg(2)+ 2*Re(wz(3)*wgC(2))); al GtM2*GsM2*wz(4)*wgC(2)= GtM2*GsM2*(-wzC(4)*wg(2)+ 2*Re(wz(4)*wgC(2))); al GtM2*GsM2*wz(5)*wgC(2)= GtM2*GsM2*(-wzC(5)*wg(2)+ 2*Re(wz(5)*wgC(2))); al GtM2*GsM2*wz(6)*wgC(2)= GtM2*GsM2*(-wzC(6)*wg(2)+ 2*Re(wz(6)*wgC(2))); al GtM2*GsM2*wz(7)*wgC(2)= GtM2*GsM2*(-wzC(7)*wg(2)+ 2*Re(wz(7)*wgC(2))); al GtM2*GsM2*wz(8)*wgC(2)= GtM2*GsM2*(-wzC(8)*wg(2)+ 2*Re(wz(8)*wgC(2))); al GtM2*GsM2*wz(9)*wgC(2)= GtM2*GsM2*(-wzC(9)*wg(2)+ 2*Re(wz(9)*wgC(2))); al GtM2*GsM2*wz(10)*wgC(2)= GtM2*GsM2*(-wzC(10)*wg(2)+ 2*Re(wz(10)*wgC(2))); * al GtM2*GsM2*wz(4)*wgC(3)= GtM2*GsM2*(-wzC(4)*wg(3)+ 2*Re(wz(4)*wgC(3))); al GtM2*GsM2*wz(5)*wgC(3)= GtM2*GsM2*(-wzC(5)*wg(3)+ 2*Re(wz(5)*wgC(3))); al GtM2*GsM2*wz(6)*wgC(3)= GtM2*GsM2*(-wzC(6)*wg(3)+ 2*Re(wz(6)*wgC(3))); al GtM2*GsM2*wz(7)*wgC(3)= GtM2*GsM2*(-wzC(7)*wg(3)+ 2*Re(wz(7)*wgC(3))); al GtM2*GsM2*wz(8)*wgC(3)= GtM2*GsM2*(-wzC(8)*wg(3)+ 2*Re(wz(8)*wgC(3))); al GtM2*GsM2*wz(9)*wgC(3)= GtM2*GsM2*(-wzC(9)*wg(3)+ 2*Re(wz(9)*wgC(3))); al GtM2*GsM2*wz(10)*wgC(3)= GtM2*GsM2*(-wzC(10)*wg(3)+ 2*Re(wz(10)*wgC(3))); * al GtM2*GsM2*wz(5)*wgC(4)= GtM2*GsM2*(-wzC(5)*wg(4)+ 2*Re(wz(5)*wgC(4))); al GtM2*GsM2*wz(6)*wgC(4)= GtM2*GsM2*(-wzC(6)*wg(4)+ 2*Re(wz(6)*wgC(4))); al GtM2*GsM2*wz(7)*wgC(4)= GtM2*GsM2*(-wzC(7)*wg(4)+ 2*Re(wz(7)*wgC(4))); al GtM2*GsM2*wz(8)*wgC(4)= GtM2*GsM2*(-wzC(8)*wg(4)+ 2*Re(wz(8)*wgC(4))); al GtM2*GsM2*wz(9)*wgC(4)= GtM2*GsM2*(-wzC(9)*wg(4)+ 2*Re(wz(9)*wgC(4))); al GtM2*GsM2*wz(10)*wgC(4)= GtM2*GsM2*(-wzC(10)*wg(4)+ 2*Re(wz(10)*wgC(4))); * al GtM2*GsM2*wz(6)*wgC(5)= GtM2*GsM2*(-wzC(6)*wg(5)+ 2*Re(wz(6)*wgC(5))); al GtM2*GsM2*wz(7)*wgC(5)= GtM2*GsM2*(-wzC(7)*wg(5)+ 2*Re(wz(7)*wgC(5))); al GtM2*GsM2*wz(8)*wgC(5)= GtM2*GsM2*(-wzC(8)*wg(5)+ 2*Re(wz(8)*wgC(5))); al GtM2*GsM2*wz(9)*wgC(5)= GtM2*GsM2*(-wzC(9)*wg(5)+ 2*Re(wz(9)*wgC(5))); al GtM2*GsM2*wz(10)*wgC(5)= GtM2*GsM2*(-wzC(10)*wg(5)+ 2*Re(wz(10)*wgC(5))); * al GtM2*GsM2*wz(7)*wgC(6)= GtM2*GsM2*(-wzC(7)*wg(6)+ 2*Re(wz(7)*wgC(6))); al GtM2*GsM2*wz(8)*wgC(6)= GtM2*GsM2*(-wzC(8)*wg(6)+ 2*Re(wz(8)*wgC(6))); al GtM2*GsM2*wz(9)*wgC(6)= GtM2*GsM2*(-wzC(9)*wg(6)+ 2*Re(wz(9)*wgC(6))); al GtM2*GsM2*wz(10)*wgC(6)= GtM2*GsM2*(-wzC(10)*wg(6)+ 2*Re(wz(10)*wgC(6))); * al GtM2*GsM2*wz(8)*wgC(7)= GtM2*GsM2*(-wzC(8)*wg(7)+ 2*Re(wz(8)*wgC(7))); al GtM2*GsM2*wz(9)*wgC(7)= GtM2*GsM2*(-wzC(9)*wg(7)+ 2*Re(wz(9)*wgC(7))); al GtM2*GsM2*wz(10)*wgC(7)= GtM2*GsM2*(-wzC(10)*wg(7)+ 2*Re(wz(10)*wgC(7))); * al GtM2*GsM2*wz(9)*wgC(8)= GtM2*GsM2*(-wzC(9)*wg(8)+ 2*Re(wz(9)*wgC(8))); al GtM2*GsM2*wz(10)*wgC(8)= GtM2*GsM2*(-wzC(10)*wg(8)+ 2*Re(wz(10)*wgC(8))); * al GtM2*GsM2*wz(10)*wgC(9)= GtM2*GsM2*(-wzC(10)*wg(9)+ 2*Re(wz(10)*wgC(9))); * id wg(n?)*Gs*GtM2*PwsC= GtM2*(-wgC(n)*GsC*Pws+2*Re(wg(n)*Gs*PwsC)); id wg(n?)*Gt*GsM2*PwtC= GsM2*(-wgC(n)*GtC*Pwt+2*Re(wg(n)*Gt*PwtC)); id wz(n?)*Gs*GtM2*PwsC= GtM2*(-wzC(n)*GsC*Pws+2*Re(wz(n)*Gs*PwsC)); id wz(n?)*Gt*GsM2*PwtC= GsM2*(-wzC(n)*GtC*Pwt+2*Re(wz(n)*Gt*PwtC)); * id GsM2*GtC*PwtC*wgC(n?)= GsM2*(-Gt*Pwt*wg(n)+2*Re(Gt*Pwt*wg(n))); id GtM2*GsC*PwsC*wgC(n?)= GtM2*(-Gs*Pws*wg(n)+2*Re(Gs*Pws*wg(n))); id GsM2*GtC*PwtC*wzC(n?)= GsM2*(-Gt*Pwt*wz(n)+2*Re(Gt*Pwt*wz(n))); id GtM2*GsC*PwsC*wzC(n?)= GtM2*(-Gs*Pws*wz(n)+2*Re(Gs*Pws*wz(n))); * id wg(n?)*Gs*GtM2= GtM2*(-GsC*wgC(n)+2*Re(Gs*wg(n))); id wg(n?)*Gt*GsM2= GsM2*(-GtC*wgC(n)+2*Re(Gt*wg(n))); id wz(n?)*Gs*GtM2= GtM2*(-GsC*wzC(n)+2*Re(Gs*wz(n))); id wz(n?)*Gt*GsM2= GsM2*(-GtC*wzC(n)+2*Re(Gt*wz(n))); * id GsM2*GtC*PwtC= GsM2*(-Gt*Pwt+2*Re(Gt*Pwt)); id GtM2*GsC*PwsC= GtM2*(-Gs*Pws+2*Re(Gs*Pws)); * id GsC*PwsC*Gt*Pwt= -Gs*Pws*GtC*PwtC+2*Re(GsC*Gt*PwsC*Pwt); *id GsC*Gt*Pws*PwtC= -Gs*GtC*PwsC*Pwt+2*Re(Gs*GtC*PwsC*Pwt); id GsC*Gt*Pws*PwtC= Re(Gs*GtC*PwsC*Pwt)+i_*Im(Gs*GtC*PwsC*Pwt); id Gs*GtC*PwsC*Pwt= Re(Gs*GtC*PwsC*Pwt)-i_*Im(Gs*GtC*PwsC*Pwt); id GsC*Gt*PwsC*PwtC= -Gs*GtC*Pws*Pwt+2*Re(Gs*GtC*Pws*Pwt); * id Gs*GtC*Pwt= -GsC*Gt*PwtC+2*Re(Gs*GtC*Pwt); * id GtM2*PwsC= GtM2*(-Pws+2*Re(Pws)); id GsM2*PwtC= GsM2*(-Pwt+2*Re(Pwt)); * id Re(wzC(2)*wg(1))= Re(wz(2)*wgC(1)); al Re(wzC(3)*wg(1))= Re(wz(3)*wgC(1)); al Re(wzC(4)*wg(1))= Re(wz(4)*wgC(1)); al Re(wzC(5)*wg(1))= Re(wz(5)*wgC(1)); al Re(wzC(6)*wg(1))= Re(wz(6)*wgC(1)); al Re(wzC(7)*wg(1))= Re(wz(7)*wgC(1)); al Re(wzC(8)*wg(1))= Re(wz(8)*wgC(1)); al Re(wzC(9)*wg(1))= Re(wz(9)*wgC(1)); al Re(wzC(10)*wg(1))= Re(wz(10)*wgC(1)); * id Re(wgC(2)*wz(1))= Re(wg(2)*wzC(1)); al Re(wgC(3)*wz(1))= Re(wg(3)*wzC(1)); al Re(wgC(4)*wz(1))= Re(wg(4)*wzC(1)); al Re(wgC(5)*wz(1))= Re(wg(5)*wzC(1)); al Re(wgC(6)*wz(1))= Re(wg(6)*wzC(1)); al Re(wgC(7)*wz(1))= Re(wg(7)*wzC(1)); al Re(wgC(8)*wz(1))= Re(wg(8)*wzC(1)); al Re(wgC(9)*wz(1))= Re(wg(9)*wzC(1)); al Re(wgC(10)*wz(1))= Re(wg(10)*wzC(1)); * id Re(wgC(3)*wz(2))= Re(wg(3)*wzC(2)); al Re(wgC(4)*wz(2))= Re(wg(4)*wzC(2)); al Re(wgC(5)*wz(2))= Re(wg(5)*wzC(2)); al Re(wgC(6)*wz(2))= Re(wg(6)*wzC(2)); al Re(wgC(7)*wz(2))= Re(wg(7)*wzC(2)); al Re(wgC(8)*wz(2))= Re(wg(8)*wzC(2)); al Re(wgC(9)*wz(2))= Re(wg(9)*wzC(2)); al Re(wgC(10)*wz(2))= Re(wg(10)*wzC(2)); * id Re(wzC(3)*wg(2))= Re(wz(3)*wgC(2)); al Re(wzC(4)*wg(2))= Re(wz(4)*wgC(2)); al Re(wzC(5)*wg(2))= Re(wz(5)*wgC(2)); al Re(wzC(6)*wg(2))= Re(wz(6)*wgC(2)); al Re(wzC(7)*wg(2))= Re(wz(7)*wgC(2)); al Re(wzC(8)*wg(2))= Re(wz(8)*wgC(2)); al Re(wzC(9)*wg(2))= Re(wz(9)*wgC(2)); al Re(wzC(10)*wg(2))= Re(wz(10)*wgC(2)); * id Re(wgC(4)*wz(3))= Re(wg(4)*wzC(3)); al Re(wgC(5)*wz(3))= Re(wg(5)*wzC(3)); al Re(wgC(6)*wz(3))= Re(wg(6)*wzC(3)); al Re(wgC(7)*wz(3))= Re(wg(7)*wzC(3)); al Re(wgC(8)*wz(3))= Re(wg(8)*wzC(3)); al Re(wgC(9)*wz(3))= Re(wg(9)*wzC(3)); al Re(wgC(10)*wz(3))= Re(wg(10)*wzC(3)); * id Re(wzC(4)*wg(3))= Re(wz(4)*wgC(3)); al Re(wzC(5)*wg(3))= Re(wz(5)*wgC(3)); al Re(wzC(6)*wg(3))= Re(wz(6)*wgC(3)); al Re(wzC(7)*wg(3))= Re(wz(7)*wgC(3)); al Re(wzC(8)*wg(3))= Re(wz(8)*wgC(3)); al Re(wzC(9)*wg(3))= Re(wz(9)*wgC(3)); al Re(wzC(10)*wg(3))= Re(wz(10)*wgC(3)); * id Re(wgC(5)*wz(4))= Re(wg(5)*wzC(4)); al Re(wgC(6)*wz(4))= Re(wg(6)*wzC(4)); al Re(wgC(7)*wz(4))= Re(wg(7)*wzC(4)); al Re(wgC(8)*wz(4))= Re(wg(8)*wzC(4)); al Re(wgC(9)*wz(4))= Re(wg(9)*wzC(4)); al Re(wgC(10)*wz(4))= Re(wg(10)*wzC(4)); * id Re(wzC(5)*wg(4))= Re(wg(5)*wzC(4)); al Re(wzC(6)*wg(4))= Re(wg(6)*wzC(4)); al Re(wzC(7)*wg(4))= Re(wg(7)*wzC(4)); al Re(wzC(8)*wg(4))= Re(wg(8)*wzC(4)); al Re(wzC(9)*wg(4))= Re(wg(9)*wzC(4)); al Re(wzC(10)*wg(4))= Re(wg(10)*wzC(4)); * id Re(wgC(6)*wz(5))= Re(wg(6)*wzC(5)); al Re(wgC(7)*wz(5))= Re(wg(7)*wzC(5)); al Re(wgC(8)*wz(5))= Re(wg(8)*wzC(5)); al Re(wgC(9)*wz(5))= Re(wg(9)*wzC(5)); al Re(wgC(10)*wz(5))= Re(wg(10)*wzC(5)); * id Re(wzC(6)*wg(5))= Re(wz(6)*wgC(5)); al Re(wzC(7)*wg(5))= Re(wz(7)*wgC(5)); al Re(wzC(8)*wg(5))= Re(wz(8)*wgC(5)); al Re(wzC(9)*wg(5))= Re(wz(9)*wgC(5)); al Re(wzC(10)*wg(5))= Re(wz(10)*wgC(5)); * id Re(wgC(7)*wz(6))= Re(wg(7)*wzC(6)); al Re(wgC(8)*wz(6))= Re(wg(8)*wzC(6)); al Re(wgC(9)*wz(6))= Re(wg(9)*wzC(6)); al Re(wgC(10)*wz(6))= Re(wg(10)*wzC(6)); * id Re(wzC(7)*wg(6))= Re(wz(7)*wgC(6)); al Re(wzC(8)*wg(6))= Re(wz(8)*wgC(6)); al Re(wzC(9)*wg(6))= Re(wz(9)*wgC(6)); al Re(wzC(10)*wg(6))= Re(wz(10)*wgC(6)); * id Re(wgC(8)*wz(7))= Re(wg(8)*wzC(7)); al Re(wgC(9)*wz(7))= Re(wg(9)*wzC(7)); al Re(wgC(10)*wz(7))= Re(wg(10)*wzC(7)); * id Re(wzC(8)*wg(7))= Re(wz(8)*wgC(7)); al Re(wzC(9)*wg(7))= Re(wz(9)*wgC(7)); al Re(wzC(10)*wg(7))= Re(wz(10)*wgC(7)); * id Re(wgC(9)*wz(8))= Re(wg(9)*wzC(8)); al Re(wgC(10)*wz(8))= Re(wg(10)*wzC(8)); * id Re(wzC(9)*wg(8))= Re(wg(9)*wgC(8)); al Re(wzC(10)*wg(8))= Re(wg(10)*wgC(8)); * id Re(wgC(10)*wz(9))= Re(wg(10)*wzC(9)); * id Re(wzC(10)*wg(9))= Re(wz(10)*wgC(9)); * id wg(n?)*wzC(n?)= -wgC(n)*wz(n)+2*Re(wg(n)*wzC(n)); * .store * #endprocedure * #procedure contractgz(dummy) * g Wd_gz= T_gz*d_(mu,nu); g Wpp_gz= T_gz*pp(mu)*pp(nu); g We_gz= T_gz*e_(mu,nu,Qm,pp); * contract; * id rs^-1= -1/(X*y*s); id pp.pp= 0; al qp.qp= 0; al Qm.Qm= QmS; id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; id qp.Qp= qp.pp; * id pp.Qm= -1/2*(QmS+y*s); al k.kb= 1/2*(-sp+ums+dms); al Qm.k= 1/2*(QmS+t-ums); al pp.kb= 1/2*(tp-dms); al Qm.kb= 1/2*(QmS+u-dms); al pp.k= 1/2*(up-ums); al pp.qp= 1/2*kp; al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(dms-zp); al k.qp= 1/2*(ums-zm); * id sp= x2*y*s; al km= z*y*s; al km^-1= 1/(z*y*s); al zm= -(x2-x1)*y*s; id QmS= X*y*s; * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id kb(mu?)= kbf(mu); id qp(mu?)= qpf(mu); al e_(kb,mu?,nu?,al?)= kbf(ab1)*e_(ab1,mu,nu,al); al e_(qp,mu?,nu?,al?)= qpf(ab2)*e_(ab2,mu,nu,al); * #call integration{dummy} * id rs^-1= -1/(X*y*s); id pp.pp= 0; al qp.qp= 0; al Qm.Qm= QmS; al kb.kb= -dms; id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; id qp.Qp= qp.pp; * id pp.Qm= -1/2*(QmS+y*s); al k.kb= 1/2*(-sp+ums+dms); al Qm.k= 1/2*(QmS+t-ums); al pp.kb= 1/2*(tp-dms); al Qm.kb= 1/2*(QmS+u-dms); al pp.k= 1/2*(up-ums); al pp.qp= 1/2*kp; al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(dms-zp); al k.qp= 1/2*(ums-zm); * id sp= x2*y*s; al km= z*y*s; al km^-1= 1/(z*y*s); al zm= -(x2-x1)*y*s; id QmS= X*y*s; * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * .store * g GZ10= (-1/2*Wd_gz-2*X/Xp^2*Wpp_gz/(y*s))*sel0*s^2; g GZ20= (-1/2*Wd_gz-6*X/Xp^2*Wpp_gz/(y*s))*sel0*s^2; g GZ30= (We_gz/(y*s)/Xp)*sel0*s^2; g GZ11= (-1/2*Wd_gz-2*X/Xp^2*Wpp_gz/(y*s))*seld*s^2; g GZ21= (-1/2*Wd_gz-6*X/Xp^2*Wpp_gz/(y*s))*seld*s^2; g GZ31= (We_gz/(y*s)/Xp)*seld*s^2; g GZ12= (-1/2*Wd_gz-2*X/Xp^2*Wpp_gz/(y*s))*seldd*s^2; g GZ22= (-1/2*Wd_gz-6*X/Xp^2*Wpp_gz/(y*s))*seldd*s^2; g GZ32= (We_gz/(y*s)/Xp)*seldd*s^2; * ratio,tdM,ups,aux; id aux= dms/(y*s); * id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * id tdM^-2*seldd= 1; id seldd= 0; * id tdM^-2*seld= 0; id tdM^-1*seld= 1; id seld= 0; * id tdM^-2*sel0= 0; id tdM^-1*sel0= 0; * id sel0= 1; * ratio,X,Xp,one; id one= 1; *-- id Re(wg(n?)*wzC(2))= 1/2*(Re(wg(n)*wzC([2+3]))+Re(wg(n)*wzC([2-3]))); id Re(wg(n?)*wzC(3))= 1/2*(Re(wg(n)*wzC([2+3]))-Re(wg(n)*wzC([2-3]))); id Re(wz(n?)*wgC(2))= 1/2*(Re(wz(n)*wgC([2+3]))+Re(wz(n)*wgC([2-3]))); id Re(wz(n?)*wgC(3))= 1/2*(Re(wz(n)*wgC([2+3]))-Re(wz(n)*wgC([2-3]))); * id Re(wg(2)*wzC(n?))= 1/2*(Re(wg([2+3])*wzC(n))+Re(wg([2-3])*wzC(n))); id Re(wg(3)*wzC(n?))= 1/2*(Re(wg([2+3])*wzC(n))-Re(wg([2-3])*wzC(n))); id Re(wz(2)*wgC(n?))= 1/2*(Re(wz([2+3])*wgC(n))+Re(wz([2-3])*wgC(n))); id Re(wz(3)*wgC(n?))= 1/2*(Re(wz([2+3])*wgC(n))-Re(wz([2-3])*wgC(n))); * id Re(wg(n?)*wzC(5))= 1/2*(Re(wg(n)*wzC([5+6]))+Re(wg(n)*wzC([5-6]))); id Re(wg(n?)*wzC(6))= 1/2*(Re(wg(n)*wzC([5+6]))-Re(wg(n)*wzC([5-6]))); id Re(wz(n?)*wgC(5))= 1/2*(Re(wz(n)*wgC([5+6]))+Re(wz(n)*wgC([5-6]))); id Re(wz(n?)*wgC(6))= 1/2*(Re(wz(n)*wgC([5+6]))-Re(wz(n)*wgC([5-6]))); * id Re(wg(5)*wzC(n?))= 1/2*(Re(wg([5+6])*wzC(n))+Re(wg([5-6])*wzC(n))); id Re(wg(6)*wzC(n?))= 1/2*(Re(wg([5+6])*wzC(n))-Re(wg([5-6])*wzC(n))); id Re(wz(5)*wgC(n?))= 1/2*(Re(wz([5+6])*wgC(n))+Re(wz([5-6])*wgC(n))); id Re(wz(6)*wgC(n?))= 1/2*(Re(wz([5+6])*wgC(n))-Re(wz([5-6])*wgC(n))); * id Re(wg(2)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)+Re(wg([2-3])*X*Y)); id Re(wg(3)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)-Re(wg([2-3])*X*Y)); id Re(wg(5)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)+Re(wg([5-6])*X*Y)); id Re(wg(6)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)-Re(wg([5-6])*X*Y)); id Re(wz(2)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)+Re(wz([2-3])*X*Y)); id Re(wz(3)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)-Re(wz([2-3])*X*Y)); id Re(wz(5)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)+Re(wz([5-6])*X*Y)); id Re(wz(6)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)-Re(wz([5-6])*X*Y)); * id Re(wg(2)*X?)= 1/2*(Re(wg([2+3])*X)+Re(wg([2-3])*X)); id Re(wg(3)*X?)= 1/2*(Re(wg([2+3])*X)-Re(wg([2-3])*X)); id Re(wg(5)*X?)= 1/2*(Re(wg([5+6])*X)+Re(wg([5-6])*X)); id Re(wg(6)*X?)= 1/2*(Re(wg([5+6])*X)-Re(wg([5-6])*X)); id Re(wz(2)*X?)= 1/2*(Re(wz([2+3])*X)+Re(wz([2-3])*X)); id Re(wz(3)*X?)= 1/2*(Re(wz([2+3])*X)-Re(wz([2-3])*X)); id Re(wz(5)*X?)= 1/2*(Re(wz([5+6])*X)+Re(wz([5-6])*X)); id Re(wz(6)*X?)= 1/2*(Re(wz([5+6])*X)-Re(wz([5-6])*X)); * id Re(wg([2+3])*wzC([2-3]))= Re(wz([2-3])*wgC([2+3])); id Re(wg([2+3])*wzC([5-6]))= Re(wz([5-6])*wgC([2+3])); id Re(wg([2-3])*wzC([5+6]))= Re(wz([5+6])*wgC([2-3])); id Re(wg([5+6])*wzC([5-6]))= Re(wz([5-6])*wgC([5+6])); id Re(wz([2+3])*wgC([2-3]))= Re(wg([2-3])*wzC([2+3])); id Re(wz([2+3])*wgC([5-6]))= Re(wg([5-6])*wzC([2+3])); id Re(wz([2-3])*wgC([5+6]))= Re(wg([5+6])*wzC([2-3])); id Re(wz([5+6])*wgC([5-6]))= Re(wg([5-6])*wzC([5+6])); * id X= [1+X]-1; id x1= 1-[1-x1]; id x2= 1-[1-x1]-[x1-x2]; id z= 1-[1-x1]-[x1-x2]-[x2-z]; id ups= [ups+X+x1]-[1+X]+[1-x1]; id Xp^-1= 1/[1+X]; ratio,X,[1+X],one; id one= 1; al one^-1= 1; * id Re(wg(X?)*wzC(Y?))= Rep(X,Y); id Re(wz(X?)*wgC(Y?))= Rep(Y,X); id GsM2*GtM2= GstM2; id [1+X] = Xp; al [1+X]^-1= Xpi; id [1-x1]= omx1; al [x1-x2]= x12; al [x2-z]= x2z; al [ups+X+x1]= ux; al zX^-1= zXi; id Xpi^2= Xpi-aux; id aux= XX; * id Xp^2= Xps; id Xpi^2= Xpis; id y^3= yc; id y^2= ys; * id omx1= O; id O^4= Oq; id O^3= Oc; id O^2= Os; id x12^3= x12c; id x12^2= x12s; id x2z^4= x2zq; id x2z^3= x2zc; id x2z^2= x2zs; id ux^2= uxs; id vv^2= vvs; id vv^3= vvc; id vvs*vv= vvc; id vvs^2= vvq; id vvc*vv= vvq; id z^-1= zi; * id Re(Gs*GtC*Pwt*PwsC)= rl1; id Re(Gs*GtC*Pwt)= rl2; id Re(Pws)= rl3; id Re(wg(n?)*Gs*PwsC)= rs1wg(n); id Re(wz(n?)*Gs*PwsC)= rs1wz(n); id Re(wg(n?)*Gs)= rs2wg(n); id Re(wz(n?)*Gs)= rs2wz(n); id Re(wg(n?)*Gt*PwtC)= rt1wg(n); id Re(wz(n?)*Gt*PwtC)= rt1wz(n); id Re(wg(n?)*Gt)= rt2wg(n); id Re(wz(n?)*Gt)= rt2wz(n); id Rep(n1?,n2?)= rwwgz(n1,n2); id Re(wg(n1?)*wzC(n2?))= rwwgz(n1,n2); id Re(wz(n1?)*wgC(n2?))= rwwgz(n2,n1); * id fun?([2+3],n?)= fun(11,n); al fun?([2-3],n?)= fun(12,n); al fun?([5+6],n?)= fun(13,n); al fun?([5-6],n?)= fun(14,n); * id fun?(n?,[2+3])= fun(n,11); al fun?(n?,[2-3])= fun(n,12); al fun?(n?,[5+6])= fun(n,13); al fun?(n?,[5-6])= fun(n,14); * id fun?([2+3])= fun(11); al fun?([2-3])= fun(12); al fun?([5+6])= fun(13); al fun?([5-6])= fun(14); *-- id ve(T)= -1/2+ss(T); al vu(T)= 1/2-2/3*ss(T); al vd(T)= -1/2+1/3*ss(T); * id ss(T)*csi(T)= csi(T)-1; id csi(T)= 1/cst; al ss(T)= sst; * .store * *----g^2-------------------- * g GZ10s= GZ10; g GZ20s= GZ20; g GZ30s= GZ30; g GZ11s= GZ11; g GZ21s= GZ21; g GZ31s= GZ31; g GZ12s= GZ12; g GZ22s= GZ22; g GZ32s= GZ32; * id gexp^4= 0; id gexp^3= 0; id gexp^2= 1; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,rs1wg,rs2wg,rt1wg,rt2wg,rwwg, rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, zXi,XX,vv,vvs,vvc,vvq,cst,sst,PZ; format fortran; print; * .store * *----g^3-------------------- * g GZ10c= GZ10; g GZ20c= GZ20; g GZ30c= GZ30; g GZ11c= GZ11; g GZ21c= GZ21; g GZ31c= GZ31; g GZ12c= GZ12; g GZ22c= GZ22; g GZ32c= GZ32; * id gexp^4= 0; id gexp^3= 1; id gexp^2= 0; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,rs1wg,rs2wg,rt1wg,rt2wg,rwwg, rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, zXi,XX,vv,vvs,vvc,vvq,cst,sst,PZ; format fortran; print; * * .store * *----g^4-------------------- * g GZ10q= GZ10; g GZ20q= GZ20; g GZ30q= GZ30; g GZ11q= GZ11; g GZ21q= GZ21; g GZ31q= GZ31; g GZ12q= GZ12; g GZ22q= GZ22; g GZ32q= GZ32; * id gexp^4= 1; id gexp^3= 0; id gexp^2= 0; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,rs1wg,rs2wg,rt1wg,rt2wg,rwwg, rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, zXi,XX,vv,vvs,vvc,vvq,cst,sst,PZ; format fortran; print; * .store * #endprocedure * #procedure totalzz(dummy) * * Now the full stuff * load wrl$grp04:[giampiero]outputampzz_B.sav; * g T_zz= PwsM2*PwtM2*wwww_zz/(GT^2*I0); *nwrite statistics; * id wz(n?)= wz(n)*rs; al wzC(n?)= wzC(n)*rs; * id wz(4)= X*wz(4); al wz(7)= X*wz(7); al wz(8)= X*wz(8); al wz(10)= X*wz(10); id wzC(4)= X*wzC(4); al wzC(7)= X*wzC(7); al wzC(8)= X*wzC(8); al wzC(10)= X*wzC(10); * id wz(n?)*wzC(n?)= wzM2(n); * id GtM2*GsM2*wz(2)*wzC(1)= GtM2*GsM2*(-wzC(2)*wz(1)+ 2*Re(wz(2)*wzC(1))); al GtM2*GsM2*wz(3)*wzC(1)= GtM2*GsM2*(-wzC(3)*wz(1)+ 2*Re(wz(3)*wzC(1))); al GtM2*GsM2*wz(4)*wzC(1)= GtM2*GsM2*(-wzC(4)*wz(1)+ 2*Re(wz(4)*wzC(1))); al GtM2*GsM2*wz(5)*wzC(1)= GtM2*GsM2*(-wzC(5)*wz(1)+ 2*Re(wz(5)*wzC(1))); al GtM2*GsM2*wz(6)*wzC(1)= GtM2*GsM2*(-wzC(6)*wz(1)+ 2*Re(wz(6)*wzC(1))); al GtM2*GsM2*wz(7)*wzC(1)= GtM2*GsM2*(-wzC(7)*wz(1)+ 2*Re(wz(7)*wzC(1))); al GtM2*GsM2*wz(8)*wzC(1)= GtM2*GsM2*(-wzC(8)*wz(1)+ 2*Re(wz(8)*wzC(1))); al GtM2*GsM2*wz(9)*wzC(1)= GtM2*GsM2*(-wzC(9)*wz(1)+ 2*Re(wz(9)*wzC(1))); al GtM2*GsM2*wz(10)*wzC(1)= GtM2*GsM2*(-wzC(10)*wz(1)+ 2*Re(wz(10)*wzC(1))); * al GtM2*GsM2*wz(3)*wzC(2)= GtM2*GsM2*(-wzC(3)*wz(2)+ 2*Re(wz(3)*wzC(2))); al GtM2*GsM2*wz(4)*wzC(2)= GtM2*GsM2*(-wzC(4)*wz(2)+ 2*Re(wz(4)*wzC(2))); al GtM2*GsM2*wz(5)*wzC(2)= GtM2*GsM2*(-wzC(5)*wz(2)+ 2*Re(wz(5)*wzC(2))); al GtM2*GsM2*wz(6)*wzC(2)= GtM2*GsM2*(-wzC(6)*wz(2)+ 2*Re(wz(6)*wzC(2))); al GtM2*GsM2*wz(7)*wzC(2)= GtM2*GsM2*(-wzC(7)*wz(2)+ 2*Re(wz(7)*wzC(2))); al GtM2*GsM2*wz(8)*wzC(2)= GtM2*GsM2*(-wzC(8)*wz(2)+ 2*Re(wz(8)*wzC(2))); al GtM2*GsM2*wz(9)*wzC(2)= GtM2*GsM2*(-wzC(9)*wz(2)+ 2*Re(wz(9)*wzC(2))); al GtM2*GsM2*wz(10)*wzC(2)= GtM2*GsM2*(-wzC(10)*wz(2)+ 2*Re(wz(10)*wzC(2))); * al GtM2*GsM2*wz(4)*wzC(3)= GtM2*GsM2*(-wzC(4)*wz(3)+ 2*Re(wz(4)*wzC(3))); al GtM2*GsM2*wz(5)*wzC(3)= GtM2*GsM2*(-wzC(5)*wz(3)+ 2*Re(wz(5)*wzC(3))); al GtM2*GsM2*wz(6)*wzC(3)= GtM2*GsM2*(-wzC(6)*wz(3)+ 2*Re(wz(6)*wzC(3))); al GtM2*GsM2*wz(7)*wzC(3)= GtM2*GsM2*(-wzC(7)*wz(3)+ 2*Re(wz(7)*wzC(3))); al GtM2*GsM2*wz(8)*wzC(3)= GtM2*GsM2*(-wzC(8)*wz(3)+ 2*Re(wz(8)*wzC(3))); al GtM2*GsM2*wz(9)*wzC(3)= GtM2*GsM2*(-wzC(9)*wz(3)+ 2*Re(wz(9)*wzC(3))); al GtM2*GsM2*wz(10)*wzC(3)= GtM2*GsM2*(-wzC(10)*wz(3)+ 2*Re(wz(10)*wzC(3))); * al GtM2*GsM2*wz(5)*wzC(4)= GtM2*GsM2*(-wzC(5)*wz(4)+ 2*Re(wz(5)*wzC(4))); al GtM2*GsM2*wz(6)*wzC(4)= GtM2*GsM2*(-wzC(6)*wz(4)+ 2*Re(wz(6)*wzC(4))); al GtM2*GsM2*wz(7)*wzC(4)= GtM2*GsM2*(-wzC(7)*wz(4)+ 2*Re(wz(7)*wzC(4))); al GtM2*GsM2*wz(8)*wzC(4)= GtM2*GsM2*(-wzC(8)*wz(4)+ 2*Re(wz(8)*wzC(4))); al GtM2*GsM2*wz(9)*wzC(4)= GtM2*GsM2*(-wzC(9)*wz(4)+ 2*Re(wz(9)*wzC(4))); al GtM2*GsM2*wz(10)*wzC(4)= GtM2*GsM2*(-wzC(10)*wz(4)+ 2*Re(wz(10)*wzC(4))); * al GtM2*GsM2*wz(6)*wzC(5)= GtM2*GsM2*(-wzC(6)*wz(5)+ 2*Re(wz(6)*wzC(5))); al GtM2*GsM2*wz(7)*wzC(5)= GtM2*GsM2*(-wzC(7)*wz(5)+ 2*Re(wz(7)*wzC(5))); al GtM2*GsM2*wz(8)*wzC(5)= GtM2*GsM2*(-wzC(8)*wz(5)+ 2*Re(wz(8)*wzC(5))); al GtM2*GsM2*wz(9)*wzC(5)= GtM2*GsM2*(-wzC(9)*wz(5)+ 2*Re(wz(9)*wzC(5))); al GtM2*GsM2*wz(10)*wzC(5)= GtM2*GsM2*(-wzC(10)*wz(5)+ 2*Re(wz(10)*wzC(5))); * al GtM2*GsM2*wz(7)*wzC(6)= GtM2*GsM2*(-wzC(7)*wz(6)+ 2*Re(wz(7)*wzC(6))); al GtM2*GsM2*wz(8)*wzC(6)= GtM2*GsM2*(-wzC(8)*wz(6)+ 2*Re(wz(8)*wzC(6))); al GtM2*GsM2*wz(9)*wzC(6)= GtM2*GsM2*(-wzC(9)*wz(6)+ 2*Re(wz(9)*wzC(6))); al GtM2*GsM2*wz(10)*wzC(6)= GtM2*GsM2*(-wzC(10)*wz(6)+ 2*Re(wz(10)*wzC(6))); * al GtM2*GsM2*wz(8)*wzC(7)= GtM2*GsM2*(-wzC(8)*wz(7)+ 2*Re(wz(8)*wzC(7))); al GtM2*GsM2*wz(9)*wzC(7)= GtM2*GsM2*(-wzC(9)*wz(7)+ 2*Re(wz(9)*wzC(7))); al GtM2*GsM2*wz(10)*wzC(7)= GtM2*GsM2*(-wzC(10)*wz(7)+ 2*Re(wz(10)*wzC(7))); * al GtM2*GsM2*wz(9)*wzC(8)= GtM2*GsM2*(-wzC(9)*wz(8)+ 2*Re(wz(9)*wzC(8))); al GtM2*GsM2*wz(10)*wzC(8)= GtM2*GsM2*(-wzC(10)*wz(8)+ 2*Re(wz(10)*wzC(8))); * al GtM2*GsM2*wz(10)*wzC(9)= GtM2*GsM2*(-wzC(10)*wz(9)+ 2*Re(wz(10)*wzC(9))); * id wz(n?)*Gs*GtM2*PwsC= GtM2*(-wzC(n)*GsC*Pws+2*Re(wz(n)*Gs*PwsC)); id wz(n?)*Gt*GsM2*PwtC= GsM2*(-wzC(n)*GtC*Pwt+2*Re(wz(n)*Gt*PwtC)); id wz(n?)*Gs*GtM2= GtM2*(-GsC*wzC(n)+2*Re(Gs*wz(n))); id wz(n?)*Gt*GsM2= GsM2*(-GtC*wzC(n)+2*Re(Gt*wz(n))); * id GsM2*GtC*PwtC*wzC(n?)= GsM2*(-Gt*Pwt*wz(n)+2*Re(Gt*Pwt*wz(n))); id GtM2*GsC*PwsC*wzC(n?)= GtM2*(-Gs*Pws*wz(n)+2*Re(Gs*Pws*wz(n))); * id GsM2*GtC*PwtC= GsM2*(-Gt*Pwt+2*Re(Gt*Pwt)); id GtM2*GsC*PwsC= GtM2*(-Gs*Pws+2*Re(Gs*Pws)); * id GsC*PwsC*Gt*Pwt= -Gs*Pws*GtC*PwtC+2*Re(GsC*Gt*PwsC*Pwt); id GsC*Gt*Pws*PwtC= Re(Gs*GtC*PwsC*Pwt)-i_*Im(Gs*GtC*PwsC*Pwt); id Gs*GtC*PwsC*Pwt= Re(Gs*GtC*PwsC*Pwt)+i_*Im(Gs*GtC*PwsC*Pwt); id GsC*Gt*PwsC*PwtC= -Gs*GtC*Pws*Pwt+2*Re(Gs*GtC*Pws*Pwt); * id Gs*GtC*Pwt= -GsC*Gt*PwtC+2*Re(Gs*GtC*Pwt); * id GtM2*PwsC= GtM2*(-Pws+2*Re(Pws)); id GsM2*PwtC= GsM2*(-Pwt+2*Re(Pwt)); * id Re(wzC(2)*wz(1))= Re(wz(2)*wzC(1)); al Re(wzC(3)*wz(1))= Re(wz(3)*wzC(1)); al Re(wzC(4)*wz(1))= Re(wz(4)*wzC(1)); al Re(wzC(5)*wz(1))= Re(wz(5)*wzC(1)); al Re(wzC(6)*wz(1))= Re(wz(6)*wzC(1)); al Re(wzC(7)*wz(1))= Re(wz(7)*wzC(1)); al Re(wzC(8)*wz(1))= Re(wz(8)*wzC(1)); al Re(wzC(9)*wz(1))= Re(wz(9)*wzC(1)); al Re(wzC(10)*wz(1))= Re(wz(10)*wzC(1)); * id Re(wzC(3)*wz(2))= Re(wz(3)*wzC(2)); al Re(wzC(4)*wz(2))= Re(wz(4)*wzC(2)); al Re(wzC(5)*wz(2))= Re(wz(5)*wzC(2)); al Re(wzC(6)*wz(2))= Re(wz(6)*wzC(2)); al Re(wzC(7)*wz(2))= Re(wz(7)*wzC(2)); al Re(wzC(8)*wz(2))= Re(wz(8)*wzC(2)); al Re(wzC(9)*wz(2))= Re(wz(9)*wzC(2)); al Re(wzC(10)*wz(2))= Re(wz(10)*wzC(2)); * id Re(wzC(4)*wz(3))= Re(wz(4)*wzC(3)); al Re(wzC(5)*wz(3))= Re(wz(5)*wzC(3)); al Re(wzC(6)*wz(3))= Re(wz(6)*wzC(3)); al Re(wzC(7)*wz(3))= Re(wz(7)*wzC(3)); al Re(wzC(8)*wz(3))= Re(wz(8)*wzC(3)); al Re(wzC(9)*wz(3))= Re(wz(9)*wzC(3)); al Re(wzC(10)*wz(3))= Re(wz(10)*wzC(3)); * id Re(wzC(5)*wz(4))= Re(wz(5)*wzC(4)); al Re(wzC(6)*wz(4))= Re(wz(6)*wzC(4)); al Re(wzC(7)*wz(4))= Re(wz(7)*wzC(4)); al Re(wzC(8)*wz(4))= Re(wz(8)*wzC(4)); al Re(wzC(9)*wz(4))= Re(wz(9)*wzC(4)); al Re(wzC(10)*wz(4))= Re(wz(10)*wzC(4)); * id Re(wzC(6)*wz(5))= Re(wz(6)*wzC(5)); al Re(wzC(7)*wz(5))= Re(wz(7)*wzC(5)); al Re(wzC(8)*wz(5))= Re(wz(8)*wzC(5)); al Re(wzC(9)*wz(5))= Re(wz(9)*wzC(5)); al Re(wzC(10)*wz(5))= Re(wz(10)*wzC(5)); * id Re(wzC(7)*wz(6))= Re(wz(7)*wzC(6)); al Re(wzC(8)*wz(6))= Re(wz(8)*wzC(6)); al Re(wzC(9)*wz(6))= Re(wz(9)*wzC(6)); al Re(wzC(10)*wz(6))= Re(wz(10)*wzC(6)); * id Re(wzC(8)*wz(7))= Re(wz(8)*wzC(7)); al Re(wzC(9)*wz(7))= Re(wz(9)*wzC(7)); al Re(wzC(10)*wz(7))= Re(wz(10)*wzC(7)); * id Re(wzC(9)*wz(8))= Re(wz(9)*wzC(8)); al Re(wzC(10)*wz(8))= Re(wz(10)*wzC(8)); * id Re(wzC(10)*wz(9))= Re(wz(10)*wzC(9)); * .store * #endprocedure * #procedure contractzz(dummy) * g Wd_zz= T_zz*d_(mu,nu); g Wpp_zz= T_zz*pp(mu)*pp(nu); g Wqq_zz= T_zz*Qm(mu)*Qm(nu); g Wpq_zz= T_zz*(pp(mu)*Qm(nu)+pp(nu)*Qm(mu)); g We_zz= T_zz*e_(mu,nu,Qm,pp); * contract; * id rs^-1= -1/(X*y*s); id pp.pp= 0; al qp.qp= 0; al Qm.Qm= QmS; al kb.kb= -dms; id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; id qp.Qp= qp.pp; * id pp.Qm= -1/2*(QmS+y*s); al k.kb= 1/2*(-sp+ums+dms); al Qm.k= 1/2*(QmS+t-ums); al pp.kb= 1/2*(tp-dms); al Qm.kb= 1/2*(QmS+u-dms); al pp.k= 1/2*(up-ums); al pp.qp= 1/2*kp; al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(dms-zp); al k.qp= 1/2*(ums-zm); * id sp= x2*y*s; al km= z*y*s; al km^-1= 1/(z*y*s); al zm= -(x2-x1)*y*s; id QmS= X*y*s; * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id kb(mu?)= kbf(mu); id qp(mu?)= qpf(mu); al e_(kb,mu?,nu?,al?)= kbf(ab1)*e_(ab1,mu,nu,al); al e_(qp,mu?,nu?,al?)= qpf(ab2)*e_(ab2,mu,nu,al); * id uM^-1= -1/(y*s)/tdM; * id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * #call integration{dummy} * id rs^-1= -1/(X*y*s); id pp.pp= 0; al qp.qp= 0; al Qm.Qm= QmS; id Qp.Qm= pp.Qm-qp.Qm; id kb.Qp= kb.pp-kb.qp; id k.Qp= k.pp-k.qp; id pp.Qp= -pp.qp; id Qp.Qp= -2*pp.qp; id qp.Qp= qp.pp; * id pp.Qm= -1/2*(QmS+y*s); al k.kb= 1/2*(-sp+ums+dms); al Qm.k= 1/2*(QmS+t-ums); al pp.kb= 1/2*(tp-dms); al Qm.kb= 1/2*(QmS+u-dms); al pp.k= 1/2*(up-ums); al pp.qp= 1/2*kp; al Qm.qp= 1/2*(QmS+km); al kb.qp= 1/2*(dms-zp); al k.qp= 1/2*(ums-zm); * id sp= x2*y*s; al km= z*y*s; al km^-1= 1/(z*y*s); al zm= -(x2-x1)*y*s; id QmS= X*y*s; * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * id t= y*s * ( - 1 - ups - z - 2*X ) + ums * ( 1 ) + dms * ( 1 ); al kp= y*s * ( - 1 + x2 - z - X ); al u= y*s * ( ups ); al up= y*s * ( 1 - x1 + ups + z + X ) + ums * ( 1 ); al tp= y*s * ( - 1 + x1 - x2 - ups - X ) + dms * ( 1 ); al zp= y*s * ( 1 - x1 ) + ums * ( 1 ) + dms * ( 1 ); * id ums= 0; id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * .store * g ZZ10= (-1/2*Wd_zz-2*X/Xp^2*Wpp_zz/(y*s)-Wpq_zz/(y*s)/Xp)*sel0*s^3; g ZZ20= (-1/2*Wd_zz-6*X/Xp^2*Wpp_zz/(y*s)-3*Wpq_zz/(y*s)/Xp- Wqq_zz/(y*s)/Xp)*sel0*s^3; g ZZ30= (We_zz/(y*s)/Xp)*sel0*s^3; g ZZ40= (-2/Xp*Wqq_zz/(y*s)-2*X/Xp^2*Wpq_zz/(y*s))*sel0*s^3; *g W5_zz_0u= (-Wqq_zz/(y*s)/X-2/Xp*Wpq_zz/(y*s))*sel0*s^3; g ZZ11= (-1/2*Wd_zz-2*X/Xp^2*Wpp_zz/(y*s)-Wpq_zz/(y*s)/Xp)*seld*s^3; g ZZ21= (-1/2*Wd_zz-6*X/Xp^2*Wpp_zz/(y*s)-3*Wpq_zz/(y*s)/Xp- Wqq_zz/(y*s)/Xp)*seld*s^3; g ZZ31= (We_zz/(y*s)/Xp)*seld*s^3; g ZZ41= (-2/Xp*Wqq_zz/(y*s)-2*X/Xp^2*Wpq_zz/(y*s))*seld*s^3; *g W5_zz_u= (-Wqq_zz/(y*s)/X-2/Xp*Wpq_zz/(y*s))*seld*s^3; g ZZ12= (-1/2*Wd_zz-2*X/Xp^2*Wpp_zz/(y*s)-Wpq_zz/(y*s)/Xp)*seldd*s^3; g ZZ22= (-1/2*Wd_zz-6*X/Xp^2*Wpp_zz/(y*s)-3*Wpq_zz/(y*s)/Xp- Wqq_zz/(y*s)/Xp)*seldd*s^3; g ZZ32= (We_zz/(y*s)/Xp)*seldd*s^3; g ZZ42= (-2/Xp*Wqq_zz/(y*s)-2*X/Xp^2*Wpq_zz/(y*s))*seldd*s^3; *g W5_zz_uu= (-Wqq_zz/(y*s)/X-2/Xp*Wpq_zz/(y*s))*seldd*s^3; * ratio,tdM,ups,aux; id aux= dms/(y*s); * id dms*tdM^-2= aux; id dms= 0; id aux= dms*tdM^-2; * id tdM^-2*seldd= 1; id seldd= 0; * id tdM^-2*seld= 0; id tdM^-1*seld= 1; id seld= 0; * id tdM^-2*sel0= 0; id tdM^-1*sel0= 0; id sel0= 1; * ratio,X,Xp,one; id one= 1; al one^-1= 1; * id cs(T)*csi(T)= 1; * *-- id Re(wz(n?)*wzC(2))= 1/2*(Re(wz(n)*wzC([2+3]))+Re(wz(n)*wzC([2-3]))); id Re(wz(n?)*wzC(3))= 1/2*(Re(wz(n)*wzC([2+3]))-Re(wz(n)*wzC([2-3]))); id Re(wz(2)*wzC(2))= wzM2(2); id Re(wz(3)*wzC(3))= wzM2(3); * id Re(wz(2)*wzC(n?))= 1/2*(Re(wz([2+3])*wzC(n))+Re(wz([2-3])*wzC(n))); id Re(wz(3)*wzC(n?))= 1/2*(Re(wz([2+3])*wzC(n))-Re(wz([2-3])*wzC(n))); id Re(wz(2)*wzC(2))= wzM2(2); id Re(wz(3)*wzC(3))= wzM2(3); * id Re(wz(n?)*wzC(5))= 1/2*(Re(wz(n)*wzC([5+6]))+Re(wz(n)*wzC([5-6]))); id Re(wz(n?)*wzC(6))= 1/2*(Re(wz(n)*wzC([5+6]))-Re(wz(n)*wzC([5-6]))); id Re(wz(5)*wzC(5))= wzM2(5); id Re(wz(6)*wzC(6))= wzM2(6); * id Re(wz(5)*wzC(n?))= 1/2*(Re(wz([5+6])*wzC(n))+Re(wz([5-6])*wzC(n))); id Re(wz(6)*wzC(n?))= 1/2*(Re(wz([5+6])*wzC(n))-Re(wz([5-6])*wzC(n))); id Re(wz(5)*wzC(5))= wzM2(5); id Re(wz(6)*wzC(6))= wzM2(6); * id Re(wz(2)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)+Re(wz([2-3])*X*Y)); id Re(wz(3)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)-Re(wz([2-3])*X*Y)); id Re(wz(5)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)+Re(wz([5-6])*X*Y)); id Re(wz(6)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)-Re(wz([5-6])*X*Y)); * id Re(wz(2)*X?)= 1/2*(Re(wz([2+3])*X)+Re(wz([2-3])*X)); id Re(wz(3)*X?)= 1/2*(Re(wz([2+3])*X)-Re(wz([2-3])*X)); id Re(wz(5)*X?)= 1/2*(Re(wz([5+6])*X)+Re(wz([5-6])*X)); id Re(wz(6)*X?)= 1/2*(Re(wz([5+6])*X)-Re(wz([5-6])*X)); * id Re(wz([2+3])*wzC([2-3]))= Re(wz([2-3])*wzC([2+3])); id Re(wz([2+3])*wzC([5-6]))= Re(wz([5-6])*wzC([2+3])); id Re(wz([2-3])*wzC([5+6]))= Re(wz([5+6])*wzC([2-3])); id Re(wz([5+6])*wzC([5-6]))= Re(wz([5-6])*wzC([5+6])); * id Re(wz([2+3])*wzC([2+3]))= wzM2(2)+wzM2(3)+2*Re(wz(2)*wzC(3)); id Re(wz([2-3])*wzC([2-3]))= wzM2(2)+wzM2(3)-2*Re(wz(2)*wzC(3)); id Re(wz([5+6])*wzC([5+6]))= wzM2(5)+wzM2(6)+2*Re(wz(5)*wzC(6)); id Re(wz([5-6])*wzC([5-6]))= wzM2(5)+wzM2(6)-2*Re(wz(5)*wzC(6)); *-- id X= [1+X]-1; id x1= 1-[1-x1]; id x2= 1-[1-x1]-[x1-x2]; id z= 1-[1-x1]-[x1-x2]-[x2-z]; id ups= [ups+X+x1]-[1+X]+[1-x1]; id Xp^-1= 1/[1+X]; ratio,X,[1+X],one; id one= 1; al one^-1= 1; * id Re(wz(X?)*wzC(Y?))= Rep(X,Y); id GsM2*GtM2= GstM2; id [1+X] = Xp; al [1+X]^-1= Xpi; id [1-x1]= omx1; al [x1-x2]= x12; al [x2-z]= x2z; al [ups+X+x1]= ux; al zX^-1= zXi; id Xpi^2= Xpi-aux; id aux= XX; * id Xp^2= Xps; id Xpi^2= Xpis; id y^3= yc; id y^2= ys; id PZ^-2= 1/PZs; * id omx1= O; id O^4= Oq; id O^3= Oc; id O^2= Os; id x12^3= x12c; id x12^2= x12s; id x2z^4= x2zq; id x2z^3= x2zc; id x2z^2= x2zs; id ux^2= uxs; id vv^2= vvs; id vv^3= vvc; id vvs*vv= vvc; id vvs^2= vvq; id vvc*vv= vvq; id csi(T)= 1/cst; id cst^-2= 1/cqt; id ve(T)= vet; al vd(T)= vdt; al vu(T)= vut; id ae(T)= aet; al ad(T)= adt; al au(T)= aut; * id vet^2= vets; al vdt^2= vdts; al vut^2= vuts; id aet^2= aets; al adt^2= adts; al aut^2= auts; * id Re(Gs*GtC*Pwt*PwsC)= rl1; id Re(Gs*GtC*Pwt)= rl2; id Re(Pws)= rl3; id Re(wz(n?)*Gs*PwsC)= rs1wz(n); id Re(wz(n?)*Gs)= rs2wz(n); id Re(wz(n?)*Gt*PwtC)= rt1wz(n); id Re(wz(n?)*Gt)= rt2wz(n); id Rep(n1?,n2?)= rwwz(n1,n2); id Re(wz(n1?)*wzC(n2?))= rwwz(n1,n2); id wzM2(n?)= rwwz(n,n); * id fun?([2+3],n?)= fun(11,n); al fun?([2-3],n?)= fun(12,n); al fun?([5+6],n?)= fun(13,n); al fun?([5-6],n?)= fun(14,n); * id fun?(n?,[2+3])= fun(n,11); al fun?(n?,[2-3])= fun(n,12); al fun?(n?,[5+6])= fun(n,13); al fun?(n?,[5-6])= fun(n,14); * id fun?([2+3])= fun(11); al fun?([2-3])= fun(12); al fun?([5+6])= fun(13); al fun?([5-6])= fun(14); * .store * *----g^2-------------------- * g ZZ10s= ZZ10; g ZZ20s= ZZ20; g ZZ30s= ZZ30; g ZZ40s= ZZ40; g ZZ11s= ZZ11; g ZZ21s= ZZ21; g ZZ31s= ZZ31; g ZZ41s= ZZ41; g ZZ12s= ZZ12; g ZZ22s= ZZ22; g ZZ32s= ZZ32; g ZZ42s= ZZ32; * id gexp^4= 0; id gexp^3= 0; id gexp^2= 1; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep,PZs, zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wz,wzC,wzM2,rs1wz,rs2wz, rt1wz,rt2wz,rwwz,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, zXi,XX,vv,vvs,vvc,vvq, cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts; * format fortran; print; * .store * *----g^3-------------------- * g ZZ10c= ZZ10; g ZZ20c= ZZ20; g ZZ30c= ZZ30; g ZZ40c= ZZ40; g ZZ11c= ZZ11; g ZZ21c= ZZ21; g ZZ31c= ZZ31; g ZZ41c= ZZ41; g ZZ12c= ZZ12; g ZZ22c= ZZ22; g ZZ32c= ZZ32; g ZZ42c= ZZ42; * id gexp^4= 0; id gexp^3= 1; id gexp^2= 0; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep,PZs, zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wz,wzC,wzM2,rs1wz,rs2wz, rt1wz,rt2wz,rwwz,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, zXi,XX,vv,vvs,vvc,vvq, cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts; format fortran; print; * .store * *----g^4-------------------- * g ZZ10q= ZZ10/GstM2; g ZZ20q= ZZ20/GstM2; g ZZ30q= ZZ30/GstM2; g ZZ40q= ZZ40/GstM2; g ZZ11q= ZZ11/GstM2; g ZZ21q= ZZ21/GstM2; g ZZ31q= ZZ31/GstM2; g ZZ41q= ZZ41/GstM2; g ZZ12q= ZZ12/GstM2; g ZZ22q= ZZ22/GstM2; g ZZ32q= ZZ32/GstM2; g ZZ42q= ZZ42/GstM2; * id gexp^4= 1; id gexp^3= 0; id gexp^2= 0; * b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Im,Mods,Rep,PZs, zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wz,wzC,wzM2,rs1wz,rs2wz, rt1wz,rt2wz,rwwz,rl1,rl2,rl3, Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis, zXi,XX,vv,vvs,vvc,vvq, cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts; format fortran; print; * #endprocedure * */////////////////////////////////////////////////////////////////////// * *#call wtensor{gg} #call totalgg{dummy} wwww_gg loaded Time = 2.22 sec Generated terms = 10539 T_gg Terms left = 4059 Bytes used = 190556 Time = 4.66 sec Generated terms = 20527 T_gg Terms left = 8174 Bytes used = 392880 Time = 5.10 sec Generated terms = 22188 T_gg Terms left = 9005 Bytes used = 435574 Time = 5.20 sec Generated terms = 22188 T_gg Terms in output = 7108 Bytes used = 338884 #call contractgg{dummy} Time = 6.72 sec Generated terms = 9622 Wd_gg Terms left = 557 Bytes used = 20428 Time = 7.41 sec Generated terms = 12635 Wd_gg Terms left = 558 Bytes used = 20494 Time = 7.41 sec Generated terms = 12635 Wd_gg Terms in output = 556 Bytes used = 20364 Time = 8.44 sec Generated terms = 7558 Wpp_gg Terms left = 481 Bytes used = 20722 Time = 9.28 sec Generated terms = 14041 Wpp_gg Terms left = 1086 Bytes used = 44114 Time = 10.43 sec Generated terms = 21107 Wpp_gg Terms left = 1710 Bytes used = 68192 Time = 11.63 sec Generated terms = 27881 Wpp_gg Terms left = 1831 Bytes used = 73486 Time = 12.10 sec Generated terms = 30228 Wpp_gg Terms left = 1913 Bytes used = 77152 Time = 12.11 sec Generated terms = 30228 Wpp_gg Terms in output = 1359 Bytes used = 55264 Time = 12.14 sec Generated terms = 556 Wd_gg Terms in output = 556 Bytes used = 20364 Time = 12.22 sec Generated terms = 1359 Wpp_gg Terms in output = 1359 Bytes used = 55264 Time = 12.25 sec Generated terms = 556 Wd_gg Terms in output = 556 Bytes used = 20364 Time = 12.33 sec Generated terms = 1359 Wpp_gg Terms in output = 1359 Bytes used = 55264 Time = 12.36 sec Generated terms = 556 Wd_gg Terms in output = 556 Bytes used = 20364 Time = 12.44 sec Generated terms = 1359 Wpp_gg Terms in output = 1359 Bytes used = 55264 Time = 12.48 sec Generated terms = 556 Wd_gg Terms in output = 556 Bytes used = 20364 Time = 12.56 sec Generated terms = 1359 Wpp_gg Terms in output = 1359 Bytes used = 55264 Time = 12.59 sec Generated terms = 556 Wd_gg Terms in output = 556 Bytes used = 20364 Time = 12.66 sec Generated terms = 1359 Wpp_gg Terms in output = 1359 Bytes used = 55264 Time = 12.77 sec Generated terms = 1147 Wd_gg Terms in output = 542 Bytes used = 21940 Time = 12.99 sec Generated terms = 2874 Wpp_gg Terms in output = 1315 Bytes used = 48362 Time = 16.88 sec Generated terms = 12640 GG10 Terms left = 373 Bytes used = 10612 Time = 20.77 sec Generated terms = 23938 GG10 Terms left = 618 Bytes used = 17810 Time = 24.49 sec Generated terms = 34911 GG10 Terms left = 869 Bytes used = 25060 Time = 24.99 sec Generated terms = 36944 GG10 Terms left = 972 Bytes used = 28110 Time = 25.01 sec Generated terms = 36944 GG10 Terms in output = 558 Bytes used = 16050 Time = 28.86 sec Generated terms = 12640 GG20 Terms left = 374 Bytes used = 10634 Time = 32.71 sec Generated terms = 23938 GG20 Terms left = 619 Bytes used = 17832 Time = 36.42 sec Generated terms = 34911 GG20 Terms left = 870 Bytes used = 25082 Time = 36.91 sec Generated terms = 36944 GG20 Terms left = 973 Bytes used = 28132 Time = 36.91 sec Generated terms = 36944 GG20 Terms in output = 561 Bytes used = 16108 Time = 40.48 sec Generated terms = 12316 GG11 Terms left = 429 Bytes used = 12336 Time = 44.38 sec Generated terms = 23263 GG11 Terms left = 698 Bytes used = 20102 Time = 48.26 sec Generated terms = 34015 GG11 Terms left = 1013 Bytes used = 28992 Time = 52.03 sec Generated terms = 44401 GG11 Terms left = 1276 Bytes used = 36300 Time = 55.76 sec Generated terms = 55323 GG11 Terms left = 1538 Bytes used = 43690 Time = 57.66 sec Generated terms = 61353 GG11 Terms left = 1693 Bytes used = 48102 Time = 57.67 sec Generated terms = 61353 GG11 Terms in output = 555 Bytes used = 15938 Time = 61.24 sec Generated terms = 12316 GG21 Terms left = 433 Bytes used = 12452 Time = 65.11 sec Generated terms = 23263 GG21 Terms left = 702 Bytes used = 20218 Time = 68.97 sec Generated terms = 34015 GG21 Terms left = 1017 Bytes used = 29108 Time = 72.72 sec Generated terms = 44401 GG21 Terms left = 1280 Bytes used = 36416 Time = 76.43 sec Generated terms = 55323 GG21 Terms left = 1542 Bytes used = 43806 Time = 78.32 sec Generated terms = 61353 GG21 Terms left = 1697 Bytes used = 48218 Time = 78.33 sec Generated terms = 61353 GG21 Terms in output = 566 Bytes used = 16280 Time = 79.17 sec Generated terms = 7567 GG12 Terms in output = 38 Bytes used = 1202 Time = 80.01 sec Generated terms = 7567 GG22 Terms in output = 38 Bytes used = 1202 Time = 80.04 sec Generated terms = 102 GG10s Terms in output = 102 Bytes used = 1996 Time = 80.06 sec Generated terms = 103 GG20s Terms in output = 103 Bytes used = 2010 Time = 80.08 sec Generated terms = 108 GG11s Terms in output = 108 Bytes used = 2220 Time = 80.10 sec Generated terms = 109 GG21s Terms in output = 109 Bytes used = 2234 Time = 80.12 sec Generated terms = 38 GG12s Terms in output = 38 Bytes used = 906 Time = 80.12 sec Generated terms = 38 GG22s Terms in output = 38 Bytes used = 906 GG10s = + + XX*rl3*vvc*GtM2*ys * ( 2*O*x12*x2z - 2./3.*O*x12*ux - 4./3.*O + *x2z*ux - 2./3.*O*x2z + 2./3.*O*x2zs + 2./3.*O*uxs - 4./3.*Os + *x12 - 4./3.*Os*x2z + 4./3.*Os*ux + 2./3.*Oc - 2./3.*x12*x2z* + ux + 4./3.*x12*x2z - 2./3.*x12*x2zs + 2./3.*x12*uxs - 2./3.* + x2z*ux + 4./3.*x2z*x12s - 2*ux*x12s + 2*x12c + 4./3.*x2zs ) GG10s = GG10s + XX*vvs*GtM2*PwsM2*y*zXi * ( 8./3.*O*x12*ux - 4./3. + *O*x12 - 4./9.*O*ux - 4./9.*O*uxs + 40./9.*Os*x12 - 16./9.*Os + *ux - 4./3.*Os - 16./9.*Oc - 4./9.*x12*ux + 4./9.*x12*uxs ) GG10s = GG10s + XX*vvs*GtM2*PwsM2*y * ( 34./9.*O*x12 + 4./9.*O* + x2z - 8./9.*O*ux - 4./3.*O - 4./3.*Os - 8./9.*x12*x2z + 10./9. + *x12*ux + 8./9.*x2z - 4./9.*ux - 2./3.*x12s ) GG10s = GG10s + Xpi*rl3*vvc*GtM2*ys * ( 5./3.*O*x12 + 1./3.*O*x2z + - 1./3.*O*ux - 2./3.*O - 1./3.*Os - 4./3.*x12*x2z + 1./3.* + x12*ux + 4./3.*x12 + 7./3.*x2z - 2./3.*ux ) GG10s = GG10s + Xpi*vvs*GtM2*PwsM2*y*zXi * ( 40./9.*O*x12 - 4./3. + *O*ux - 20./9.*O - 8./3.*Os + 4./3.*x12*ux - 8./9.*x12 - 4./9. + *ux ) GG10s = GG10s + Xpi*vvs*GtM2*PwsM2*y * ( - 4./9.*O ) GG10s = GG10s + rl1*vvs*y*Xp * ( 1./3. - 1./3.*O - x12 ) GG10s = GG10s + rl1*vvs*y * ( 1./3.*O*x2z - 1./3.*O*ux - 1./3.*O + + 2./3.*x12*x2z - 1./3.*x12*ux + 1./3.*x12 - 1./3.*x2z + 2./ + 3.*x12s ) GG10s = GG10s + rl3*vvc*GtM2*Xps*ys * ( - 2./3. ) GG10s = GG10s + rl3*vvc*GtM2*ys*Xp * ( 2 - 5./3.*O - 1./3.*x12 + + 4./3.*x2z - 2./3.*ux ) GG10s = GG10s + rl3*vvc*GtM2*ys * ( - 1 - 2./3.*O*x12 + O*x2z + + 1./3.*O*ux + 1./3.*O + 2./3.*Os + 2./3.*x12*x2z - x12 + 2./3. + *x2z*ux - 10./3.*x2z + 2./3.*ux - 2./3.*x2zs ) GG10s = GG10s + vvs*GtM2*PwsM2*y*Xp*zXi * ( - 2./3. - 14./9.*O + + 2./3.*x12 ) GG10s = GG10s + vvs*GtM2*PwsM2*y*Xp * ( 1./9. ) GG10s = GG10s + vvs*GtM2*PwsM2*y*zXi * ( 8./9. - 34./9.*O*x12 + + 14./9.*O*ux + 28./9.*O + 10./3.*Os - 10./9.*x12*ux - 4./9.* + x12 + 4./9.*ux ) GG10s = GG10s + vvs*GtM2*PwsM2*y * ( - 2./9. + O + 1./3.*x12 - 1. + /9.*x2z ) GG20s = + + XX*rl3*vvc*GtM2*ys * ( 6*O*x12*x2z - 2*O*x12*ux - 4*O*x2z*ux + - 2*O*x2z + 2*O*x2zs + 2*O*uxs - 4*Os*x12 - 4*Os*x2z + 4*Os* + ux + 2*Oc - 2*x12*x2z*ux + 4*x12*x2z - 2*x12*x2zs + 2*x12*uxs + - 2*x2z*ux + 4*x2z*x12s - 6*ux*x12s + 6*x12c + 4*x2zs ) GG20s = GG20s + XX*vvs*GtM2*PwsM2*y*zXi * ( 8*O*x12*ux - 4*O*x12 + - 4./3.*O*ux - 4./3.*O*uxs + 40./3.*Os*x12 - 16./3.*Os*ux - + 4*Os - 16./3.*Oc - 4./3.*x12*ux + 4./3.*x12*uxs ) GG20s = GG20s + XX*vvs*GtM2*PwsM2*y * ( 34./3.*O*x12 + 4./3.*O* + x2z - 8./3.*O*ux - 4*O - 4*Os - 8./3.*x12*x2z + 10./3.*x12*ux + + 8./3.*x2z - 4./3.*ux - 2*x12s ) GG20s = GG20s + Xpi*rl3*vvc*GtM2*ys * ( 5*O*x12 + O*x2z - O*ux - + 2*O - Os - 4*x12*x2z + x12*ux + 4*x12 + 7*x2z - 2*ux ) GG20s = GG20s + Xpi*vvs*GtM2*PwsM2*y*zXi * ( 40./3.*O*x12 - 4*O* + ux - 20./3.*O - 8*Os + 4*x12*ux - 8./3.*x12 - 4./3.*ux ) GG20s = GG20s + Xpi*vvs*GtM2*PwsM2*y * ( - 4./3.*O ) GG20s = GG20s + rl1*vvs*y*Xp * ( 1./3. - 1./3.*O - x12 ) GG20s = GG20s + rl1*vvs*y * ( 1./3.*O*x2z - 1./3.*O*ux - 1./3.*O + + 2./3.*x12*x2z - 1./3.*x12*ux + 1./3.*x12 - 1./3.*x2z + 2./ + 3.*x12s ) GG20s = GG20s + rl3*vvc*GtM2*Xps*ys * ( - 2./3. ) GG20s = GG20s + rl3*vvc*GtM2*ys*Xp * ( 4 - 5./3.*O - 5./3.*x12 + + 4./3.*x2z - 2./3.*ux ) GG20s = GG20s + rl3*vvc*GtM2*ys * ( - 3 - 4*O*x12 + 1./3.*O*x2z + + O*ux + 5./3.*O + 4./3.*Os + 10./3.*x12*x2z - 2./3.*x12*ux + - 7./3.*x12 + 2./3.*x2z*ux - 8*x2z + 2*ux - 2./3.*x2zs ) GG20s = GG20s + vvs*GtM2*PwsM2*y*Xp*zXi * ( - 22./9. - 10./3.*O + + 22./9.*x12 ) GG20s = GG20s + vvs*GtM2*PwsM2*y*Xp * ( 1./9. ) GG20s = GG20s + vvs*GtM2*PwsM2*y*zXi * ( 8./3. - 38./3.*O*x12 + + 38./9.*O*ux + 28./3.*O + 26./3.*Os - 34./9.*x12*ux - 4./9.* + x12 + 4./3.*ux ) GG20s = GG20s + vvs*GtM2*PwsM2*y * ( - 2./9. + 17./9.*O + 1./3.* + x12 - 1./9.*x2z ) GG11s = + + XX*rl3*vvc*GtM2*ys * ( - 2./3.*O*x12*x2zs - 4./3.*O*x2z*x12s + - 2./3.*O*x12c - 2./3.*x12*x12c + 2./3.*x12*x2zc - 4./3.*x12 + *x2zs - 2./3.*x2z*x12s - 2./3.*x2z*x12c + 2./3.*x12s*x2zs - 2. + /3.*x2zc ) GG11s = GG11s + XX*vvs*GtM2*PwsM2*y*zXi * ( - 4./9.*Os*x12 + 8./ + 9.*Oc*x12 - 4./9.*Oc ) GG11s = GG11s + XX*vvs*GtM2*PwsM2*y * ( - 8./9.*O*x12*x2z + 4./9. + *O*x2z - 8./9.*O*x12s + 8./9.*Os*x12 - 4./9.*Os - 4./9.*x12* + x2z + 2./9.*x12*x2zs - 2./9.*x12c - 4./9.*x2zs ) GG11s = GG11s + Xpi*rl3*vvc*GtM2*ys * ( - 5./3.*O*x12*x2z - 5./3. + *O*x12s - 3*x12*x2z + 5./3.*x12*x2zs + 2./3.*x2z*x12s - 2./3. + *x12s - x12c - 7./3.*x2zs ) GG11s = GG11s + Xpi*vvs*GtM2*PwsM2*y*zXi * ( - 4./9.*O*x12 + 8./ + 9.*Os*x12 - 8./9.*Os ) GG11s = GG11s + Xpi*vvs*GtM2*PwsM2*y * ( - 8./9.*O*x12 + 4./9.* + x12*x2z - 4./9.*x12 - 8./9.*x2z ) GG11s = GG11s + rl1*vvs*Xps*y * ( 2./3. - 2./3.*O - 4./3.*x12 ) GG11s = GG11s + rl1*vvs*y*Xp * ( - 1./3. + 2./3.*O*x2z + 1./3.*O + + 4./3.*x12*x2z + 1./3.*x12 - x2z + 4./3.*x12s ) GG11s = GG11s + rl1*vvs*y * ( 1./3.*O*x12 - 1./3.*O*x2z - 1./3.* + Os*x12 - 1./3.*x12*x2zs - 2./3.*x2z*x12s + 1./3.*x2z - 1./3.* + x12s - 1./3.*x12c + 1./3.*x2zs ) GG11s = GG11s + rl3*vvc*GtM2*Xps*ys*Xp * ( - 2./3. ) GG11s = GG11s + rl3*vvc*GtM2*Xps*ys * ( 7./3. - 4./3.*O - 2./3.* + x12 + 4./3.*x2z ) GG11s = GG11s + rl3*vvc*GtM2*ys*Xp * ( - 4./3. + 1./3.*O*x12 + 2 + *O*x2z + 1./3.*O + 5./3.*x12*x2z - 2*x12 - 5*x2z + x12s - 2./ + 3.*x2zs ) GG11s = GG11s + rl3*vvc*GtM2*ys * ( 1./3.*O*x12*x2z + 4./3.*O*x12 + - 1./3.*O*x2z + O*x12s - 2./3.*O*x2zs - 1./3.*Os*x12 + 8./3. + *x12*x2z - 2*x12*x2zs + 5./3.*x12 - 4./3.*x2z*x12s + 3*x2z + + 2./3.*x12s + 2./3.*x12c + 10./3.*x2zs ) GG11s = GG11s + vvs*GtM2*PwsM2*Xps*y * ( 2./9. ) GG11s = GG11s + vvs*GtM2*PwsM2*y*Xp*zXi * ( - 2./9. + 4./9.*O* + x12 - 4./9.*Os ) GG11s = GG11s + vvs*GtM2*PwsM2*y*Xp * ( - 7./9. + 1./3.*O + 1./3. + *x12 - 1./3.*x2z ) GG11s = GG11s + vvs*GtM2*PwsM2*y*zXi * ( - 4./9.*O*x12 + 2./3.*O + - 4./9.*Os*x12 + 4./9.*Os + 4./9.*Oc + 2./9.*x12 ) GG11s = GG11s + vvs*GtM2*PwsM2*y * ( 2./3. + 7./9.*O*x12 - 1./3.* + O*x2z - 4./9.*O + 4./9.*Os - 4./9.*x12*x2z + 1./9.*x12 + 11./ + 9.*x2z - 1./9.*x12s + 1./9.*x2zs ) GG21s = + + XX*rl3*vvc*GtM2*ys * ( - 2*O*x12*x2zs - 4*O*x2z*x12s - 2*O* + x12c - 2*x12*x12c + 2*x12*x2zc - 4*x12*x2zs - 2*x2z*x12s - 2* + x2z*x12c + 2*x12s*x2zs - 2*x2zc ) GG21s = GG21s + XX*vvs*GtM2*PwsM2*y*zXi * ( - 4./3.*Os*x12 + 8./ + 3.*Oc*x12 - 4./3.*Oc ) GG21s = GG21s + XX*vvs*GtM2*PwsM2*y * ( - 8./3.*O*x12*x2z + 4./3. + *O*x2z - 8./3.*O*x12s + 8./3.*Os*x12 - 4./3.*Os - 4./3.*x12* + x2z + 2./3.*x12*x2zs - 2./3.*x12c - 4./3.*x2zs ) GG21s = GG21s + Xpi*rl3*vvc*GtM2*ys * ( - 5*O*x12*x2z - 5*O*x12s + - 9*x12*x2z + 5*x12*x2zs + 2*x2z*x12s - 2*x12s - 3*x12c - 7* + x2zs ) GG21s = GG21s + Xpi*vvs*GtM2*PwsM2*y*zXi * ( - 4./3.*O*x12 + 8./ + 3.*Os*x12 - 8./3.*Os ) GG21s = GG21s + Xpi*vvs*GtM2*PwsM2*y * ( - 8./3.*O*x12 + 4./3.* + x12*x2z - 4./3.*x12 - 8./3.*x2z ) GG21s = GG21s + rl1*vvs*Xps*y * ( 2./3. - 2./3.*O - 4./3.*x12 ) GG21s = GG21s + rl1*vvs*y*Xp * ( - 1./3. + 2./3.*O*x2z + 1./3.*O + + 4./3.*x12*x2z + 1./3.*x12 - x2z + 4./3.*x12s ) GG21s = GG21s + rl1*vvs*y * ( 1./3.*O*x12 - 1./3.*O*x2z - 1./3.* + Os*x12 - 1./3.*x12*x2zs - 2./3.*x2z*x12s + 1./3.*x2z - 1./3.* + x12s - 1./3.*x12c + 1./3.*x2zs ) GG21s = GG21s + rl3*vvc*GtM2*Xps*ys*Xp * ( - 2./3. ) GG21s = GG21s + rl3*vvc*GtM2*Xps*ys * ( 13./3. - 4./3.*O - 2*x12 + + 4./3.*x2z ) GG21s = GG21s + rl3*vvc*GtM2*ys*Xp * ( - 10./3. - 5./3.*O*x12 + + 2*O*x2z + 1./3.*O + 5*x12*x2z - 4*x12 - 31./3.*x2z + 5./3.* + x12s - 2./3.*x2zs ) GG21s = GG21s + rl3*vvc*GtM2*ys * ( 11./3.*O*x12*x2z + 10./3.*O* + x12 - 1./3.*O*x2z + 13./3.*O*x12s - 2./3.*O*x2zs - 1./3.*Os* + x12 + 16./3.*x12*x2z - 16./3.*x12*x2zs + 5*x12 - 8./3.*x2z* + x12s + 25./3.*x2z + 4./3.*x12s + 8./3.*x12c + 8*x2zs ) GG21s = GG21s + vvs*GtM2*PwsM2*Xps*y * ( 2./9. ) GG21s = GG21s + vvs*GtM2*PwsM2*y*Xp*zXi * ( - 2./9. + 4./9.*O* + x12 - 8./9.*O - 4./9.*Os ) GG21s = GG21s + vvs*GtM2*PwsM2*y*Xp * ( - 5./3. + 1./3.*O + 7./9. + *x12 - 1./3.*x2z ) GG21s = GG21s + vvs*GtM2*PwsM2*y*zXi * ( 4./9.*O*x12 + 14./9.*O + - 20./9.*Os*x12 + 20./9.*Os + 4./9.*Oc + 2./9.*x12 ) GG21s = GG21s + vvs*GtM2*PwsM2*y * ( 14./9. + 23./9.*O*x12 - 1./3. + *O*x2z - 4./9.*O + 4./9.*Os - 4./3.*x12*x2z + 5./9.*x12 + 3* + x2z - 1./9.*x12s + 1./9.*x2zs ) GG12s = + + XX*vvs*dms*GtM2*PwsM2*s**(-1) * ( - 2./9.*x12*x2zs - 8./9.* + x2z*x12s - 2./3.*x12c ) GG12s = GG12s + XX*vvs*GtM2*PwsM2*y * ( 2./9.*x12*x12c + 4./9.* + x2z*x12c + 2./9.*x12s*x2zs ) GG12s = GG12s + Xpi*vvs*dms*GtM2*PwsM2*s**(-1) * ( - 8./9.*x12* + x2z - 4./3.*x12s ) GG12s = GG12s + Xpi*vvs*GtM2*PwsM2*y * ( 2./9.*x12*x2zs + 8./9.* + x2z*x12s + 2./3.*x12c ) GG12s = GG12s + vvs*dms*GtM2*PwsM2*s**(-1)*Xps * ( 2./9. ) GG12s = GG12s + vvs*dms*GtM2*PwsM2*s**(-1)*Xp * ( - 1./9. - 1./9. + *O - x12 - 1./3.*x2z ) GG12s = GG12s + vvs*dms*GtM2*PwsM2*s**(-1) * ( 1./9.*O*x12 + 1./9. + *O*x2z + 10./9.*x12*x2z + 7./9.*x12 + 1./9.*x2z + 13./9.*x12s + + 1./9.*x2zs ) GG12s = GG12s + vvs*GtM2*PwsM2*Xps*y*Xp * ( 1./9. ) GG12s = GG12s + vvs*GtM2*PwsM2*Xps*y * ( - 1./9. - 4./9.*x12 - 2. + /9.*x2z ) GG12s = GG12s + vvs*GtM2*PwsM2*y*Xp * ( 2./3.*x12*x2z + 4./9.*x12 + + 2./9.*x2z + 7./9.*x12s + 1./9.*x2zs ) GG12s = GG12s + vvs*GtM2*PwsM2*y * ( - 2./3.*x12*x2z - 2./9.*x12 + *x2zs - 8./9.*x2z*x12s - 7./9.*x12s - 2./3.*x12c - 1./9.*x2zs + ) GG22s = + + XX*vvs*dms*GtM2*PwsM2*s**(-1) * ( - 2./3.*x12*x2zs - 8./3.* + x2z*x12s - 2*x12c ) GG22s = GG22s + XX*vvs*GtM2*PwsM2*y * ( 2./3.*x12*x12c + 4./3.* + x2z*x12c + 2./3.*x12s*x2zs ) GG22s = GG22s + Xpi*vvs*dms*GtM2*PwsM2*s**(-1) * ( - 8./3.*x12* + x2z - 4*x12s ) GG22s = GG22s + Xpi*vvs*GtM2*PwsM2*y * ( 2./3.*x12*x2zs + 8./3.* + x2z*x12s + 2*x12c ) GG22s = GG22s + vvs*dms*GtM2*PwsM2*s**(-1)*Xps * ( 2./9. ) GG22s = GG22s + vvs*dms*GtM2*PwsM2*s**(-1)*Xp * ( - 1./9. - 1./9. + *O - 7./3.*x12 - 1./3.*x2z ) GG22s = GG22s + vvs*dms*GtM2*PwsM2*s**(-1) * ( 1./9.*O*x12 + 1./9. + *O*x2z + 26./9.*x12*x2z + 19./9.*x12 + 1./9.*x2z + 37./9.* + x12s + 1./9.*x2zs ) GG22s = GG22s + vvs*GtM2*PwsM2*Xps*y*Xp * ( 1./9. ) GG22s = GG22s + vvs*GtM2*PwsM2*Xps*y * ( - 1./9. - 8./9.*x12 - 2. + /9.*x2z ) GG22s = GG22s + vvs*GtM2*PwsM2*y*Xp * ( 14./9.*x12*x2z + 8./9.* + x12 + 2./9.*x2z + 19./9.*x12s + 1./9.*x2zs ) GG22s = GG22s + vvs*GtM2*PwsM2*y * ( - 14./9.*x12*x2z - 2./3.* + x12*x2zs - 8./3.*x2z*x12s - 19./9.*x12s - 2*x12c - 1./9.*x2zs + ) Time = 80.18 sec Generated terms = 456 GG10c Terms in output = 456 Bytes used = 8536 Time = 80.21 sec Generated terms = 458 GG20c Terms in output = 458 Bytes used = 8564 Time = 80.24 sec Generated terms = 447 GG11c Terms in output = 447 Bytes used = 8670 Time = 80.27 sec Generated terms = 457 GG21c Terms in output = 457 Bytes used = 8850 Time = 80.27 sec Generated terms = 0 GG12c Terms in output = 0 Bytes used = 2 Time = 80.27 sec Generated terms = 0 GG22c Terms in output = 0 Bytes used = 2 GG10c = + + rs1wg(1)*XX*vvc*GtM2*ys * ( 4*O*x12*x2z - 4./3.*O*x12*ux - 8./ + 3.*O*x2z*ux - 4./3.*O*x2z + 4./3.*O*x2zs + 4./3.*O*uxs - 8./3. + *Os*x12 - 8./3.*Os*x2z + 8./3.*Os*ux + 4./3.*Oc - 4./3.*x12* + x2z*ux + 8./3.*x12*x2z - 4./3.*x12*x2zs + 4./3.*x12*uxs - 4./ + 3.*x2z*ux + 8./3.*x2z*x12s - 4*ux*x12s + 4*x12c + 8./3.*x2zs + ) GG10c = GG10c + rs1wg(1)*Xpi*vvc*GtM2*ys * ( 16./3.*O*x12 + 4./3. + *O*x2z - 4./3.*O*ux - 4./3.*O - 4./3.*Os - 4*x12*x2z + 4./3.* + x12*ux + 8./3.*x12 + 16./3.*x2z - 4./3.*ux - 4./3.*x12s ) GG10c = GG10c + rs1wg(1)*vvc*GtM2*Xps*ys * ( - 4./3. ) GG10c = GG10c + rs1wg(1)*vvc*GtM2*ys*Xp * ( 14./3. - 10./3.*O - 2 + *x12 + 8./3.*x2z - 4./3.*ux ) GG10c = GG10c + rs1wg(1)*vvc*GtM2*ys * ( - 8./3. - 4*O*x12 + 2./ + 3.*O*x2z + 2*O*ux + 2./3.*O + 8./3.*Os + 8./3.*x12*x2z - 2./3. + *x12*ux - 2./3.*x12 + 4./3.*x2z*ux - 22./3.*x2z + 4./3.*ux + + 4./3.*x12s - 4./3.*x2zs ) GG10c = GG10c + rs1wg(4)*XX*vvc*GtM2*ys * ( 2./3.*O*x12*x2z*ux - + 1./3.*O*x12*x2zs - 2./3.*O*x2z*x12s - 1./3.*O*x2z*uxs + 1./3. + *O*ux*x2zs - 1./3.*O*x2zs + 1./3.*Os*x12*x2z - 1./3.*Os*x2z* + ux + 1./3.*Os*x2zs - 1./3.*x12*x2z*uxs + 2./3.*x12*ux*x2zs - + 1./3.*x12*x2zc - 1./3.*x12*x2zs + x2z*ux*x12s - x2z*x12c + 1./ + 3.*ux*x2zs - 4./3.*x12s*x2zs - 1./3.*x2zc ) GG10c = GG10c + rs1wg(4)*Xpi*vvc*GtM2*ys * ( - 2./3.*O*x12*x2z + + 2./3.*O*x12*ux + 1./3.*O*x2z*ux - 2./3.*O*x2z - 2./3.*O* + x12s - 1./3.*O*uxs + 1./3.*Os*x12 + 2./3.*Os*x2z - 1./3.*Os* + ux + 2./3.*x12*x2z*ux - 2./3.*x12*x2z - 2./3.*x12*x2zs - 1./3. + *x12*uxs + 2./3.*x2z*ux - 2*x2z*x12s + ux*x12s - x12c - x2zs + ) GG10c = GG10c + rs1wg(4)*vvc*GtM2*Xps*ys * ( 2./3. - 1./3.*O - 1./ + 3.*x12 ) GG10c = GG10c + rs1wg(4)*vvc*GtM2*ys*Xp * ( - 1./3. + 1./6.*O* + x12 + 1./2.*O*x2z - 1./6.*O*ux - 1./2.*O + 1./2.*Os + 1./6.* + x12*x2z - 1./6.*x12*ux - 2./3.*x12 - 5./3.*x2z + 1./2.*ux - 1. + /3.*x12s ) GG10c = GG10c + rs1wg(4)*vvc*GtM2*ys * ( 1./3.*O*x12*x2z - 5./6.* + O*x12*ux + 1./3.*O*x12 - 1./3.*O*x2z*ux + 5./6.*O*x2z + 2./3. + *O*x12s - 1./6.*O*x2zs + 1./2.*O*uxs + 1./3.*O - 1./2.*Os*x12 + - 5./6.*Os*x2z + 1./2.*Os*ux - 1./3.*Os - 5./6.*x12*x2z*ux + + 4./3.*x12*x2z + 1./2.*x12*x2zs + 1./2.*x12*uxs + 1./3.*x12 + - 5./6.*x2z*ux + 2*x2z*x12s + x2z - 4./3.*ux*x12s - 1./3.*ux + + 2./3.*x12s + 7./6.*x12c + 4./3.*x2zs ) GG10c = GG10c + rs1wg(7)*XX*vvc*GtM2*ys * ( 1./6.*O*x12*x2zs - 1./ + 6.*O*x2z*x12s + 1./6.*O*x2z*uxs - 1./3.*O*ux*x2zs + 1./6.*O* + x2zc - 1./3.*O*x2zs - 1./6.*Os*x12*x2z + 1./3.*Os*x2z*ux - 1./ + 6.*Os*x2zs + 1./6.*Oc*x2z + 1./6.*x12*x2z*uxs - 1./6.*x12*ux* + x2zs - 1./6.*x12*x2zc + 2./3.*x12*x2zs - 1./3.*x2z*ux*x12s + + 1./6.*x2z*x12c - 1./6.*ux*x2zs + 1./3.*x2zc ) GG10c = GG10c + rs1wg(7)*Xpi*vvc*GtM2*ys * ( - 1./6.*O*x12*x2z + - 1./3.*O*x2z + 1./6.*Os*x2z + 2./3.*x12*x2z - 1./6.*x12* + x2zs - 1./6.*x2z*ux - 1./3.*x2z*x12s + 1./2.*x2zs ) GG10c = GG10c + rs1wg(7)*vvc*GtM2*Xps*ys * ( - 1./6. + 1./6.*O + + 1./6.*x12 ) GG10c = GG10c + rs1wg(7)*vvc*GtM2*ys*Xp * ( - 1./2.*O*x12 - 1./2. + *O*x2z + 1./6.*O*ux - 1./6.*O + 1./12.*Os - 1./2.*x12*x2z + 1. + /6.*x12*ux + 1./2.*x12 + 2./3.*x2z - 1./6.*ux - 7./12.*x12s ) GG10c = GG10c + rs1wg(7)*vvc*GtM2*ys * ( 2./3.*O*x12*x2z - 1./6.* + O*x12*ux - 1./6.*O*x2z*ux + 1./2.*O*x2z + 1./4.*O*x12s + 1./3. + *O*x2zs - 1./6.*Os*x2z - 1./12.*Os*ux - 1./12.*Oc - 1./6.*x12 + *x2z*ux - 7./6.*x12*x2z + 1./2.*x12*x2zs + 1./3.*x2z*ux + 5./ + 6.*x2z*x12s - 1./6.*x2z - 1./12.*ux*x12s + 1./6.*x12c - 5./6. + *x2zs ) GG10c = GG10c + rs1wg(8)*Xpi*vvc*GtM2*ys * ( 1./3.*O*x12*x2z - 2./ + 3.*O*x12*ux - 1./3.*O*x2z*ux + 1./3.*O*x2z + 2./3.*O*x12s + 1. + /3.*O*uxs - 1./3.*Os*x12 - 1./3.*Os*x2z + 1./3.*Os*ux - 2./3. + *x12*x2z*ux + 1./3.*x12*x2z + 1./3.*x12*x2zs + 1./3.*x12*uxs + - 1./3.*x2z*ux + 4./3.*x2z*x12s - ux*x12s + x12c + 1./3.* + x2zs ) GG10c = GG10c + rs1wg(8)*vvc*GtM2*Xps*ys * ( - 1./3. ) GG10c = GG10c + rs1wg(8)*vvc*GtM2*ys*Xp * ( 1./3. + 1./6.*O*x12 + + 1./6.*O*x2z - 1./6.*O*ux + 2./3.*O - 2./3.*Os + 1./2.*x12* + x2z - 1./6.*x12*ux + 1./6.*x12 + 1./2.*x2z - 1./6.*ux + 5./6. + *x12s ) GG10c = GG10c + rs1wg(8)*vvc*GtM2*ys * ( - 1./3.*O*x12*x2z + 5./ + 6.*O*x12*ux - 1./3.*O*x12 + 1./2.*O*x2z*ux - 2./3.*O*x2z - 2./ + 3.*O*x12s - 1./2.*O*uxs - 1./3.*O + 1./2.*Os*x12 + 1./2.*Os* + x2z - 1./2.*Os*ux + 1./3.*Os + x12*x2z*ux - 1./2.*x12*x2z - 1. + /2.*x12*x2zs - 1./2.*x12*uxs - 1./3.*x12 + 1./6.*x2z*ux - 5./ + 3.*x2z*x12s - 2./3.*x2z + 4./3.*ux*x12s + 1./3.*ux - 2./3.* + x12s - 7./6.*x12c - 1./6.*x2zs ) GG10c = GG10c + rs1wg(9)*Xpi*vvc*GtM2*ys * ( - 2*O*x12 - 2./3.*O + *x2z + 2./3.*O*ux + 2./3.*Os + 4./3.*x12*x2z - 2./3.*x12*ux + - 2./3.*x2z + 4./3.*x12s ) GG10c = GG10c + rs1wg(9)*vvc*GtM2*ys*Xp * ( - 2./3. + 4./3.*x12 + ) GG10c = GG10c + rs1wg(9)*vvc*GtM2*ys * ( 2./3. + 8./3.*O*x12 + 4./ + 3.*O*x2z - 4./3.*O*ux - 4./3.*Os - 4./3.*x12*x2z + 2./3.*x12* + ux - 4./3.*x12 + 2./3.*x2z - 4./3.*x12s ) GG10c = GG10c + rs1wg(10)*XX*vvc*GtM2*ys * ( - 1./6.*O*x12*x2zs + + 1./6.*O*x2z*x12s - 1./6.*O*x2z*uxs + 1./3.*O*ux*x2zs - 1./ + 6.*O*x2zc + 1./3.*O*x2zs + 1./6.*Os*x12*x2z - 1./3.*Os*x2z*ux + + 1./6.*Os*x2zs - 1./6.*Oc*x2z - 1./6.*x12*x2z*uxs + 1./6.* + x12*ux*x2zs + 1./6.*x12*x2zc - 2./3.*x12*x2zs + 1./3.*x2z*ux* + x12s - 1./6.*x2z*x12c + 1./6.*ux*x2zs - 1./3.*x2zc ) GG10c = GG10c + rs1wg(10)*Xpi*vvc*GtM2*ys * ( 1./3.*O*x2z*ux + 2./ + 3.*O*x2z + 1./6.*O*x12s - 1./6.*O*x2zs - 1./6.*O*uxs + 1./6.* + Os*x12 - 1./3.*Os*ux - 1./6.*Oc + 1./6.*x12*x2z*ux - 4./3.* + x12*x2z + 1./3.*x12*x2zs - 1./6.*x12*uxs + 1./3.*x2z*ux + 1./ + 3.*x2z*x12s + 1./3.*ux*x12s - 1./6.*x12c - 5./6.*x2zs ) GG10c = GG10c + rs1wg(10)*vvc*GtM2*Xps*ys * ( 1./3. - 1./6.*O - 1. + /6.*x12 ) GG10c = GG10c + rs1wg(10)*vvc*GtM2*ys*Xp * ( - 1./6. + 7./12.*O* + x12 + 5./12.*O*x2z - 1./12.*O*ux + 1./2.*O - 1./4.*Os + 7./12. + *x12*x2z - 1./12.*x12*ux - 13./12.*x12 - 13./12.*x2z + 1./4.* + ux + 5./6.*x12s ) GG10c = GG10c + rs1wg(10)*vvc*GtM2*ys * ( - 1./3.*O*x12*x2z + 1./ + 12.*O*x12*ux - 1./6.*O*x12 - 5./12.*O*x2z*ux - 5./6.*O*x2z - + 5./12.*O*x12s + 1./4.*O*uxs - 1./3.*O - 1./4.*Os*x12 - 1./12. + *Os*x2z + 1./2.*Os*ux + 1./6.*Os + 1./4.*Oc - 1./6.*x12*x2z* + ux + 19./12.*x12*x2z - 7./12.*x12*x2zs + 1./4.*x12*uxs + 2./3. + *x12 - 5./12.*x2z*ux - 2./3.*x2z*x12s + 2./3.*x2z - 5./12.*ux + *x12s - 1./6.*ux - 1./3.*x12s + 1./12.*x12c + 13./12.*x2zs ) GG10c = GG10c + rs1wg(11)*vvc*GtM2*Xps*ys * ( - 2./3. ) GG10c = GG10c + rs1wg(11)*vvc*GtM2*ys*Xp * ( - 1./3.*O + 5./3.* + x12 + 4./3.*x2z - 2./3.*ux ) GG10c = GG10c + rs1wg(11)*vvc*GtM2*ys * ( 1./3.*O*x12 + 1./3.*O* + x2z - 5./3.*x12*x2z + 2./3.*x12*ux + 2./3.*x2z*ux - x12s - 2./ + 3.*x2zs ) GG10c = GG10c + rs1wg(12)*XX*vvc*GtM2*ys * ( - 2./3.*O*x12*ux - + 2./3.*O*x2z*ux - 2./3.*O*x2z + 2./3.*O*uxs - 4./3.*Os*x12 - 2. + /3.*Os*x2z + 4./3.*Os*ux + 2./3.*Oc - 4./3.*x12*x2z*ux + 4./3. + *x12*x2z + 2./3.*x12*x2zs + 2./3.*x12*uxs - 2./3.*x2z*ux + 8./ + 3.*x2z*x12s - 2*ux*x12s + 2*x12c + 2./3.*x2zs ) GG10c = GG10c + rs1wg(12)*Xpi*vvc*GtM2*ys * ( 2./3.*O*x12 - 2./3. + *O + 2./3.*x12*x2z + 4./3.*x12 + 4./3.*x2z - 2./3.*ux + 2./3. + *x12s ) GG10c = GG10c + rs1wg(12)*vvc*GtM2*ys*Xp * ( 4./3. - 2./3.*O - 2./ + 3.*x12 ) GG10c = GG10c + rs1wg(12)*vvc*GtM2*ys * ( - 2./3. + 1./3.*O*x2z + + 1./3.*O + 1./3.*Os - 1./3.*x12*x2z - 5./3.*x12 - 2*x2z + 2. + /3.*ux - 1./3.*x12s ) GG10c = GG10c + rs1wg(13)*XX*vvc*GtM2*ys * ( O*x12*x2z - 1./3.*O* + x12*ux - 2./3.*O*x2z*ux - 1./3.*O*x2z + 1./3.*O*x2zs + 1./3.* + O*uxs - 2./3.*Os*x12 - 2./3.*Os*x2z + 2./3.*Os*ux + 1./3.*Oc + - 1./3.*x12*x2z*ux + 2./3.*x12*x2z - 1./3.*x12*x2zs + 1./3.* + x12*uxs - 1./3.*x2z*ux + 2./3.*x2z*x12s - ux*x12s + x12c + 2./ + 3.*x2zs ) GG10c = GG10c + rs1wg(13)*Xpi*vvc*GtM2*ys * ( 1./3.*O*x12 - 1./3. + *O - 1./3.*x12*x2z + 2./3.*x12 + x2z - 1./3.*ux + 1./3.*x12s + ) GG10c = GG10c + rs1wg(13)*vvc*GtM2*ys*Xp * ( 1./2. ) GG10c = GG10c + rs1wg(13)*vvc*GtM2*ys * ( - 1./3. - 1./6.*O*x12 + + 1./6.*O*ux + 1./2.*O + 1./2.*x12*x2z + 1./6.*x12*ux - 7./6. + *x12 - 7./6.*x2z + 1./3.*ux - 1./6.*x12s ) GG10c = GG10c + rs1wg(14)*XX*vvc*GtM2*ys * ( - O*x12*x2z + 1./3. + *O*x2z*ux - 1./3.*O*x2zs + 1./3.*Os*x2z - 1./3.*x12*x2z*ux + + 2./3.*x12*x2zs + 2./3.*x2z*x12s - 1./3.*x2zs ) GG10c = GG10c + rs1wg(14)*Xpi*vvc*GtM2*ys * ( - O*x12 - 1./3.*O* + x2z + 1./3.*O*ux + 1./3.*Os + 4./3.*x12*x2z - 1./3.*x12*ux - + 2./3.*x2z + 2./3.*x12s ) GG10c = GG10c + rs1wg(14)*vvc*GtM2*Xps*ys * ( - 1./3. ) GG10c = GG10c + rs1wg(14)*vvc*GtM2*ys*Xp * ( - 1./6. - 1./6.*O + + 3./2.*x12 + 2./3.*x2z - 1./3.*ux ) GG10c = GG10c + rs1wg(14)*vvc*GtM2*ys * ( 1./3. + 2./3.*O*x12 + 1. + /3.*O*x2z - 1./6.*O*ux - 1./6.*Os - 13./6.*x12*x2z + 1./2.* + x12*ux - 2./3.*x12 + 1./3.*x2z*ux + 1./2.*x2z - 7./6.*x12s - + 1./3.*x2zs ) GG20c = + + rs1wg(1)*XX*vvc*GtM2*ys * ( 12*O*x12*x2z - 4*O*x12*ux - 8*O* + x2z*ux - 4*O*x2z + 4*O*x2zs + 4*O*uxs - 8*Os*x12 - 8*Os*x2z + + 8*Os*ux + 4*Oc - 4*x12*x2z*ux + 8*x12*x2z - 4*x12*x2zs + 4 + *x12*uxs - 4*x2z*ux + 8*x2z*x12s - 12*ux*x12s + 12*x12c + 8* + x2zs ) GG20c = GG20c + rs1wg(1)*Xpi*vvc*GtM2*ys * ( 16*O*x12 + 4*O*x2z + - 4*O*ux - 4*O - 4*Os - 12*x12*x2z + 4*x12*ux + 8*x12 + 16* + x2z - 4*ux - 4*x12s ) GG20c = GG20c + rs1wg(1)*vvc*GtM2*Xps*ys * ( - 4./3. ) GG20c = GG20c + rs1wg(1)*vvc*GtM2*ys*Xp * ( 10 - 10./3.*O - 22./3. + *x12 + 8./3.*x2z - 4./3.*ux ) GG20c = GG20c + rs1wg(1)*vvc*GtM2*ys * ( - 8 - 44./3.*O*x12 - 2* + O*x2z + 14./3.*O*ux + 10./3.*O + 16./3.*Os + 32./3.*x12*x2z + - 10./3.*x12*ux - 2./3.*x12 + 4./3.*x2z*ux - 18*x2z + 4*ux + + 4*x12s - 4./3.*x2zs ) GG20c = GG20c + rs1wg(4)*XX*vvc*GtM2*ys * ( 2*O*x12*x2z*ux - O* + x12*x2zs - 2*O*x2z*x12s - O*x2z*uxs + O*ux*x2zs - O*x2zs + Os + *x12*x2z - Os*x2z*ux + Os*x2zs - x12*x2z*uxs + 2*x12*ux*x2zs + - x12*x2zc - x12*x2zs + 3*x2z*ux*x12s - 3*x2z*x12c + ux*x2zs + - 4*x12s*x2zs - x2zc ) GG20c = GG20c + rs1wg(4)*Xpi*vvc*GtM2*ys * ( - 2*O*x12*x2z + 2*O + *x12*ux + O*x2z*ux - 2*O*x2z - 2*O*x12s - O*uxs + Os*x12 + 2* + Os*x2z - Os*ux + 2*x12*x2z*ux - 2*x12*x2z - 2*x12*x2zs - x12* + uxs + 2*x2z*ux - 6*x2z*x12s + 3*ux*x12s - 3*x12c - 3*x2zs ) GG20c = GG20c + rs1wg(4)*vvc*GtM2*Xps*ys * ( 4./3. - 1./3.*O - 1./ + 3.*x12 ) GG20c = GG20c + rs1wg(4)*vvc*GtM2*ys*Xp * ( - 1 - 1./2.*O*x12 + + 1./2.*O*x2z - 1./6.*O*ux - 7./6.*O + 7./6.*Os - 1./2.*x12*x2z + - 1./6.*x12*ux - 4./3.*x12 - 11./3.*x2z + 7./6.*ux - 5./3.* + x12s ) GG20c = GG20c + rs1wg(4)*vvc*GtM2*ys * ( 5./3.*O*x12*x2z - 13./6. + *O*x12*ux + O*x12 - O*x2z*ux + 13./6.*O*x2z + 2*O*x12s - 1./6. + *O*x2zs + 7./6.*O*uxs + O - 7./6.*Os*x12 - 13./6.*Os*x2z + 7./ + 6.*Os*ux - Os - 13./6.*x12*x2z*ux + 10./3.*x12*x2z + 11./6.* + x12*x2zs + 7./6.*x12*uxs + x12 - 13./6.*x2z*ux + 6*x2z*x12s + + 3*x2z - 10./3.*ux*x12s - ux + 2*x12s + 19./6.*x12c + 10./3. + *x2zs ) GG20c = GG20c + rs1wg(7)*XX*vvc*GtM2*ys * ( 1./2.*O*x12*x2zs - 1./ + 2.*O*x2z*x12s + 1./2.*O*x2z*uxs - O*ux*x2zs + 1./2.*O*x2zc - + O*x2zs - 1./2.*Os*x12*x2z + Os*x2z*ux - 1./2.*Os*x2zs + 1./2. + *Oc*x2z + 1./2.*x12*x2z*uxs - 1./2.*x12*ux*x2zs - 1./2.*x12* + x2zc + 2*x12*x2zs - x2z*ux*x12s + 1./2.*x2z*x12c - 1./2.*ux* + x2zs + x2zc ) GG20c = GG20c + rs1wg(7)*Xpi*vvc*GtM2*ys * ( - 1./2.*O*x12*x2z + - O*x2z + 1./2.*Os*x2z + 2*x12*x2z - 1./2.*x12*x2zs - 1./2.* + x2z*ux - x2z*x12s + 3./2.*x2zs ) GG20c = GG20c + rs1wg(7)*vvc*GtM2*Xps*ys * ( - 1./6. + 1./6.*O + + 1./6.*x12 ) GG20c = GG20c + rs1wg(7)*vvc*GtM2*ys*Xp * ( - 1./2.*O*x12 - 1./2. + *O*x2z + 1./6.*O*ux - 1./6.*O + 1./12.*Os - 1./2.*x12*x2z + 1. + /6.*x12*ux + 1./2.*x12 + x2z - 1./6.*ux - 7./12.*x12s ) GG20c = GG20c + rs1wg(7)*vvc*GtM2*ys * ( O*x12*x2z - 1./6.*O*x12* + ux - 1./6.*O*x2z*ux + 7./6.*O*x2z + 1./4.*O*x12s + 1./3.*O* + x2zs - 1./2.*Os*x2z - 1./12.*Os*ux - 1./12.*Oc - 1./6.*x12* + x2z*ux - 5./2.*x12*x2z + 5./6.*x12*x2zs + 2./3.*x2z*ux + 3./2. + *x2z*x12s - 1./2.*x2z - 1./12.*ux*x12s + 1./6.*x12c - 11./6.* + x2zs ) GG20c = GG20c + rs1wg(8)*Xpi*vvc*GtM2*ys * ( O*x12*x2z - 2*O*x12* + ux - O*x2z*ux + O*x2z + 2*O*x12s + O*uxs - Os*x12 - Os*x2z + + Os*ux - 2*x12*x2z*ux + x12*x2z + x12*x2zs + x12*uxs - x2z*ux + + 4*x2z*x12s - 3*ux*x12s + 3*x12c + x2zs ) GG20c = GG20c + rs1wg(8)*vvc*GtM2*Xps*ys * ( - 1 ) GG20c = GG20c + rs1wg(8)*vvc*GtM2*ys*Xp * ( 1 + 5./6.*O*x12 + 1./ + 6.*O*x2z - 1./6.*O*ux + 4./3.*O - 4./3.*Os + 7./6.*x12*x2z - + 1./6.*x12*ux + 5./6.*x12 + 11./6.*x2z - 5./6.*ux + 13./6.* + x12s ) GG20c = GG20c + rs1wg(8)*vvc*GtM2*ys * ( - O*x12*x2z + 13./6.*O* + x12*ux - O*x12 + 7./6.*O*x2z*ux - 4./3.*O*x2z - 2*O*x12s - 7./ + 6.*O*uxs - O + 7./6.*Os*x12 + 7./6.*Os*x2z - 7./6.*Os*ux + Os + + 7./3.*x12*x2z*ux - 11./6.*x12*x2z - 7./6.*x12*x2zs - 7./6. + *x12*uxs - x12 + 5./6.*x2z*ux - 13./3.*x2z*x12s - 2*x2z + 10./ + 3.*ux*x12s + ux - 2*x12s - 19./6.*x12c - 5./6.*x2zs ) GG20c = GG20c + rs1wg(9)*Xpi*vvc*GtM2*ys * ( - 6*O*x12 - 2*O*x2z + + 2*O*ux + 2*Os + 4*x12*x2z - 2*x12*ux - 2*x2z + 4*x12s ) GG20c = GG20c + rs1wg(9)*vvc*GtM2*ys*Xp * ( - 2 + 4*x12 ) GG20c = GG20c + rs1wg(9)*vvc*GtM2*ys * ( 2 + 20./3.*O*x12 + 8./3. + *O*x2z - 8./3.*O*ux - 8./3.*Os - 4*x12*x2z + 2*x12*ux - 4*x12 + + 2*x2z - 4*x12s ) GG20c = GG20c + rs1wg(10)*XX*vvc*GtM2*ys * ( - 1./2.*O*x12*x2zs + + 1./2.*O*x2z*x12s - 1./2.*O*x2z*uxs + O*ux*x2zs - 1./2.*O* + x2zc + O*x2zs + 1./2.*Os*x12*x2z - Os*x2z*ux + 1./2.*Os*x2zs + - 1./2.*Oc*x2z - 1./2.*x12*x2z*uxs + 1./2.*x12*ux*x2zs + 1./ + 2.*x12*x2zc - 2*x12*x2zs + x2z*ux*x12s - 1./2.*x2z*x12c + 1./ + 2.*ux*x2zs - x2zc ) GG20c = GG20c + rs1wg(10)*Xpi*vvc*GtM2*ys * ( O*x2z*ux + 2*O*x2z + + 1./2.*O*x12s - 1./2.*O*x2zs - 1./2.*O*uxs + 1./2.*Os*x12 + - Os*ux - 1./2.*Oc + 1./2.*x12*x2z*ux - 4*x12*x2z + x12*x2zs + - 1./2.*x12*uxs + x2z*ux + x2z*x12s + ux*x12s - 1./2.*x12c + - 5./2.*x2zs ) GG20c = GG20c + rs1wg(10)*vvc*GtM2*Xps*ys * ( 2./3. - 1./6.*O - 1. + /6.*x12 ) GG20c = GG20c + rs1wg(10)*vvc*GtM2*ys*Xp * ( - 1./2. + 11./12.*O + *x12 + 5./12.*O*x2z - 1./12.*O*ux + 7./6.*O - 7./12.*Os + 11./ + 12.*x12*x2z - 1./12.*x12*ux - 29./12.*x12 - 29./12.*x2z + 7./ + 12.*ux + 3./2.*x12s ) GG20c = GG20c + rs1wg(10)*vvc*GtM2*ys * ( - 1./3.*O*x12*x2z + 1./ + 12.*O*x12*ux - 1./2.*O*x12 - 13./12.*O*x2z*ux - 13./6.*O*x2z + - 3./4.*O*x12s + 1./3.*O*x2zs + 7./12.*O*uxs - O - 7./12.*Os + *x12 - 1./12.*Os*x2z + 7./6.*Os*ux + 1./2.*Os + 7./12.*Oc - 1. + /2.*x12*x2z*ux + 47./12.*x12*x2z - 5./4.*x12*x2zs + 7./12.* + x12*uxs + 2*x12 - 13./12.*x2z*ux - 4./3.*x2z*x12s + 2*x2z - + 13./12.*ux*x12s - 1./2.*ux - x12s + 5./12.*x12c + 11./4.*x2zs + ) GG20c = GG20c + rs1wg(11)*vvc*GtM2*Xps*ys * ( - 2./3. ) GG20c = GG20c + rs1wg(11)*vvc*GtM2*ys*Xp * ( - 1./3.*O + 5./3.* + x12 + 4./3.*x2z - 2./3.*ux ) GG20c = GG20c + rs1wg(11)*vvc*GtM2*ys * ( 1./3.*O*x12 + 1./3.*O* + x2z - 5./3.*x12*x2z + 2./3.*x12*ux + 2./3.*x2z*ux - x12s - 2./ + 3.*x2zs ) GG20c = GG20c + rs1wg(12)*XX*vvc*GtM2*ys * ( - 2*O*x12*ux - 2*O* + x2z*ux - 2*O*x2z + 2*O*uxs - 4*Os*x12 - 2*Os*x2z + 4*Os*ux + + 2*Oc - 4*x12*x2z*ux + 4*x12*x2z + 2*x12*x2zs + 2*x12*uxs - 2* + x2z*ux + 8*x2z*x12s - 6*ux*x12s + 6*x12c + 2*x2zs ) GG20c = GG20c + rs1wg(12)*Xpi*vvc*GtM2*ys * ( 2*O*x12 - 2*O + 2* + x12*x2z + 4*x12 + 4*x2z - 2*ux + 2*x12s ) GG20c = GG20c + rs1wg(12)*vvc*GtM2*ys*Xp * ( 8./3. - 2./3.*O - 2./ + 3.*x12 ) GG20c = GG20c + rs1wg(12)*vvc*GtM2*ys * ( - 2 - 4./3.*O*x12 + 1./ + 3.*O*x2z + 5./3.*O + 1./3.*Os - 5./3.*x12*x2z - 13./3.*x12 - + 14./3.*x2z + 2*ux - 5./3.*x12s ) GG20c = GG20c + rs1wg(13)*XX*vvc*GtM2*ys * ( 3*O*x12*x2z - O*x12* + ux - 2*O*x2z*ux - O*x2z + O*x2zs + O*uxs - 2*Os*x12 - 2*Os* + x2z + 2*Os*ux + Oc - x12*x2z*ux + 2*x12*x2z - x12*x2zs + x12* + uxs - x2z*ux + 2*x2z*x12s - 3*ux*x12s + 3*x12c + 2*x2zs ) GG20c = GG20c + rs1wg(13)*Xpi*vvc*GtM2*ys * ( O*x12 - O - x12*x2z + + 2*x12 + 3*x2z - ux + x12s ) GG20c = GG20c + rs1wg(13)*vvc*GtM2*ys*Xp * ( 7./6. ) GG20c = GG20c + rs1wg(13)*vvc*GtM2*ys * ( - 1 - 5./6.*O*x12 + 1./ + 6.*O*ux + 7./6.*O + 7./6.*x12*x2z + 1./6.*x12*ux - 5./2.*x12 + - 19./6.*x2z + ux - 5./6.*x12s ) GG20c = GG20c + rs1wg(14)*XX*vvc*GtM2*ys * ( - 3*O*x12*x2z + O* + x2z*ux - O*x2zs + Os*x2z - x12*x2z*ux + 2*x12*x2zs + 2*x2z* + x12s - x2zs ) GG20c = GG20c + rs1wg(14)*Xpi*vvc*GtM2*ys * ( - 3*O*x12 - O*x2z + + O*ux + Os + 4*x12*x2z - x12*ux - 2*x2z + 2*x12s ) GG20c = GG20c + rs1wg(14)*vvc*GtM2*Xps*ys * ( - 1./3. ) GG20c = GG20c + rs1wg(14)*vvc*GtM2*ys*Xp * ( - 5./6. - 1./6.*O + + 17./6.*x12 + 2./3.*x2z - 1./3.*ux ) GG20c = GG20c + rs1wg(14)*vvc*GtM2*ys * ( 1 + 8./3.*O*x12 + O*x2z + - 5./6.*O*ux - 5./6.*Os - 29./6.*x12*x2z + 7./6.*x12*ux - 2* + x12 + 1./3.*x2z*ux + 11./6.*x2z - 5./2.*x12s - 1./3.*x2zs ) GG11c = + + rs1wg(1)*XX*vvc*GtM2*ys * ( - 4./3.*O*x12*x2zs - 8./3.*O*x2z + *x12s - 4./3.*O*x12c - 4./3.*x12*x12c + 4./3.*x12*x2zc - 8./3. + *x12*x2zs - 4./3.*x2z*x12s - 4./3.*x2z*x12c + 4./3.*x12s*x2zs + - 4./3.*x2zc ) GG11c = GG11c + rs1wg(1)*Xpi*vvc*GtM2*ys * ( - 4*O*x12*x2z - 4*O + *x12s - 20./3.*x12*x2z + 4*x12*x2zs + 8./3.*x2z*x12s - 4./3.* + x12s - 4./3.*x12c - 16./3.*x2zs ) GG11c = GG11c + rs1wg(1)*vvc*GtM2*Xps*ys*Xp * ( - 4./3. ) GG11c = GG11c + rs1wg(1)*vvc*GtM2*Xps*ys * ( 16./3. - 8./3.*O - 8. + /3.*x12 + 8./3.*x2z ) GG11c = GG11c + rs1wg(1)*vvc*GtM2*ys*Xp * ( - 10./3. + 4*O*x2z + + 2./3.*O + 16./3.*x12*x2z - 10./3.*x12 - 34./3.*x2z + 4* + x12s - 4./3.*x2zs ) GG11c = GG11c + rs1wg(1)*vvc*GtM2*ys * ( 4./3.*O*x12*x2z + 10./3. + *O*x12 - 2./3.*O*x2z + 8./3.*O*x12s - 4./3.*O*x2zs - 2./3.*Os + *x12 + 4*x12*x2z - 14./3.*x12*x2zs + 4*x12 - 4*x2z*x12s + 22./ + 3.*x2z - 2./3.*x12s + 2./3.*x12c + 22./3.*x2zs ) GG11c = GG11c + rs1wg(4)*XX*vvc*GtM2*ys * ( 2./3.*O*x12*x2zs + 2./ + 3.*O*x2z*x12c + 2./3.*O*x12s*x2zs + 1./3.*Os*x2z*x12s + 1./3. + *x12*x2z*x12c + 1./3.*x12s*x2zc + 2./3.*x12c*x2zs + 1./3.* + x2zc ) GG11c = GG11c + rs1wg(4)*Xpi*vvc*GtM2*ys * ( 4./3.*O*x12*x2z + O* + x12*x2zs + 5./3.*O*x2z*x12s + 2./3.*O*x12c + 1./3.*Os*x12s + + 1./3.*x12*x12c + 1./3.*x12*x2zc - 1./3.*x12*x2zs + 5./3.*x2z* + x12c + 5./3.*x12s*x2zs + 1./3.*x2zc + x2zs ) GG11c = GG11c + rs1wg(4)*vvc*GtM2*Xps*ys*Xp * ( 2./3. - 1./3.*O + - 1./3.*x12 ) GG11c = GG11c + rs1wg(4)*vvc*GtM2*Xps*ys * ( - 5./6. - 2./3.*O* + x12 + 2./3.*O*x2z + 1./6.*O + 1./3.*x12*x2z + 1./6.*x12 - 11./ + 6.*x2z - 2./3.*x12s ) GG11c = GG11c + rs1wg(4)*vvc*GtM2*ys*Xp * ( 1./3. + 4./3.*O*x12* + x2z + 11./6.*O*x12 - 1./3.*O*x2z + 2./3.*O*x12s - 1./3.*O* + x2zs - 1./6.*Os*x12 + 1./6.*x12*x2z + 1./6.*x12*x2zs - 1./3.* + x12 + 5./3.*x2z*x12s + 7./3.*x2z + 7./6.*x12s + 5./6.*x12c + + 5./3.*x2zs ) GG11c = GG11c + rs1wg(4)*vvc*GtM2*ys * ( - 7./2.*O*x12*x2z - 2./ + 3.*O*x12*x2zs - 2./3.*O*x12 - 4./3.*O*x2z*x12s - O*x12s - 2./ + 3.*O*x12c + 1./6.*O*x2zs + 1./6.*Os*x12*x2z - 1./3.*Os*x12s + + 2./3.*x12*x2z - 1./3.*x12*x12c - 1./6.*x12*x2zc - 2./3.* + x12*x2zs - 5./2.*x2z*x12s - 3./2.*x2z*x12c - x2z - 4./3.*x12s + *x2zs - x12c - 1./2.*x2zc - 13./6.*x2zs ) GG11c = GG11c + rs1wg(7)*XX*vvc*GtM2*ys * ( 1./6.*O*x12*x2zs - 1./ + 6.*O*x2zc - 1./6.*Os*x12*x2zs + 1./6.*x12*x2zq - 2./3.*x12* + x2zc + 1./3.*x12s*x2zc - 1./2.*x12s*x2zs + 1./6.*x12c*x2zs - + 1./6.*x2zq + 1./6.*x2zc ) GG11c = GG11c + rs1wg(7)*Xpi*vvc*GtM2*ys * ( 1./6.*O*x12*x2z + 1./ + 2.*O*x12*x2zs + 1./6.*O*x2z*x12s - 1./3.*O*x2zs - 1./6.*Os* + x12*x2z + 1./3.*x12*x2zc - 5./3.*x12*x2zs - 1./2.*x2z*x12s + + 1./3.*x2z*x12c + 2./3.*x12s*x2zs - 1./2.*x2zc + 1./3.*x2zs ) GG11c = GG11c + rs1wg(7)*vvc*GtM2*Xps*ys*Xp * ( - 1./6. + 1./6.* + O + 1./6.*x12 ) GG11c = GG11c + rs1wg(7)*vvc*GtM2*Xps*ys * ( - 1./12. - 5./6.*O* + x12 - 1./2.*O*x2z + 1./6.*O - 1./6.*Os - 1./2.*x12*x2z + 2./3. + *x12 + 2./3.*x2z - 2./3.*x12s ) GG11c = GG11c + rs1wg(7)*vvc*GtM2*ys*Xp * ( 11./6.*O*x12*x2z - 1./ + 2.*O*x2z + 2./3.*O*x12s + 1./2.*O*x2zs + 1./6.*Os*x12 + 1./3. + *Os*x2z - 13./6.*x12*x2z + 2./3.*x12*x2zs + 3./2.*x2z*x12s + + 1./6.*x2z - 1./3.*x12s + 1./2.*x12c - x2zs ) GG11c = GG11c + rs1wg(7)*vvc*GtM2*ys * ( - 2./3.*O*x12*x2z - O* + x12*x2zs - 2./3.*O*x2z*x12s + 1./6.*O*x2z - 1./6.*O*x12c - 1./ + 6.*O*x2zc + 1./2.*O*x2zs - 1./12.*Os*x12s - 1./6.*Os*x2zs + + x12*x2z - 1./12.*x12*x12c - 1./2.*x12*x2zc + 2*x12*x2zs + 1./ + 3.*x2z*x12s - 2./3.*x2z*x12c - 1./6.*x2z - 13./12.*x12s*x2zs + + 2./3.*x2zc + 1./12.*x2zs ) GG11c = GG11c + rs1wg(8)*Xpi*vvc*GtM2*ys * ( - 2./3.*O*x12*x2z + - 2./3.*O*x2z*x12s - 2./3.*O*x12c - 1./3.*Os*x12s - 1./3.* + x12*x12c - 2./3.*x2z*x12c - 1./3.*x12s*x2zs - 1./3.*x2zs ) GG11c = GG11c + rs1wg(8)*vvc*GtM2*Xps*ys*Xp * ( - 1./3. ) GG11c = GG11c + rs1wg(8)*vvc*GtM2*Xps*ys * ( 2./3. + O*x12 + 1./3. + *x12*x2z - 1./3.*x12 + 2./3.*x2z + x12s ) GG11c = GG11c + rs1wg(8)*vvc*GtM2*ys*Xp * ( - 1./3. - O*x12*x2z + - 5./3.*O*x12 - O*x12s - 1./3.*x12*x2zs + 1./3.*x12 - 4./3.* + x2z*x12s - 4./3.*x2z - x12s - x12c - 1./3.*x2zs ) GG11c = GG11c + rs1wg(8)*vvc*GtM2*ys * ( 5./3.*O*x12*x2z + 2./3.* + O*x12 + 2./3.*O*x2z*x12s + O*x12s + 2./3.*O*x12c + 1./3.*Os* + x12s - 1./3.*x12*x2z + 1./3.*x12*x12c + 1./3.*x12*x2zs + 4./3. + *x2z*x12s + 2./3.*x2z*x12c + 2./3.*x2z + 1./3.*x12s*x2zs + + x12c + 2./3.*x2zs ) GG11c = GG11c + rs1wg(9)*Xpi*vvc*GtM2*ys * ( 2./3.*O*x12*x2z + 2./ + 3.*O*x12s + 2./3.*x12*x2z - 2./3.*x12*x2zs - 4./3.*x2z*x12s + - 2./3.*x12c + 2./3.*x2zs ) GG11c = GG11c + rs1wg(9)*vvc*GtM2*Xps*ys * ( - 2./3. + 4./3.*x12 + ) GG11c = GG11c + rs1wg(9)*vvc*GtM2*ys*Xp * ( 2./3. + 2./3.*O*x12 + - 2*x12*x2z - 2./3.*x12 + 4./3.*x2z - 2*x12s ) GG11c = GG11c + rs1wg(9)*vvc*GtM2*ys * ( - 2./3.*O*x12*x2z - 2./ + 3.*O*x12 - 2./3.*O*x12s + 4./3.*x12*x2z + 2./3.*x12*x2zs - 2./ + 3.*x12 + 4./3.*x2z*x12s - 4./3.*x2z + 2*x12s + 2./3.*x12c - 2. + /3.*x2zs ) GG11c = GG11c + rs1wg(10)*XX*vvc*GtM2*ys * ( - 1./6.*O*x12*x2zs + + 1./6.*O*x2zc + 1./6.*Os*x12*x2zs - 1./6.*x12*x2zq + 2./3.* + x12*x2zc - 1./3.*x12s*x2zc + 1./2.*x12s*x2zs - 1./6.*x12c* + x2zs + 1./6.*x2zq - 1./6.*x2zc ) GG11c = GG11c + rs1wg(10)*Xpi*vvc*GtM2*ys * ( - 1./3.*O*x12*x2z + - 1./2.*O*x12*x2zs - 1./6.*O*x2z*x12s + 1./2.*O*x2zs + 1./3. + *Os*x12*x2z - 1./2.*x12*x2zc + 7./3.*x12*x2zs + x2z*x12s - 1./ + 2.*x2z*x12c - x12s*x2zs + 2./3.*x2zc - 1./2.*x2zs ) GG11c = GG11c + rs1wg(10)*vvc*GtM2*Xps*ys*Xp * ( 1./3. - 1./6.*O + - 1./6.*x12 ) GG11c = GG11c + rs1wg(10)*vvc*GtM2*Xps*ys * ( 1./12. + 4./3.*O* + x12 + 1./2.*O*x2z - 1./3.*O + 1./6.*Os + 2./3.*x12*x2z - 5./3. + *x12 - 7./6.*x2z + 7./6.*x12s ) GG11c = GG11c + rs1wg(10)*vvc*GtM2*ys*Xp * ( - 1./6. - 7./3.*O* + x12*x2z - 2./3.*O*x12 + 5./6.*O*x2z - 5./6.*O*x12s - 1./2.*O* + x2zs + 1./6.*O - 1./3.*Os*x2z + 11./3.*x12*x2z - x12*x2zs + + x12 - 13./6.*x2z*x12s + 1./3.*x12s - 5./6.*x12c + 3./2.*x2zs + ) GG11c = GG11c + rs1wg(10)*vvc*GtM2*ys * ( 4./3.*O*x12*x2z + O*x12 + *x2zs + 1./6.*O*x12 + 2./3.*O*x2z*x12s - 1./2.*O*x2z + 1./6.* + O*x12s + 1./6.*O*x12c + 1./6.*O*x2zc - 2./3.*O*x2zs - 1./6.* + Os*x12*x2z - 1./6.*Os*x12 + 1./12.*Os*x12s + 1./6.*Os*x2zs - + 8./3.*x12*x2z + 1./12.*x12*x12c + 2./3.*x12*x2zc - 7./3.*x12* + x2zs - 1./6.*x2z*x12s + 5./6.*x2z*x12c + 1./2.*x2z + 17./12.* + x12s*x2zs - 1./2.*x12s + 1./3.*x12c - 5./6.*x2zc - 5./12.* + x2zs ) GG11c = GG11c + rs1wg(11)*vvc*GtM2*Xps*ys*Xp * ( - 2./3. ) GG11c = GG11c + rs1wg(11)*vvc*GtM2*Xps*ys * ( 1./3. + 4./3.*x12 + + 4./3.*x2z ) GG11c = GG11c + rs1wg(11)*vvc*GtM2*ys*Xp * ( - 1./3.*O*x12 - 5./ + 3.*x12*x2z - 1./3.*x12 - 2./3.*x2z - x12s - 2./3.*x2zs ) GG11c = GG11c + rs1wg(11)*vvc*GtM2*ys * ( 1./3.*O*x12*x2z + 1./3. + *O*x12s + 1./3.*x12*x2z + 1./3.*x12*x2zs + 2./3.*x2z*x12s + 1. + /3.*x12c + 1./3.*x2zs ) GG11c = GG11c + rs1wg(12)*XX*vvc*GtM2*ys * ( - 2./3.*O*x2z*x12s + - 2./3.*O*x12c - 2./3.*x12*x12c - 2./3.*x12*x2zs - 2./3.*x2z + *x12s - 4./3.*x2z*x12c - 2./3.*x12s*x2zs ) GG11c = GG11c + rs1wg(12)*Xpi*vvc*GtM2*ys * ( - 2./3.*O*x12*x2z + - 4./3.*O*x12s - 2*x12*x2z - 2./3.*x12*x2zs - 2*x2z*x12s - 2. + /3.*x12s - 4./3.*x12c - 2./3.*x2zs ) GG11c = GG11c + rs1wg(12)*vvc*GtM2*Xps*ys * ( 4./3. - 2./3.*O - 2. + /3.*x12 ) GG11c = GG11c + rs1wg(12)*vvc*GtM2*ys*Xp * ( - 1 + 2./3.*O*x2z + + 1./3.*O - 5./3.*x12 - 7./3.*x2z ) GG11c = GG11c + rs1wg(12)*vvc*GtM2*ys * ( O*x12 - 1./3.*O*x2z + 2. + /3.*O*x12s - 1./3.*Os*x12 + 10./3.*x12*x2z + 1./3.*x12*x2zs + + 4./3.*x12 + 4./3.*x2z*x12s + 5./3.*x2z + 5./3.*x12s + x12c + + x2zs ) GG11c = GG11c + rs1wg(13)*XX*vvc*GtM2*ys * ( - 1./3.*O*x12*x2zs + - 2./3.*O*x2z*x12s - 1./3.*O*x12c - 1./3.*x12*x12c + 1./3.* + x12*x2zc - 2./3.*x12*x2zs - 1./3.*x2z*x12s - 1./3.*x2z*x12c + + 1./3.*x12s*x2zs - 1./3.*x2zc ) GG11c = GG11c + rs1wg(13)*Xpi*vvc*GtM2*ys * ( - 2./3.*O*x12*x2z + - 2./3.*O*x12s - 4./3.*x12*x2z + 2./3.*x12*x2zs - 1./3.*x12s + - 2./3.*x12c - x2zs ) GG11c = GG11c + rs1wg(13)*vvc*GtM2*Xps*ys * ( 1./2. ) GG11c = GG11c + rs1wg(13)*vvc*GtM2*ys*Xp * ( - 1./6. + 1./6.*O* + x12 + 1./3.*O*x2z - 1./6.*O + 5./6.*x12*x2z - 5./3.*x12 - 3./ + 2.*x2z + 1./6.*x12s ) GG11c = GG11c + rs1wg(13)*vvc*GtM2*ys * ( 1./6.*O*x12*x2z + 1./6. + *O*x12 + 1./6.*O*x2z + 1./2.*O*x12s - 1./3.*O*x2zs + 1./6.*Os + *x12 + 11./6.*x12*x2z - x12*x2zs + 2./3.*x12 - 2./3.*x2z*x12s + + 5./6.*x2z + 7./6.*x12s + 1./3.*x12c + 4./3.*x2zs ) GG11c = GG11c + rs1wg(14)*XX*vvc*GtM2*ys * ( 1./3.*O*x12*x2zs + 1. + /3.*O*x2z*x12s - 1./3.*x12*x2zc + 1./3.*x12*x2zs - 1./3.*x2z* + x12c - 2./3.*x12s*x2zs + 1./3.*x2zc ) GG11c = GG11c + rs1wg(14)*Xpi*vvc*GtM2*ys * ( 2./3.*O*x12*x2z + 1. + /3.*O*x12s + 2./3.*x12*x2z - 4./3.*x12*x2zs - 5./3.*x2z*x12s + - 1./3.*x12c + x2zs ) GG11c = GG11c + rs1wg(14)*vvc*GtM2*Xps*ys*Xp * ( - 1./3. ) GG11c = GG11c + rs1wg(14)*vvc*GtM2*Xps*ys * ( - 1./3. + 1./3.*O + + 5./3.*x12 + 2./3.*x2z ) GG11c = GG11c + rs1wg(14)*vvc*GtM2*ys*Xp * ( 1./3. - 1./3.*O*x12 + - 2./3.*O*x2z - 3*x12*x2z - 1./3.*x12 + x2z - 2*x12s - 1./3. + *x2zs ) GG11c = GG11c + rs1wg(14)*vvc*GtM2*ys * ( - 1./3.*O*x12 + 1./3.* + O*x2zs + x12*x2z + 5./3.*x12*x2zs - 1./3.*x12 + 7./3.*x2z* + x12s - x2z + x12s + 2./3.*x12c - x2zs ) GG21c = + + rs1wg(1)*XX*vvc*GtM2*ys * ( - 4*O*x12*x2zs - 8*O*x2z*x12s - + 4*O*x12c - 4*x12*x12c + 4*x12*x2zc - 8*x12*x2zs - 4*x2z*x12s + - 4*x2z*x12c + 4*x12s*x2zs - 4*x2zc ) GG21c = GG21c + rs1wg(1)*Xpi*vvc*GtM2*ys * ( - 12*O*x12*x2z - 12 + *O*x12s - 20*x12*x2z + 12*x12*x2zs + 8*x2z*x12s - 4*x12s - 4* + x12c - 16*x2zs ) GG21c = GG21c + rs1wg(1)*vvc*GtM2*Xps*ys*Xp * ( - 4./3. ) GG21c = GG21c + rs1wg(1)*vvc*GtM2*Xps*ys * ( 32./3. - 8./3.*O - 8 + *x12 + 8./3.*x2z ) GG21c = GG21c + rs1wg(1)*vvc*GtM2*ys*Xp * ( - 26./3. - 16./3.*O* + x12 + 4*O*x2z + 2./3.*O + 16*x12*x2z - 6*x12 - 74./3.*x2z + + 28./3.*x12s - 4./3.*x2zs ) GG21c = GG21c + rs1wg(1)*vvc*GtM2*ys * ( 28./3.*O*x12*x2z + 26./3. + *O*x12 - 2./3.*O*x2z + 32./3.*O*x12s - 4./3.*O*x2zs - 2./3.* + Os*x12 + 20./3.*x12*x2z - 38./3.*x12*x2zs + 12*x12 - 28./3.* + x2z*x12s + 62./3.*x2z - 10./3.*x12s + 10./3.*x12c + 18*x2zs ) GG21c = GG21c + rs1wg(4)*XX*vvc*GtM2*ys * ( 2*O*x12*x2zs + 2*O* + x2z*x12c + 2*O*x12s*x2zs + Os*x2z*x12s + x12*x2z*x12c + x12s* + x2zc + 2*x12c*x2zs + x2zc ) GG21c = GG21c + rs1wg(4)*Xpi*vvc*GtM2*ys * ( 4*O*x12*x2z + 3*O* + x12*x2zs + 5*O*x2z*x12s + 2*O*x12c + Os*x12s + x12*x12c + x12 + *x2zc - x12*x2zs + 5*x2z*x12c + 5*x12s*x2zs + x2zc + 3*x2zs ) GG21c = GG21c + rs1wg(4)*vvc*GtM2*Xps*ys*Xp * ( 4./3. - 1./3.*O + - 1./3.*x12 ) GG21c = GG21c + rs1wg(4)*vvc*GtM2*Xps*ys * ( - 13./6. - 8./3.*O* + x12 + 2./3.*O*x2z + 1./6.*O - 1./3.*x12*x2z + 5./6.*x12 - 23./ + 6.*x2z - 8./3.*x12s ) GG21c = GG21c + rs1wg(4)*vvc*GtM2*ys*Xp * ( 1 + 16./3.*O*x12*x2z + + 31./6.*O*x12 - 1./3.*O*x2z + 8./3.*O*x12s - 1./3.*O*x2zs + - 1./6.*Os*x12 - 1./2.*x12*x2z + 3./2.*x12*x2zs - x12 + 19./ + 3.*x2z*x12s + 19./3.*x2z + 19./6.*x12s + 17./6.*x12c + 11./3. + *x2zs ) GG21c = GG21c + rs1wg(4)*vvc*GtM2*ys * ( - 61./6.*O*x12*x2z - 8./ + 3.*O*x12*x2zs - 2*O*x12 - 14./3.*O*x2z*x12s - 3*O*x12s - 2*O* + x12c + 1./6.*O*x2zs + 1./6.*Os*x12*x2z - Os*x12s + 2*x12*x2z + - x12*x12c - 5./6.*x12*x2zc - 4./3.*x12*x2zs - 43./6.*x2z* + x12s - 29./6.*x2z*x12c - 3*x2z - 14./3.*x12s*x2zs - 3*x12c - + 7./6.*x2zc - 37./6.*x2zs ) GG21c = GG21c + rs1wg(7)*XX*vvc*GtM2*ys * ( 1./2.*O*x12*x2zs - 1./ + 2.*O*x2zc - 1./2.*Os*x12*x2zs + 1./2.*x12*x2zq - 2*x12*x2zc + + x12s*x2zc - 3./2.*x12s*x2zs + 1./2.*x12c*x2zs - 1./2.*x2zq + + 1./2.*x2zc ) GG21c = GG21c + rs1wg(7)*Xpi*vvc*GtM2*ys * ( 1./2.*O*x12*x2z + 3./ + 2.*O*x12*x2zs + 1./2.*O*x2z*x12s - O*x2zs - 1./2.*Os*x12*x2z + + x12*x2zc - 5*x12*x2zs - 3./2.*x2z*x12s + x2z*x12c + 2*x12s + *x2zs - 3./2.*x2zc + x2zs ) GG21c = GG21c + rs1wg(7)*vvc*GtM2*Xps*ys*Xp * ( - 1./6. + 1./6.* + O + 1./6.*x12 ) GG21c = GG21c + rs1wg(7)*vvc*GtM2*Xps*ys * ( - 1./12. - 5./6.*O* + x12 - 1./2.*O*x2z + 1./6.*O - 1./6.*Os - 1./2.*x12*x2z + 2./3. + *x12 + x2z - 2./3.*x12s ) GG21c = GG21c + rs1wg(7)*vvc*GtM2*ys*Xp * ( 17./6.*O*x12*x2z - 5./ + 6.*O*x2z + 2./3.*O*x12s + 1./2.*O*x2zs + 1./6.*Os*x12 + 1./3. + *Os*x2z - 25./6.*x12*x2z + x12*x2zs + 5./2.*x2z*x12s + 1./6.* + x2z - 1./3.*x12s + 1./2.*x12c - 2*x2zs ) GG21c = GG21c + rs1wg(7)*vvc*GtM2*ys * ( - 2*O*x12*x2z - 2*O*x12 + *x2zs - O*x2z*x12s + 1./2.*O*x2z - 1./6.*O*x12c - 1./6.*O* + x2zc + 7./6.*O*x2zs + 1./3.*Os*x12*x2z - 1./12.*Os*x12s - 1./ + 6.*Os*x2zs + 3*x12*x2z - 1./12.*x12*x12c - 7./6.*x12*x2zc + 5 + *x12*x2zs + 1./3.*x2z*x12s - 4./3.*x2z*x12c - 1./2.*x2z - 29./ + 12.*x12s*x2zs + 5./3.*x2zc + 5./12.*x2zs ) GG21c = GG21c + rs1wg(8)*Xpi*vvc*GtM2*ys * ( - 2*O*x12*x2z - 2*O + *x2z*x12s - 2*O*x12c - Os*x12s - x12*x12c - 2*x2z*x12c - x12s + *x2zs - x2zs ) GG21c = GG21c + rs1wg(8)*vvc*GtM2*Xps*ys*Xp * ( - 1 ) GG21c = GG21c + rs1wg(8)*vvc*GtM2*Xps*ys * ( 2 + 3*O*x12 + x12* + x2z - x12 + 2*x2z + 3*x12s ) GG21c = GG21c + rs1wg(8)*vvc*GtM2*ys*Xp * ( - 1 - 3*O*x12*x2z - + 5*O*x12 - 3*O*x12s - x12*x2zs + x12 - 4*x2z*x12s - 4*x2z - 3* + x12s - 3*x12c - x2zs ) GG21c = GG21c + rs1wg(8)*vvc*GtM2*ys * ( 5*O*x12*x2z + 2*O*x12 + + 2*O*x2z*x12s + 3*O*x12s + 2*O*x12c + Os*x12s - x12*x2z + x12* + x12c + x12*x2zs + 4*x2z*x12s + 2*x2z*x12c + 2*x2z + x12s*x2zs + + 3*x12c + 2*x2zs ) GG21c = GG21c + rs1wg(9)*Xpi*vvc*GtM2*ys * ( 2*O*x12*x2z + 2*O* + x12s + 2*x12*x2z - 2*x12*x2zs - 4*x2z*x12s - 2*x12c + 2*x2zs + ) GG21c = GG21c + rs1wg(9)*vvc*GtM2*Xps*ys * ( - 2 + 4*x12 ) GG21c = GG21c + rs1wg(9)*vvc*GtM2*ys*Xp * ( 2 + 2*O*x12 - 6*x12* + x2z - 2*x12 + 4*x2z - 6*x12s ) GG21c = GG21c + rs1wg(9)*vvc*GtM2*ys * ( - 2*O*x12*x2z - 2*O*x12 + - 2*O*x12s + 4*x12*x2z + 2*x12*x2zs - 2*x12 + 4*x2z*x12s - 4 + *x2z + 6*x12s + 2*x12c - 2*x2zs ) GG21c = GG21c + rs1wg(10)*XX*vvc*GtM2*ys * ( - 1./2.*O*x12*x2zs + + 1./2.*O*x2zc + 1./2.*Os*x12*x2zs - 1./2.*x12*x2zq + 2*x12* + x2zc - x12s*x2zc + 3./2.*x12s*x2zs - 1./2.*x12c*x2zs + 1./2.* + x2zq - 1./2.*x2zc ) GG21c = GG21c + rs1wg(10)*Xpi*vvc*GtM2*ys * ( - O*x12*x2z - 3./2. + *O*x12*x2zs - 1./2.*O*x2z*x12s + 3./2.*O*x2zs + Os*x12*x2z - + 3./2.*x12*x2zc + 7*x12*x2zs + 3*x2z*x12s - 3./2.*x2z*x12c - 3 + *x12s*x2zs + 2*x2zc - 3./2.*x2zs ) GG21c = GG21c + rs1wg(10)*vvc*GtM2*Xps*ys*Xp * ( 2./3. - 1./6.*O + - 1./6.*x12 ) GG21c = GG21c + rs1wg(10)*vvc*GtM2*Xps*ys * ( 1./12. + 7./3.*O* + x12 + 1./2.*O*x2z - 2./3.*O + 1./6.*Os + x12*x2z - 11./3.*x12 + - 5./2.*x2z + 13./6.*x12s ) GG21c = GG21c + rs1wg(10)*vvc*GtM2*ys*Xp * ( - 1./2. - 13./3.*O* + x12*x2z - 2*O*x12 + 11./6.*O*x2z - 7./6.*O*x12s - 1./2.*O* + x2zs + 1./2.*O + 1./3.*Os*x12 - 1./3.*Os*x2z + 26./3.*x12*x2z + - 2*x12*x2zs + 3*x12 - 9./2.*x2z*x12s + 1./3.*x2z + 1./3.* + x12s - 3./2.*x12c + 7./2.*x2zs ) GG21c = GG21c + rs1wg(10)*vvc*GtM2*ys * ( 4*O*x12*x2z + 2*O*x12* + x2zs + 1./2.*O*x12 + O*x2z*x12s - 3./2.*O*x2z + 1./2.*O*x12s + + 1./6.*O*x12c + 1./6.*O*x2zc - 5./3.*O*x2zs - 5./6.*Os*x12* + x2z - 1./2.*Os*x12 + 1./12.*Os*x12s + 1./6.*Os*x2zs - 8*x12* + x2z + 1./12.*x12*x12c + 5./3.*x12*x2zc - 6*x12*x2zs + 1./6.* + x2z*x12s + 11./6.*x2z*x12c + 3./2.*x2z + 41./12.*x12s*x2zs - + 3./2.*x12s + x12c - 13./6.*x2zc - 17./12.*x2zs ) GG21c = GG21c + rs1wg(11)*vvc*GtM2*Xps*ys*Xp * ( - 2./3. ) GG21c = GG21c + rs1wg(11)*vvc*GtM2*Xps*ys * ( 1./3. + 4./3.*x12 + + 4./3.*x2z ) GG21c = GG21c + rs1wg(11)*vvc*GtM2*ys*Xp * ( - 1./3.*O*x12 - 5./ + 3.*x12*x2z - 1./3.*x12 - 2./3.*x2z - x12s - 2./3.*x2zs ) GG21c = GG21c + rs1wg(11)*vvc*GtM2*ys * ( 1./3.*O*x12*x2z + 1./3. + *O*x12s + 1./3.*x12*x2z + 1./3.*x12*x2zs + 2./3.*x2z*x12s + 1. + /3.*x12c + 1./3.*x2zs ) GG21c = GG21c + rs1wg(12)*XX*vvc*GtM2*ys * ( - 2*O*x2z*x12s - 2* + O*x12c - 2*x12*x12c - 2*x12*x2zs - 2*x2z*x12s - 4*x2z*x12c - + 2*x12s*x2zs ) GG21c = GG21c + rs1wg(12)*Xpi*vvc*GtM2*ys * ( - 2*O*x12*x2z - 4* + O*x12s - 6*x12*x2z - 2*x12*x2zs - 6*x2z*x12s - 2*x12s - 4* + x12c - 2*x2zs ) GG21c = GG21c + rs1wg(12)*vvc*GtM2*Xps*ys * ( 8./3. - 2./3.*O - 2. + /3.*x12 ) GG21c = GG21c + rs1wg(12)*vvc*GtM2*ys*Xp * ( - 7./3. - 4./3.*O* + x12 + 2./3.*O*x2z + 1./3.*O - 4./3.*x12*x2z - 13./3.*x12 - 5* + x2z - 4./3.*x12s ) GG21c = GG21c + rs1wg(12)*vvc*GtM2*ys * ( 4./3.*O*x12*x2z + 7./3. + *O*x12 - 1./3.*O*x2z + 10./3.*O*x12s - 1./3.*Os*x12 + 26./3.* + x12*x2z + 5./3.*x12*x2zs + 4*x12 + 16./3.*x2z*x12s + 13./3.* + x2z + 13./3.*x12s + 11./3.*x12c + 7./3.*x2zs ) GG21c = GG21c + rs1wg(13)*XX*vvc*GtM2*ys * ( - O*x12*x2zs - 2*O* + x2z*x12s - O*x12c - x12*x12c + x12*x2zc - 2*x12*x2zs - x2z* + x12s - x2z*x12c + x12s*x2zs - x2zc ) GG21c = GG21c + rs1wg(13)*Xpi*vvc*GtM2*ys * ( - 2*O*x12*x2z - 2* + O*x12s - 4*x12*x2z + 2*x12*x2zs - x12s - 2*x12c - 3*x2zs ) GG21c = GG21c + rs1wg(13)*vvc*GtM2*Xps*ys * ( 7./6. ) GG21c = GG21c + rs1wg(13)*vvc*GtM2*ys*Xp * ( - 5./6. - 1./2.*O* + x12 + 1./3.*O*x2z - 1./6.*O + 3./2.*x12*x2z - 3*x12 - 7./2.* + x2z - 1./2.*x12s ) GG21c = GG21c + rs1wg(13)*vvc*GtM2*ys * ( 3./2.*O*x12*x2z + 5./6. + *O*x12 + 1./6.*O*x2z + 11./6.*O*x12s - 1./3.*O*x2zs + 1./6.* + Os*x12 + 23./6.*x12*x2z - 7./3.*x12*x2zs + 2*x12 - 2./3.*x2z* + x12s + 17./6.*x2z + 5./2.*x12s + 5./3.*x12c + 10./3.*x2zs ) GG21c = GG21c + rs1wg(14)*XX*vvc*GtM2*ys * ( O*x12*x2zs + O*x2z* + x12s - x12*x2zc + x12*x2zs - x2z*x12c - 2*x12s*x2zs + x2zc ) GG21c = GG21c + rs1wg(14)*Xpi*vvc*GtM2*ys * ( 2*O*x12*x2z + O* + x12s + 2*x12*x2z - 4*x12*x2zs - 5*x2z*x12s - x12c + 3*x2zs ) GG21c = GG21c + rs1wg(14)*vvc*GtM2*Xps*ys*Xp * ( - 1./3. ) GG21c = GG21c + rs1wg(14)*vvc*GtM2*Xps*ys * ( - 1 + 1./3.*O + 3* + x12 + 2./3.*x2z ) GG21c = GG21c + rs1wg(14)*vvc*GtM2*ys*Xp * ( 1 + 1./3.*O*x12 - 2./ + 3.*O*x2z - 19./3.*x12*x2z - x12 + 3*x2z - 4*x12s - 1./3.*x2zs + ) GG21c = GG21c + rs1wg(14)*vvc*GtM2*ys * ( - 4./3.*O*x12*x2z - O* + x12 - 2./3.*O*x12s + 1./3.*O*x2zs + 3*x12*x2z + 13./3.*x12* + x2zs - x12 + 17./3.*x2z*x12s - 3*x2z + 3*x12s + 4./3.*x12c - + 3*x2zs ) GG12c = 0 GG22c = 0 *---------------------------------------------------------------- *#call wtensor{gz} *#call totalgz{dummy} *#call contractgz{dummy} *---------------------------------------------------------------- *#call wtensor{zz} *#call totalzz{dummy} *#call contractzz{dummy} * .end Time = 80.35 sec Generated terms = 0 GG10q Terms in output = 0 Bytes used = 2 Time = 80.37 sec Generated terms = 0 GG20q Terms in output = 0 Bytes used = 2 Time = 80.38 sec Generated terms = 0 GG11q Terms in output = 0 Bytes used = 2 Time = 80.40 sec Generated terms = 0 GG21q Terms in output = 0 Bytes used = 2 Time = 80.40 sec Generated terms = 0 GG12q Terms in output = 0 Bytes used = 2 Time = 80.40 sec Generated terms = 0 GG22q Terms in output = 0 Bytes used = 2 GG10q = 0 GG20q = 0 GG11q = 0 GG21q = 0 GG12q = 0 GG22q = 0 GIAMPIERO job terminated at 2-DEC-1999 15:59:06.63 Accounting information: Buffered I/O count: 828 Peak working set size: 7104 Direct I/O count: 1640 Peak page file size: 100560 Page faults: 582 Mounted volumes: 0 Charged CPU time: 0 00:01:21.20 Elapsed time: 0 00:01:39.32