vtkoutb()

void vtkoutb	(	int	icounter,
		COMPLEX	cfreq
	)

Routine vtkoutb(icounter,CFREQ). Produces the file 'runname.surf.fHz.vtk', which can bee used to visualize the solution on the boundary. If dirneu==1 the file contains dphi/dn on S. If dirneu==2 the file contains phi on S. Arguments: 'icounter' is the frequency index (equal to ifreq); 'cfreq' is the complex frequency. Data format is 'DATASET POLYDATA'. Scalar fields included are:

• Re Im and Abs of the incident potential
• Re Im and Abs of the scattered potential
• Re Im and Abs of the total potential

Parameters

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

                                         {

  int i,j,iv;
  FILE *f;
  char fname[MAX_PATH];
  Vec gPhinc = NULL;
  Vec gPsinc = NULL;
  COMPLEX val;
  COMPLEX vali;

  int locRows;
  int gRows;
  int mm,nn;
  int descl[9];
  int descg[9];
  int izero=0;
  int ione=1;
  int nvrt=4;
  int info,ive1,ive2,ive3,ive4;

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

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

  // Gathering vectors on node 0 
  mm = gRows;// * runinfo->nprows;
  nn = runinfo->npcols;
  descinit_(descl   , &gRows  ,&ione, &runinfo->row_block_size ,&ione,&izero,&izero,&runinfo->ctxt,&locRows  ,&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);

  if (0 != rank) return;
  get_filename(fname,"%s-surf-%.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 POLYDATA\n");
  fprintf(f,"POINTS %d float \n",modgeominfo->nelmb*4);
  for(i=0;i<modgeominfo->nelmb;i++){
    fprintf(f,"  %f   %f   %f \n",(double)geometry->elements[i].x[0].x,(double)geometry->elements[i].x[0].y,(double)geometry->elements[i].x[0].z);
    fprintf(f,"  %f   %f   %f \n",(double)geometry->elements[i].x[1].x,(double)geometry->elements[i].x[1].y,(double)geometry->elements[i].x[1].z);
    fprintf(f,"  %f   %f   %f \n",(double)geometry->elements[i].x[2].x,(double)geometry->elements[i].x[2].y,(double)geometry->elements[i].x[2].z);
    fprintf(f,"  %f   %f   %f \n",(double)geometry->elements[i].x[3].x,(double)geometry->elements[i].x[3].y,(double)geometry->elements[i].x[3].z);
  }
  fprintf(f,"POLYGONS %d %d \n",modgeominfo->nelmb,5*(modgeominfo->nelmb));
  iv= 0;
  for(i=0;i<modgeominfo->nelmb;i++){
    ive1 = iv++;
    ive2 = iv++;
    ive3 = iv++;
    ive4 = iv++;
                fprintf(f,"%d   %d   %d   %d   %d\n",nvrt,ive1,ive2,ive3,ive4);
  }
  fprintf(f,"CELL_DATA %d\n", modgeominfo->nelmb);
  fprintf(f,"SCALARS Re(Incident_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    (1 == modsolinfo->dirneu) ? (val = gPsinc[j]) : (val = gPhinc[j]);
                    fprintf(f,f1F,CREAL(val));
          }
  }
  fprintf(f,"SCALARS Im(Incident_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    (1 == modsolinfo->dirneu) ? (val = gPsinc[j]) : (val = gPhinc[j]);
                    fprintf(f,f1F,CIMAG(val));
          }
  }
  fprintf(f,"SCALARS Abs(Incident_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    (1 == modsolinfo->dirneu) ? (val = gPsinc[j]) : (val = gPhinc[j]);
                    fprintf(f,f1F,CABS(val));
          }
  }
  fprintf(f,"SCALARS Re(Scattered_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    val = solution->sol[j];
                    (1 == modsolinfo->dirneu) ? (vali = gPsinc[j]) : (vali = gPhinc[j]);
                    if(1 == modsolinfo->knw)val -= vali;
                    fprintf(f,f1F,CREAL(val));
          }
  }
  fprintf(f,"SCALARS Im(Scattered_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    val = solution->sol[j];
                    (1 == modsolinfo->dirneu) ? (vali = gPsinc[j]) : (vali = gPhinc[j]);
                    if(1 == modsolinfo->knw)val -= vali;
                    fprintf(f,f1F,CIMAG(val));
          }
  }
  fprintf(f,"SCALARS Abs(Scattered_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    val = solution->sol[j];
                    (1 == modsolinfo->dirneu) ? (vali = gPsinc[j]) : (vali = gPhinc[j]);
                    if(1 == modsolinfo->knw)val -= vali;
                    fprintf(f,f1F,CABS(val));
          }
  }
  fprintf(f,"SCALARS Re(Total_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    val = solution->sol[j];
                    (1 == modsolinfo->dirneu) ? (vali = gPsinc[j]) : (vali = gPhinc[j]);
                    if(2 == modsolinfo->knw)val += vali;
                    fprintf(f,f1F,CREAL(val));
          }
  }
  fprintf(f,"SCALARS Im(Total_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    val = solution->sol[j];
                    (1 == modsolinfo->dirneu) ? (vali = gPsinc[j]) : (vali = gPhinc[j]);
                    if(2 == modsolinfo->knw)val += vali;
                    fprintf(f,f1F,CIMAG(val));
          }
  }
  fprintf(f,"SCALARS Abs(Total_Field) float\n");
  fprintf(f,"LOOKUP_TABLE default\n");
  for(i=0;i<runinfo->nsymm;i++){
          for(j=0;j<modgeominfo->ncntr;j++){
                    val = solution->sol[j];
                    (1 == modsolinfo->dirneu) ? (vali = gPsinc[j]) : (vali = gPhinc[j]);
                    if(2 == modsolinfo->knw)val += vali;
                    fprintf(f,f1F,CABS(val));
          }
  }

  fclose(f);

  VecDestroy(gPhinc);
  VecDestroy(gPsinc);

}