if(g_mu == 0.){
if(g_use_clover_flag == 1){
void deri(double q_off,double q_off2) {  int j,jmax=1;  int i,mu;  double qo;  for(i=0;i<(VOLUME+RAND);i++){     for(mu=0;mu<4;mu++){       _zero_su3adj(df0[i][mu]);     }  }  if(g_mu == 0.){    /* All the clover stuff */    for(i=0;i<(VOLUME+RAND);i++){       for(mu=0;mu<4;mu++){ 	_zero_su3adj(dclover[i][mu]);      }    }        for(i=0;i<VOLUME;i++){       for(mu=0;mu<4;mu++){ 	_su3_zero(swm[i][mu]); 	_su3_zero(swp[i][mu]);       }    }    if(q_off==0.){      jmax=1;    }     else{      jmax=2;    }    if(q_off2>0.){      jmax=3;    }  }  for(j=0;j<jmax;j++){     if(j==0){      /*contributions from field 0 */      gamma5(DUM_DERI,0);      /* Invert first Q_+ */      count00+=bicg(DUM_DERI,0,q_off,EPS_SQ1);      gamma5(DUM_DERI+1,DUM_DERI);      /* Now Q_- */      g_mu = -g_mu;      count01+=bicg(DUM_DERI+1,DUM_DERI,q_off,EPS_SQ1);      g_mu = -g_mu;    }    if(j==1){      /* contributions from field 1 */      gamma5(DUM_DERI,1);      qo=q_off-q_off2;      count10+=bicg(DUM_DERI,1,q_off2,EPS_SQ2/qo);      deri_linalg(DUM_DERI+2,qo*qo,DUM_DERI,qo,1);      gamma5(DUM_DERI+1,DUM_DERI+2);      count11+=bicg(DUM_DERI+1,DUM_DERI+2,q_off2,EPS_SQ2*qo);    }    if(j==2){      /* contributions from field 2 (stored on 4) */      gamma5(DUM_DERI,4);      count20+=bicg(DUM_DERI,4,0.,EPS_SQ3/q_off2);      deri_linalg(DUM_DERI+2,q_off2*q_off2,DUM_DERI,q_off2,4);      gamma5(DUM_DERI+1,DUM_DERI+2);      count21+=bicg(DUM_DERI+1,DUM_DERI+2,0.,EPS_SQ3*q_off2);    }    if(g_mu == 0.){      /* apply H_eo to  Q^{-2} phi */      H_eo_psi(1,DUM_DERI+2,DUM_DERI+1);      /* result resides on odd sites */            deriv_Sb(0,DUM_DERI,DUM_DERI+2);            /* add the other contibution */      H_eo_psi(1,DUM_DERI+3,DUM_DERI);      /* includes (1+T_oo)^{-1} now */      /* apply (1+T_oo)^{-1} to the result !!!! */      deriv_Sb(1,DUM_DERI+3,DUM_DERI+1);            /* add the contribution from inside */      gamma5(DUM_DERI+2,DUM_DERI+2);      sw_spinor(1,DUM_DERI+2,DUM_DERI+3);            /* compute the contribution for the det-part */      gamma5(DUM_DERI,DUM_DERI);      sw_spinor(0,DUM_DERI,DUM_DERI+1);    }    else{      /* apply Hopping Matrix */      /* to get the even sites of X */      /* Check the signs here ! */      H_eo_tm_inv_psi(DUM_DERI+2, DUM_DERI+1, 1, -1.);      deriv_Sb(0, DUM_DERI, DUM_DERI+2);      /* to get the odd sites of X */      H_eo_tm_inv_psi(DUM_DERI+3, DUM_DERI, 1, +1);      deriv_Sb(1, DUM_DERI+3, DUM_DERI+1);    }  }}