if(g_mu == 0.){
if(g_use_clover_flag == 1){
int solve_cg(int k,int l, double q_off, double eps_sq) {  static double normsq,pro,err,alpha_cg,g_beta_cg;  int iteration;    /* initialize residue r and search vector p */    if(g_mu == 0.){    Q_psi(DUM_SOLVER,k,q_off);    Q_psi(DUM_SOLVER,DUM_SOLVER,q_off);  }  else{    Qtm_pm_psi(DUM_SOLVER, k);  }  diff_field(DUM_SOLVER+1,l,DUM_SOLVER);  assign_field(DUM_SOLVER+2,DUM_SOLVER+1);  normsq=square_norm(DUM_SOLVER+1);/* main loop */   for(iteration=1;iteration<=ITER_MAX;iteration++) {     if(g_mu == 0.){       Q_psi(DUM_SOLVER,DUM_SOLVER+2,q_off);       Q_psi(DUM_SOLVER,DUM_SOLVER,q_off);     }     else {       Qtm_pm_psi(DUM_SOLVER, DUM_SOLVER+2);     }     pro=vprod(DUM_SOLVER+2,DUM_SOLVER);     alpha_cg=normsq/pro;     add_assign_field(k, alpha_cg, DUM_SOLVER+2);     add_assign_field2(DUM_SOLVER, -alpha_cg, DUM_SOLVER+1);     err=square_norm(DUM_SOLVER);     if (err <= eps_sq){       break;     }     g_beta_cg = err/normsq;     add_assign_field2(DUM_SOLVER+2, g_beta_cg, DUM_SOLVER);     assign_field(DUM_SOLVER+1, DUM_SOLVER);     normsq=err;   }   return iteration;}