fermiqcd_fermilab_coefficients.h Source File

00001 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 // include clover + c1,c2,c3,c4,c5,c10
00014 
00015  mdp_matrix M(4,4), SE[4], SB[4];
00016 
00017 
00018  M=1.0+(2.0*42*c3+2.0*18*c4)*kappat+Gamma[0]*2.0*kappat*(mac1-2.0*mac2+6.0*mac4);
00019  for(a=0;a<4; a++)
00020    for(b=0; b<4; b++) 
00021      c_id[4*a+b]=M(b,a);
00022  
00023  M=kappat*Gamma[0]-rt*kappat+Gamma[0]*2.0*kappat*(mac2-mac3-4.0*mac4);
00024  for(a=0;a<4; a++)
00025    for(b=0; b<4; b++) 
00026      c_u[16*a+4*b]=M(b,a);
00027  
00028  M=-kappat*Gamma[0]-rt*kappat+Gamma[0]*2.0*kappat*(mac2+mac3-4.0*mac4);
00029  for(a=0;a<4; a++)
00030    for(b=0; b<4; b++) 
00031      c_d[16*a+4*b]=M(b,a);
00032 
00033  M=Gamma[0]*2.0*kappat*(0.5*mac3+mac4);
00034  for(a=0;a<4; a++)
00035    for(b=0; b<4; b++) 
00036      c_uu[64*a+16*b+4*0+0]=M(b,a);
00037 
00038  M=Gamma[0]*2.0*kappat*(-0.5*mac3+mac4);
00039  for(a=0;a<4; a++)
00040    for(b=0; b<4; b++) 
00041      c_dd[64*a+16*b+4*0+0]=M(b,a);
00042 
00043  for(i=1; i<4; i++) {
00044    switch(i) {
00045    case 1: j=2; k=3; break;
00046    case 2: j=1; k=3; break;
00047    case 3: j=1; k=2; break;
00048    }   
00049    SE[i]=(-I)*Sigma[0][i];
00050    SB[i]=(-I)*Sigma[j][k];
00051 
00052    M=(zeta-12.0*c1-2.0*c2)*kappat*Gamma[i]+(-rs*zeta+(-24.0)*c3+(-8.0)*c4)*kappat;
00053    for(a=0;a<4; a++)
00054      for(b=0; b<4; b++)
00055        c_u[16*a+4*b+i]=M(b,a);
00056    
00057    M=(-zeta+12.0*c1+2.0*c2)*kappat*Gamma[i]+(-rs*zeta+(-24.0)*c3+(-8.0)*c4)*kappat;
00058    for(a=0;a<4; a++)
00059      for(b=0; b<4; b++) 
00060        c_d[16*a+4*b+i]=M(b,a);
00061        
00062    M=-cE*cSW*zeta*kappat*SE[i];
00063      for(a=0;a<4; a++)
00064        for(b=0; b<4; b++) 
00065          c_e[16*a+4*b+i]=M(b,a);
00066    
00067    M=-cB*cSW*zeta*kappat*SB[i]-16*c10*kappat*SB[i];
00068    for(a=0;a<4; a++)
00069      for(b=0; b<4; b++) 
00070        c_b[16*a+4*b+i]=M(b,a);
00071  } 
00072 
00073  for(i=1; i<4; i++)
00074    for(j=1; j<4; j++) {
00075 
00076      M=c1*kappat*(Gamma[i]+Gamma[j])
00077        +c2*kappat*Gamma[i]*delta(i,j)
00078        +2.0*c3*kappat
00079        +2.0*c4*kappat*delta(i,j);
00080      for(a=0;a<4; a++)
00081        for(b=0; b<4; b++) 
00082          c_uu[64*a+16*b+4*i+j]=M(b,a);
00083 
00084      M=c1*kappat*(Gamma[i]-Gamma[j])
00085        +c2*kappat*Gamma[i]*delta(i,j)
00086        +2.0*c3*kappat
00087        +2.0*c4*kappat*delta(i,j);
00088      for(a=0;a<4; a++)
00089        for(b=0; b<4; b++) 
00090          c_ud[64*a+16*b+4*i+j]=M(b,a);
00091 
00092      M=-c1*kappat*(Gamma[i]-Gamma[j])
00093        -c2*kappat*Gamma[i]*delta(i,j)
00094        +2.0*c3*kappat
00095        +2.0*c4*kappat*delta(i,j);
00096      for(a=0;a<4; a++)
00097        for(b=0; b<4; b++) 
00098          c_du[64*a+16*b+4*i+j]=M(b,a);
00099 
00100      M=-c1*kappat*(Gamma[i]+Gamma[j])
00101        -c2*kappat*Gamma[i]*delta(i,j)
00102        +2.0*c3*kappat
00103        +2.0*c4*kappat*delta(i,j);
00104      for(a=0;a<4; a++)
00105        for(b=0; b<4; b++) 
00106          c_dd[64*a+16*b+4*i+j]=M(b,a);
00107 
00108      M=c5*kappat*SB[i]*Gamma[j]+(1.0-delta(i,j))*2.0*c10*SB[i];
00109      for(a=0;a<4; a++)
00110        for(b=0; b<4; b++) 
00111          c_bu[64*a+16*b+4*i+j]=M(b,a);
00112 
00113      M=-c5*kappat*SB[i]*Gamma[j]+(1.0-delta(i,j))*2.0*c10*SB[i];
00114      for(a=0;a<4; a++)
00115        for(b=0; b<4; b++) 
00116          c_bd[64*a+16*b+4*i+j]=M(b,a);
00117 
00118      M=c5*kappat*Gamma[i]*SB[j]+(1.0-delta(i,j))*2.0*c10*SB[j];
00119      for(a=0;a<4; a++)
00120        for(b=0; b<4; b++) 
00121          c_ub[64*a+16*b+4*i+j]=M(b,a);
00122 
00123      M=-c5*kappat*Gamma[i]*SB[j]+(1.0-delta(i,j))*2.0*c10*SB[j];
00124      for(a=0;a<4; a++)
00125        for(b=0; b<4; b++) 
00126          c_bu[64*a+16*b+4*i+j]=M(b,a);
00127    }
00128  
00129