printsol_on_master()

void printsol_on_master	(	int	icounter,
		COMPLEX	cfreq
	)

Prints solution on Master node Output file name is in the form: <run title>-surf-<frequency in Hertz>Hz.out.

Parameters

[in]	cfreq	COMPLEX frequency
[in]	icounter	frequency index

                                                   {
  FILE *file;
  char fname[MAX_PATH];
  int mpi_ncntr;
  COMPLEX phiinc,phisc,phitot;
  COMPLEX psiinc,psisc,psitot;
  DOUBLE rs,is;
  //Vec locSol;
  int nstart,nend;
  int i = 0;

  Vec gPhinc = NULL;
  Vec gPsinc = NULL;
  int gRows,mm,nn;
  int descl[9];
  int descg[9];
  int izero = 0;
  int ione = 1;
  int info;

  mpi_ncntr  = modgeominfo->ncntr + modgeominfo->nchief;

  calculate_bounds(rank,mpi_ncntr,size, &nstart,&nend);


  if (0 == rank) gPhinc  = calloc(modgeominfo->ncntr,sizeof(COMPLEX));
  if (0 == rank) gPsinc  = calloc(modgeominfo->ncntr,sizeof(COMPLEX));
  gRows = modgeominfo->ncntr;

  mm = gRows;
  nn = runinfo->npcols;

  // Gathering vectors on node 0 
  descinit_(descl   , &mm     ,&ione, &runinfo->row_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->vecPhinc ,&ione, &ione,descl, gPhinc ,&ione,&ione,descg,&runinfo->ctxt);
  pzgemr2d_(&gRows,&ione,solution->vecPsinc ,&ione, &ione,descl, gPsinc ,&ione,&ione,descg,&runinfo->ctxt);

  file = NULL;
  rs = 0.0;
  is = 0.0;
  if(0==rank){
    get_filename(fname,"%s-surf-%.4fHz.out",rundetails->title,(double)HZFREQ(cfreq));
    file = fopen(fname,"w");
    fprintf(file, "# Run name  = %s\n",rundetails->title);
    fprintf(file, "# Run owner = %s\n",rundetails->owner);
    fprintf(file, "# Cfg file  = %s\n",rundetails->cfg_filename);
    fprintf(file, "# Frequency = %f + i*%f (%8.4f Hz)\n",(double)CREAL(cfreq),(double)CIMAG(cfreq), (double)HZFREQ(cfreq));
    fprintf(file, "# \n");
    fprintf(file, "# 1    2    3    4    5       6       7        8        9        10        11      12      13\n");
    fprintf(file, "# inod xnod ynod znod re(inc) im(inc) abs(inc) re(scat) im(scat) abs(scat) re(tot) im(tot) abs(tot)\n");
    for(i=0;i<mpi_ncntr;i++){
        if (modsolinfo->dirneu == 1){
                          psiinc = gPsinc[i];
                          psisc  = (modsolinfo->knw == 1) ? solution->sol[i]-psiinc : solution->sol[i];
                          psitot = (modsolinfo->knw == 1) ? solution->sol[i]        : solution->sol[i] + psiinc;
                          fprintf(file,f1I12F,
                                   i,   // CONTROL POINT INFO
                                   geometry->cntr[i].x,
                                   geometry->cntr[i].y,
                                   geometry->cntr[i].z,

                                   CREAL(psiinc), // INCIDENT
                                   CIMAG(psiinc),
                                   CABS(psiinc),

                                   CREAL(psisc), // SCATTERING
                                   CIMAG(psisc),
                                   CABS(psisc),

                                   CREAL(psitot), // TOTAL
                                   CIMAG(psitot),
                                   CABS(psitot)
                                   );
        }else if(modsolinfo->dirneu == 2){
              phiinc = gPhinc[i];
              phisc  = (modsolinfo->knw == 1) ? solution->sol[i]-phiinc : solution->sol[i];
              phitot = (modsolinfo->knw == 1) ? solution->sol[i]        : solution->sol[i] + phiinc;
              fprintf(file,f1I12F,
                   i,   // CONTROL POINT INFO
                   geometry->cntr[i].x,
                   geometry->cntr[i].y,
                   geometry->cntr[i].z,

                   CREAL(phiinc), // INCIDENT
                   CIMAG(phiinc),
                   CABS(phiinc),

                   CREAL(phisc), // SCATTERING
                   CIMAG(phisc),
                   CABS(phisc),

                   CREAL(phitot), // TOTAL
                   CIMAG(phitot),
                   CABS(phitot)
                   );
        }
    }
    
    fclose(file);
    free(gPhinc);
    free(gPsinc);
  }
}