if (g_proc_id == 0) {
if (g_proc_coords[0] == 0 && g_proc_coords[1]==0) {
void source_spinor_field(spinor * const P, spinor * const Q, int is, int ic) {  int ix;  spinor * s;  zero_spinor_field(P);  zero_spinor_field(Q);  if (g_proc_id == 0) {    s = P;        /* put source to 1.0 */    if (is==0){      if      (ic==0) (*s).s0.c0.re=1.0;      else if (ic==1) (*s).s0.c1.re=1.0;      else if (ic==2) (*s).s0.c2.re=1.0;    }    else if (is==1){      if      (ic==0) (*s).s1.c0.re=1.0;      else if (ic==1) (*s).s1.c1.re=1.0;      else if (ic==2) (*s).s1.c2.re=1.0;    }    else if (is==2){      if      (ic==0) (*s).s2.c0.re=1.0;      else if (ic==1) (*s).s2.c1.re=1.0;      else if (ic==2) (*s).s2.c2.re=1.0;    }    else if (is==3){      if      (ic==0) (*s).s3.c0.re=1.0;      else if (ic==1) (*s).s3.c1.re=1.0;      else if (ic==2) (*s).s3.c2.re=1.0;    }  }}