fr_print_pre_mic()

void fr_print_pre_mic	(	FILE *	file,
		int	imic,
		int	icounter,
		COMPLEX	cfreq
	)

                                                                     {

  Vec gPhiscm = NULL;
  Vec gPhincm = NULL;
  int locRows;
  int gRows;
  int mm,nn;
  int descl[9];
  int descg[9];
  int izero=0;
  int ione=1;
  int info;
  COMPLEX rho_s;


  logger(LOG_DEBUG,"Collecting solution at microphones to print frequency response file...\n");
  rho_s = cfreq*modsolinfo->rho;

  if (0 == rank) gPhiscm  = calloc(modgeominfo->nmics,sizeof(COMPLEX));
  if (0 == rank) gPhincm  = calloc(modgeominfo->nmics,sizeof(COMPLEX));

  gRows = modgeominfo->nmics;
  locRows  = numroc_(&gRows  ,&runinfo->mics_block_size,&runinfo->myrow,&izero,&runinfo->nprows);


  // Gathering vectors on node 0 
  mm = gRows;// * runinfo->nprows;
  nn = runinfo->npcols;
  descinit_(descl   , &gRows  ,&ione, &runinfo->mics_block_size ,&ione,&izero,&izero,&runinfo->ctxt,&gRows  ,&info);
  descinit_(descg   , &mm     ,&ione, &gRows                    ,&ione,&izero,&izero,&runinfo->ctxt,&gRows  ,&info);

  pzgemr2d_(&gRows,&ione,solution->vecPhium  ,&ione, &ione,descl, gPhiscm  ,&ione,&ione,descg,&runinfo->ctxt);
  pzgemr2d_(&gRows,&ione,solution->vecPhincm ,&ione, &ione,descl, gPhincm ,&ione,&ione,descg,&runinfo->ctxt);

  if(0 == rank){
  fprintf(file,"%d %10.8E %10.8E %10.8E %10.8E %10.8E %10.8E %10.8E %10.8E %10.8E %10.8E %10.8E\n",
        icounter,
        CIMAG(cfreq),
        HZFREQ(cfreq),
        CREAL(gPhincm[imic]),
        CIMAG(gPhincm[imic]),
        CABS (gPhincm[imic]),
        CREAL(gPhiscm[imic]),
        CIMAG(gPhiscm[imic]),
        CABS (gPhiscm[imic]),
        CREAL(gPhincm[imic]+gPhiscm[imic]),
        CIMAG(gPhincm[imic]+gPhiscm[imic]),
        CABS (gPhincm[imic]+gPhiscm[imic])
        );
  }
  logger(LOG_DEBUG,"Done.\n");

  if(0 == rank && NULL != gPhiscm) free(gPhiscm);
  if(0 == rank && NULL != gPhincm) free(gPhincm);
}