*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Sun Mar 6 15:34:04 2005 * * This FORM code processes the result produced by the (PQint) code: * * Simplifications and sums over sea and valence quark indices are * * implemented. Includes the possiblilty of producing output intended for * * fortran functions and LaTeX. * * ************************************************************************** * * comp 1 : Simplification of the pseudoscalar meson mass to NNLO * * comp 2 : Simplification of the pseudoscalar meson decay constant to NNLO * * #define comp "1" * #define nfl "2" * * * The parameters (mode) and (type) select the values of (dval) and (dsea), * * according to the tables below: * * * nfl = 2: * * dval + dsea Mode Type * ***************************** * * 1 + 1 1 1 * * 1 + 2 2 1 * * 2 + 1 4 1 * * 2 + 2 5 1 * * * nfl = 3: * * dval + dsea Mode Type * ***************************** * * 1 + 1 2 1 * * 1 + 2 4 1 * * 1 + 3 5 1 * * 2 + 1 2 2 * * 2 + 2 4 2 * * 2 + 3 5 2 * * #define mode "1" * #define type "1" * * #define pout "0" * #define nsim "0" * * #- * Sigma2 loaded * Sigma4 loaded * Sigma6 loaded * *Time = 0.04 sec Generated terms = 1 * M0p2 Terms in output = 1 * Bytes used = 24 * *Time = 0.04 sec Generated terms = 6 * M0p4 Terms in output = 6 * Bytes used = 226 * *Time = 0.04 sec Generated terms = 100 * M0p6 Terms in output = 100 * Bytes used = 3660 L M0p2timo11 = + X11 ; L M0p4timo11 = + F0^-2*L4r * ( - 16*X11*X33 ) + F0^-2*L5r * ( - 8*X11^2 ) + F0^-2*L6r * ( + 32*X11*X33 ) + F0^-2*L8r * ( + 16*X11^2 ) + Ab(X11)*F0^-2 * ( - 1/2*X11 ) + Bb(X11,X11,0)*F0^-2 * ( - 1/2*RD(X11)*X11 ); L M0p6timo11 = + F0^-4*pi16*L0r * ( + 2*X11*X33^2 + 8*X11^2*X33 - 4*X11^3 ) + F0^-4*pi16*L1r * ( + 4*X11^3 ) + F0^-4*pi16*L2r * ( + 6*X11*X33^2 + 2*X11^3 ) + F0^-4*pi16*L3r * ( + X11*X33^2 + 6*X11^2*X33 - 5*X11^3 ) + F0^-4*pi16^2 * ( + 47/96*X11*X33^2 + 35/144*X11^2*X33 - 17/288*X11^3 ) + F0^-4*L4r*L5r * ( + 256*X11^2*X33 ) + F0^-4*L4r*L6r * ( - 512*X11*X33^2 ) + F0^-4*L4r*L8r * ( - 256*X11^2*X33 ) + F0^-4*L4r*L11r * ( + 128*X11^2*X33 ) + F0^-4*L4r^2 * ( + 256*X11*X33^2 ) + F0^-4*L5r*L6r * ( - 256*X11^2*X33 ) + F0^-4*L5r*L8r * ( - 128*X11^3 ) + F0^-4*L5r*L11r * ( + 64*X11^3 ) + F0^-4*L5r^2 * ( + 64*X11^3 ) + F0^-4*L6r*L11r * ( - 256*X11^2*X33 ) + F0^-4*L8r*L11r * ( - 128*X11^3 ) + F0^-4*K17r * ( - 32*X11^3 ) + F0^-4*K18r * ( - 64*X11^2*X33 ) + F0^-4*K19r * ( - 16*X11^3 ) + F0^-4*K20r * ( - 32*X11^2*X33 ) + F0^-4*K21r * ( - 32*X11*X33^2 ) + F0^-4*K22r * ( - 64*X11*X33^2 ) + F0^-4*K23r * ( - 16*X11^3 ) + F0^-4*K25r * ( + 48*X11^3 ) + F0^-4*K26r * ( + 32*X11*X33^2 + 64*X11^2*X33 ) + F0^-4*K27r * ( + 192*X11*X33^2 ) + F0^-4*K39r * ( + 32*X11^3 ) + F0^-4*K40r * ( + 64*X11^2*X33 ) + Ab(X11)*F0^-4*L0r * ( - 12*X11^2 - 12*RD(X11)*X11 ) + Ab(X11)*F0^-4*L1r * ( + 8*X11^2 ) + Ab(X11)*F0^-4*L2r * ( + 20*X11^2 ) + Ab(X11)*F0^-4*L3r * ( - 12*X11^2 - 12*RD(X11)*X11 ) + Ab(X11)*F0^-4*L4r * ( + 16*X11*X33 ) + Ab(X11)*F0^-4*L5r * ( + 16*X11^2 + 8*RD(X11)*X11 ) + Ab(X11)*F0^-4*L6r * ( + 16*X11^2 + 16*RD(X11)*X11 ) + Ab(X11)*F0^-4*L7r * ( + 32*RD(X11)*X11 ) + Ab(X11)*F0^-4*L8r * ( - 32*X11^2 ) + Ab(X11)^2*F0^-4 * ( + 5/8*X11 ) + Ab(X11)*Bb(X11,X11,0)*F0^-4 * ( + 3/2*X11^2 + 1/4*RD(X11)*X11 ) + Ab(X11)*Cb(X11,X11,X11,0)*F0^-4 * ( + 1/2*RD(X11)*X11^2 ) + Ab(X11,eps)*F0^-4*pi16 * ( - 7/4*X11^2 + 1/4*RD(X11)*X11 ) + Ab(X13)*F0^-4*pi16 * ( + 4/3*X11*X33 ) + Ab(X13)*F0^-4*L0r * ( + 16*X11*X13 ) + Ab(X13)*F0^-4*L3r * ( + 40*X11*X13 ) + Ab(X13)*F0^-4*L5r * ( - 32*X11*X13 ) + Ab(X13)*F0^-4*L8r * ( + 64*X11*X13 ) + Ab(X13)^2*F0^-4 * ( - X11 ) + Ab(X13)*Bb(X11,X11,0)*F0^-4 * ( - 2*X11*X33 ) + Ab(X13,eps)*F0^-4*pi16 * ( - 2*X11*X33 ) + Ab(X33)*F0^-4*L1r * ( + 48*X11*X33 ) + Ab(X33)*F0^-4*L2r * ( + 12*X11*X33 ) + Ab(X33)*F0^-4*L4r * ( - 48*X11*X33 ) + Ab(X33)*F0^-4*L6r * ( + 48*X11*X33 ) + Ab(X33)*Bb(X11,X11,0)*F0^-4 * ( + 3/4*X11*X33 ) + Ab(X33,eps)*F0^-4*pi16 * ( - 3/4*X11*X33 ) + Bb(X11,X11,0)*F0^-4*L0r * ( - 12*RD(X11)*X11^2 ) + Bb(X11,X11,0)*F0^-4*L3r * ( - 12*RD(X11)*X11^2 ) + Bb(X11,X11,0)*F0^-4*L4r * ( - 24*X11*X33^2 + 32*X11^2*X33 ) + Bb(X11,X11,0)*F0^-4*L5r * ( + 4*X11^3 + 24*RD(X11)*X11^2 ) + Bb(X11,X11,0)*F0^-4*L6r * ( + 32*X11*X33^2 - 48*X11^2*X33 ) + Bb(X11,X11,0)*F0^-4*L7r * ( + 16*RD(X11)^2*X11 ) + Bb(X11,X11,0)*F0^-4*L8r * ( - 8*X11^3 - 48*RD(X11)*X11^2 + 8*RD(X11)^2*X11 ) + Bb(X11,X11,0)^2*F0^-4 * ( + 1/2*RD(X11)*X11^2 + 1/8*RD(X11)^2*X11 ) + Bb(X11,X11,0)*Cb(X11,X11,X11,0)*F0^-4 * ( + 1/2*RD(X11)^2*X11^2 ) + Bb(X11,X11,0,eps)*F0^-4*pi16 * ( + 1/4*RD(X11)*X11^2 ) + Cb(X11,X11,X11,0)*F0^-4*L4r * ( + 16*RD(X11)*X11^2*X33 ) + Cb(X11,X11,X11,0)*F0^-4*L5r * ( + 8*RD(X11)*X11^3 ) + Cb(X11,X11,X11,0)*F0^-4*L6r * ( - 32*RD(X11)*X11^2*X33 ) + Cb(X11,X11,X11,0)*F0^-4*L8r * ( - 16*RD(X11)*X11^3 ) + Hbb(1,X11,X11,X11,X11)*F0^-4 * ( + 5/6*X11^2 ) + Hbb(1,X11,X13,X13,X11)*F0^-4 * ( + 1/4*X11*X33 - X11^2 ) + Hbb(1,X13,X13,X33,X11)*F0^-4 * ( + 3/4*X11*X33 ) + Hbb(2,X11,X11,X11,X11)*F0^-4 * ( + RD(X11)*X11^2 ) + Hbb(2,X11,X13,X13,X11)*F0^-4 * ( + 3/4*RD(X11)*X11^2 ) + Hbb(5,X11,X11,X11,X11)*F0^-4 * ( + 1/2*RD(X11)^2*X11^2 ) + H1bb(3,X13,X11,X13,X11)*F0^-4 * ( - 4*RD(X11)*X11^2 ) + H21bb(1,X11,X13,X13,X11)*F0^-4 * ( + 3/4*X11^2 ) + H21bb(1,X33,X13,X13,X11)*F0^-4 * ( + 9/4*X11^2 ) + H21bb(2,X11,X13,X13,X11)*F0^-4 * ( - 3/4*RD(X11)*X11^2 );