print_premic_on_master()

void print_premic_on_master

(

COMPLEX

cfreq

)

{
  FILE *f;
  char fname[MAX_PATH];

  int imics;
  COMPLEX phiscmval,phincmval;

  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;

  rho_s = cfreq*modsolinfo->rho;
  if (0 == rank) {
    gPhiscm = calloc(modgeominfo->nmics,sizeof(COMPLEX));
    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);

  /* Printing pressures */
  if(0== rank){

    // defining output file name
    get_filename(fname,"%s-mics-%.4fHz.out",rundetails->title,(double)HZFREQ(cfreq));
    f = fopen(fname,"w");

    fprintf(f, "# Run name  = %s\n",rundetails->title);
    fprintf(f, "# Run owner = %s\n",rundetails->owner);
    fprintf(f, "# Cfg file  = %s\n",rundetails->cfg_filename);
    fprintf(f, "# Frequency = %f + i*%f (%8.4f Hz)\n",(double)CREAL(cfreq),(double)CIMAG(cfreq), (double)HZFREQ(cfreq));
    fprintf(f, "#\n");
    fprintf(f, "# 1    2  3  4  5       6       7        8        9        10        11      12      13    \n");
    fprintf(f, "# imic xm ym zm re(inc) im(inc) abs(inc) re(scat) im(scat) abs(scat) re(tot) im(tot) abs(tot)\n");
    for(imics=0;imics<modgeominfo->nmics;imics++){
      phiscmval  = gPhiscm[imics];
      phincmval  = gPhincm[imics];
  
       fprintf(f,f1I12F,
           /*1*/imics,        // MIC NODES INFO
           /*2*/geometry->mics[imics].x,
           /*3*/geometry->mics[imics].y,
           /*4*/geometry->mics[imics].z,


           /*5*/CREAL(phincmval),  // INCIDENT
           /*6*/CIMAG(phincmval),  
           /*7*/CABS(phincmval) ,

           /*8*/CREAL( phiscmval ),  // SCATTERING 
           /*9*/CIMAG( phiscmval ),  
           /*10*/CABS( phiscmval ),

           /*11*/CREAL(phincmval + phiscmval),  // TOTAL 
           /*12*/CIMAG(phincmval + phiscmval),  
           /*13*/CABS (phincmval + phiscmval) 

           );
    }
    fclose(f);

    VecDestroy(gPhincm);
    VecDestroy(gPhiscm);
  }
}