vtkoutmu2()

void vtkoutmu2	(	int	icounter,
		COMPLEX	cfreq
	)

Routine vtkoutmu2(icounter, cfreq). Produces the file 'runname.mics.fHz.vtk', which can bee used to visualize the solution at the microphones array. Arguments: 'icounter' is the frequency index (equal to ifreq); 'cfreq' is the complex frequency (the Laplace variable). Data format is 'DATASET UNSTRUCTURED_GRID'. AN UNSTRUCTURED GRID IS BUILT FROM THE MICS VECTOR. IT USES ONE SINGLE POLY-VERTEX CELL Scalar fields included are:

• Re Im and Abs of the scattered field
• Re Im and Abs of the total field
• Total pressure in Db
• Insertion Loss = 20 Log [(total field)/(incident field)]

Parameters

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

                                           {

  int i;
  FILE *f;
  char fname[MAX_PATH];
  COMPLEX val;
  Vec gPhiscm = NULL;
  Vec gPhincm = NULL;
  DOUBLE iloss,prat,spl;

  int locRows;
  int gRows;
  int mm,nn;
  int descl[9];
  int descg[9];
  int izero=0;
  int ione=1;
  int info, celltype, cellvrtx;

  celltype=2;
  cellvrtx=modgeominfo->nmics;

  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);

  if (0 != rank) return;
  get_filename(fname,"%s-mics-%.4fHz.vtk",rundetails->title,(double)HZFREQ(cfreq));
  f = fopen(fname,"w");

  fprintf(f,"# vtk DataFile Version 2.0\n");
  fprintf(f,"Solution at Mics\n");
  fprintf(f,"ASCII\n");
  fprintf(f,"DATASET UNSTRUCTURED_GRID\n");
  fprintf(f,"POINTS %d float \n",modgeominfo->nmics);
  for(i=0;i<modgeominfo->nmics;i++){
            fprintf(f,"  %f   %f   %f \n",(double)geometry->mics[i].x,(double)geometry->mics[i].y,(double)geometry->mics[i].z);
  }

  fprintf(f,"CELLS %d %d\n", 1, modgeominfo->nmics+1);
  fprintf(f,"%d",cellvrtx);
  for(i=0;i<modgeominfo->nmics;i++){
            fprintf(f," %d",i);
  }
  fprintf(f,"\nCELL_TYPES %d\n", 1);
  fprintf(f,"%d\n",celltype);

  fprintf(f,"POINT_DATA %d\n", gRows-0);

  fprintf(f,"SCALARS Re(Inc_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<modgeominfo->nmics;i++){
    val = gPhincm[i];
    fprintf(f,f1F,CREAL(val));
  }
  fprintf(f,"SCALARS Im(Inc_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<modgeominfo->nmics;i++){
    val = gPhincm[i];
    fprintf(f,f1F,CIMAG(val));
  }
  fprintf(f,"SCALARS Abs(Inc_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<modgeominfo->nmics;i++){
    val = gPhincm[i];
    fprintf(f,f1F,CABS(val));
  }

  fprintf(f,"SCALARS Re(Scat_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<modgeominfo->nmics;i++){
    val = gPhiscm[i];
    fprintf(f,f1F,CREAL(val));
  }
  fprintf(f,"SCALARS Im(Scat_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<modgeominfo->nmics;i++){
    val = gPhiscm[i];
    fprintf(f,f1F,CIMAG(val));
  }
  fprintf(f,"SCALARS Abs(Scat_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<modgeominfo->nmics;i++){
    val = gPhiscm[i];
    fprintf(f,f1F,CABS(val));
  }

  fprintf(f,"SCALARS Re(Total_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<gRows;i++){
    val = gPhiscm[i] + gPhincm[i];
    fprintf(f,f1F,CREAL(val));
  }
  fprintf(f,"SCALARS Im(Total_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<gRows;i++){
    val = gPhiscm[i] + gPhincm[i];
    fprintf(f,f1F,CIMAG(val));
  }
  fprintf(f,"SCALARS Abs(Total_Solution) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<gRows;i++){
    val = gPhiscm[i] + gPhincm[i];
    fprintf(f,f1F,CABS(val));
  }
  fprintf(f,"SCALARS Insertion_Loss float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<gRows;i++){
          prat  = CABS(gPhincm[i]+gPhiscm[i])/CABS(gPhincm[i]);
          iloss = 20*log10(1./prat);
          fprintf(f,f1F,iloss);
  }
  fprintf(f,"SCALARS SPL float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<gRows;i++){
          prat  = 0.70710678*CABS(gPhiscm[i]+gPhincm[i])/0.00002;
          spl = 20*log10(prat);
          fprintf(f,f1F,spl);
  }

  fclose(f);

  VecDestroy(gPhincm);
  VecDestroy(gPhiscm);

}