*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Mon Oct 4 00:51:02 2004 * * This program loads the result for the finite part of the meson mass * * produced by PQint.frm and simplifies it as much as possible. Should * * also include a facility to print out the result it fortran format. * * * NB! Definition of mode must be as in the save file: * * * Mode 1 : full 3+3 theory, all quarks nondegenerate * * Mode 2 : 1+1 theory, all valence and all sea quarks degenerate * * Mode 3 : unquenched SU(2) symmetric limit * * Mode 4 : SU(2) + 2 theory, sea quarks 4 and 6 nondegenerate * * Mode 5 : SU(2) + 3 theory, all sea quarks nondegenerate * * * Choice of quark indices to be assigned to the external legs: * * The choices available are (in the notation ijkl): * * * Type 1 : pi^+/- meson - 1221 * * Type 2 : K^+/- meson - 1331 * * #define comp "2" * #define mode "5" * #define type "1" * * #- * Sigma2 loaded * Sigma4 loaded * Sigma6 loaded * *Time = 0.17 sec Generated terms = 1 * F0p2 Terms in output = 1 * Bytes used = 24 * *Time = 0.17 sec Generated terms = 3 * F0p4 Terms in output = 3 * Bytes used = 98 * *Time = 0.17 sec Generated terms = 108 * F0p6 Terms in output = 108 * Bytes used = 4516 L F0p2timo13 = + F0 * ( + 1 ); L F0p4timo13 = + F0^-1*L4r * ( + 12*Xb(1) ) + F0^-1*L5r * ( + 4*X11 ) + Ab(X1s)*F0^-1 * ( + 1/2 ); L F0p6timo13 = + F0^-3*pi16*L0r * ( + 4/9*Xp0*Xe0 + 1/2*X11^2 - 13/3*Xb(1)*X11 - 35/18*Xb(2) ) + F0^-3*pi16*L1r * ( - 2*X11^2 ) + F0^-3*pi16*L2r * ( - 11/3*Xp0*Xe0 - X11^2 - 13/3*Xb(2) ) + F0^-3*pi16*L3r * ( + 4/9*Xp0*Xe0 + 5/4*X11^2 - 17/6*Xb(1)*X11 - 43/36*Xb(2) ) + F0^-3*pi16^2 * ( - 15/64*Xp0*Xe0 + 1/64*X11^2 - 1/2*Xb(1)*X11 - 43/128*Xb(2) ) + F0^-3*L4r*L5r * ( - 48*Xb(1)*X11 ) + F0^-3*L4r^2 * ( - 72*Xb(1)^2 ) + F0^-3*L5r^2 * ( - 8*X11^2 ) + F0^-3*K19r * ( + 8*X11^2 ) + F0^-3*K20r * ( + 24*Xb(1)*X11 ) + F0^-3*K21r * ( + 24*Xb(2) ) + F0^-3*K22r * ( + 72*Xb(1)^2 ) + F0^-3*K23r * ( + 8*X11^2 ) + Ab(Xn1)*F0^-3*L0r * ( + 4*R(Xn1,Xn2,X11,X11)*Xn1 ) + Ab(Xn1)*F0^-3*L1r * ( + 8*Xb(Xn1,Xn2,0)*Xn1 ) + Ab(Xn1)*F0^-3*L2r * ( + 2*Xb(Xn1,Xn2,0)*Xn1 ) + Ab(Xn1)*F0^-3*L3r * ( + 4*R(Xn1,Xn2,X11,X11)*Xn1 ) + Ab(Xn1)*F0^-3*L4r * ( - 4*Xb(Xn1,Xn2,1) ) + Ab(Xn1)*F0^-3*L5r * ( - 4/3*R(Xn1,Xn2,X11,X11)*X11 ) + Ab(Xn1)*Bb(X1s,X1s,0)*F0^-3 * ( - 1/6*R(Xn1,Xn2,X11,Xss)*X1s ) + Ab(Xn1,eps)*F0^-3*pi16 * ( - 7/24*Xb(Xn1,Xn2,0)*Xn1 + 1/6*Xb(Xn1,Xn2,1) - 1/8*R(Xn1,Xn2,X11,X11)*Xn1 ) + Ab(X11)*F0^-3*L0r * ( + 4*RC(X11)*X11 + 4*RD(X11) ) + Ab(X11)*F0^-3*L1r * ( - 4*X11 ) + Ab(X11)*F0^-3*L2r * ( - 10*X11 ) + Ab(X11)*F0^-3*L3r * ( + 4*RC(X11)*X11 + 4*RD(X11) ) + Ab(X11)*F0^-3*L5r * ( - 4/3*RC(X11)*X11 ) + Ab(X11)*Bb(X1s,X1s,0)*F0^-3 * ( - 1/6*R(X11,Xss,Xp0,Xe0)*X1s ) + Ab(X11,eps)*F0^-3*pi16 * ( + 1/4*X11 - 1/8*RC(X11)*X11 - 1/8*RD(X11) ) + Ab(X1s)*F0^-3*pi16 * ( - 1/4*X1s - 3/8*Xb(1) ) + Ab(X1s)*F0^-3*L0r * ( - 4*X1s ) + Ab(X1s)*F0^-3*L3r * ( - 10*X1s ) + Ab(X1s)*F0^-3*L4r * ( - 6*Xb(1) ) + Ab(X1s)*F0^-3*L5r * ( + 2*X11 ) + Ab(X1s,eps)*F0^-3*pi16 * ( + 3/8*Xss + 3/8*Xb(1) ) + Ab(Xss)*F0^-3*L1r * ( - 8*Xss + 8/3*RC(Xss)*Xss ) + Ab(Xss)*F0^-3*L2r * ( - 2*Xss + 2/3*RC(Xss)*Xss ) + Ab(Xss)*F0^-3*L4r * ( + 4*Xss - 4/3*RC(Xss)*Xss ) + Ab(Xss,eps)*F0^-3*pi16 * ( + 1/8*Xss - 1/24*RC(Xss)*Xss ) + Ab(Xst)*F0^-3*L1r * ( - 16*Xst ) + Ab(Xst)*F0^-3*L2r * ( - 4*Xst ) + Ab(Xst)*F0^-3*L4r * ( + 8*Xst ) + Ab(Xst,eps)*F0^-3*pi16 * ( + 1/4*Xst ) + Bb(X11,X11,0)*F0^-3*L0r * ( + 4*RD(X11)*X11 ) + Bb(X11,X11,0)*F0^-3*L3r * ( + 4*RD(X11)*X11 ) + Bb(X11,X11,0)*F0^-3*L5r * ( - 4/3*RD(X11)*X11 ) + Bb(X11,X11,0,eps)*F0^-3*pi16 * ( - 1/8*RD(X11)*X11 ) + Bb(X1s,X1s,0)*F0^-3*L4r * ( - 12*Xb(1)*X1s ) + Bb(X1s,X1s,0)*F0^-3*L5r * ( - 4*X1s^2 ) + Bb(X1s,X1s,0)*F0^-3*L6r * ( + 24*Xb(1)*X1s ) + Bb(X1s,X1s,0)*F0^-3*L8r * ( + 8*X1s^2 ) + Hbb(1,Xn1,X1s,X1s,X11)*F0^-3 * ( + 1/24*RV(Xn1,Xn2,X11,Xss)*Xn1 ) + Hbb(1,X11,X1s,X1s,X11)*F0^-3 * ( - 1/12*R(X11,Xss,Xp0,Xe0)*X11 + 1/24*RC(X11)*X11 + 1/24*RD(X11) ) + Hbb(1,X1s,X1s,Xss,X11)*F0^-3 * ( - 1/8*Xss + 1/24*RC(Xss)*Xss ) + Hbb(1,X1s,X1t,Xst,X11)*F0^-3 * ( - 1/4*Xst ) + Hbb(2,X11,X1s,X1s,X11)*F0^-3 * ( + 1/24*RD(X11)*X11 ) + dHbb(1,Xn1,Xn1,X11,X11)*F0^-3 * ( + 1/9*R(Xn1,Xn2,X11,X11)^2*X11^2 ) + dHbb(1,Xn1,X11,X11,X11)*F0^-3 * ( + 2/9*R(Xn1,Xn2,X11,X11)*RC(X11)*X11^2 ) + dHbb(1,Xn1,X1s,X1s,X11)*F0^-3 * ( - 1/6*R(Xn1,Xn2,X11,X11)*X11^2 - 1/24*RV(Xn1,Xn2,X11,Xss)*Xn1*X11 ) + dHbb(1,Xp0,Xe0,X11,X11)*F0^-3 * ( + 2/9*R(Xp0,Xe0,X11,X11)*R(Xe0,Xp0,X11,X11)*X11^2 ) + dHbb(1,X11,X11,X11,X11)*F0^-3 * ( + 1/6*X11^2 + 1/9*RC(X11)^2*X11^2 ) + dHbb(1,X11,X1s,X1s,X11)*F0^-3 * ( - 1/4*R(X11,Xss,Xp0,Xe0)*X11^2 + 1/8*RC(X11)*X11^2 - 1/24*RD(X11)*X11 ) + dHbb(1,X1s,X1s,Xss,X11)*F0^-3 * ( + 1/8*X11*Xss - 1/24*RC(Xss)*X11*Xss ) + dHbb(1,X1s,X1t,Xst,X11)*F0^-3 * ( + 1/4*X11*Xst ) + dHbb(2,X11,Xn1,X11,X11)*F0^-3 * ( + 2/9*R(Xn1,Xn2,X11,X11)*RD(X11)*X11^2 ) + dHbb(2,X11,X11,X11,X11)*F0^-3 * ( + 2/9*RC(X11)*RD(X11)*X11^2 ) + dHbb(2,X11,X1s,X1s,X11)*F0^-3 * ( + 1/8*RD(X11)*X11^2 ) + dHbb(5,X11,X11,X11,X11)*F0^-3 * ( + 1/9*RD(X11)^2*X11^2 ) + dH1bb(1,Xn1,X1s,X1s,X11)*F0^-3 * ( - 1/3*R(Xn1,Xn2,X11,Xss)*RS(X11,Xss,Xn1)*X11^2 ) + dH1bb(1,X1s,X1s,X11,X11)*F0^-3 * ( - 2/3*R(Xn1,Xn2,X11,Xss)*RS(X11,Xss,Xn1)*X11^2 ) + dH1bb(3,X1s,X11,X1s,X11)*F0^-3 * ( - 2/3*RD(X11)*X11^2 ) + dH21bb(1,Xn1,X1s,X1s,X11)*F0^-3 * ( - 1/8*RV(Xn1,Xn2,X11,Xss)*X11^2 ) + dH21bb(1,X11,X1s,X1s,X11)*F0^-3 * ( + 1/4*R(X11,Xss,Xp0,Xe0)*X11^2 - 1/8*RC(X11)*X11^2 ) + dH21bb(1,Xss,X1s,X1s,X11)*F0^-3 * ( + 3/8*X11^2 - 1/8*RC(Xss)*X11^2 ) + dH21bb(1,Xst,X1s,X1t,X11)*F0^-3 * ( + 3/4*X11^2 ) + dH21bb(2,X11,X1s,X1s,X11)*F0^-3 * ( - 1/8*RD(X11)*X11^2 );