vtkoutmmesh()

void vtkoutmmesh	(	int	icounter,
		COMPLEX	cfreq
	)

Routine vtkoutmmesh(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 POLYDATA'. THE SUBROUDTINE CAN BE USED ONLY IF A MESH TOPOLOGY IS PROVIDED FOR MICROPHONES. COMPATIBLE WITH GMSH MICS FORMAT. 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,i1,i2,i3,i4;
  FILE *f;
//  FILE *fb;
//  FILE *ft;
  char fname[MAX_PATH];
  //char fnameb[MAX_PATH];
  //char fnamet[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=7;
  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 on the Boundary\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", modgeominfo->nemics, 5*modgeominfo->nemics);
  for(i=0;i<modgeominfo->nemics;i++){
            i1  = geometry->jnodm[0 + i*4];
            i2  = geometry->jnodm[1 + i*4];
            i3  = geometry->jnodm[2 + i*4];
            i4  = geometry->jnodm[3 + i*4];
            fprintf(f," %d %d %d %d %d\n",4,i1,i2,i3,i4);
  }
  fprintf(f,"\nCELL_TYPES %d\n", modgeominfo->nemics);
  for(i=0;i<modgeominfo->nemics;i++){
          if (2 == geometry->kelmm[i] || 3 == geometry->kelmm[i]) celltype = 7;
          if (4 == geometry->kelmm[i]) celltype = 10;
          fprintf(f,"%d\n",celltype);
  }

  fprintf(f,"POINT_DATA %d\n",modgeominfo->nmics);

  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]);
          spl = 20*log10(prat) + 94.0;
          fprintf(f,f1F,spl);
  }

  fclose(f);

  VecDestroy(gPhincm);
  VecDestroy(gPhiscm);

}